most likely either 'lower' or 'upper' is NA. put
options(error = recover)
in your script to stop on the error and examine the value. you need to learn
"debugging 101" to help yourself out.
Sent from my iPad
On Mar 5, 2013, at 16:00, li li <hannah....@gmail.com> wrote:
> Dear all,
> I got an error message when running the following code.
> Can anyone give any suggestions on fixing this type of error?
> Thank you very much in advance.
> Hanna
>
>
>
>> integrand <- function(x, rho, a, b, z){
> + x1 <- x[1]
> + x2 <- x[2]
> + Sigma <- matrix(c(1, rho, rho, 1), 2,2)
> + mu <- rep(0,2)
> + f <- pmnorm(c((z-a*x1)/b, (z-a*x2)/b), mu,
> Sigma)*dmnorm(c(0,0), mu, diag(2))
> + f
> + }
>>
>> adaptIntegrate(integrand, lower=rep(-Inf, 2), upper=c(2,2),
> + rho=0.1, a=0.6, b=0.3, z=3, maxEval=10000)
> Error in if (any(lower > upper)) stop("lower>upper integration limits") :
> missing value where TRUE/FALSE needed
>
> [[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.
______________________________________________
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.
