Displaying 1 result from an estimated 1 matches for "variable_id".
Did you mean:
variable_1
2007 Apr 06
0
translating sas proc mixed to lme()
...ng to translate a proc mixed into a lme() syntax. It seems that I was
able to do it for part of the model, but a few things are still different.
It is a 2-level bivariate model (some call it a pseudo-3-level model).
PROC MIXED DATA=psdata.bivar COVTEST METHOD = ml;
CLASS cluster_ID individual_id variable_id ;
MODEL y = Dp Dq / SOLUTION NOINT;
RANDOM Dp Dq / SUBJECT = cluster_ID TYPE=UN G GCORR;
REPEATED variable_id / SUBJECT = individual_ID(cluster_ID) TYPE=UN R RCORR;
RUN;
Here is my try:
dta = sqlQuery(odbcConnect("sasodbc", believeNRows=FALSE), "select * from
psdata.bivar")
v...