[R] Bug of pmvnorm?

Daniel Yang yungjui at gate.sinica.edu.tw
Sat May 27 09:31:00 CEST 2006


Dear R-help,

Is the following a bug of pmvnorm?

R<-matrix(c(1,.5,.5,1),2,2)
pmvnorm(lower=c(-Inf,-Inf),upper=c(Inf,Inf),corr=R)

It returns,

[1] 0
attr(,"error")
[1] 0
attr(,"msg")
[1] "Normal Completion"

But I think it should return 1.

Yung-jui



More information about the R-help mailing list