DATA EXAMPLE4_12;
INPUT g1 g2 m;
CARDS;
1	1	0.671 
1	1	0.947 
1	1	1.241 
1	1	0.953 
1	1	0.711 
1	1	0.977 
1	2	1.191 
1	2	1.613 
1	2	1.585 
1	2	2.112 
1	2	1.956 
1	2	2.181 
2	1	0.270 
2	1	0.423 
2	1	0.507 
2	1	0.399 
2	1	0.454 
2	1	0.535 
2	2	0.705 
2	2	0.589 
2	2	0.751 
2	2	0.804 
2	2	0.725 
2	2	0.652 
;
PROC GLM;
CLASS g1 g2;                    /*CLASS语句指明有两个分组因素g1和g2*/
MODEL m=g1 g2 g1*g2;           
/*MODEL语句指明要分析g1、g2及g1和g2的乘积即交互项的效应*/
LSMEANS g1*g2/slice=g1;       
/*LSMEANS语句输出最小二乘均数，g1*g2交叉组成四组，故输出四组的均数。SLICE语句表示控制g1比较g2的差异，即分别给出g1=1和g1=2时g2的两组比较结果*/
LSMEANS g1*g2/slice=g2;
/* SLICE语句表示控制g2比较g1的差异，即分别给出g2=1和g2=2时g1的两组比较结果*/
LSMEANS g1*g2/tdiff adjust=bon;
/* 该语句表示对g1*g2交叉形成的四组采用bonferroni法进行两两比较*/
RUN;

