# Algorithm underlying omxMnor

5 posts / 0 new
Offline
Joined: 03/19/2022 - 09:48
Algorithm underlying omxMnor

Dear all,

I am currently using omxMnor to calculate the multivariate normal integral in my paper. Wonder what is the underlying algorithm implemented by omxMnor. Any references would be greatly appreciated!

Cheers,
Ou

Offline
Joined: 01/24/2014 - 12:15

The algorithm underlying omxMnor() is a FORTRAN subroutine called 'SADMVN'. It starts on line 558 of src/sadmvn.f in the OpenMx source repository. SADMVN's code contains the following comment, which provides a citation for the algorithm:

*     A subroutine for computing multivariate normal probabilities.
*     This subroutine uses an algorithm given in the paper
*     "Numerical Computation of Multivariate Normal Probabilities", in
*     J. of Computational and Graphical Stat., 1(1992), pp. 141-149, by
*          Alan Genz
*          Department of Mathematics
*          Washington State University
*          Pullman, WA 99164-3113
*          Email : alangenz@wsu.edu

Professor Genz has a PDF of that paper on his website here.

Offline
Joined: 03/19/2022 - 09:48

This is really helpful! Thanks heaps!!!

Offline
Joined: 01/24/2014 - 12:15

You're welcome. Glad to be of help.

Offline
Joined: 07/31/2009 - 15:26
Reference in help page

I've added the Genz reference to the help page of the omxMnor() function (i.e., man page obtained from ?omxMnor).