Question: 17. Consider three approaches for sampling from a N(0, 1) target density, p(x) exp(1 2x2): a standard Metropolis algorithm using a Gaussian proposal density,

17. Consider three approaches for sampling from a N(0, 1) target density, p(x) ∝ exp(−1 2x2):

a standard Metropolis algorithm using a Gaussian proposal density, x∗|x(t−1) ∼ N(x(t−1), σ2), the simple slice sampler (3.20), and a Langevin-Hastings algorithm (3.18).

(a) Evaluate and compare the Metropolis and Langevin-Hastings acceptance ratios, the latter of which is given by (3.19).

(b) Write an R program to run single chains of 1000 iterations for each of these samplers, starting both chains at x(0) = 0. For which values of the Metropolis proposal variance σ2 do the various samplers produce the smallest lag 1 autocorrelations (i.e., faster convergence)?

Step by Step Solution

There are 3 Steps involved in it

1 Expert Approved Answer
Step: 1 Unlock blur-text-image
Question Has Been Solved by an Expert!

Get step-by-step solutions from verified subject matter experts

Step: 2 Unlock
Step: 3 Unlock

Students Have Also Explored These Related Categorical Data Analysis Questions!