/* avgain is the mean average daily gain of all calves in a pen. */
/* Average daily gain is (final weight minus initial weight)/(number of days) on study. */

data weightgain;
input farm pen diet $ avgain;
cards;
1 1 C 2.19
1 2 D 2.44
1 3 B 3.02
1 4 A 2.66
2 5 D 2.36
2 6 C 2.60
2 7 A 2.85
2 8 B 3.37
3 9 B 2.01
3 10 C 1.30
3 11 D 1.57
3 12 A 1.88
4 13 A 2.63
4 14 B 2.99
4 15 D 2.45
4 16 C 2.17
5 17 D 2.43
5 18 C 2.18
5 19 B 2.82
5 20 A 2.55
6 21 A 1.73
6 22 D 1.49
6 23 B 1.96
6 24 C 1.33
;

/* Without Blcok Effect*/
proc glm data=weightgain;
class diet;
model avgain=diet;
lsmeans diet / pdiff;
run;


/* with block effect */
proc glm data=weightgain;
class farm diet;
model avgain=farm diet / p clm;
lsmeans diet / pdiff cl adjust=tukey;
output out=OutData p=predict r=residual;
run;