

Author: Brodtkorb Per
Publisher: Springer Publishing Company
ISSN: 1387-5841
Source: Methodology And Computing In Applied Probability, Vol.8, Iss.1, 2006-03, pp. : 65-91
Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.
Abstract
The numerical computation of expectations for (nearly) singular multivariate normal distribution is a difficult problem, which frequently occurs in widely varying statistical contexts. In this article we discuss several strategies to improve the algorithm proposed by Genz and Kwong (2000) when either a moderate accuracy is requested, the correlation structure is strong, and, most importantly, the dimension of the integral is large. Test results for typical problems show an average speedup of 10 using the modified algorithm, but even more is gained as the dimension of the problem increases. We apply the modified algorithm to compute long-run distributions of Gaussian wave characteristics, a difficult problem where previous algorithms fail to compute accurate values in reasonable time.
Related content


The chi squared and multinormal distributions
JOURNAL OF CHEMOMETRICS (ELECTRONIC), Vol. 29, Iss. 1, 2015-01 ,pp. :


A new transformation technique for evaluating nearly singular integrals
By Ye Wenjing
Computational Mechanics, Vol. 42, Iss. 3, 2008-08 ,pp. :


STUDIES IN APPLIED MATHEMATICS, Vol. 97-1421, Iss. 1, 1938-04 ,pp. :


STUDIES IN APPLIED MATHEMATICS, Vol. 17, Iss. 1, 1938-04 ,pp. :


Singular deformations of nearly R2 inflation potentials
Classical and Quantum Gravity, Vol. 32, Iss. 23, 2015-12 ,pp. :