I couldn't find a simple R code for random-walk Metropolis sampling (the symmetric proposal version of Metropolis Hastings sampling) from a multivariate target distribution in arbitrary dimensions, so I wrote one. This is also my first R code. It requires the package MASS to sample from the multivariate normal proposal distribution using the mvrnorm function. …

Continue reading R code for multivariate random-walk Metropolis sampling