Message-ID: <AANLkTimav-+PPP_cDH5udwAWAZa1gNMRO_6zuEKk86q9@mail.gmail.com>
Date: 2010-10-25T15:55:25Z
From: Carlos J. Gil Bellosta
Subject: [R-es] Ayuda con una función
In-Reply-To: <AANLkTinMyZoyMV9=eoVTKbMFgFn7NGPG2yX1MimsO--y@mail.gmail.com>
Hola, ¿qué tal?
Utiliza un valor mayor (¿1000?) para el parámetro "subdivisions" de "integrate".
Alternativamente, no uses integración numérica. Puedes reescribir tu
integral como una suma de unos cuantos términos (parte entera de y más
uno) habida cuenta de que ppois es una función de escalonada (con
escalones en los naturales).
Un saludo,
Carlos J. Gil Bellosta
http://www.datanalytics.com
El dÃa 25 de octubre de 2010 17:32, Cesar Escalante
<c.escalante.c en gmail.com> escribió:
> Reciban un cordial saludo.
>
> De manera muy comedida les pido ayuda con la función indicada abajo. La usan
> en un artÃculo de inventarios. No veo dónde estoy errado.
>
> Gracias de antemano.
>
> César Escalante C.
>
>> lamb<-50
>> L<-1
>> h.<-10
>> p<-25
>> K<-5
>> integ<-function(y){integrate(function(x)ppois(x,lamb*L),0,y)$value}
>> G<-function(y){(h.+p)*integ(y) + p*(lamb*L-y)}
>> G(0:43)
> Â [1] 1250 1225 1200 1175 1150 1125 1100 1075 1050 1025 1000 Â 975 Â 950 Â 925
> 900 Â 875 Â 850 Â 825 Â 800 Â 775 Â 750 Â 725 Â 700 Â 675 Â 650 Â 625 Â 600 Â 575 Â 550
> 525
> [31] Â 500 Â 475 Â 450 Â 425 Â 400 Â 375 Â 350 Â 325 Â 300 Â 275 Â 250 Â 225 Â 200 Â 175
>> G(44)
> Error en integrate(function(x) ppois(x, lamb * L), 0, y) :
> Â maximum number of subdivisions reached
>> G(500)
> [1] 4499.785
>
> Â Â Â Â [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>