Hi when using lme or nlme, I could use intervals() to obtain 95%CI for sd estimates of random effects. However intervals does not work with lmer. is there any equivalent function? or is there any way to obtain those 95%CI with lmer? Thanks Julien