powered by
This function runs a simple Metropolis sampler with standard normal target distribution and uniform innovations.
metrop(n, alpha)
An R vector containing the output of the sampler.
The number of iterations of the Metropolis sampler.
The tuning parameter of the sampler. The innovations of the sampler are of the form U(-alpha,alpha).
normgibbs
normvec=metrop(1000,1) op=par(mfrow=c(2,1)) plot(ts(normvec)) hist(normvec,20) par(op)
Run the code above in your browser using DataLab