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]]

