第十四章不完全区组设计和统计分析课件

第十四章不完全区组设计和统计分析第一节不完全区组设计的主要类型第二节重复内分组和分组内重复设计的统计分析第三节简单格子设计的统计分析第四节平衡不完全区组设计的统计分析第十四章不完全区组设计和统计分析第一节不完全区组设1第一节不完全区组设计的主要类型一、田间试验常用设计的归类二、重复内分组和分组内重复设计三、格子设计四、平衡不完全区组设计第一节不完全区组设计的主要类型一、田间试验常用设计的归类2一、田间试验常用设计的归类完全区组(completeblock)：每一区组包含全套处理。不完全区组(incompleteblock)：即一套处理分成几个区组，或一个区组并不包含全部处理，但同样要通过区组实施地区控制。

一、田间试验常用设计的归类完全区组(completeblo3二、重复内分组和分组内重复设计重复内分组设计(blockinreplication)：将供试品种分为几个组，看作为主区，每个组内包含的各个品种看作为副区，重复若干次，主副区都按随机区组布置的设计。例如20个品种，分为4组，每组包含5个品种，若重复3次，则田间布置可设计如下图：

二、重复内分组和分组内重复设计重复内分组设计(blocki4

重复内分组设计的田间布置该例中重复内分组设计的自由度分析如下：重复Ⅰ重复Ⅱ重复Ⅲ区组(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)420111017751591912331815816621381713121913918811271615251612620103146201441171471994111018115重复Ⅰ重复Ⅱ重复Ⅲ区组(1)(2)(3)(4)(5)(6)(5变异来源DF重复2组间3误差(Ea)6组内品种间16误差(Eb)32总59组内品种间比较的误差将为：；变异来源DF6各组平均数间比较的误差将为：；不同组品种间比较的误差(仿照裂区的情况)将为：。由于Ea与Eb常取不同数值，Ea往往大于Eb，例如

=3，若如此，则：组内品种间比较的误差将为：不同组品种间比较的误差将为：各组平均数间比较的误差将为：7两者比值为：即不同组品种间比较的方差将比组内品种间比较的方差大40%，因而像这种不完全区组设计的方法，并不能保证任何两个品种间比较具有相近的精确度。分组内重复设计(replicationinblock)：将供试材料分组后放在连片土地上的几组随机区组试验，通过土地连片而进行联合分析与比较。

两者比值为：8

分组内重复设计

分组1分组2分组3分组4区组(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)191618131511541898181917121114352710716201915121313510692017161114152149710171820141312423686分组1分组2分组3分组4区组(1)(2)(3)(4)(5)(9三、格子设计格子设计(latticedesign)：为了克服重复内分组设计中组间品种比较和组内品种比较精确度悬殊的问题，对品种分组的方法可考虑从固定的分组改进为不固定的分组，使一个品种有机会和许多其他品种，甚至其他各个品种都在同一区组中相遇过。三、格子设计格子设计(latticedesign)：为了10(一)格子设计的类别平方格子设计(squaredlattice)：供试品种数为区组内品种数的平方，区组内品种数为p，供试品种数为p2；立方格子设计(cubiclattice)：供试品种数为区组内品种数的立方，区组内品种数为p，供试品种数为p3；矩形格子设计：区组内品种数为p，供试品种数为p(p+1)。

(一)格子设计的类别11(二)平方格子设计1.仿照随机区组式的设计按品种分组方法的变换次数有：(1)简单格子设计(simplelattice)品种分组方法为二种，试验重复次数为2或2的倍数。

重复I重复Ⅱ(1)123(4)147区组(2)456(5)258(3)789(6)369(二)平方格子设计重复I重复Ⅱ(1)123(4)12(2)三重格子设计(triplelattice)：品种分组方法为三种，即在简单格子设计二种分组方法的基础上再增加对角线分组一种，重复次数为3或3的倍数。

(3)四重格子设计(quadruplelattice)：在三重格子设计的基础上，再增加对角线一组，

重复I重复Ⅱ重复III(1)123(4)147(7)159区组(2)456(5)258(8)267(3)789(6)369(9)348(2)三重格子设计(triplelattice)：品种分13

(4)平衡格子设计(balancedlattice)：品种分组方法增加到使每一对品种都能在同一区组中相遇一次。分组法X分组法Y分组法Z分组法L区组(1)12345(6)16111621(11)17131925(16)18151724(2)678910(7)27121722(12)28142021(17)29111825(3)1112131415(8)38131823(13)39151622(18)310121921(4)1617181920(9)49141924(14)410111723(19)46132022(5)2122232425(10)510152025(15)56121824(20)571416235×5四重格子设计方法分组法X分组法Y分组法Z分组法L区组(1)12314

2.仿照拉丁方的格子设计(1)平衡格子方设计(balancedlatticesquare)

①重复数r=(p+1)/2，每对品种在行或列区组中共相遇一次；

重复Ⅰ重复Ⅱ重复Ⅲ重复Ⅳ(1)123(4)147(7)159(10)168区组(2)456(5)258(8)267(11)249(3)789(6)369(9)348(12)3573×3平衡格子设计重复Ⅰ重复Ⅱ重复Ⅲ重复Ⅳ(1)123(4)1415

3×3平衡格子方设计[在行或列中相遇一次，r=(p+1)/2]ⅠⅡ123168456924789573ⅠⅡ12316845692416②重复数r=(p+1)，每对品种在行及列区组中均相遇一次，亦即共相遇二次。

ⅠⅡⅢ159131234111166261014658712251537111511129101483948121616151413713104ⅣⅤ17121411015882131192716101635136312159645141144×4平衡格子方设计[在行及列中共相遇二次，r=(p+1)]

②重复数r=(p+1)，每对品种在行及列区组中均相遇一次，亦17(2)部分平衡格子方设计(partiallybalancedlatticesquare)：重复次数少于最小平衡重复数。与三重、四重格子设计类似，不一定每一对品种都在行或列区组中相遇。格子设计的优点是：考虑了供试品种间平衡比较的问题。但由于供试品种数多，这常只能实施部分平衡，而事实上很难实施完全平衡，因为完全平衡所需的重复次数导致试验规模过大。(2)部分平衡格子方设计(partiallybalanc18育种工作中产量比较在早、中期阶段，因供试材料多需要考虑适合大量处理的设计，但这时每份材料的种子数少，一般不可能进行小区较大的精确试验，因而实际应用中部分平衡的格子设计已可满足要求。育种工作中产量比较在早、中期阶段，因供试材料多需要考虑适合大19四、平衡不完全区组设计平衡不完全区组设计(balancedincompleteblockdesign)：设计的供试处理数不多，不须按格子设计那样每一重复包含有区组大小为k的k个区组，而可将各重复寓于全部区组之中，区组数与区组大小不一定相等，即全试验包括大小为k的区组共t(处理数)或t倍个。

四、平衡不完全区组设计平衡不完全区组设计(balanced20图14.7一种平衡不完全区组设计例如品尝试验，对于一个人的味觉来说，品尝的对象增加太多时鉴别差异的灵敏度便下降，因而每个人只能品尝一部分。图14.7的情况，若有7个水果品种供鉴评，每人品尝3个，请7位品尝家作鉴评，便共品尝21次，每个品种品尝3次。此处每位专家区组（1）（2）（3）（4）（5）（6）（7）123456723456714567123区组（1）（2）（3）（4）（5）（6）（7）123456721便是一个区组，每区组包含3个品种。这时尽管每人并未将7个品种全部鉴评过，但因是均衡的，每个品种至少和其他6个品种比较过1次。这一试验可增加至14位专家则每对品种相遇2次，21位专家则相遇3次。因而可以请许多专家作出综合评判。便是一个区组，每区组包含3个品种。这时尽管每人并未将7个22第二节重复内分组和分组内重复设计的统计分析一、重复内分组设计的统计分析二、分组内重复设计的统计分析第二节重复内分组和分组内重复设计的统计分析一、重复内分组23一、重复内分组设计的统计分析重复内分组用于品种(系)试验时有二种情况：一是大量品种(系)间的比较目的在于选拔高产优系（固定模型试验）；另一是从一个群体内随机抽出大量家系进行试验，通过供试的样本推论总体的情况（随机模型试验）。

一、重复内分组设计的统计分析重复内分组用于品种(系)试验时有24假定重复内分组设计的供试品种为m=a×b个，分a组，每组有b个品种(系)，重复r次，则重复内分组设计的线性模型为：

(14·1)

固定模型时：，，～，～；随机模型时：Ak～，Bkl～～，～。

假定重复内分组设计的供试品种为m=a×b个，分a组，每组有b25重复内分组设计的自由度及期望均方变异来源DFMSEMS固定模型随机模型重复r-1MS1分组(区组，主区)a-1MS2重复×分组(Ea)(r-1)(a-1)MS3分组内品种(系)a(b-1)MS4重复×分组内品种(系)(Eb)a(b-1)(r-1)MS5变异来源DFMSEMS固定模型随机模型重复26固定模型时分组间差异的测验，F=MS2/MS3；分组内品种(系)间差异的测验F=MS4/MS5。重复内分组设计着重在分组内品种间的比较，其分组间比较，其

(14·3)(14·2)固定模型时分组间差异的测验，F=MS2/MS3；(1427不同组品种间比较，其

(14·4)

随机模型时分组间变异的测验:

(14·5)分组内变异的测验：

F=MS4/MS5

(14·6)

不同组品种间比较，其28F=(MS2+MS5)/(MS3+MS4)时，其有效自由度可用Satterthwaite公式计算：

(14·7)(14·7)中fi为各均方对应的自由度。由(14·5)及(14·6)的关系可分别估计出及。

F=(MS2+MS5)/(MS3+MS4)时，其有效自由度可29二、分组内重复设计的统计分析分组内重复的设计的线性模型为：

(14·8)

固定模型时：，，～；随机模型时，Ak～，Bkl～，二、分组内重复设计的统计分析分组内重复的设计的线性模型为：30～。分组内重复设计的自由度及期望均方变异来源DFMSEMS固定模型随机模型分组a-1MS1分组内品种a(b-1)MS2分组内重复(区组)a(r-1)MS3重复×组内品种(E)a(b-1)(r-1)MS4～。变异来源D31固定模型时分组间差异的测验，F=MS1/MS4；分组内品种(系)间差异的测验F=MS2/MS4。分组内重复设计着重在分组内品种间的比较，其

(14·9)

分组间可以比较，其

(14·10)固定模型时分组间差异的测验，F=MS1/MS4；32不同组品种间的比较，其

(14·11)随机模型时分组间差异的测验:

(14·12)其有效自由度按Satterthwaite公式。分组内品种间差异测验：

F=MS2/MS4

(14·13)不同组品种间的比较，其33由(14·12)及(14·13)测验及。在各分组品种(系)均为总体一随机样本的前题下，可假定分组平均数相等，从而对品种(系)平均数作统一调整。重复内分组和分组内重复是目前品系产量早期比较试验较常用的设计，并常用于遗传参数的估计，尤其前者更为常用。由(14·12)及(14·13)测验及34第三节简单格子设计的统计分析一、简单格子设计分析的基本原理二、简单格子设计的例题第三节简单格子设计的统计分析一、简单格子设计分析的基本原35一、简单格子设计分析的基本原理设有9个品种，重复2次的简单格子设计试验，这9个品种分别给以二位数的代号如下：品种按横行、纵行分组，分别设置为一个重复，则其分组安排如下：123111213456212223789313233一、简单格子设计分析的基本原理设有9个品种，重复2次的简单格36由重复Ⅰ所得产量以x表示，重复Ⅱ以y表示，各品种总和以t表示，则可以将试验结果整理如表14.3的形式(虚线表示区组)。

重复Ⅰ111213212223313233重复Ⅱ112131122232132333x11x12x13X1·y11y12y13Y1·t11t12t13T1·x21x22x23X2·y21y22y23Y2·t21t22t23T2·x31x32x33X3·y31y32y33Y3·t31t32t33T3·X·1X·2X·3X··Y·1Y·2Y·3Y··T·1T·2T·3T··X组Y组品种总和简单格子设计试验结果符号表重复Ⅰ111213212223313237横行总和作为试验因子A(X分组)的效应，纵列为B(Y分组)的效应。此试验可看作为每个因子各具3个级别的二因子试验，其自由度为：由于重复Ⅰ中A因子的效应和区组效应混杂，重复Ⅱ中B因子与区组混杂，整个试验相当于一个虚拟的二因子部分混杂试验，其混杂的效应是A与B主效。DFA2B2A×B4总8横行总和作为试验因子A(X分组)的效应，纵列为B(Y分组)的38若将重复当作区组，那么本试验可按随机区组的方法进行方差分析，其自由度为（左图）现在每一重复又划分为区组，要把区组的变异从误差中扣去以减小试验误差，故其自由度分析将为（右图）DF重复1品种8误差8总17DF重复1区组(Eb)4品种8区组内误差(Ei)4总17若将重复当作区组，那么本试验可按随机区组的方法进行方差分析，39由t11、t12、…、t33计算品种平方和中包含有区组的效应，夸大了品种的效应；由X1·、X2·、X3·，Y·1、Y·2、Y·3计算区组平方和则又包含了品种的效应，夸大了区组的效应。关键：从品种效应中扣去区组部分，得到可以共同比较的调整的品种平均数及品种平方和；估计出除去品种效应的区组间变异，得到一个无偏的试验误差估计，进行合理的统计推断。由t11、t12、…、t33计算品种平方和中包含有区组的效应40(一)品种调整平均数的计算

1·=T1·/6为A因子第一级别的未调整平均数；

·1=T·1/6为B因子第一级别的未调整平均数。如品种12的未调整平均数为v12，则:

(14·14)其中，m为全试验总平均数。(一)品种调整平均数的计算41(14·14)说明任一品种总的离均差为横行离均差、纵行离均差以及横行×纵行互作效应三部分之和。令：Ai表示不包含区组效应A因子效应估计值；

Bi表示不包含区组效应B因子效应估计值。则：A因子第一个级别的估计值，

B因子第一个级别的估计值(14·14)说明任一品种总的离均差为横行离均差、纵42又令Ab表示与区组混杂的A因子效应估计值，Bb表示与区组混杂的B因子效应估计值则A因子第一个级别的估计值，

B因子第一个级别的估计值

若A0，B0分别表示X组及Y组综合在一起未调整的A因子及B因子效应，则：又令Ab表示与区组混杂的A因子效应估计值，43

求A及B的调整值比较合理的方法是以Ai、Bi及Ab、Bb各分组所获得结果的可靠程度进行加权平均，这里Ai、Bi效应没有区组效应在内，可用衡量其可靠程度，其中代表区组内误差的理论方差。Ab、Bb效应混有区组效应，区组效应越大，Ab、Bb估计A及B的可靠程度越小，可用衡量其可靠程度，代表重复内区组间的理论方差(以小区为单位)。(14·15)(14·15)44

(14·16)当区组间没有真实差异时，，Ai、Bi和Ab、Bb同等重要，故：第十四章不完全区组设计和统计分析课件45得到A及B的估计值后，可得：

(14·17)因未调整的(v0-A0-B0+m)与调整后的(v-A-B+m)应是相等的，两者相减

v-v0=(A-A0)+(B-B0)

(14·18)表示调整的品种平均数可由v0、(A-A0)及(B-B0)三部分计算。得到A及B的估计值后，可得：46由(14·16)及(14·15)可得：令则

(14·19)由(14·16)及(14·15)可得：(14·19)47

以品种11为例，需求出A及B各第一级别的A0、Ab、

B0及Bb，其中以品种11为例，需求出A及B各第一级别的A0、Ab、48若令以上二矫正数分别以及代表，则：

(14·20)

其中vef中的ef代表以二位数字表示的某品种，在具有二个重复参试材料为p2的简单格子设计中及的通式可写为：

第十四章不完全区组设计和统计分析课件49

如果简单格子设计，每种分组重复二次，全试验共有四次重复，则：

(14·21)(14·22)(14·21)(14·22)50

在品种平均数的横行及纵行旁求出，求

出，就可计算出各个品种的调整平均数。但为便于计算，一般直接在品种总和表旁求出品种总和的矫正数，计算出各个品种的调整总和，再求调整平均数。2次重复时调整品种总和为：

(14·23)

在品种平均数的横行及纵行旁求出，51(二)与及w与的估计上述品种调整平均数的计算需按，

进行调整。

可以由区组内均方Ei直接估计，主要需估计出。区组间均方的计算需由二部分平方和合并，要了解清楚这二部分平方和的计算，从一个四次重复的试验比较容易说明。(二)与及w与的估计52

表14.4四次重复简单格子设计试验结果符号表X分组法Y分组法ⅠⅡⅢⅣ111213g11111213g12111213111213212223g21212223g22212223212223313233g31313233g32313233313233G1G2g13g23g33G3g14g24g34G4x11x12x13X1·y11y12y13Y1·t11t12t13T1·x21x22x23X2·y21y22y23Y2·t21t22t23T2·x31x32x33X3·y31y32y33Y3·t31t32t33T3·X·1X·2X·3X··Y·1Y·2Y·3Y··T·1T·2T·3T··X分组法Y分组53在X、Y两种分组各有重复时，从相同品种组的区组两次重复间的差异的效应扣去整个重复间差异的效应，可以估计出区组效应。其计算方法为(14·24)二式之和。

(14·24)在X、Y两种分组各有重复时，从相同品种组的区组两次重复间的54这部分平方和相当于A因子与重复的互作和B因子与重复的互作之和，称为成分(a)。两种分组方法各对应X1·与Y1·之间差异的效应扣去整个分组方法总差异间的效应，也将属于区组的效应，其计算方法为(14·25)二式之和。

(14·25)这部分平方和相当于A因子与重复的互作和B因子与(14·25)55这部分平方和相当于A因子与分组方法的互作和B因子与分组方法的互作之和，称为成分(b)。因T1·-2X1·=(X1·+Y1·-2X1·)=Y1·-X1·故成分(b)也可写为：

(14·26)

这部分平方和相当于A因子与分组方法的互作和B因56在3×3简单格子设计具有4个重复时，成分(a)具有2+2=4个自由度，成分(b)也具有2+2=4个自由度，(a)与(b)两者相加共有8个区组自由度。在只有2个重复时，显然成分(a)无从计算，因此仅由成分(b)代表区组的平方和。不过(14·26)中分母将相应改变为2×3及2×9。

在3×3简单格子设计具有4个重复时，成分(a)具有57分析成分(a)均方所估计的方差分量为，其中为区组内误差，为区组间的方差。成分(b)均方所估计的方差分量为，这是因为成分(b)的两部分是从同一材料计算来的，所以只估计了。当只有二个重复时，只能由成分(b)计得区组的均方()，但是由方差分析原理，正常的区组项均方应由组成。所以对区组的理论方差的估计要作适当调整。

分析成分(a)均方所估计的方差分量为，其中58所以，

(14·27)

当有四次重复时，成分(a)与(b)综合的均方所估计的分量为，即第十四章不完全区组设计和统计分析课件59所以，

(14·28)

(三)品种平均数间比较的误差计算同区组内品种间比较：第十四章不完全区组设计和统计分析课件60

异区组品种间比较：

不论区组异同，品种间相互比较：(14·29)(14·30)(14·29)(14·30)61

若由成分(a)单独估计，则，。当Eb≤Ei时，，上列各公式均变为，这就类似随机区组时的公式。当Eb很大时，接近于1，(14·29)、(14·30)、(14·31)三公式相应变为：

(14·31)(14·31)62

，和这种情况下，A与B的效应相当于由Ai及Bi单独估计，Ab及Bb对A、B均未提供信息。(四)品种平方和的调整直接按格子设计进行测验，则要对品种平方和进行调整，对于简单格子设计，其矫正数为：

63

(14·32)其中，Ku为未调整的成分(b)平方和，Kb为调整的成分(b)平方和。Kb由(14·25)计算，表14.3中的Ku可由下式计算：

(14·33)(14·33)64

表14.5简单格子设计方差分析表变异来源DF重复r-1区组(调整的)r(p-1)2(p-1)2(p-1)品种(未调整的)p2-1区组内误差(Ei)(p-1)(rp-p-1)总rp2-1变异来源DF重复65(五)期望均方简单格子设计用于单因素试验，其期望均方和随机区组的情况一样，区组内误差估计了，调整的品种均方估计了(随机模型)或(固定模型)。二、简单格子设计的例题(一)二次重复简单格子设计的例题(五)期望均方66[例14.1]表14.6为一个5×5大豆品种重复二次简单格子设计的试验结果。其田间排列是随机的。随机的步骤：①在每一重复内分别独立地随机安排区组；②在每一区组内分别独立地随机安排品种代号；③将各品种随机决定品种代号。[例14.1]表14.6为一个5×5大豆品种重复二次简67表14.65×5大豆品种简单格子设计的产量试验结果(r=2，kg/区)第十四章不完全区组设计和统计分析课件68分析步骤如下：1.从表14.6计算各区组总和(这里即Xe·及Y·f)，重复总和(这里即X··及Y··)各品种(未调整)总和(tef)以及Te·、T·f值。并按随机区组进行方差分析。结果列于表14.7。随机区组方差分析结果品种间无显著差异。进一步再按格子设计分析。分析步骤如下：69表14.7随机区组方差分析表2.计算消去品种效应的区组平方和。由成分(b)单独估计。按(14·25)，r=2时为：变异来源DFSSMSF重复1212.18品种24559.2823.30＜1误差24720.3230.01总491491.78表14.7随机区组方差分析表变异来源DFSSMSF重70在表14.6上分别计算Te·-2Xe·及T·f-2Y·f值，代进上式得：3.列出分解有区组变异的方差分析表(表14.8)。=501.84表14.85×5简单格子设计(r=2)方差分析表变异来源DFSSMSF重复1212.18品种(未调整)24559.2823.30重复内区组(调整)8501.8462.73(Eb)4.59**区组内误差16218.4813.66(Ei)总491491.78在表14.6上分别计算Te·-2Xe·及T·f-2Y·f值71调整后重复内区组间的变异很显著，说明将区组划出是很必要的。4.计算调整的品种总和()。由(14·23)，在简单格子设计两个重复时：=0.7820=0.1564

调整后重复内区组间的变异很显著，说明将区组划出是很必要的。=72调整品种总和=在表14.6中分别计算及然后计算各品种调整的总和，以品种(1)为例：=30+9.5-1.4=38.1。其余类推，全部结果列于表14.6的末端。5.计算品种平均数间比较的误差。同区组品种平均数间比较：调整品种总和=73

异区组品种平均数间比较：全试验品种平均数相互比较：一般用2.93作标准误进行品种间比较即可。异区组品种平均数间比较：全试验品种平均数相互比较：一般用2746.计算调整的品种平方和再进一步测验品种差异的显著性.按(14·32)品种平方和的矫正数为:其中Ku仿(14·33)为:6.计算调整的品种平方和再进一步测验品种差异的显著性.其中75

Kb为调整的区组成分(b)平方和，即表14.9中的501.84。

w=1/Ei=1/13.66=0.073211/(2Eb-Ei)=1/(2×62.73-13.66)=0.008945=559.28+85.30=644.58故调整品种平方和Kb为调整的区组成分(b)平方和，即表14.9中的501.76

调整的品种均方及F测验如下：按照简单格子设计的分析结果调整以后的品种均方比未调整时增大了，误差比随机区组时降低了，因而提高了试验的精确性。它与随机区组设计相比较，所提高的效率可估计如下：

变异来源自由度平方和均方F品种(调整)24644.5826.861.97区组内误差16218.4813.66调整的品种均方及F测验如下：按照简单格子设计的分析结果77即提高了74％。本试验品种间无显著差异，所以不必进一步再做品种平均数间的比较。第十四章不完全区组设计和统计分析课件78(二)四次重复简单格子设计的例题[例14.2]上例5×5大豆试验，原为一个四次重复的简单格子设计，若表14.6中的是第一重复及第三重复，今将第二重复，第四重复的结果补充列在表14.9中，重复Ⅱ与重复Ⅰ属同一种分组，重复Ⅳ与重复Ⅲ属另一种分组。分析步骤如下：1.从表14.6及14.9计算各重复各区组的总和g，重复总和G，同品种的两个区组总和Xe·及Y·f，各品种

(二)四次重复简单格子设计的例题79表14.95×5大豆品种简单格子设计Ⅱ、Ⅳ重复的产量结果(r=2，kg/区)第十四章不完全区组设计和统计分析课件80(未调整)总和tef以及Te·、T·f值。按随机区组预先进行方差分析(表14.10)。随机区组方差分析结果品种间无显著差异，进一步按格子设计分析。表14.10随机区组方差分析表变异来源DFSSMSF重复r-1=4-1=3226.19品种p2-1=25-1=24791.2432.961.53误差(r-1)(p2-1)=721547.5621.49总992564.99(未调整)总和tef以及Te·、T·f值。变异来源DFS812.计算消去品种效应的区组平方和。这里包括成分(a)及成分(b)两部分。成分(a)的计算：2.计算消去品种效应的区组平方和。82成分(a)的另一种计算方法可适用于更多次重复的分析。即由相同分组方法内品种组与二次重复的交互作用项计算。

区组平方和(区组总SS)成分(a)的另一种计算方法可适用于更多次重复的分析。83

重复间平方和(重复SS)品种组间平方和(品种组SS)重复间平方和(重复SS)品种组间平方和(品种组SS)84

成分(a)=区组总SS-重复SS-品种组SS=602.18-309.28-128.14=164.72计算结果与前相同。成分(b)r=4时，为：3.列出分解有区组变异的方差分析表(表14.11)。成分(a)=区组总SS-重复SS-品种组SS=60285表14.115×5简单格子设计（r=4）方差分析表4.计算调整的品种总和。

变异来源DFSSMS重复r-1=3226.19品种(未调整)p2-1=24791.2432.96重复内区组间r(p-1)=16786.0049.12(Eb)2(p-1)=82(p-1)=8164.72621.28区组内误差(p-1)(rp-p-1)=56761.5613.60(Ei)总rp2-1=992564.99表14.115×5简单格子设计（r=4）方差分析表4.86

调整品种总和调整品种总和87

在表14.9中分别计算出及然后计算各品种调整的总和，方法同上例。如品种15=72+(+8.8)+(-61)=74.7，余类推。全部计算结果列于表14.9的末端。5.计算品种平均数间比较的误差。同区组品种在表14.9中分别计算出及然后计算各品种调整的总和，方法同88异区组品种

全试验品种6.计算调整品种平方和并进一步测验品种差异的显著性。异区组品种全试验品种6.计算调整品种平89

此即计算成分(a)时的品种组间平方和一项。

调整品种平方和

=791.24+154.33=945.57此即计算成分(a)时的品种组间平方和一项。90调整的品种均方及F测验如下：按格子设计分析，扣除了重复内区组间的变异，降低了试验误差，使品种间的变异呈现出显著性。7.进一步可以计算出调整的平均数，并由全试验品种SE计算LSD进行品种间的比较。方法同随机区组，此处从略。变异来源自由度平方和均方F品种(调整)24945.5739.402.90**区组内误差56761.5613.60调整的品种均方及F测验如下：变异来源自由度平方和均方F品种91第四节平衡不完全区组设计的统计分析[例14.3]设若对某种水果7个品种进行风味品尝，请7位专家评分，每位专家按图14.7的计划鉴评3个品种，其第1号为对照品种，评分范围为最低0分，最高5分，结果列于表14.12。该试验具有处理数t=7，区组数k=3，重复数r=k=3，两两品种在同一区组相遇1次。这一设计的线性模型为:

(14·34)

第四节平衡不完全区组设计的统计分析[例14.3]设若92表14.12七个品种风味的专家评分结果(平衡不完全区组设计)区组(专家)品种与评分yij区组总和B(1)①3.5②3.8④4.111.4(2)②3.4③4.0⑤3.310.7(3)③4.1④4.3⑥4.613.0(4)④4.3⑤4.2⑦4.613.1(5)⑤3.7⑥4.6①3.912.2(6)⑥4.0⑦4.8②3.712.5(7)⑦4.9①4.0③4.513.4G=86.3区组(专家)品种与评分yij区组总和B(1)①3.5②93其分析步骤如下：1.在表14.12中计算未调整的区组总和(B)及全试验总和(G)。计算未调整的品种总和(Tt)列于表14.13；同时计算出品种所在区组各区组总和的和数(Bt)，如品种1为11.4+12.2+13.4=37.0等，列于表14.13。应与kG相等，可用以验算数据。2.计算各品种的W值。W=(t-k)T-(t-1)Bt+(k-1)G=4T-6Bt+2G(本例情况)。按(14·34)将各小区的线性组成相加、减，可以发现不同品种的W值只包含区组效应，因而W值间的变异表示了调整后区组间的变异，其总和ΣW应为0。其分析步骤如下：94表14.13平衡不完全区组设计数据分析表调整处理平均数品种TtBtW调整处理总和Tc=Tt+wW1(CK)11.437.0-3.811.263.75210.934.68.611.223.74312.637.10.412.614.20412.737.5-1.612.644.21511.236.01.411.253.75613.237.7-0.813.174.39714.339.0-4.214.144.7186.3258.90.086.29调整处理平均品种TtBtW调整处理总和Tc=Tt+wW953.进行方差分析。全试验21个小区的总变异中包含有品种间纯变异、区组间纯变异、由于区组不完全而导致的品种与区组相混杂的一部分变异、以及区组内的误差四部分。其中品种与区组相混杂的一部分变异包含在处理总和(T)间的变异中，也包含在区组总和(B)间的变异中。因混杂的这一部分变异不论在前者还是在后者是同一个成分，因此在方差分析中只须考虑一个方面便可。由W值计算调整的区组间平方和的公式为：3.进行方差分析。96

(14·35)本例中为·3·4·2)=0.6629未调整的品种平方和:全试验总平方和区组内平方和=4.0981-3.0114-0.6629=0.4238

97表14.14平衡不完全区组设计的方差分析表

此处所获的Ee，实际上只是一个初步估计值，并不立即用于进行F测验，而需作进一步调整。4.计算加权因子w，并调整处理总和及平方和。

(14·36)

变异来源DFSSMS调整MSF品种(未调整)t

-1=63.01140.5020.42787.035**

F0.05=3.58,F0.01=6.37区组(已调整)b-1=60.66290.110(Eb)区组内误差

tr-2t+1=80.42380.053(Ee)0.0608总tr-1=204.0981表14.14平衡不完全区组设计的方差分析表变异来源98

按(Tt+wW)计算调整的品种总和(Tc)，如品种1(CK)为11.4+(-3.8)(0.0370)=11.26等，填入表14.13。

(14·37)

本例中本例中相应的均方为2.5665/6=0.4278。本例中本例中相应的均方为2.5665/6=0.4278995.计算有效误差并作进一步方差分析。有效误差E=Ee[1+(t-k)w]本例中E=0.053[1+(7-3)0.0370]=0.0608将调整的品种均方和有效误差填入表14.14右端，这时可进行F测验。F测验的结果表明品种间风味评价上有很显著的差异。必须说明平衡不完全区组设计的方差分析中根据加权因子w调整的处理均方和误差均方都是近似的，包括w值本身也有抽样波动，所以这一F测验也是一种近似的测验。5.计算有效误差并作进一步方差分析。1006.处理间的比较。处理平均数间比较可用LSD法，此例中已经F测验证实品种间有显著差异，故实际上已用了Fisher保护最小显著差数法(FPLSD)。FPLSD0.05=测验结果如下：6.处理间的比较。FPLSD0.05=测验结果如下：101

品种2153467评分3.743.753.754.024.214.394.71显著性比较结果，品种2、5、3与对照间无显著差异，品种4、6、7的风味评价均优于对照，尤其品种7最佳，优于品种3、4。品种2153467评分3.743.753.754.02102第十四章不完全区组设计和统计分析第一节不完全区组设计的主要类型第二节重复内分组和分组内重复设计的统计分析第三节简单格子设计的统计分析第四节平衡不完全区组设计的统计分析第十四章不完全区组设计和统计分析第一节不完全区组设103第一节不完全区组设计的主要类型一、田间试验常用设计的归类二、重复内分组和分组内重复设计三、格子设计四、平衡不完全区组设计第一节不完全区组设计的主要类型一、田间试验常用设计的归类104一、田间试验常用设计的归类完全区组(completeblock)：每一区组包含全套处理。不完全区组(incompleteblock)：即一套处理分成几个区组，或一个区组并不包含全部处理，但同样要通过区组实施地区控制。

一、田间试验常用设计的归类完全区组(completeblo105二、重复内分组和分组内重复设计重复内分组设计(blockinreplication)：将供试品种分为几个组，看作为主区，每个组内包含的各个品种看作为副区，重复若干次，主副区都按随机区组布置的设计。例如20个品种，分为4组，每组包含5个品种，若重复3次，则田间布置可设计如下图：

二、重复内分组和分组内重复设计重复内分组设计(blocki106

重复内分组设计的田间布置该例中重复内分组设计的自由度分析如下：重复Ⅰ重复Ⅱ重复Ⅲ区组(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)420111017751591912331815816621381713121913918811271615251612620103146201441171471994111018115重复Ⅰ重复Ⅱ重复Ⅲ区组(1)(2)(3)(4)(5)(6)(107变异来源DF重复2组间3误差(Ea)6组内品种间16误差(Eb)32总59组内品种间比较的误差将为：；变异来源DF108各组平均数间比较的误差将为：；不同组品种间比较的误差(仿照裂区的情况)将为：。由于Ea与Eb常取不同数值，Ea往往大于Eb，例如

=3，若如此，则：组内品种间比较的误差将为：不同组品种间比较的误差将为：各组平均数间比较的误差将为：109两者比值为：即不同组品种间比较的方差将比组内品种间比较的方差大40%，因而像这种不完全区组设计的方法，并不能保证任何两个品种间比较具有相近的精确度。分组内重复设计(replicationinblock)：将供试材料分组后放在连片土地上的几组随机区组试验，通过土地连片而进行联合分析与比较。

两者比值为：110

分组内重复设计

分组1分组2分组3分组4区组(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)191618131511541898181917121114352710716201915121313510692017161114152149710171820141312423686分组1分组2分组3分组4区组(1)(2)(3)(4)(5)(111三、格子设计格子设计(latticedesign)：为了克服重复内分组设计中组间品种比较和组内品种比较精确度悬殊的问题，对品种分组的方法可考虑从固定的分组改进为不固定的分组，使一个品种有机会和许多其他品种，甚至其他各个品种都在同一区组中相遇过。三、格子设计格子设计(latticedesign)：为了112(一)格子设计的类别平方格子设计(squaredlattice)：供试品种数为区组内品种数的平方，区组内品种数为p，供试品种数为p2；立方格子设计(cubiclattice)：供试品种数为区组内品种数的立方，区组内品种数为p，供试品种数为p3；矩形格子设计：区组内品种数为p，供试品种数为p(p+1)。

(一)格子设计的类别113(二)平方格子设计1.仿照随机区组式的设计按品种分组方法的变换次数有：(1)简单格子设计(simplelattice)品种分组方法为二种，试验重复次数为2或2的倍数。

重复I重复Ⅱ(1)123(4)147区组(2)456(5)258(3)789(6)369(二)平方格子设计重复I重复Ⅱ(1)123(4)114(2)三重格子设计(triplelattice)：品种分组方法为三种，即在简单格子设计二种分组方法的基础上再增加对角线分组一种，重复次数为3或3的倍数。

(3)四重格子设计(quadruplelattice)：在三重格子设计的基础上，再增加对角线一组，

重复I重复Ⅱ重复III(1)123(4)147(7)159区组(2)456(5)258(8)267(3)789(6)369(9)348(2)三重格子设计(triplelattice)：品种分115

(4)平衡格子设计(balancedlattice)：品种分组方法增加到使每一对品种都能在同一区组中相遇一次。分组法X分组法Y分组法Z分组法L区组(1)12345(6)16111621(11)17131925(16)18151724(2)678910(7)27121722(12)28142021(17)29111825(3)1112131415(8)38131823(13)39151622(18)310121921(4)1617181920(9)49141924(14)410111723(19)46132022(5)2122232425(10)510152025(15)56121824(20)571416235×5四重格子设计方法分组法X分组法Y分组法Z分组法L区组(1)123116

2.仿照拉丁方的格子设计(1)平衡格子方设计(balancedlatticesquare)

①重复数r=(p+1)/2，每对品种在行或列区组中共相遇一次；

重复Ⅰ重复Ⅱ重复Ⅲ重复Ⅳ(1)123(4)147(7)159(10)168区组(2)456(5)258(8)267(11)249(3)789(6)369(9)348(12)3573×3平衡格子设计重复Ⅰ重复Ⅱ重复Ⅲ重复Ⅳ(1)123(4)14117

3×3平衡格子方设计[在行或列中相遇一次，r=(p+1)/2]ⅠⅡ123168456924789573ⅠⅡ123168456924118②重复数r=(p+1)，每对品种在行及列区组中均相遇一次，亦即共相遇二次。

ⅠⅡⅢ159131234111166261014658712251537111511129101483948121616151413713104ⅣⅤ17121411015882131192716101635136312159645141144×4平衡格子方设计[在行及列中共相遇二次，r=(p+1)]

②重复数r=(p+1)，每对品种在行及列区组中均相遇一次，亦119(2)部分平衡格子方设计(partiallybalancedlatticesquare)：重复次数少于最小平衡重复数。与三重、四重格子设计类似，不一定每一对品种都在行或列区组中相遇。格子设计的优点是：考虑了供试品种间平衡比较的问题。但由于供试品种数多，这常只能实施部分平衡，而事实上很难实施完全平衡，因为完全平衡所需的重复次数导致试验规模过大。(2)部分平衡格子方设计(partiallybalanc120育种工作中产量比较在早、中期阶段，因供试材料多需要考虑适合大量处理的设计，但这时每份材料的种子数少，一般不可能进行小区较大的精确试验，因而实际应用中部分平衡的格子设计已可满足要求。育种工作中产量比较在早、中期阶段，因供试材料多需要考虑适合大121四、平衡不完全区组设计平衡不完全区组设计(balancedincompleteblockdesign)：设计的供试处理数不多，不须按格子设计那样每一重复包含有区组大小为k的k个区组，而可将各重复寓于全部区组之中，区组数与区组大小不一定相等，即全试验包括大小为k的区组共t(处理数)或t倍个。

四、平衡不完全区组设计平衡不完全区组设计(balanced122图14.7一种平衡不完全区组设计例如品尝试验，对于一个人的味觉来说，品尝的对象增加太多时鉴别差异的灵敏度便下降，因而每个人只能品尝一部分。图14.7的情况，若有7个水果品种供鉴评，每人品尝3个，请7位品尝家作鉴评，便共品尝21次，每个品种品尝3次。此处每位专家区组（1）（2）（3）（4）（5）（6）（7）123456723456714567123区组（1）（2）（3）（4）（5）（6）（7）1234567123便是一个区组，每区组包含3个品种。这时尽管每人并未将7个品种全部鉴评过，但因是均衡的，每个品种至少和其他6个品种比较过1次。这一试验可增加至14位专家则每对品种相遇2次，21位专家则相遇3次。因而可以请许多专家作出综合评判。便是一个区组，每区组包含3个品种。这时尽管每人并未将7个124第二节重复内分组和分组内重复设计的统计分析一、重复内分组设计的统计分析二、分组内重复设计的统计分析第二节重复内分组和分组内重复设计的统计分析一、重复内分组125一、重复内分组设计的统计分析重复内分组用于品种(系)试验时有二种情况：一是大量品种(系)间的比较目的在于选拔高产优系（固定模型试验）；另一是从一个群体内随机抽出大量家系进行试验，通过供试的样本推论总体的情况（随机模型试验）。

一、重复内分组设计的统计分析重复内分组用于品种(系)试验时有126假定重复内分组设计的供试品种为m=a×b个，分a组，每组有b个品种(系)，重复r次，则重复内分组设计的线