Hello, I need to obtain the residuals and then use them as the outcome in a separate model. In the past, I was able to extract the residuals using proc reg, as shown below: proc reg data=abstract.pedometer_speedquint2; model stepsperday=stepmin_avg_valid ; OUTPUT out=b r=res; run; proc print data=abstract.pedometer_speedquint2; var r; /*residuals per participant*/ run; However, I have not been able to figure out how to do the same using proc mixed. This is the code I am using for the analysis from which I want to extract the residuals: proc mixed data=trial4; class idno time; model CASI= time WMFA WMH ICV TOTAL_GM_VOL TOTAL_WM_VOL AGE RACE GENDER/s chisq; repeated time/type=un subject=idno r rcorr; run; Thanks!
... View more