options nodate pagesize=54 linesize=80;

/* Read in the input file   */

data ex0601;
* infile '../../../data/biostat/ex0601.dat';
infile 'C:\biostat\data\ex0601.dat';
input trt cd4;
if trt=1 then treatment='SZZ';
if trt=2 then treatment='SZ';
if trt=3 then treatment='ZZ';
run;

/* Run the ANOVA and compare means based on Tukey's and Bonferroni's methods  */

proc glm;
class treatment;
model cd4 = treatment;
means treatment / bon tukey;
run;

quit;