topic Re: How to run proc GLIMMIX for repeated measures in Statistical Procedures

How to run proc GLIMMIX for repeated measures

Therain — Sat, 18 Mar 2017 18:10:42 GMT

Hello,

I use SAS 9.4 and I have a data set like this:

ID pcs year disease

1 45 1 1

1 47 2 1

1 50 3 1

2 42 1 0

2 43 2 0

2 46 3 0

3 40 1 0

3 40 2 0

3 49 3 0

...

I want to see how the presence or absence of disease affect the PCS scores rate of change from year 1 to year 3 for each subject. I run the following codes:

proc glimmix data=have;
class disease (ref='0');
model pcs = disease year disease*year / cl;
run;

Then I got the results but I do not think it is correct, since the estimate for intercept is around 12, but the average pcs at year 1 is over 40 instead of 12. Does anyone can help me figure out the problem?

Thank you!

Re: How to run proc GLIMMIX for repeated measures

Ksharp — Sun, 19 Mar 2017 02:30:44 GMT

In GLIMMIX , use RANDOM to get repeated measures ( R- side random effect).

proc glimmix data=have;
class disease (ref='0')  id year ;
model pcs = disease year disease*year / cl;
random year/subject=id type=ar(1) residual;
run;

Re: How to run proc GLIMMIX for repeated measures

Therain — Sun, 19 Mar 2017 17:41:14 GMT

Thank you for the reply! When I run your codes I met another problem.

I can get parameter estimates for intercept, disease, year, and disease*year. But when I add some other covariates into the model, i found that the estimates for intercept and disease changed but the estimate for disease*year remained the same. I am not quite sure if this is correct? (I think the estimate for disease*year should change and the estimates for intercept and disease should remain the same).