×

An accept-reject algorithm for the positive multivariate normal distribution. (English) Zbl 1306.65032

Summary: The need to simulate from a positive multivariate normal distribution arises in several settings, specifically in Bayesian analysis. A variety of algorithms can be used to sample from this distribution, but most of these algorithms involve Gibbs sampling. Since the sample is generated from a Markov chain, the user has to account for the fact that sequential draws in the sample depend on one another and that the sample generated only follows a positive multivariate normal distribution asymptotically. The user would not have to account for such issues if the sample generated was i.i.d. In this paper, an accept-reject algorithm is introduced in which variates from a positive multivariate normal distribution are proposed from a multivariate skew-normal distribution. This new algorithm generates an i.i.d. sample and is shown, under certain conditions, to be very efficient.

MSC:

62-08 Computational methods for problems pertaining to statistics

Software:

QSIMVN
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Albert JH, Chib S (1993) Bayesian analysis of binary and polychotomous response data. J Am Stat Assoc 88:669-679 · Zbl 0774.62031 · doi:10.1080/01621459.1993.10476321
[2] Chen M, Deely J (1996) Bayesian analysis for a constrained linear multiple regression problem for predicting the new crop of apples. J Agric Biol Environ Stat 1:467-489 · doi:10.2307/1400440
[3] Chib S (1992) Bayes inference in the Tobit censored regression model. J Econom 51:79-99 · Zbl 0742.62033 · doi:10.1016/0304-4076(92)90030-U
[4] Chib S, Greenberg E (1998) Analysis of multivariate probit models. Biometrika 85:347-361 · Zbl 0938.62020 · doi:10.1093/biomet/85.2.347
[5] Damien P, Walker SG (2001) Sampling truncated normal, beta, and gamma densities. J Comput Graph Stat 10:206-215 · doi:10.1198/10618600152627906
[6] Gelfand AE, Smith AFM, Lee TM (1992) Bayesian analysis of constrained parameter and truncated data problems using Gibbs sampling. J Am Stat Assoc 87:523-532 · doi:10.1080/01621459.1992.10475235
[7] Genz A (1992) Numerical computation of multivariate normal probabilities. J Comput Graph Stat 1:141-150
[8] Genz A (1993) Comparison of methods for the computation of multivariate normal probabilities. Comput Sci Stat 25:400-405
[9] Geweke J (1991) Efficient Simulation From the Multivariate Normal and Student t-distributions Subject to Linear Constraints. Computer Science and Statistics: Proceedings of the 23rd Symposium Interface, pp 571-577 · Zbl 0742.62033
[10] Gupta AK, González-Farías G, Domínguez-Molina JA (2004) A multivariate skew normal distribution. J Multivar Anal 89:181-190 · Zbl 1036.62043
[11] Hajivassiliou VA, McFadden D, Ruud PA (1996) Simulation of multivariate normal rectangle probabilities and their derivatives: theoretical and computational results. J Econom 72:85-134 · Zbl 0849.62064 · doi:10.1016/0304-4076(94)01716-6
[12] Linardakis M, Dellaportas P (2003) Assessment of Athen’s metro passenger behavior via a multiranked probit model. J R Stat Soc Ser C 52:185-200 · Zbl 1111.62393 · doi:10.1111/1467-9876.00397
[13] Liu X, Daniels MJ, Marcus B (2009) Joint models for the association of longitudinal binary and continuous processes with application to a smooking cessation trial. J Am Stat Assoc 104:429-438 · Zbl 1388.62328 · doi:10.1198/016214508000000904
[14] Mira A, Moller J, Roberts GO (2001) Perfect slice samplers. J R Stat Soc Ser B 63:593-606 · Zbl 0993.65015 · doi:10.1111/1467-9868.00301
[15] Phillipe A, Robert CA (2003) Perfect simulation of positive Gaussian distributions. Stat Comput 13:179-186 · doi:10.1023/A:1023264710933
[16] Pitt M, Chan D, Kohn R (2006) Efficient Bayesian inference for Gaussian copula regression models. Biometrika 93:537-554 · Zbl 1108.62027 · doi:10.1093/biomet/93.3.537
[17] Robert CP (1995) Simulation of truncated normal variables. Stat Comput 5:121-125 · Zbl 0866.90030 · doi:10.1007/BF00143942
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.