These days I use SAS mixed model to compute the variances.( For example, proc mixed data=phe convf=1e-8 maxiter=50; )Through google, I know some R packages, like lme4 and nlme. But I can not get the results like SAS. Now I want to know the corresponding R package having the similar function with SAS mixed model. Thank you very much!