The following SAS code is from http://www.ats.ucla.edu/stat/sas/faq/anovmix1.htm
PROC MIXED DATA=long;
CLASS sub group trial;
MODEL dv = group trial group*trial;
REPEATED trial / SUBJECT=sub TYPE=CS;
run;
I was just wondering how you would use PROC GLIMMIX to write this model. Would it be:
PROC GLIMMIX DATA =long;
class sub group trial;
MODEL dv = group;
RANDOM trial/ SUBJECT = sub TYPE = CS;
RANDOM trial sub trial*sub trial*group;
Also out of curiosity, would the model be as follows:
DV = GROUP + TRIAL + SUB + TRIAL*SUB + TRIAL*GROUP where the terms after GROUP are all random?
No, your glimmix code is fitting a different model. Just take your MIXED code and change your REPEATED statement to:
RANDOM TRIAL / SUBJECT=SUB RESIDUAL TYPE=CS;
If you do not have missing values for trial (i.e., same trials for each subject), you can use:
RANDOM _RESIDUAL_ / SUBJECT=SUB TYPE=CS;
What does the RESIDUAL option do?
It sets up the R side matrix (repeated).
Steve Denham
As Steve says, the statements in GLIMMIX
random _residual_ / sub=___ type=____;
or
random ___ / sub=____ type=____ RESIDUAL;
do the same thing as the REPEATED statement in MIXED. Since you had a model in MIXED that you wanted to fit in GLIMMIX, I showed you the appropriate code.
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.