[Major] remove dependency on MCMCpack (the MCMCmetrop1R function use fixed seed for random number generation, which is problematic for Metropolis-within-Gibbs sampling) and replace with hand written MH algorithm