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.