[R] maximum number of subdivisions reached
. .
xkziloj at gmail.com
Mon Aug 29 17:16:08 CEST 2011
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.
More information about the R-help
mailing list