[R-es] Ayuda con una función
Carlos J. Gil Bellosta
cgb en datanalytics.com
Lun Oct 25 17:55:25 CEST 2010
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
>
>
Más información sobre la lista de distribución R-help-es