plmm function

Fit a linear mixed model via penalized maximum likelihood.