adapump {mcsm} | R Documentation |
This function constructs an update of the location and scale matrix of a normal proposal in a Metropolis-Hastings algorithm, based on earlier simulations in order to show the danger of online adaptivity.
adapump(T = 10^2, MM = 10^2)
T |
Number of steps between updates |
MM |
Number of updates, leading to a chain of length MM*T |
The function simply plots the sequence of beta's along iterations, which should collapse, as well as the range of the variability of the proposed values.
Christian P. Robert and George Casella
Chapter 8 of EnteR Monte Carlo Statistical Methods
adapump(T=50,MM=50)