SAS方差分析

《统计计算》第二次上机练习

10 统计1 何斌

一、由SAS函数产生正态分布N（3，9）的观察值500个，并画直方图。 SAS程序如下所示 data norm_rand; do i=1 to 500; k=rannor(0); r=3*k+3; output; end;

proc print data=norm_rand; proc gchart data=norm_rand; vbar r; run;

FREQUENCY1009080706050403020100-6.0-4.5-3.0-1.50.01.53.04.56.07.59.010.5r MIDPOINT

二、考察的因子有水稻品种A和施肥量B；考察的指标为水稻的产量Y。设因子A有三个水平：A1（窄叶青），A2（珍珠矮）和A3（江二矮）；因子B有四个水平：B1（无肥），B2（低肥），B3（中肥）和B4（高肥）。对这12种搭配的每一种，在两块试验田上做试验。每一块试验田分为12块面积相同小田，随机地安排12种搭配条件进行试验。得数据如下表（小田的产量）：

因子B 因子A A1（窄叶青） A2（珍珠矮） A3（江二矮） B1（无肥） 19.3, 19.2 21.7, 22.6 20.0, 20.1 B2（低肥） 24.0, 27.3 27.5, 30.3 24.2, 27.3 B3（中肥） 26.0, 28.5 29.0, 28.7 24.5, 27.1 B4（高肥） 27.8, 28.5 30.2, 29.8 28.1, 27.7 1、 进行单因子方差分析，要求完成：

（1） 不同的施肥量对产量是否有明显的影响；

（2） 若不同的施肥量对产量有明显影响，请进行多重比较(SNK)，指出哪种施肥量有

显著的不同；

2、 进行两因子方差分析，要求完成：

（1） 稻种和施肥量之间有无交互作用；

（2） 稻种和施肥量对产量是否有显著性影响； （3） 使产量达最高的生产条件是什么；

答：

1. 单因子方差分析 SAS程序如下所示 data rice;

do a=1 to 3; do b=1 to 4; do i=1 to 2;

input production @@; output; end; end; end; cards;

19.3 19.2 24.0 27.3 26.0 28.5 27.8 28.5 21.7 22.6 27.5 30.3 29.0 28.7 30.2 29.8 20.0 20.1 24.2 27.3 24.5 27.1 28.1 27.7 ;

proc print;run; proc anova; class b;

model production=b;

means b/snk alpha=0.05; run;

得到结果如下所示

从方差分析表的输出结果来看，由于P值<0.0001，故可以认为模型显著有效，由于施肥程度的P值<0.0001，可以认为在显著性水平为0.05的条件下施肥量对产量有显著性的影响。

由SNK法得到的多重比较的结果来看，B2（低肥），B3（中肥）和B4（高肥）这三种施肥量对产量的影响没有差异，而B1（无肥）与前三种施肥量对产量的影响有显著不同。

2. 两因子方差分析 SAS程序如下所示

proc anova data=rice; class a b;

model production=a|b; run;

得到结果如下所示

从方差分析表的输出结果来看，由于P值<0.0001，故可以认为模型显著有效，由于A*B的交互作用的P值=0.8995>0.05，所以在显著性水平为0.05的条件下可以认为稻种和施肥量之间的交互作用对水稻的产量没有显著影响。

重新建立SAS程序如下所示 proc anova data=rice; class a b;

model production=a b; means a b; run;

得到结果如下所示

从方差分析表的输出结果来看，由于P值<0.0001，故可以认为模型显著有效，由于去除稻种与施肥量的交互作用后在显著性水平为0.05的条件下，稻种的检验P值和施肥量的检验P值均小于0.05，故可以认为稻种和施肥量对产量均有显著性影响。

因素A的水平为2时，其所对应的产量的均值最大，对于因素B，则水平4所对应的产量的均值最大，因此可以选出使产量达最高的生产条件是A=2 B=4，即水稻品种为珍珠矮和施肥量为高肥时能使产量达到最高。

三、有四种药物，每种药物有三个水平：A药（A1=2%，A2=4%，A3=5%），B药（B1=0%，B2=1%，B3=2%），C药（C1=0%，C2=1%，C3=3%）及D药（D1=0%，D2=1%，D3=3%）。问四种药物对淡色库蚊的50%击倒时间有无差别？并进行多重比较，得出结论。采用正交表L9（34），相同试验条件下均做4次重复，其数据如下：

本文来源：https://www.bwwdw.com/article/0ggh.html