I have a Bayesian MCMC in R, and I have the code below:
RWM = function(Niter,Y,X){ p = ncol(X) alpha = 0.7 beta = matrix(0,ncol=1,nrow=3) beta