Hi, I ran into some problems with adding a continous covariate in repeated measures using proc mixed, and there is not much information on the internet.
The situation is like measuring weight gain over time and I want to add initial weights as covariate, then compare weight gain between time points. Analysis ran well without covariate, however, when I added in the covariate, lsmeans of weight gain at different time points were not estimatable.
here are the codes I used
data survey;
input id gain time blk wt;
cards;
1 6.00 1 2 90
1 4.80 1 1 83
2 5.80 1 2 74
2 4.40 1 1 87
1 5.25 2 2 90
1 4.00 2 1 83
2 5.00 2 2 74
2 4.70 2 1 87
1 5.25 3 2 90
1 4.50 3 1 83
2 4.25 3 2 74
2 5.10 3 1 87
;
proc mixed data=survey;
class id time blk;
model gain = time blk id(blk) wt;
repeated time / type = ar(1) subject = id(blk);
lsmeans time / pdiff adjust = bon;
run;
Can anyone help? Thanks.
First, recognize that blk is a random effect, and that you have to have a unique record for each subject.
Try:
proc mixed data=survey;
class id time blk;
model gain = time wt;
random intercept id/subject=blk;
repeated time / type = ar(1) subject = id(blk);
lsmeans time / pdiff adjust = bon;
run;
This gives an AR+RE model for the response variable.
HOWEVER, modeling a change variable with a covariate that is one of the elements of the change is looking for trouble. See the following link
http://biostat.mc.vanderbilt.edu/wiki/Main/MeasureChange2
and especially the section:
In a two-group parallel design, analysis of change is not recommended at all. The response variable should be the final measurement and the baseline measurement should be adjusted for as a covariate using analysis of covariance, with treatment assigned as one of the other variables. Besides the issues listed above, change scores are affected by regression to the mean. The slope of the baseline value may not be 1.0.
What this means is that you are adjusting for the initial body weight twice in your proposed analysis.
Steve Denham
Don't miss out on SAS Innovate - Register now for the FREE Livestream!
Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.