[R] SAS LMM to lme syntax

Pedro Mardones mardones.p at gmail.com
Sat Dec 20 22:28:23 CET 2008


Dear all;

I have the following SAS code:

proc mixed data=mydata covtest asycov update;
class rep cross;
 model ht = rep;
 random col1-col5 / type=toep(1);
 random cross/type=vc;
ods output covparms=_varcomp asycov=_cov;
run;

where col1-col5 are columns of 0 and 1 (Z matrix of dummy variables)

what I'm looking for is some hint of how to replicate this analysis
using lme in R, especifically how to translate the variance structures
used for the random effects: the commands type=toep(1) (banded
Toeplitz in SAS terminology) and type=vc (variance components in SAS
terminology).

please let me know if you have any suggestion

thanks in advance

PM



More information about the R-help mailing list