Page 1 of 1

[astro-ph/0405462] Fast and reliable MCMC for cosmological p

Posted: September 24 2004
by Antony Lewis
The result for scaling the proposal density with dimension given in the paper appears to be generally useful and correct. Two comments:

1. There are some potential problems using the autocorrelation function to assess chain efficiency: it is quite possible for the autocorrelation to be small but the separation between effectively independence samples to be large. Jo has confirmed that similar problems can happen with the spectral estimator.

2. Results are restricted to using N-D Gaussians as the proposal density. In high-N (N>2) this is probably rather un-robust, and it may be better to use a lower dimensional proposal.

I've written up these comments in detail in the cosmomc notes

http://cosmologist.info/notes/cosmomc.ps.gz

Similar comments apply to [arxiv]astro-ph/0310723[/arxiv]

[astro-ph/0405462] Fast and reliable MCMC for cosmological p

Posted: November 22 2004
by Michael Doran
The scaling is indeed useful (and easy to implement). Funny enough, the adaptive
stepsize I used in AnlalyzeThis! provided roughly the same acceptance rate.
But as the adaptive stepping in AnalyzeThis! stops after "sensing" the distribution (i.e. some extended burn-in), the formula is still useful for the step-proposal after "burn-in".