Error message
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.hlx at 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 at 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.