/* Section 1.3.4 (p. 12) Treatment x age combination example */

data interact;
input obs age old new diff;
cards;
1 1 57 72 -15
2 1 53 27 26
3 1 28 26 2
4 1 60 71 -11
5 1 40 60 -20
6 1 48 45 3
7 2 70 52 18
8 2 85 26 59
9 2 50 46 4
10 2 61 52 9
11 2 83 53 30
12 2 51 58 -7
13 3 55 83 -28
14 3 36 65 -29
15 3 31 40 -9
16 3 28 66 -38
17 3 41 50 -9
18 3 32 44 -12
19 4 18 40 -22
20 4 39 55 -16
21 4 53 70 -17
22 4 44 55 -11
23 4 63 61 2
24 4 80 60 20
25 5 76 60 16
26 5 67 37 30
27 5 75 45 30
28 5 78 58 20
29 5 67 54 13
30 5 80 69 11
;
proc glm;
class age;
model diff=age;
lsmeans age/cl alpha=.010206218;

/* illustration of the product inequality method */

run;