Optimal direction Gibbs sampler for truncated multivariate normal distributions

J. Andrés Christen, Colin Fox, Mario Santana-Cibrian
2015 Communications in statistics. Simulation and computation  
Generalized Gibbs samplers simulate from any direction, not necessarily limited to the coordinate directions of the parameters of the objective function. We study how to optimally choose such directions in a random scan Gibbs sampler setting. We consider that optimal directions will be those that minimize the Kullback-Leibler divergence of two Markov chain Monte Carlo steps. Two distributions over direction are proposed for the multivariate Normal objective function. The resulting algorithms
more » ... used to simulate from a truncated multivariate Normal distribution, and the performance of our algorithms is compared with the performance of two algorithms based on the Gibbs sampler.
doi:10.1080/03610918.2015.1053926 fatcat:kjqd4xdk6fhevbpkkfqne6zcny