Page 1 of 1

Mode jumping MCMC

PostPosted: Sat Feb 01, 2014 5:42 pm
by ste_s
Hi,
I am struggling in getting estimated parameters with mode jumping MCMC. As a proposal density I tried either an identity matrix and the prior variance but there seems to be a scaling problem (see priors&posteriors). I am attaching log files under both proposal densities. The model is baseline Smets & Wouters and the estimation works regularly using the default Hessian at a posterior mode.
According to my little understanding, it looks like the MCMC sampler is having hard time to cross regions where the density is very low but the reason why one may use jumping MCMC is to address bimodality (e.g. the cases discussed here: http://www.dynare.org/phpBB3/viewtopic.php?f=1&t=5174). Is that correct?
Can you help me to get the jumping MCMC working?
Do you have any advice in how to uncover genuine multimodality in the model?
Thank you very much

Re: Mode jumping MCMC

PostPosted: Wed Feb 05, 2014 9:07 am
by jpfeifer
Hi, could it be that you are confusing something? The
Code: Select all
mcmc_jumping_covariance

option is for specifying a the covariance for the standard Random-Walk Metropolis-Hastings sampler. The reason I used that command was that it allowed me to scale the proposal density to a size that allowed jumping to a differerent mode. Also I knew where the second mode was so that I had a feeling how large the scaling needed to be. While this wide proposal density allowed crossing regions, the sampler became hugely inefficient. For higher-dimensional problems, this is completely impractical.

As should be clear, this procedure is not a mode-jumping MCMC. While we would like to have something like that in Dynare, it is not yet there.

Regarding finding multimodality: my only advice would be to start mode-finding from different starting points and see where it converges to. Slices through the likelihood then provide an indication whether that is actually a mode.