SAS学习第7章：方差分析之两因素系统分组次级样本

5.两因素系统次级样本含量相等资料方差分析

例：为测定3种不同的来源的鱼粉的蛋白质消化率，在不含蛋白质的饲料中按一定比例分别加入不同鱼粉A1，A2，A3制成饲料，各喂给3头实验动物B，分别从每头动物排泄物中各取两份样品做化学分析，结果如下图，试分析不同来源鱼粉的蛋白质消化率是否有差异。

一级因素 鱼粉A	二级因素 个体B	测定结果
A1	B11	82.5	82.4
	B12	87.1	86.5
	B13	84	83.9
A2	B21	86.6	85.8
	B22	86.2	85.7
	B23	87	87.6
A3	B31	82	81.5
	B32	80	80.5
	B33	79.5	80.3

data fish;
do a=1 to 3;
do b=1 to 3;
do c=1 to 2;
input d@@;
output;
end;
end;
end;
drop c;
cards;
82.5 82.4 87.1 86.5 84 83.9
86.6 85.8 86.2 85.7 87 87.6
82 81.5 80 80.5 79.5 80.3
;
proc anova;
class a b;
model d=a b (a);
means a/snk;
run;

 

结论：3种鱼粉消化率差异显著，A2最高，A3最低。

 6.两因素系统次级样本含量不相等资料方差分析

例：某品种3头公猪和8头母猪所生仔猪的35日龄断奶重(kg)资料如下表，分析不同公猪与母猪对仔猪断奶重的影响是否有差异？

公猪A	配种母猪B	仔猪断奶重（kg）
A1	B11	10.5	8.3	8.8	9.8	10	9.5	8.8	9.3	7.3
	B12	7	7.8	8.3	9	8	7.5	9.3
A2	B21	12	11.3	12	10	11	11.5	11	11.3
	B22	9.5	9.8	10	11.8	9.5	10.5	8.3
	B23	8	8	7.8	10.3	7	8.8	7.3	7.8	9.5
A3	B31	7.5	6.5	6.8	6.3	8.3	6.8	8	8.8
	B32	9.5	10.5	10.8	9.5	7.8	10.5	10.8
	B33	11.3	10.5	10.8	9.5	7.3	10	11.8	11

data pig;
input a b c d@@;
cards;
1 1 1 10.5 1 1 2 8.3 1 1 3 8.8 1 1  4 9.8 1 1 5 10 1 1 6 9.5 1 1 7 8.8 1 1 8 9.3 1 1 9 7.3
1 2 1 7 1 2 2 7.8 1 2 3 8.3 1 2 4 9 1 2 5 8  1 2 6 7.5 1 2 7 9.3
2 1 1 12 2 1 2 11.3 2 1 3 12 2 1 4 10 2 1 5 11 2 1 6 11.5 2 1 7 11 2 1 8 11.3
2 2 1 9.5 2 2 2 9.8 2 2 3 10 2 2 4 11.8 2 2 5 9.5 2 2 6 10.5 2 2 7 8.3
2 3 1 8 2 3 2 8 2 3 3 7.8 2 3 4 10.3 2 3 5 7 2 3 6 8.8 2 3 7 7.3 2 3 8 7.8 2 3 9 9.5
3 1 1 7.5 3 1 2 6.5 3 1 3 6.8 3 1 4 6.3 3 1 5 8.3 3 1 6 6.8 3 1 7 8 3 1 8 8.8
3 2 1 9.5 3 2 2 10.5 3 2 3 10.8 3 2 4 9.5 3 2 5 7.8 3 2 6 10.5 3 2 7 10.8
3 3 1 11.3 3 3 2 10.5 3 3 3 10.8 3 3 4 9.5 3 3 5 7.3 3 3 6 10 3 3 7 11.8 3 3 8 11
;
proc glm;
class a b c;
model d=a b (a);
means a/snk;
run;

 答：公猪母猪的品种对仔猪断奶重影响均有显著差异

注：侧重分析谁，谁在model后要再加（）予以标注，如侧重分析a：model c=a b (a);

means 后也可只写 a ：means a/snk;