working on new gibbs sampler for multiple rho modle