sas编程习题与实例应用要点

input a$ b@@; cards; run;

proc ttest data=test2 alpha=0.01; class a; var b;

run;

因为0.6189>0.01,所以方差相等，看第一个equal,0.0024<0.01 拒绝原假设，两个独立样本有显著差异。 四、方差分析

1、生产配方问题。生产某种产品需要两种原材料： A和B；考察的指标为产品的产量Y。设因子A有三个品牌：A1，A2，A3；因子B有四个品牌：B1，B2，B3，B4。对这12中配搭的每一种，安排两次实验，得数据如下表（每次实验的产量）：

B1 B2 B3 26,28.5 B4 27.8，28.5 30.2,29.8 A1 19.3，19.2 24，27.3 A2 21.7,22.6 27.5，30.3 29，28.7 A3 20,20.1 24.2,27.3 24.5，27.1 28.1，27.7 （1）、原材料A的不同品牌对产量是否有显著的差别？哪个品牌更好些（??0.05）？

9

（2）、原材料B的不同品牌对产量是否有显著的影响？ （3）、A和B对产量的影响哪个更大些？ （4）、A和B是否有交互作用？

（5）、使产量达到最大的生产条件是什么？

data test3; do a=1 to 3; do b=1 to 4; do c=1 to 2; input x @@; output; end; end; end; cards; run;

proc glm data=test3; class a b; model x=a b a*b;

means a b a*b/t; /*多重比较设置*/ run;

因为0.0036<0.05 拒绝原假设，A品牌对产量有显著影响。

10

0.001<0.05 拒绝原假设，A品牌对产量有显著影响。 0.8995>0.05 接受原假设，A、B相互没有作用。

在A品牌中分成两组，a2和a1、a3.因为a2均值最大所以最好。 （如何判断看a1a3可分成一组看上面那个表格第五行，只要差值小于1.4596，即可看成一组，两个效果差不多。）

同理，b4最好。

若A、B相互不影响，那么使产量最大的组合可直接选择A、B分

11

别均值最大那个，即a2*b4.

若A、B相互影响，那么看这个表格均值最大那个。

2、考虑合成纤维收缩率（因子A）和总拉伸倍数（因子B）对纤维弹性Y的影响。设收缩率（A）有4个水平：A1，A2，A3，A4；总拉伸倍数（B）也有4个水平B1，B2，B3，B4。在每个组合下重复二次试验，弹性的数据如下表：

A1（=0） A2（=4） A3（=8） A4（=12） B1（=460） B2（=520） B3（=560） B4（=640） 71,73 73,75 76,73 75,73 72,73 76,74 79,77 73,72 75,73 78,77 74,75 70,71 77,75 74,74 74,73 69,69 （1）、收缩率对弹性Y是否有显著的影响？哪个收缩率影响更大些？（??0.05）

（2）、不同拉伸倍数对产量是否有显著的影响？ （3）、收缩率和拉伸倍数是否有交互作用？ （4）、使弹性达到最大的生产条件是什么？

12