options ps=60 ls=80 pageno=1 formdlim='_';

data cholest;
input drug $ plant $ conc;
cards;
d1 p1 102
d1 p1 104
d1 p2 103
d1 p2 104
d2 p1 108
d2 p1 110
d2 p2 109
d2 p2 108
d3 p1 104
d3 p1 106
d3 p2 105
d3 p2 107
;

proc print;
title 'cholesterol data (subsamples)';
title2 'ANOVA using proc glm';
run;

proc glm;
class drug plant;
model conc = drug plant(drug);
random plant(drug);
test h=drug e=plant(drug);
lsmeans drug/pdiff cl adjust=tukey e=plant(drug);
run;