Adaptive radial-based direction sampling; Some flexible and robust Monte Carlo integration methods
Adaptive radial-based direction sampling (ARDS) algorithms are specified for Bayesian analysis of models with nonelliptical, possibly, multimodal target distributions. A key step is a radial-based transformation to directions and distances. After the transformations a Metropolis-Hastings method or, alternatively, an importance sampling method is applied to evaluate generated directions. Next, distances are generated from the exact target distribution by means of the numerical inverse transformation method. An adaptive procedure is applied to update the initial location and covariance matrix in order to sample directions in an efficient way. Tested on a set of canonical mixture models that feature multimodality, strong correlation, and skewness, the ARDS algorithms compare favourably with the standard Metropolis-Hastings and importance samplers in terms of flexibility and robustness. The empirical examples include a regression model with scale contamination and a mixture model for economic growth of the USA.
|Keywords||Markov chain Monte Carlo, importance sampling, radial coordinates|
|Sponsor||HCM grant ERBCHRXCT 940514 of EC|
Bauwens, L., Bos, C.S., van Dijk, H.K., & van Oest, R.D.. (2003). Adaptive radial-based direction sampling; Some flexible and robust Monte Carlo integration methods (No. EI 2003-22). Retrieved from http://hdl.handle.net/1765/1722