Fast simulation of truncated Gaussian distributions

  • Authors:
  • Nicolas Chopin

  • Affiliations:
  • ENSAE-CREST, Malakoff Cedex, France 92245

  • Venue:
  • Statistics and Computing
  • Year:
  • 2011

Quantified Score

Hi-index 0.00

Visualization

Abstract

We consider the problem of simulating a Gaussian vector X, conditional on the fact that each component of X belongs to a finite interval [a i ,b i ], or a semi-finite interval [a i ,+驴). In the one-dimensional case, we design a table-based algorithm that is computationally faster than alternative algorithms. In the two-dimensional case, we design an accept-reject algorithm. According to our calculations and numerical studies, the acceptance rate of this algorithm is bounded from below by 0.5 for semi-finite truncation intervals, and by 0.47 for finite intervals. Extension to three or more dimensions is discussed.