第六讲 方差分析(三)

65 62 56 61 59 60
71 75 85 67 78 89
区分是否是系统分组设计的关键是看因素之间的地 位是否平等。应把一级因素不同水平中的次级因素同一 水平看作是不同水平。分析侧重于一级因素。
变异来源：
❖ A因素（一级因素）：A因素各水平间的变； ❖ B(A)因素：A因素各水平内B因素各水平间的变异； ❖ 随机误差：B因素各水平内重复观测值间的变异。
产量
A1 B1 B2 B3 85 72 70 89 70 67
A2 B4 B5 B6 82 91 85 84 88 83
A3 B7 B8 B9 65 59 60 61 62 56
(三)拉丁方设计的ANOVA
1.拉丁方设计的概念：
从横向和纵向两个方向对试验环境条件进行局部控制， 试验精确性比随机区组设计高。由于重复数与处理数必须相 等，缺乏灵活性。
5×5 阶标准拉丁方设计 ABCDE BAECD CDAEB DEBAC ECDBA
n×n 阶方阵，各处理在每一横行与纵列都只出现一次。
（1）列随机（按14325）
设计步骤：
14325
1 ADCB E
2 BCEAD
3 CEADB
4 DAB EC
5 EBDCA
（2）行随机（按53124）
14325 5 EBDCA 3 CEADB 1 ADCB E 2 BCEAD 4 DAB EC
根据局部控制的原理，将试验的所有供试单元先按重复 划分成非处理条件相对一致的若干单元组（block）。
再同一区组内各处理单元的排列顺序随机而定，这样的 区组即为随机区组（random block）。
例6-6：某作物品种比较试验，有8个品种（含对照），
设3次重复。考虑到不同试验田土壤肥沃程度存在差异，使 用随机区组设计。
叶片间存在极显著差异。
课堂SPSS演示：
例题6-14 随机选取3株植物，在每一株内随机选取两片叶子，
用取样器从每一片叶子上选取同样面积的两个样品，称取湿重 （g），结果见表。试分析不同植株和同一植株上的不同叶片间湿 重是否有显著差异。
植株A A1 A2 A3
叶片B B11 B12 B21 B22 B31 B32
0.44
例题6-10：玉米品种（A）与施肥量（B）两因素试验，因素A
有A1、A2、A3、A4四个水平，因素B有B1、B2两个水平，共有8个
水平组合即处理，3次重复（r=3），小区计产面积20 m2，田间排列
和产量（kg/20 m2）见表，采用随机区组设计，试对试验结果进行方
差分析。
处理
区组
I
II
III
（3）处理随机（按4=A， 1=B，2=C，3=D，5=E）
51324 25431 43215 12543 34152
例题6-16：有一冬小麦施氮肥时期试验，5个处理
为：A不施氮肥（CK）；B播种期（10月29日）施氮；C越冬期 （12月13日）施氮；D拔节期（3月17日）施氮；E抽穗期（5月1 日）施氮。采用5*5拉丁方设计，小区计产面积32m2，其田间排 列和产量（k/32m2）见表。对试验结果进行方差分析。
例题 空气污染状况调查
因素A （地区）
因素B （季节）
重复
广州
春季 冬季
123123
北京 春季 冬季
123123
例题 不同催化剂在不同温度下对某化合物转化率（%）影响
实
催化剂A1
验
批
温度℃
次 70 80 90
催化剂A2 温度℃
55 65 75
催化剂A3 温度℃
90 95 100
1 82 91 85 2 84 88 83
12.8
14.4
14.4
14.7
14.5
23.1
23.4
28.1
28.8
SPSS分析演示： AnalyzeGeneral Linear ModelUnivariate…
第一步：各项平方和与自由度的计算（略） 第二步：建立方差分析表（略）
植株p < 0.05，差异显著；叶片（植株），p < 0.01，差异极显著。 不同植株的不同叶片湿重存在显著差异；同一植株上不同
SPSS分析演示： AnalyzeGeneral Linear ModelUnivariate… 注意：“model”选择“main effect”
第一步：各项平方和与自由度的计算（略） 第二步：建立方差分析表（略）
处理p < 0.01，差异极显著；行、列区组，p >0.05，差异不显著。
第三步：多重比较
混合模型（Mixed model）：既有固定因素，又有随机因素。
例题6-10：玉米品种（A）与施肥量（B）两因素试验，因素A
有A1、A2、A3、A4四个水平，因素B有B1、B2两个水平，共有wenku.baidu.com个
水平组合即处理，3次重复（r=3），小区计产面积20 m2，田间排列
和产量（kg/20 m2）见表，采用随机区组设计，试对试验结果进行方
C 10.1 A 7.0 E 7.6 D 10.5 B 8.9
A 7.9 D 10.0 C 9.7 B 9.6 E 8.9
B 9.8 E 7.0 D 10.0 C 9.8 A 8.6
E 7.1 C 9.7 B 9.3 A 6.6 D 10.6
D 9.6 B 9.1 A 6.8 E 7.9 C 10.1
第三步：多重比较
Tukey法多重比较
处理
mean±SE
A1
0.22 ±0.01 Cc
A2
0.34 ±0.01 Bb
A3
0.39 ±0.01 Aa
第四步：写出统计结论
不同提取时间对山楂色素的提取结果有极显著影响，其中A3处 理（4 h）的结果最优。乙醇浓度，提取时间与乙醇浓度的互作，及 重复效应对提取结果的影响不显著。
B因素的不同水平分别嵌套在A因素的不同水平下，这种 设计称为嵌套设计。也称作系统分组设计、巢式设计、套设 计等。农业试验中常见的试验设计方法。
（1）观察对象本身具有分组再分组的各种分组因素，处 理（即最终的试验条件）是各因素各水平的全面组合，且因 素之间在专业上有主次之分。
例题：随机取若干地块，每地块随机取若干个样点，
第六讲 方差分析（三）
内容框架
(一)随机区组设计的ANOVA
随机区组设计的概念 ANOVA基本步骤
(二)嵌套设计的ANOVA
嵌套设计的概念 ANOVA基本步骤
(三)拉丁方设计的ANOVA
拉丁方设计的概念 ANOVA基本步骤
(四)ANOVA的基本假定及数据转换
(一)随机区组设计的ANOVA
1.随机区组设计的概念：
忽略，则需采取随机区组设计。
害虫数量
少
多
产品效果试验
例6-8：在蛋糕加工工艺研究中，欲考察不同食品添加剂（A1、
A2、A3、A4）对不同配方（B1、B2、B3）制成的蛋糕质量的影响。 因试验用烤箱容量不大，不能一次性将全部试验蛋糕烘烤完，只能分 次烘烤，故选用随机区组设计进行试验。3次烘烤分为区组I、II、III。
例题6-14 随机选取3株植物，在每一株内随机选取两片叶子，
用取样器从每一片叶子上选取同样面积的两个样品，称取湿重 （g），结果见表。试分析不同植株和同一植株上的不同叶片间湿 重是否有显著差异。
植株A A1 A2 A3
叶片B B11 B12 B21 B22 B31 B32
湿重
12.1
12.1
12.8
SPSS分析演示： AnalyzeGeneral Linear ModelUnivariate… 注意：“model”选择“interaction” 例题6-9 输入：
第一步：各项平方和与自由度的计算（略） 第二步：建立方差分析表（略）
A因素（提取时间），p < 0.01，差异极显著；B因素（乙醇 浓度）、区组和A*B互作，p > 0.05，差异不显著。
差分析。
处理
区组
I
II
III
A1
B1
12.0
13.0
13.0
田
B2
11.0
10.0
13.0
A2
B1
19.0
16.0
12.0
间
B2
20.0
19.0
17.0
试
A3
B1
19.0
18.0
16.0
B2
10.0
8.0
7.0
验
A4
B1
17.0
16.0
15.0
B2
11.0
9.0
8.0
(二)嵌套设计的ANOVA
嵌套设计的概念：
A1
B1
12.0
13.0
13.0
田
B2
11.0
10.0
13.0
A2
B1
19.0
16.0
12.0
间
B2
20.0
19.0
17.0
试
A3
B1
19.0
18.0
16.0
B2
10.0
8.0
7.0
验
A4
B1
17.0
16.0
15.0
B2
11.0
9.0
8.0
小知识：
固定模型（Fiexed model）：表示你打算比较的就是你现在 选中的这几个处理。例如，比较3种药物的疗效的差别；这3种 药不是从很多种药中抽样出来的，不想推广到其他的药物，结 论仅限于这3种药。 随机模型（Random model）：表示你打算比较的不仅是你设 计中的几种处理，而是想通过这几种处理的比较，推广到他们 所能代表的总体。例如，想知道名牌大学的就业率是否高于普 通大学，你选择了清华、浙大、浙江工商大学、浙江农林大学4 所学校进行比较。
每一样点的土样又作数次分析所获得的资料。
因素A （地块）
因素B （土样）
样点
地块1
上层
下层
123123
地块2
上层
下层
123123
例题： 高校大学生新生入学成绩调查
因素A （大学）
因素B （专业）
学生
某重点大学
农学 机械
123123
某非重点大学 农学 机械
123123
（2）受试对象本身并非具有分组再分组的各种分组因素， 处理（即最终的试验条件）不是各因素各水平的全面组 合，而是各因素按其隶属关系系统分组，且因素之间在 专业上有主次之分。
区组
蛋糕在烘箱中的位置
I 6 12 3 5 1 7 11 2 8 4 10 9 A2B A4B A1B A2B A1B A3B A4B A1B A3B A2B A4B A3B
II 8 1 4 9 10 6 3 12 2 5 7 11 A3B A1B A2B A3B A4B A2B A1B A4B A1B A2B A3B A4B
Tukey法多重比较
处理
mean±SE
A不施氮肥
7.38 ±0.84 Bb
B播种期施氮
9.34 ±0.36 Aa
C越冬期施氮
9.88 ±0.20 Aa
D拔节期施氮
10.14 ±0.41 Aa
E抽穗期施氮
7.70 ±0.76 Bb
第四步：写出统计结论
不同时期施氮处理存在极显著差异， 其中拔节期、播种期和越冬期施氮处理的 产量极显著高于抽穗期施氮和CK。
SPSS分析演示：
例题6-16：有一冬小麦施氮肥时期试验，5个处理
为：A不施氮肥（CK）；B播种期（10月29日）施氮；C越冬期 （12月13日）施氮；D拔节期（3月17日）施氮；E抽穗期（5月1 日）施氮。采用5*5拉丁方设计，小区计产面积32m2，其田间排 列和产量（k/32m2）见表。对试验结果进行方差分析。
湿重
12.1
12.1
12.8
12.8
14.4
14.4
14.7
14.5
23.1
23.4
28.1
28.8
课堂SPSS操作：
例题6-15 比较不同公鱼的产鱼效应。随机选取3条公鱼，
每条公鱼随机与3条同品种的母鱼交配受精后，所生小鱼各分 两池养殖，长大为成鱼后检测各池产鱼量，试作方差分析。
公鱼（A） 母鱼（B）
重复（区组）
B1
B2
B3
A1
0.22
0.18
0.25
工
I
A2
0.33
A3
0.39
0.35 0.42
0.36 0.35
艺
A1
0.18
0.22
0.22
优
II
A2
0.32
0.30
0.37
A3
0.37
0.40
0.38
化
A1
0.24
0.20
0.27
III
A2
0.35
0.32
0.38
A3
0.41
0.37
0.44
贫瘠 区组Ⅰ 区组Ⅱ
肥沃 区组 Ⅲ
E
C
F
A
对 照
D
BG
G
B
对 照
D
F
C
E
A
D
A
E
G
B
F
对 照
C
育种试验&推广试验
例6-7：比较4种杀虫剂（A1,A2,A3,A4）对棉花田的害虫杀虫效
果，特选了20块田，设重复5次。不同试验田害虫多少、植物长势、
土壤肥沃程度等因素都可能对杀虫效果产生影响。如果此种影响不可
III 10 7 2 11 4 8 5 9 1 12 6 3 A4B A3B A1B A4B A2B A3B A2B A3B A1B A4B A2B A1B
产品配方&研发
例题6-9：为研究山楂色素的最佳提取条件，选取提取时间(A)
和乙醇浓度(B)为试验因素。提取时间（A/h）取2、3、4三个水平， 乙醇浓度（B/%）取55、75、95三个水平，每个水平组合重复3次， 试验结果如表所示。现以重复为区组，对试验结果进行统计分析。
课堂SPSS操作：
重复（区组）
B1
B2
B3
A1
0.22
0.18
0.25
I
A2
0.33
0.35
0.36
A3
0.39
0.42
0.35
A1
0.18
0.22
0.22
II
A2
0.32
0.30
0.37
A3
0.37
0.40
0.38
A1
0.24
0.20
0.27
III
A2
0.35
0.32
0.38
A3
0.41
0.37