Description Usage Arguments Value References See Also Examples
A Robbins-Munro stochastic approximation update is used to adapt the tuning parameter of the proposal kernel. The idea is to update the tuning parameter at each iteration of the sampler:
h^{(i+1)} = h^{(i)} + η^{(i+1)}(α^{(i)} - α_{opt}),
where h^{(i)} and α^{(i)} are the tuning parameter and acceptance probability at iteration i and α_{opt} is a target acceptance probability. For Gaussian targets, and in the limit as the dimension of the problem tends to infinity, an appropriate target acceptance probability for MALA algorithms is 0.574. The sequence {η^{(i)}} is chosen so that ∑_{i=0}^∞η^{(i)} is infinite whilst ∑_{i=0}^∞≤ft(η^{(i)}\right)^{1+ε} is finite for ε>0. These two conditions ensure that any value of h can be reached, but in a way that maintains the ergodic behaviour of the chain. One class of sequences with this property is,
η^{(i)} = \frac{C}{i^α},
where α\in(0,1] and C>0.The scheme is set via
the mcmcpars
function.
1 | andrieuthomsh(inith, alpha, C, targetacceptance = 0.574)
|
inith |
initial h |
alpha |
parameter α |
C |
parameter C |
targetacceptance |
target acceptance probability |
an object of class andrieuthomsh
Andrieu C, Thoms J (2008). A tutorial on adaptive MCMC. Statistics and Computing, 18(4), 343-373.
Robbins H, Munro S (1951). A Stochastic Approximation Methods. The Annals of Mathematical Statistics, 22(3), 400-407.
Roberts G, Rosenthal J (2001). Optimal Scaling for Various Metropolis-Hastings Algorithms. Statistical Science, 16(4), 351-367.
mcmcpars, lgcpPredict
1 | andrieuthomsh(inith=1,alpha=0.5,C=1,targetacceptance=0.574)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.