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.weyla...@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, . . <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. > > ______________________________________________ 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.