Have you considered using PROC GLIMMIX, as in the following:
proc glimmix data=MAST_long_VAS;
class PPID Time;
model BC_stress = Time Time*Time / solution;
random _residual_ / subject=PPID type=ar(1);
run;
In large part, this is adapted from the example - Example 22.10 Correlation and Covariance Matrices in Using the Output Delivery System documentation.
SteveDenham
... View more
