ACCAspace_sitemap
PPclass_sitemap
sitemap_google
sitemap_baidu
CFA Forums
返回列表 发帖

!!挑战高手!!来看看这个程序!!

谁能看懂,求解释!!

AIDS Clinical Trial Group (ACTG) 193A Study

 

 

Linear Mixed Effects Model (Random Intercept and Slopes)

 

 

data cd4;
     infile 'cd4.dat';
     input id group age sex week logcd4;

****************************************************************;

*   Create new variable combining groups 1, 2, and 3   *;

****************************************************************;

trt=0;
if (group=4) then trt=1;

week_16=max(week - 16, 0);
w16=week_16 - week;


title1 Mixed Effects Model for log CD4;
title2 AIDS Clinical Trial Group (ACTG) 193A Study;

proc mixed method=reml noclprint=10 covtest;
     class id;
     model logcd4 = week week_16 trt*week trt*week_16 / s chisq;
     random intercept week week_16  / subject=id type=un g gcorr;
run;

<Selected Output>

 

title1 Mixed Effects Model for log CD4 adjusting for Age and Gender;
title2 AIDS Clinical Trial Group (ACTG) 193A Study;

ods output solutionr=bluptable;

proc mixed method=reml noclprint=10 covtest;
     class id;
     model logcd4 = age sex week week_16 trt*w16 / s chisq outpred=yhat;
     random intercept week week_16  / subject=id type=un g gcorr solution;
run;


<Selected Output>

 

 

 

Linear Mixed Effects Model (Random Intercept and Slopes): Predicted Means

 

 

******************************************;

*   Print first 60 observations only   *;

******************************************;

proc print data=yhat (obs=60);
     var id trt week logcd4 pred;
run;


<Selected Output>

 

 

 

 

Linear Mixed Effects Model (Random Intercept and Slopes): Empirical BLUPs

 

 

******************************************;

*   Print first 60 observations only   *;

******************************************;

proc print data=bluptable (obs=60);
run;

<Selected Output>

以下连接为CD4.DAT

http://www.biostat.harvard.edu/~fitzmaur/ala/cd4.txt

返回列表