I noticed a message when i try variance analysis with more groups that the GRR analysis switch to using bayesian estimates because of negative REML variance. Was wondering if there is code or formula of how to calculate the variance using the bayesian variance component analysis?