12 第十二章 单因素试验的统计分析.ppt

第1页：第十二章单因素试验的统计分析封面

第十二章单因素试验的统计分析

• 第一节对比法和间比法试验的统计分析
• 第二节完全随机和随机区组试验的统计分析
• 第三节拉丁方试验的统计分析
• 第四节试验处理的合并比较

第2页：第一节对比法和间比法试验的统计分析

第一节对比法和间比法试验的统计分析

一、对比法试验结果的统计分析

二、间比法试验结果的统计分析

第3页：对比法试验结果的统计分析

一、对比法试验结果的统计分析

百分比法：设对照(CK)的产量(或其它性状)为100，然后将各处理产量和对照相比较，求出其百分数。但因处理作顺序排列不能估计无偏的试验误差，难以进行假设测验和统计推断

第4页：例12.1 有A、B、C、D、E、F 6个

例12.1 有A、B、C、D、E、F 6个

玉米品种的比较试验，设标准品种CK，采用3次重复的对比设计，田间排列在表12.1第1列基础上作阶梯式更替，此处图形从略。小区计产面积40m2，所得产量结果列于表12.1，试作分析

表12.1首先将各品种在各重复中的小区产量相加，得面积上的产量总和。然后，将各个除以重复次数得各小区平均产量。再计算各品种产量对邻近CK产量的百分数

第5页：例12.1 有A、B、C、D、E、F （续）

例12.1 有A、B、C、D、E、F （续）

例如或其余品种皆类推

第6页：表12.1 玉米品比试验

表12.1 玉米品比试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

对比法)的产量结果分析总和平均品种名称各重复小区产量(kg)对邻近CK的%ⅠⅡⅢCK

第7页：例12.1 有A、B、C、D、E、F （续）

例12.1 有A、B、C、D、E、F （续）

计算各品种对邻近CK的百分数是为得到一个比较精确的、表示各品种相对生产力的指标。相对生产力＞100%的品种，其相对生产力愈高，就愈可能显著地优于对照品种。但是决不能认为相对生产力＞100%，所有品种都是显著地优于对照的

由于误差的存在，一般田间试验很难察觉处理间差异在5%以下的显著性。对于对比法(以及后面的间比法)的试验结果，要判断某品种的生产力确优于对照，其相对生产力一般至少应超过对照10%以上

第8页：相对生产力仅超过对照5%左右的品种

相对生产力仅超过对照5%左右的品种

相对生产力仅超过对照5%左右的品种，宜继续试验，再作结论。当然，由于不同试验的误差大小不同，上述标准仅具有参考性质

在本例，B品种产量最高，超过对照19.3%；C品种占第二位，超过对照11.7%；大体上可以认为它们确是优于对照。D品种占第三位，仅超过对照6.7%；再查看各重复的产量，有两个重复(Ⅰ和Ⅲ)D超过CK，一个重复(Ⅱ)D低于CK；因而显然不能作出D品种确优于对照的结论

第9页：作物产量习惯于用每亩产量表示

作物产量习惯于用每亩产量表示

作物产量习惯于用每亩产量表示。其折算方法，先算得对照区的总产量。然后将对照区总产量乘以化对照区总产量为亩产量的改算系数cf，得到对照的亩产量

(12·1)中的A是小区计产面积，以m2为单位；n是小区数目。最后可用各品种的相对生产力乘对照的亩产量，即得各品种的亩产量。(12·1)

第10页：例12.1 有A、B、C、D、E、F （续）

例12.1 有A、B、C、D、E、F （续）

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

如本例，由表12.1可算得对照区总产量=109.9+91.8+91.2+98.0=390.0(kg)，cf= (12×40)=1.3889，所以

对照种亩产量= =541.7(kg)

A品种亩产量= =532.5(kg)

……，依此类推

本例题的田间排列方法也可以按第二章第五节所提排列，即A，CK，B，C，CK，D，E，CK，F，这样可以减少一个对照小区，分析方法相同

第11页：间比法试验结果的统计分析

二、间比法试验结果的统计分析

1：计算前后两个对照产量的平均数

2：计算各品系产量相对应产量的百分数，即得各品系的相对生产力

间比法设计中，采用推广良种作为对照计算肥力指数调整供试家系产量，所以在参试家系数目较多时一般常用两个或两个以上的对照品种

第12页：例12.2 有12个小麦新品系鉴定试验

例12.2 有12个小麦新品系鉴定试验

另加一推广品种CK，采用5次重复间比法设计，田间排列在表12.2第1列基础上按阶梯式更替，小区计产面积70m2，每隔4个品系设一个CK，所得产量结果列于表12.2，试作分析

第13页：表12.2 小麦品系鉴定试验

表12.2 小麦品系鉴定试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

间比法)的产量结果与分析总和平均对照对品系各重复小区产量(kg)的%ⅠⅡⅢⅣⅤCK

第14页：例12.2 有12个小麦新品系鉴定试验（续）

例12.2 有12个小麦新品系鉴定试验（续）

首先，计算前后两个对照产量的平均数如A、B、C、D4品系的=(33.1+33.5)/2=33.3(kg)……，然后，计算各品系产量相对应产量的百分数，即得各品系的相对生产力。如品系A的相对生产力(%)=×100=109.9，……等

结果表明，相对生产力超过对照10%以上的品系有K、B、D、E、J、G6个，其中K品系增产幅度最大，达15.4%

间比法设计中，采用推广良种作为对照计算肥力指数调整供试家系产量，所以在参试家系数目较多时一般常用两个或两个以上的对照品种

第15页：第二节完全随机和随机区组试验的统计分析

第二节完全随机和随机区组试验的统计分析

一、完全随机试验设计的统计分析

二、随机区组试验结果的分析示例

三、随机区组的线性模型与期望均方

四、随机区组试验的缺区估计和结果分析

第16页：完全随机试验设计的统计分析

一、完全随机试验设计的统计分析

完全随机试验设计是指每一供试单位都有同等机会(等概率)接受所有可能处理的试验设计方法，没有局部控制，但要求在尽可能一致的环境中进行试验。它用于估算试验误差的自由度最多，统计显著性要求的F 值最小

第17页：随机区组试验结果的分析示例

二、随机区组试验结果的分析示例

随机区组试验结果的统计分析，可应用第六章所述两向分组单个观察值资料的方差分析法

这里可将处理看作A因素，区组看作B因素，其剩余部分则为试验误差。设试验有个处理，个区组，则其自由度和平方和的分解式如下

第18页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

总自由度=区组自由度+处理自由度+误差自由度(12·3)(12·2)

第19页：y表示各小区产量

y表示各小区产量

y表示各小区产量(或其他性状)，表示区组平均数，表示处理平均数，表示全试验平均数。总平方和=区组平方和+处理平方和+试验误差平方和

[例12.3]有一小麦品比试验，共有A、B、C、D、E、F、G、H8个品种(k=8)，其中A是标准品种，采用随机区组设计，重复3次(n=3)，小区计产面积25m2，其产量结果列于表12.3，试作分析

第20页：表12.3 小麦品比试验

表12.3 小麦品比试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

随机区组)的产量结果(kg)品种区组ⅠⅡⅢA = =11.6

第21页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

(1)自由度和平方和的分解

①自由度的分解

区组

品种

误差

②平方和的分解

矫正数

第22页：总，区组，品种，误差==22.97

总，区组，品种，误差==22.97

区组

品种

误差= =22.97

第23页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

(2)F测验将上述计算结果列入表12.4，算得各变异来源的MS值。表12.4 表12.3结果的方差分析变异来源DFSSMSFF0.05区组间 品种间 误差 总变异

第24页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

对区组间MS作F测验，在此有H0：，HA：、、不全相等( 、、分别代表区组Ⅰ、Ⅱ、Ⅲ的总体平均数)

得F = =8.40＞F0.05，所以H0应予否定，说明3个区组间的土壤肥力有显著差别。在这个试验中，区组作为局部控制的一项手段，对于减少误差是相当有效的(一般区组间的F测验可以不必进行，因为试验目的不是研究区组效应)

第25页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

对品种间MS作F测验，有H0：，HA：、、…、不全相等(、、…、分别代表品种A、B、…、H的总体平均数)，得F= =2.97＞F0.05，所以H0应予否定，说明8个供试品种的总体平均数有显著差异。需进一步作多重比较

(3)品种间平均数的多重比较

①最小显著差数法(LSD法)本例目的是要测验各供试品种是否与标准品种A有显著差异，宜应用LSD

第26页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

法。首先应算得品种间平均数或总和数差数的标准误。在以各品种的小区平均产量作比较时，差数标准误为：(12·4) 从而(12·5) 如果以各品种的小区总产量作比较，则因总产量大n 倍，故差数标准误为

第27页：12·4)—(12·7)中

12·4)—(12·7)中

(12·4)—(12·7)中，为方差分析表中的误差项均方MS；t值的，即误差项自由度。凡品种与对照的差异达到或超过者为显著，达到或超过者为极显著

如果试验结果需以亩产量表示，只要将总产量和总产量的LSD皆乘以cf即可。(12·6) 并有：(12·7)

第28页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

在此，如以各品种的小区平均产量(即表12.3的)进行比较，则

由于时，=2.145，=2.977，故

如对各品种的三个小区总产量(表12.3的)进行比较，则(kg)

第29页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

如以亩产量表示试验结果，则可算得化各品种总产量为亩产量的改算系数

因此，品种A的亩产量=(kg)(kg) (kg) (kg)

第30页：品种B的亩产量=

品种B的亩产量=

品种B的亩产量=(kg)

……，余类推

并且有

亩产量(kg)

亩产量(kg)

上述结果皆列于表12.5不论哪一种比较，结果都完全一样，只有E品种与对照有极显著的差异，其余品种都和对照没有显著差异

第31页：表12.5 表12.3资料各品种产量和对照相

表12.5 表12.3资料各品种产量和对照相比的

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

差异显著性品种的比较的比较亩产量的比较差异差异kg/亩差异E (CK)

第32页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

②新复极差测验(LSR法)如果我们不仅要测验各品种和对照相比的差异显著性，而且要测验各品种相互比较的差异显著性，则宜应用LSR法。首先，应算得品种的标准误SE

在小区平均数的比较时为(12·8)在小区总数的比较时为(12·9)

第33页：随机区组试验结果的分析示例（续）

随机区组试验结果的分析示例（续）

在亩产量的比较时为

然后，查附表8当时，自2至的SSR0.05和SSR0.01值，进而算得LSR0.05和LSR0.01值

本例如以小区平均数为比较标准，则有

查附表8，得到自由度、不同显著水平和秩次距p下的SSR值，进而算得LSR值(表12.6)。品种平均数差(12·10)(kg)

第34页：表12.7

表12.7

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

表12.6 表12.3资料新复极差测验的最小显著极差p

第35页：表12.7 表12.3资料的新复极差测验结果

表12.7 表12.3资料的新复极差测验结果产量

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第36页：结果表明

结果表明

E品种与H、C、F、A、D5个品种有5%水平上的差异显著性，E品种与D品种有1%水平上的差异显著性，其余各品种之间都没有显著差异

以上是以各品种的小区平均产量为比较标准。如以各品种总产量或亩产量为比较标准，则只要应用由(12·9)或(12·10)算出的SE值即可，方法类同，不再赘述

用时，仅需选择上述3种比较的任一种

第37页：随机区组的线性模型与期望均方

三、随机区组的线性模型与期望均方

(一) 线性模型

一个随机区组的试验结果，若以I 代横行(处理)，则i=1，2，…，k ；以j 代纵行(区组)，则j =1，2，…，n，整个资料共有k行n列。所以，在第i行j列的方格可以ij表示之(参见表12.3)。如果每一方格仅有一个观察值yij，则其线性模型为

第38页：上式中，为总体平均，为行的效应或处理效应

上式中，为总体平均，为行的效应或处理效应

上式中，为总体平均，为行的效应或处理效应，可为固定模型或随机模型，在固定模型中，假定，在随机模型中，假定～N(0，）；为列的效应或区组效应，一般为随机模型，假定～N(0，），若为固定模型则假定；而则为相互独立的随机误差，服从N(0，）。(12·11)

第39页：期望均方

（二）期望均方

随机区组的各种效应一般有3种模型(参见第六章)，即固定模型(称模型Ⅰ)，随机模型(称模型Ⅱ)和混合模型。这3种模型的期望均方(EMS)列于表12.8

第40页：表12.8 随机区组设计的期望均方变异来源D

表12.8 随机区组设计的期望均方变异来源DFM

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第41页：固定模型

固定模型

随机区组中仅有两个或较少的处理或品种时是适用

随机区组试验是随机模型，则表示处理(或品种)和区组都是从处理(或品种)总体和区组总体中随机抽取的。试验结论则推断到有关处理(或品种)和区组总体，而不是仅涉及某一特定处理(或特定品种)

第42页：随机区组试验的缺区估计和结果分析

四、随机区组试验的缺区估计和结果分析

缺区估计可采用最小二乘法

移项可得缺区估计值为：(12·13)

第43页：随机区组试验缺一个小区产量的结果分析

（一）随机区组试验缺一个小区产量的结果分析

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

[例12.4] 有一玉米栽培试验，缺失一区产量ye(kg)，其结果如表12.9，试作分析。表12.9 玉米随机区组试验缺一区产量(kg)的试验结果处理区组ⅠⅡⅢⅣA +yeC +ye +ye

第44页：首先，应估计出缺值ye

首先，应估计出缺值ye

首先，应估计出缺值ye。根据(12·12)可得

如将表12.9中有关数值代入(12·13)，也同样可得：即所以，(kg)

第45页：然后，将该置入表12.9ye的位置

然后，将该置入表12.9ye的位置

然后，将该置入表12.9 ye的位置，得表12.10。表12.10的形式和表12.3完全一样，因此可同样进行平方和的分解；但在分解自由度时需注意：因为是一个没有误差的理论值，它不占有自由度，所以误差项和总变异项的自由度都要比常规的少1个。由此得到的方差分析表如表12.11

第46页：表12.10 玉米随机区组试验结果处理区组Ⅰ

表12.10 玉米随机区组试验结果处理区组ⅠⅡⅢⅣ

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第47页：表12.11 玉米栽培试验

表12.11 玉米栽培试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

缺一区)的方差分析变异来源DFSSMSFF0.05区组 处理 误差 总变异 在进行处理间的比较时，一般用t 测验

对于非缺区处理间的比较，其仍由(12·4)式算出(假定以小区平均产量作为比较标准)，对于缺区处理和非缺区处理间的比较

第48页：12·14)上式的MSe为误差项均方

12·14)上式的MSe为误差项均方

(12·14)上式的MSe为误差项均方，n为区组数，k为处理数。在本例可求得：(kg)

第49页：随机区组试验缺二个小区产量的结果分析

（二）随机区组试验缺二个小区产量的结果分析

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

[例12.5]有一水稻栽培试验，假定缺失两区产量(yc和ya)，其结果如表12.12，试分析。表12.12 水稻随机区组试验缺两区产量(kg/小区)的试验结果处理区组ⅠⅡⅢⅣⅤⅥA +yaB +yc +yc22+ya187+yc+ya

第50页：首先，应估计出缺值yc和ya

首先，应估计出缺值yc和ya

首先，应估计出缺值yc和ya。采用解方程法，根据(12·12)，对yc有方程对ya有方程整理成二元一次联立方程

第51页：解之得

解之得

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

将yc≈18kg，ya≈10kg置入表12.12中，即得表12.13。然后进行方差分析，得到表12.14。由于表12.13中有两个缺区估计值，它们不占有自由度，故表12.14中误差项和总变异项的自由度均比通常的少2个。表12.13 水稻随机区组试验结果处理区组ⅠⅡⅢⅣⅤⅥA (10) (18)

第52页：表12.14 水稻随机区组试验

表12.14 水稻随机区组试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

缺二区)的方差分析

在进行处理间比较时，非缺区处理间比较的差数标准误仍由(12·4)给出(当以各处理小区平均数相比较时)

若相互比较的处理中有缺区的，则其平均数差数的标准误为变异来源DFSSMSFF0.05区组574.28处理 误差 总变异

第53页：上式中的为误差项均方

上式中的为误差项均方

上式中的为误差项均方，和分别表示两个相比较处理的有效重复数，其计算方法是：在同一区组内，若两处理都不缺区，则各记为1；若一处理缺区，另一处理不缺区，则缺区处理记0，不缺区处理记，其中k为试验的处理数目。例如，本试验在A和B比较时，(12·15)

第54页：解之得（续）

解之得（续）

A的有效重复数B的有效重复数在A和C比较时，A的有效重复数C的有效重复数故故(kg)(kg)

第55页：第三节拉丁方试验的统计分析

第三节拉丁方试验的统计分析

一、拉丁方试验结果的分析示例

二、拉丁方的线性模型与期望均方

三、拉丁方试验的缺区估计和结果分析

第56页：拉丁方试验结果的分析示例

一、拉丁方试验结果的分析示例

拉丁方试验在纵横两个方向都应用了局部控制，使得纵横两向皆成区组

设有k个处理(或品种)作拉丁方试验，则必有横行区组和纵行区组各k个，其自由度和平方和的分解式为：(12·16)

第57页：拉丁方试验结果的分析示例（续）

拉丁方试验结果的分析示例（续）

总自由度=横行自由度+纵行自由度+处理自由度+ 误差自由度(12·17)

总平方和= 横行平方和+纵行平方和+处理平方和+误差平方和

第58页：例12.6 有A、B、C、D、E 5个水稻

例12.6 有A、B、C、D、E 5个水稻

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

品种作比较试验，其中E为标准品种，采用5×5拉丁方设计，其田间排列和产量结果见表12.15，试作分析。

表12.15 水稻品比5×5拉丁方试验的产量结果(kg)横行区组纵行区组ⅠⅡⅢⅣⅤⅠD(37)A(38)C(38)B(44)E(38)195ⅡB(48)E(40)D(36)C(32)A(35)191ⅢC(27)B(32)A(32)E(30)D(26)147ⅣE(28)D(37)B(43)A(38)C(41)187ⅤA(34)C(30)E(27)D(30)B(41) =882

第59页：表12.15 算得各横行区组总和和各纵

表12.15 算得各横行区组总和和各纵

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

行区组总和，并得全试验总和

再在表12.16算得各品种的总和和小区平均产量

然后进入以下步骤：表12.16 表12.15资料的品种总和和品种平均数品种A38+35+32+38+34= +48+32+43+41= +32+27+41+30= +36+26+37+30= +40+30+28+27=

第60页：例12.6 有A、B、C、D、E 5个（续）

例12.6 有A、B、C、D、E 5个（续）

(1)自由度和平方和的分解

①自由度的分解由(12·16)可得

横行

纵行

品种

误差

②平方和的分解由(12·17)可得

第61页：矫正数，总，横行区组，纵行区组

矫正数，总，横行区组，纵行区组

矫正数

横行区组

纵行区组

第62页：表12.17 算得各变异来源的MS值

表12.17 算得各变异来源的MS值

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第63页：表12.17 表12.15资料的方差分析

表12.17 表12.15资料的方差分析

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

对品种间作F 测验，，、…、不全相等( 、、…、)分别代表A、B、…、E品种的总体平均数)得，变异来源DFSSMSFF0.05横行区组 纵行区组 品种 试验误差 总变异

第64页：F0.05=3.26，所以应被否定

F0.05=3.26，所以应被否定

F0.05=3.26，所以应被否定，即各供试品种的产量有显著差异。(3)品种平均数间的比较①最小显著差数法(LSD法)应用(12·4)，得当=12时，，，应用(12·5)(kg)(kg)(kg)

第65页：例12.6 有A、B、C、D、E 5个（续）

例12.6 有A、B、C、D、E 5个（续）

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

以之为尺度，在表12.18测验各品种对标准品种(E)的差异显著性。结果只有B品种的产量极显著地高于对照，其余品种皆与对照无显著差异。表12.18 表12.15资料各品种与标准品种相比的差异显著性品种小区平均产量(kg)差异B (CK)32.6**表示达1%显著水平

第66页：例12.6 有A、B、C、D、E 5个（续）

例12.6 有A、B、C、D、E 5个（续）

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

②新复极差测验(LSR 法) 根据(12·8)计算，求得

再根据=12时的和的值算得=2，3，4，5时的和的值于表12.19。根据表12.19的和的尺度，测验各品种小区平均产量的差异显著性于表12.20。表12.19 表12.15资料各品种小区平均产量( )互比时的LSR 值(kg)

第67页：表12.20 水稻品比试验的新复极差测验

表12.20 水稻品比试验的新复极差测验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

由表12.20可见，B品种与其他各品种的差异都达到=0.05水平，而B品种与D、E品种的差异达到=0.01水平，A、C、D、E 4品种之间则无显著差异

小区平均产量(品种)差异显著性

第68页：拉丁方的线性模型与期望均方

二、拉丁方的线性模型与期望均方

用代表拉丁方的横行、纵行的交叉观察值，再以代表处理，则拉丁方试验的线性模型为：(12·18)

上式中为总体平均数；为横行效应，为纵行效应，若两者为固定模型，有

第69页：拉丁方的线性模型与期望均方（续）

拉丁方的线性模型与期望均方（续）

若两者均为随机模型，有～N(0，)，～N(0，)；为处理效应，固定模型有，随机模型时～N(0，)；相互独立的随机误差～N(0，)

拉丁方设计的固定模型和随机模型的期望均方如下表

第70页：拉丁方的线性模型与期望均方（续）

拉丁方的线性模型与期望均方（续）

拉丁方设计的期望均方变异来源DF固定模型随机模型横行间纵行间处理间试验误差

第71页：拉丁方试验的缺区估计和结果分析

三、拉丁方试验的缺区估计和结果分析

缺值ye的估计公式为：(12·19)

移项可得：(12·20)

第72页：拉丁方试验的缺区估计和结果分析（续）

拉丁方试验的缺区估计和结果分析（续）

当仅有一个缺区时，可由(12·19)或(12·20) 直接解得ye 值

当有多个缺区时，可由(12·19)建立联立方程组，解出各个缺区估计值

[例12.7] 有一甘蔗品比试验，采用5×5拉丁方设计，缺失一区产量，其结果见表12.22，试求该缺区估计值ye并作分析

第73页：表12.22 5×5甘蔗试验缺失一区产量的试

表12.22 5×5甘蔗试验缺失一区产量的试验结果

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第74页：首先求缺区估计值ye

首先求缺区估计值ye

首先求缺区估计值ye。将表12.21的有关数值代入(12·19)可得

同样，代入(12·20)得：(100kg)

第75页：拉丁方试验的缺区估计和结果分析（续）

拉丁方试验的缺区估计和结果分析（续）

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

将置入表12.22的ye位置，得表12.23。表12.23 5×5甘蔗试验具有一个估计值的试验结果(100kg/区)横行区组纵行区组ⅠⅡⅢⅣⅤⅠA 14E 22D 20C 18B 2599ⅡD 19B 21A 16E 23C 1897ⅢB 23A 15C 20D 18E 2399ⅣC 21D (18)E 24B 21A 23C 16B 23A 17D =495A=79B=113C=93D=95E=

第76页：表12.23 可按没有缺区的拉丁方资料作出方

表12.23 可按没有缺区的拉丁方资料作出方差分析

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

仅误差项和总变异项的自由度比没有缺区的拉丁方资料少一个，因为有一个缺区估计值，它不占有自由度

由此所得的结果列于表12.24

表12.24 甘蔗5×5拉丁方试验(缺一区)的方差分析变异来源DFSSMSFF0.05横行 纵行 品种 误差 总变异

第77页：拉丁方试验的缺区估计和结果分析（续）

拉丁方试验的缺区估计和结果分析（续）

在对各品种的小区平均数作t测验时，没有缺区品种间的比较仍用(12·4)；但当缺区品种与非缺区品种比较时，其差数标准误应为：(12·21) 如在本例中，(100kg)

第78页：拉丁方试验的缺区估计和结果分析（续）

拉丁方试验的缺区估计和结果分析（续）

以上是有一个缺区的拉丁方试验的分析。如果拉丁方试验有几个缺区，则首先应算得各个缺区的估计值。这些估计值可由(12·19)建立联立方程解出。算得各缺区估计值后，可按正常(没有缺区的)拉丁方资料计算各变异来源的平方和，但误差项和总变异项的自由度要比正常的少l个（l为缺区数目）。在对各处理小区平均数作t测验时，没有缺区的处理间比较的差数标准误仍由(12·4)给出；若相互比较的处理中有缺区存在，则其平均数差数的标准误为

第79页：12·22)中的为误差项均方

12·22)中的为误差项均方

(12·22)中的为误差项均方，和分别代表两个相互比较的处理的有效重复数，其计算方法是

1. 若相互比较的甲、乙二处理在横行和纵行皆不缺区，则分别记1
2. 若甲处理不缺区，而其所在的横行或纵行的乙处理缺一区，则甲记2/3
3. 若甲处理不缺区，而其所在的横行和纵行的乙处理皆缺区，则甲记1/3
4. 若甲处理本身为缺区，则记0。例如，有一5×5拉丁方试验为：(12·22)

第80页：拉丁方试验的缺区估计和结果分析（续）

拉丁方试验的缺区估计和结果分析（续）

以上有圆圈者表示缺区。则在B、C、E处理间比较时，其用(12·4)，其余各处理相互比较都要先计算有效重复数，再代入(12·22)计算

如：A与E比较时

A的有效重复数AEDCBDBAECBACDEC

第81页：E的有效重复数故而A与D比较时

E的有效重复数故而A与D比较时

E的有效重复数故而A与D比较时，A的有效重复数D的有效重复数故其余类推

第82页：第四节试验处理的合并比较

第四节试验处理的合并比较

正交性(orthogonality)：任两比较间的系数乘积之和为0

正交对比或正交比较(orthogonal comparison)：若在所有对比(比较)中，两两比较间的比较系数( )乘积之和都为0的比较。这种比较系数称为正交系数(orthogonal coefficient)

第83页：第四节试验处理的合并比较（续）

第四节试验处理的合并比较（续）

例如，处理为氨水1(A)、氨水2 (B)、碳酸氢铵(C)、尿素(D)和不施肥(E)按完全随机设计进行试验的例6.10，就预先安排了以下几个特定的比较

①施肥与不施肥，即A+B+C+D与E

②液态氮与固态氮，即A+B与C+D

③液态氮之间的比较，即A与B

④固态氮之间的比较，即C与D

可以看出，比较①和②是将处理合并后进行比较的，因此，这种比较称为试验处理的合并比较。在

第84页：比较②、③、④中

比较②、③、④中

比较②、③、④中，相比较的处理个数是相同的，这时在同一比较的所有处理都占相同比例，其系数均为1，但是，在比较①中，用A+B+C+D 4处理与E 1个处理比较，其处理数是不同的，显然不是均衡比较，缺乏可比性，应该将4处理的合并值与E处理的4倍进行比较，即A、B、C、D各用1份而E用4份才有可比性，换句话说，前者的系数为1，后者的系数为4。

同时，在同一比较(或对比)中，一方系数为正，另一方系数为负，如比较②，若A与B的系数为正，则C与D的系数为负。按此方法将4种比较的系数填于表12.25

第85页：表12.25 表6.11资料单一自由度比较

表12.25 表6.11资料单一自由度比较的正交系

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

数(Ci)和计算处理ABCDESSQ(MS) 比较①A+B+C+D对E ②A+B对C+D ③A对B ④C对D 总和301.20

第86页：第四节试验处理的合并比较（续）

第四节试验处理的合并比较（续）

可以发现，该比较①～④中，任两比较间的系数乘积之和为0，如①与②对比

这称为正交性(orthogonality)。若在所有对比(比较)中，两两比较间的比较系数()乘积之和都为0，则称这种对比为正交对比或正交比较(orthogonal comparison)，这种比较系数称为正交系数(orthogonal coefficient)

第87页：例12.8 在采用完全随机设计

例12.8 在采用完全随机设计的

表6.11资料中，已事先确定要研究以下4种比较的差异显著性；①施肥对不施肥，②施液体肥与施固体肥，③施氨水1对施氨水2，④施碳酸氢铵对施尿素。试作比较

分析比较步骤

1. 将资料各处理的总产量列于表12.25(为便于计算，不用平均产量，但后面所得结果仍是关于平均产量的)
2. 写出各个预定比较的正交系数Ci(见上表)

按上述方法获得正交系数Ci后，可以计算每一比较的差数

第88页：例12.8 在采用完全随机设计（续）

例12.8 在采用完全随机设计（续）

(12·23)因而，由进一步计算每一比较的SS(即MS，因为每一比较的自由度都是1）：(12·24)

第89页：例如比较①的同理

例如比较①的同理

第90页：例如比较①的同理（续）

例如比较①的同理（续）

这里可注意，正是表6.12中的。也就是我们已将表6.12处理间具4个自由度的平方和再分解为属于4个独立比较的平方和，各具，因此，将这种比较也称为单一自由度的独立比较(independentcomparisonofsingledegreeoffreedom)，其方差分析表于12.26(与表6.12对照）

第91页：表12.26 表6.11资料单一自由度的方差

表12.26 表6.11资料单一自由度的方差分析变

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

异来源DFQSSQMSQF施肥对不施肥 施固体肥对施液体肥 施氨水1对施氨水 施尿素对施碳酸氢铵 试验误差

第92页：例如比较①的同理（续）

例如比较①的同理（续）

将表中各个MSQ与MSe比，得到F值，查F表当，时，，，结果表明该试验预定的4个比较中，施肥对不施肥、施固体肥对施液体肥的差异极显著，其余两种比较的差异不显著

如果要计算平均数，可由下式求得：(12·25)式中表示比较中取正值的正交系数，如(g)

第93页：例如比较①的同理（续）

例如比较①的同理（续）

(g)即表示施肥比不施肥平均每盆增产7.875克，施固体肥比液体肥平均每盆增产4.250克，皆为极显著。由上分析可知，处理的合并比较十分简便，值得推广应用。正确进行处理合并比较的关键是正确确定比较的内容和正确写出比较的正交系数。为此，须满足下列3个条件

1. 比较的数目必须为k-1，以使每一比较占有而且仅占有1个自由度
2. 每一独立比较的正交系数之和必须为0，即

第94页：例如比较①的同理（续）

例如比较①的同理（续）

以使每一比较都是均衡的

(3)任何两个独立比较的相应正交系数乘积之和必须为0，即，以保证却好分解为个。例如表12.25比较①和②的比较①和③的……等，如果出现，则比较就不再是独立的，其后果是各个比较的之和必不等于

第95页：例如比较①的同理（续）

例如比较①的同理（续）

这一点在试验设计时需特别注意。如有相等的3个处理，其总和数分别为、、，则2-1-1或 等皆为两个独立的比较，而1-10或

第96页：等则为不独立的，因为前者的，后者的

等则为不独立的，因为前者的，后者的

等则为不独立的，因为前者的，后者的

在具体写出各个独立比较的正交系数时，可按下列规则进行

(1) 若被比较的两个组的处理数目相等，则给一个组的各处理以系数+1，另一个组的各处理都是系数-1。到底哪个组取+1，哪个组取-1，可以随便，一般以“在前”的组(如施肥对不施肥的比较，施肥在前)取“+”号

第97页：例如比较①的同理（续）

例如比较①的同理（续）

(2) 若被比较的两个组的处理数目不相等，则第一组的系数为第二组的处理数，第二组的系数为第一组的处理数，例如2个处理(第一组)与3个处理的一个比较，其写作3，3，-2，-2，-2

1. 如果写出的有公约数，则应将其约为最小的整数，例如4个处理与2个处理的一个比较，按规则
2. 其为2，2，2，2，-4，-4，应简化成1，1，1，1，-2，-2

(4)如果某一处理已经和所有其余处理作过一次比较(如表12.25的处理)，则该处理不能再参加其余比较，否则一定破坏了

第98页：例如比较①的同理（续）

例如比较①的同理（续）

因子式试验在供试因子增加时，处理组合数迅速增加，这给试验带来了巨大的工作量和试验误差的增大。若试验目的并不在于研究各因子间的交互作用而是主要了解各因子的主效。那么，有些试验可以删去一些次要的处理组合。这时，因子间不是正交的关系，而是一种分枝式关系。如例12.8和例12.9

第99页：例12.9 现有一肥料比较的盆钵试验

例12.9 现有一肥料比较的盆钵试验

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

每处理含氮量相等，重复六个盆钵，每盆栽5棵稻苗，随机区组排列，成熟后收获计产，结果列于表12.27。试进行处理的合并比较。表12.27 水稻肥料盆钵试验产量结果(单位：10g)变异来源DFSSMSFF0.01区组 处理 误差 总

第100页：例12.9 现有一肥料比较的盆钵试验（续）

例12.9 现有一肥料比较的盆钵试验（续）

然而，表12.28的结果并未研究因子的主效。为此，应该进行处理的合并比较

①施与未施，即A与B+C+D+E+F+G+H

②化肥与有机肥，即B+C与D+E+F+G+H

③化肥间，即B与C

④畜粪尿与堆肥间，即D+E与F+G+H

⑤畜粪尿间，即D与E

⑥堆肥老方法与新方法间，即F与G+H

第101页：表12.29

表12.29

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

对于比较①，施肥处理合并与未施肥相比其效应的总和Q1与其平方和分别为

其余比较同理可得，并将结果一并填入表12.29

并将与误差均方比较，进行单一自由度的F 测验( ，)，其结果也见表12.29

第102页：表12.29 水稻肥料盆钵试验各种处理合并

表12.29 水稻肥料盆钵试验各种处理合并的计算和

本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。

第103页：例12.9 现有一肥料比较的盆钵试验（续）

例12.9 现有一肥料比较的盆钵试验（续）

由表12.29可知，比较①、②、④、⑤、⑥均存在极显著差异，比较⑦差异显著，只有比较③不显著