NPML estimation for random effect models in R