The Normal Law Under Linear Restrictions: Simulation and Estimation via Minimax Tilting

Abstract

Simulation from the truncated multivariate normal distribution in high dimensions is a recurrent problem in statistical computing and is typically only feasible by using approximate Markov chain Monte Carlo sampling. We propose a minimax tilting method for exact independently and identically distributed data simulation from the truncated multivariate normal distribution. The new methodology provides both a method for simulation and an efficient estimator to hitherto intractable Gaussian integrals. We prove that the estimator has a rare vanishing relative error asymptotic property. Numerical experiments suggest that the scheme proposed is accurate in a wide range of set-ups for which competing estimation schemes fail. We give an application to exact independently and identically distributed data simulation from the Bayesian posterior of the probit regression model.

The next hundred years of growth and convergence

“…The uniform priors on θ c and θ d lead to truncated normal posteriors for θ c and θ d . Somewhat surprisingly, the draws from the truncated normal sometimes run into numerical difficulties, which are solved by using the algorithm given in Botev () as well as code Botev makes publicly available. The priors for the precision are diffuse, with the prior mean for each country set to the unconditional variance of and the prior standard deviations set to four times the prior means.…”

Optical Jamming Enhances the Secrecy Performance of the Generalized Space-Shift-Keying-Aided Visible-Light Downlink

“…is a time-varying jamming signal vector, whose entries are from a real-valued truncated Gaussian distribution [42] confined to the interval of [−…”

Efficient Simulation for Dependent Rare Events with Applications to Extremes

“…More effective estimators can be constructed if we use βn to estimate terms from α n−1 (10). We label the random terms in αn as…”

“…Sampling from X i | X i > γ can be easily done by acceptance-rejection with shifted exponential proposals [35] (or by inverse transform sampling [7,Remark 2.4], though this can be problematic using only double precision arithmetic). To simulate (X i , X j ) | min{X i , X j } > γ we use Botev's Matlab library [10], but also remark that a Gibb's sampler is a commonly used alternative [11,35].…”