Re: [R] Computing multivariate normal probabilities. Was: Re: Problem with Numerical derivatives (numDeriv) and mvtnorm

2009-11-23 Thread SL
Hi Torsten, Thanks for you comment. If you have some free time to spare, partial derivatives with respect to bounds and correlation coefficients would be great for pmvnorm! In complex problems, optim is not very good at estimating the hessian numerically and first order derivatives help to build

Re: [R] Computing multivariate normal probabilities. Was: Re: Problem with Numerical derivatives (numDeriv) and mvtnorm

2009-11-23 Thread Torsten Hothorn
On Sun, 22 Nov 2009, Ravi Varadhan wrote: Hi Torsten, Hi Ravi, It would be useful to "warn" the users that the multivariate normal probability calculated by "pmvnorm" using the GenzBretz algorithm is "random", i.e. the result can vary between repeated executions of the function. only i

Re: [R] Computing multivariate normal probabilities. Was: Re: Problem with Numerical derivatives (numDeriv) and mvtnorm

2009-11-22 Thread stephane Luchini
I'm now making some trials with sadmvn which provides results similar to pmvnorm for optimization but I know compute my OPG estimator of the covariance matrix with sadmvn (by the way Ravi, when I was refering to "exist in theory" I was refering to the theory not to the computation - would an approp

[R] Computing multivariate normal probabilities. Was: Re: Problem with Numerical derivatives (numDeriv) and mvtnorm

2009-11-22 Thread Ravi Varadhan
Hi Torsten, It would be useful to "warn" the users that the multivariate normal probability calculated by "pmvnorm" using the GenzBretz algorithm is "random", i.e. the result can vary between repeated executions of the function. This would prevent inappropriate use of pmvnorm such as computin