options ps=54 ls=76;

data caffeine;
* infile '../../data/caffeine.dat';
infile 'C:\caffeine.dat';
input subject dose endurance;
run;

proc plot;
plot endurance*subject=dose;
run;

proc glm;
class subject dose;
model endurance = dose subject;
means dose / bon tukey;
output out=caffout r=caffres p=caffpred;
run;

proc chart;
vbar caffpred;
run;

quit;