Dear All, Herewith, is a sample of mixed model command from the help menu. fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age fm1 I noticed that there is no standard error for the fixed effects. Could someone advice me the way to extract it? Thank you Fir