The application of Gibbs sampling to nested variance components models with heterogeneous within-group variance