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.
