options ls=75; data one; input junk CP AMBW21D; cards; 1 21.8 942.50 2 21.8 985.00 3 21.8 977.50 4 21.8 978.75 5 21.8 1000.00 6 21.8 956.25 7 21.8 1018.75 8 21.8 1090.00 9 21.8 1001.43 10 21.8 1016.25 11 21.8 1013.75 12 21.8 971.25 13 21.8 930.00 14 21.8 1025.00 15 21.8 1040.00 16 21.8 973.75 17 21.8 962.50 18 21.8 992.50 19 23.5 1031.43 20 23.5 1005.00 21 23.5 997.50 22 23.5 992.50 23 23.5 941.25 24 23.5 1025.00 25 23.5 946.25 26 23.5 985.00 27 23.5 1031.25 28 23.5 1012.50 29 23.5 1003.75 30 23.5 976.25 31 23.5 1006.25 32 23.5 1003.75 33 23.5 1002.50 34 23.5 1050.00 35 23.5 1028.75 36 23.5 1016.25 37 25.2 982.50 38 25.2 1045.00 39 25.2 1040.00 40 25.2 1020.00 41 25.2 965.00 42 25.2 1037.50 43 25.2 960.00 44 25.2 968.75 45 25.2 1106.25 46 25.2 1006.25 47 25.2 1047.50 48 25.2 1011.25 49 25.2 1041.25 50 25.2 1065.00 51 25.2 1050.00 52 25.2 992.50 53 25.2 1001.43 54 25.2 1051.25 55 26.9 1020.00 56 26.9 1091.25 57 26.9 1025.00 58 26.9 1056.25 59 26.9 1008.75 60 26.9 1062.50 61 26.9 1003.75 62 26.9 1037.50 63 26.9 1053.75 64 26.9 1046.25 65 26.9 1127.50 66 26.9 1088.75 67 26.9 1015.00 68 26.9 1063.75 69 26.9 1078.75 70 26.9 1041.25 71 26.9 1037.50 72 26.9 1050.00 ; proc glm; class cp; model AMBW21D = cp; contrast "lof of linear model" cp 0 0 0 0, cp 1 -1 -1 1, cp -1 3 -3 1; *means cp/tukey; *lsmeans cp/cl pdiff adj=tukey; run;