SAS分析两因素(A,B)的作用
SAS分析得到两因素(A,B)的作用,因素A主效应显著,因素B主效应不显著,但交互作用显著,如何进行下一步分析,得到具体A,B交互作用的具体情况
B1 B2 B3
21.4 19.6 17.6
A1 21.2 18.8 16.6
20.1 16.4 17.5
12 13 13.3
A2 14.2 13.7 14
12.1 12 13.9
12.8 14.2 12
A3 13.8 13.6 14.6
13.7 13.3 14
这是我的SAS程序,运行结果显示因素A主效应显著,因素B主效应不显著,但交互作用显著,然后显示了因素A、B之间显著情况,我现在想知道接下来怎么看AB交互作用具体显著情况
data new;
do a = 1 to 3;
do i = 1 to 3;
do b = 1 to 3;
input y@@;
output;
end;
end;
end;
drop i;
cards;
21.4 19.6 17.6
21.2 18.8 16.6
20.1 16.4 17.5
12.0 13 13.3
14.2 13.7 14
12.1 12 13.9
12.8 14.2 12
13.8 13.6 14.6
13.7 13.3 14.0
;
proc anova;
class a b;
model y = a b a*b;
means a/duncan alpha = 0.05;
means a/duncan alpha = 0.01;
means b/duncan alpha = 0.05;
means b/duncan alpha = 0.01;
run;
返回小木虫查看更多
京公网安备 11010802022153号
采用两两比较就可以知道了。具体程序如下:
data new;
do a = 1 to 3;
do i = 1 to 3;
do b = 1 to 3;
input y@@;
output;
end;
end;
end;
drop i;
cards;
21.4 19.6 17.6
21.2 18.8 16.6
20.1 16.4 17.5
12.0 13 13.3
14.2 13.7 14
12.1 12 13.9
12.8 14.2 12
13.8 13.6 14.6
13.7 13.3 14.0
;
PROC GLM;
CLASS a b;
MODEL y = a b a*b;
LSMEANS a b a*b/tdiff adjust=bon;
/*LSMEANS语句表明按a b a*b分组输出每组的最小二乘均数,tdiff选项表示输出
两两比较的t值和P值,adjust选项指定两两比较方法采用bon,即bonferroni法。
其它方法还有bootstrap和permutation法,nsample= 指定后两种法重复抽样的
次数*/
RUN;
结果中,看两两比较后的p值,就可以知道具体是谁跟谁差异有统计学意义,希望能帮到你,
谢谢,输出结果不是太会看,能帮我解释一下吗
结果怎么看:比如你要看1组和3组两两比较的结果,就看第一行和第三列交叉处,有两行数字,上面的数字是两组比较的t值,下面的为两组比较后t值对应的p值,主要看下面这一行。凡是p值小于0.05,则表明在0.05水平两组差异有统计学意义,否则认为两组差异无统计学意义