Implementation of metropolis function, p. 375.
metropolis(theta.matrix,reps,I.mat)
theta.matrix
reps
I.mat