data(wafers)
m1 <- HLfit(y ~X1+X2+(1|batch),
resid.formula = ~ 1 ,data=wafers,HLmethod="ML")
fixef(m1)
vcov(m1)
ranef(m1)
## see 'predict' for an example with predictionCoeffs and getDistMat
Run the code above in your browser using DataLab