Chapter 5, Section 5.4
Analysis of Response Profiles of Blood Lead Levels
data tlc;
infile 'tlc.dat';
input id group $ lead0 lead1 lead4 lead6;
y=lead0; time=0; output;
y=lead1; time=1; output;
y=lead4; time=4; output;
y=lead6; time=6; output;
drop lead0 lead1 lead4 lead6;
title1 Analysis of Response Profiles of data on Blood Lead Levels;
title2 Treatment of Lead Exposed Children (TLC) Trial;
**************************************************************************;
* Uncomment to change reference levels for group and time *;
**************************************************************************;
*proc sort;
*by group descending time;
proc mixed noclprint=10 order=data;
class id group time;
model y = group time group*time / s chisq;
repeated time / type=un subject=id r;
run;