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, . . <xkzi...@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@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. > [[alternative HTML version deleted]] ______________________________________________ R-help@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.