options nodate ps=54 ls=76; data ex0603; * infile '../../../data/biostat/ex0603.dat'; infile 'C:\biostat\data\theoph.dat'; input subject intagnt thcl; length fintagnt $8.; if intagnt=1 then fintagnt='placebo'; if intagnt=2 then fintagnt='pepcid'; if intagnt=3 then fintagnt='tagamet'; run; /* Get mean theophylline clearance by subject */ proc means mean; class subject; var thcl; run; /* Get mean theophylline clearance by interacting agent */ proc means mean; class fintagnt; var thcl; run; /* Run ANOVA for RBD and use Tukey's method to compare agent effects */ proc glm; class subject fintagnt; model thcl = fintagnt subject; means fintagnt / tukey; run; quit;