some changes in lme4 and correction on G marix calculation in MCMC