Chapter 5, Section 5.4
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;