Title

Numerical computation of multivariate normal and multivariate-t probabilities over convex regions

Authors

Authors

P. N. Somerville

Comments

Authors: contact us about adding a copy of your work at STARS@ucf.edu

Abbreviated Journal Title

J. Comput. Graph. Stat.

Keywords

convex region; integral equation; Monte Carlo; multiple comparisons; multivariate normal integral; multivariate-t integral; INTEGRALS; Statistics & Probability

Abstract

A methodology has been developed and Fortran 90 programs have been written to evaluate multivariate normal and multivariate-t integrals over convex regions. The Cholesky transformation is used to transform the integrand into a product of standard normal or spherically symmetric t variables. For any random direction from the origin, an unbiased estimate of the value of the integral is Pr[chi(2) less than or equal to r(2)] (multivariate normal) or Pr[F less than or equal to r(2)/k] (multivariate-t), where r is the distance from the origin to the boundary in a randomly chosen direction, and k is the dimension of the integral. Two Fortran 90 programs have been written. MVI uses the average of many estimates. MVIB uses a binning procedure to obtain an empirical distribution of the distance from the origin to the boundary. Gauss-Legendre quadrature is then used to estimate the value of the integral. The running time for MVIB is modestly smaller than that for MVI. However, in solving certain integral equations (e.g., using an iterative procedure to find the percentage point of a statistic); using MVIB usually requires no Monte Carlo sampling after the first iteration, and is considerably more efficient. MVIB and MVI are highly efficient for the evaluation of integrals whose value is large. "Naive" Monte Carlo (MC) may be competitive with MVI or MVIB only if the value of the probability integral is small or the shape of the region is "extreme."

Journal Title

Journal of Computational and Graphical Statistics

Volume

7

Issue/Number

4

Publication Date

1-1-1998

Document Type

Article

Language

English

First Page

529

Last Page

544

WOS Identifier

WOS:000077362100008

ISSN

1061-8600

Share

COinS