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

    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