[R] maximum number of subdivisions reached

[. .]

Ooops,

sorry!

The problem occurs when

func(1:2,0.1,0.1,sad=Exp)

On Mon, Aug 29, 2011 at 12:27 PM, R. Michael Weylandt
<michael.weylandt at gmail.com> wrote:
> Can't help, code runs fine on my machine once you change "valu" to "value."
> Are you sure it fails in a vanilla run of R and isn't caused by any other
> choices you have made along the way?
>
> Michael
>
> PS -- Here's the code
>
> func <- function(y, a, rate, sad){
>     f3 <- function(z){
>         f1 <- function(y,a,n){
>             dpois(y,a*n)
>         }
>
>         f2 <- function(n,rate){
>             dexp(n,rate)
>         }
>
>         f <- function(n){
>             f1(y,a,n)*f2(n,rate)
>         }
>
>         r <- 0
>         r1 <- 1
>         x1 <- 0
>         dx <- 20
>
>         while(r1 > 10e-500){
>             r1 <- integrate(f,x1,x1+dx)$value
>             r <- r + r1
>             x1 <- x1 + dx
>         }
>
>         r + integrate(f,x1,Inf)$value
>     }
>  sapply(y,f3)
> }
>
> V = func(200,0.1,0.1,sad=Exp)
>
> On Mon, Aug 29, 2011 at 11:16 AM, . . <xkziloj at gmail.com> wrote:
>>
>> Why I am getting
>>
>> Error in integrate(f, x1, x1 + dx) :
>>  maximum number of subdivisions reached
>>
>> and can I avoid this?
>>
>> func <- function(y, a, rate, sad){
>> f3 <- function(z){
>>  f1 <- function(y,a,n){
>>    dpois(y,a*n)
>>  }
>>  f2 <- function(n,rate){
>>    dexp(n,rate)
>>  }
>>  f <- function(n){
>>    f1(y,a,n)*f2(n,rate)
>>  }
>>  r <- 0
>>  r1 <- 1
>>  x1 <- 0
>>  dx <- 20
>>  while(r1 > 10e-500){
>>  r1 <- integrate(f,x1,x1+dx)$value
>>    r <- r + r1
>>    x1 <- x1 + dx
>>  }
>>  r + integrate(f,x1,Inf)$valu
>> }
>>  sapply(y,f3)
>> }
>> func(200,0.1,0.1,sad=Exp)
>>
>> Thanks in advance.
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
>