NPML estimation for random effect models in R8 years ago