随机区组设计试验结果统计分析.docx

1、随机区组设计试验结果统计分析第9章 随机区组设计试验结果统计分析第一节 单因素随机区组设计试验结果统计分析随机区组试验设计是把试验各处理随机排列在一个区组中，区组内条件基本上是一致的，区组间可以有适当的差异，它相当于一个试验因素，可以从误差中扣除出来，因而降低试验误差以提高试验结果分析的精确性。分析这类资料时，可应用两向分组资料的方差分析方法进行分析。在单因素随机区组试验结果分析时，处理看作A因素，区组看作B因素，其剩余部分则为试验误差。设有A个处理，n个区组，则其平方和与自由度分解 =+ (9-1)总平方和=区组平方和+处理平方和+试验误差平方和 nk1=（n1）+（k1）+（n1）（k1）

2、 (9-2)总自由度=区组自由度+处理自由度+误差项自由度 例9.1有一包括A、B、C、D、E、F、G7个小麦品种的品种比较试验，G为对照品种，随机区组设计，重复3次，小区计产面积30m2，其产量(kg30m2)结果见图9-1，试作分析。1、资料整理2、平方和及自由度的分解根据9-1式和9-2式计算平方和及自由度。总平方和 总自由度 dfT=nk-1处理平方和 处理自由度 dfs=k-1区组平方和 区组自由度 dfb=n-1误差平方和 误差自由度 dfe=(k-1)(n-1)平方和计算如下：C=T2/nk=341.32/(37)=5546.94总=14.02+16.12+17.12C=137.

3、5区组 =(103.22+112.92+125.22)/7=34.73品种 =(49.82+51.92+45.22)/3-C=66.60=137.5-34.73-66.60=35.92将以上结果填入表9-2，并将自由度计算结果直接填入表9-2。表9-2 资料的方差分析表变异来源自由度平方和方差F值F0.05F0.01区组间品种间误差261234.7366.6035.9217.3711.102.995.81*3.71*3.893.006.934.82总和20137.253、F测验 列方差分析表，并算得各变异来源的s2值，进而进行F测验。对区组间MS作F测验，结果表明3个区组间的土壤肥力有显著差异

4、，区组控制是有效的，如果不进行区蛆设计，没有局部控制，误差就大。在试验中，区组作为减少误差的手段，一般不作F测验。对品种间MS作F测验，结果表明7个总体平均数间有显著差异。需进一步作多重比较，以明了哪些品种间有显著差异，哪些品种间没有显著差异。 4、多重比较(1)最小显著差数法(LSD法) 根据品种比较试验要求，各供试品种应分别与对照品种进行比较，宜应用LSD法。首先应算得品种间平均数差数的标准误。在利用LSD法进行比较时，平均数差数标准误为：从而得到各品种与对照品种(G)的差数及显著性，并列于表9-3。表9-3 图9-1资料各品种与对照产量相比较差异显著性 品种小区平均产量与G(CK)差异E

5、19.84.7*B17.32.2A16.61.5C15.90.8F15.10G（CK15.10D13.9-1.2从表9-3可以看出，仅有品种正比对照品种增产极显著，其余品种与对照比较均没有显著差异。(2)最小显著相差法(LSR法) 如果不仅要测验品种和对照相比的差异显著性，而且要测验品种间相互比较的差异显著性，则应该应用LSR法。用这种方法比较，首先应算得品种的标准误SE。其标准误为查SSR值表，当v=12时得k=2、37的SSRa值，并根据公式LSRa=SESSRa，算得LSRa值列于表9-4，然后用字母标记法以表9-4的LSRa衡量不同品种间产量差异显著性将比较结果列于表9-5。表9-4

6、图9-1资料最小显著极差法测验LSRa值234567SSR0.053.083.233.333.363.403.42SSR0.014.324.554.684.764.844.92LSR0.053.083.233.333.363.403.42LSR0.014.324.554.684.764.844.92表9-5 图9-1资料最小显著极差法测验结果品种差异显著性5%1%E19.8aA*B17.3 abABA16.6 abcABC15.9 bcABF15.1 bcABG（CK15.1 bcABD13.9 c B结果表明：E品种与C、F、G、D4个品种和B与D品种达到5水平的显著差异；E品种与D品种达到

7、1水平的显著差异，其余品种间均无显著差异。第二节 多因素随机区组设计试验结果统计分析有两个以上试验因素的试验称为多因素试验。这里重点说明两因素随机区组试验结果的统计分析方法。 设有A和B两个因素，各具有a和b个水平，则有ab个处理组合(处理)。采用随机区组设计，重复r次，共有rab个观察值。由于处理项是由A和B两个因素不同水平的组合，因此处理间差异又可分解为A因素水平间差异(A)、B因素水平间差异(B)和A与B的交互作用(AB)三部分。 总平方和SST=SSr区组平方和+SSt处理平方和+SSe误差项平方和 （9-3） (rab1)=(r1)十(ab1)十(r1)(ab1) （9-4）总自由度

8、DFT=DFr区组间自由度+DFt处理自由度+DFe误差自由度其中处理项平方和及自由度可进一步分解：SSt处理平方和=SSA +SSB+SSAB （9-5）(ab1)= (a1)+ (b1)+ (a1) (b1) （9-6）DFt处理自由度=DFA+DFB+DFAB例9 2现对苹果新品种进行氮肥经济用量试验。苹果新品种为A，有3个水平，A1为甲品种，A2为乙品种，A3为丙品种，氮肥用量为B，有4个水平，B1为不增氮肥，B2低氮(增氮低水平)，B3为中氮，B4为高氮，全试验共有12个处理组合，随机区组设计，采用单株小区，重复3次，调查第三年果实产量(kg株)如表9-7，试作分析。表9-7 苹果新

9、品种施肥试验区组与处理两向表处理区组TtA1B140413912040A1B243444212943A1B346474413745.67A1B443424613143.67A2B136363811036.67A2B248444213444.67A2B344494914247.33A2B446414012742.33A3B130343810234A3B240425013244A3B364526017658.67A3B444443612441.33Tr524516524T=1564=43.441、资料整理 将调查结果填入表9-7中，计算出各区组总和Tr，各处理总和Tt及平均数。然后根据各处理总和数

10、作品种和施肥量两向分组整理成表9-8。表9-8 表9-7资料品种（A)和施肥量(日)两向表 品种施肥量B1B2B3B4TAA112012913713151743.08A211013414212751342.75A310213217612453444.5TB332395455382T=156436.88943.88950.55642.444 2、平方和与自由度的分解根据表9-6计算各变异来源的平方和及自由度，平方和计算如下：C=T2/rab=15642/(334)=67947.11总=402+412+362C=1500.89区组 =(5242+5162+5242)/（34）=3.56 =(120

11、2+1292+1242)/3-C=1219.56=1500.89-3.56-1219.56=277.77将以上结果填人表9-9中，井将自由度计算直接填入表9-9。3、F测验 列方差分析表9-9，井计算各项MS值，进而进行F测验。表9-9 苹果品种与施氮量二因素试验的方差分析变异来源DFSSMSFF0.05F0.01区组间23.561.781处理间111219.56110.878.783.263.18品 种220.7210.3613.445.72施肥量3852.67284.2222.503.004.82品种施氮量6346.1757.704.572.553.75误 差22277.7712.63总变

12、异351500.89F测验结果说明：区组间、品种间差异不显著，而处理间、施氮量间、品种施氮量间的差异极显著。由此说明，不同的施肥用量对苹果产量影响不同，而不同苹果品种对氮肥用量又有不同的要求，所以需进一步作氮肥用量间与品种施氮量间差异显著性测验。处理间差异体现在各个试验因素和因素的交互作用，一般不必另行对处理间进行多重比较。4、多重比较(1)施肥量间比较 以各施氮量的小区平均数进行最小显著极差法（LSR法）测验查SSR值表，当v=22，k=2、3、4时的SSR值，并根据公式LSR=SESSR计算各LSR值列于表9-10。根据表9-10的LSR值，对四种氮施量间差异显著性进行测验，其结果列于表9

13、-11。表9-10 不同施氮量间比较的LSR值k234SSR0.05293308317SSR0.01309417428LSR0.05346363374LSR0.01471492505表9-11 不同施氮量间平均产量的差异显著性施氮量差异显著性5%1%B350.56aAB243.81bBB442.44bBB136.89cC 表9-11说明，三种施氮量处理的品种产量均极显著高于不施氮的产量；在增施氮肥的过程中，中氮肥的产量极显著高于低氮量和高氮量的产量，高氮量和低氮量之间的产量差异不显著。 (2)品种施氮量的互作 F测验说明，品种与施氮量的交互作用极显著，表明不同品种适应不同的施氮肥量。进行这种相

14、互比较，可以比较全部处理(水平组合)的平均数；也可将每一个品种的四种施氮肥用量之间进行比较。现将这两种比较方法分别介绍如下：各品种不同施氮肥量间的比较 查SSR值表，当v=22时，得k=2、3、4的SSR值，根据公式LSR=SESSR计算出各LSR值，并列于表9-12，用此标准对每一个品种的四种施氮用量之间进行比较，其结果列于表9-13。表9-12 表9-7资料LSR值表k234567891012SSR0.052933083173.243.293.323.353.373.393.42SSR0.013094174284.364.424.484.534.574.604.65LSR0.056.016

15、.136.506.646.746.816.876.916.957.01LSR0.018.188.558.778.949.069.189.299.379.439.53表9-13 各品种在不同施氮量的差异显著性A1品种 施氮量差异显著性5%1%B345.67aAB443.63aAB243.00aAB140aAA2品种 施氮量差异显著性5%1%B347.33aAB244.67aAB442.33aAB134BBA3品种 施氮量差异显著性5%1%B358.67aAB244.00bBB441.33bBCB134cC 表9-13结果表明：Al品种对施氮肥量的大小不敏感，四种施氮肥量间无显著差异。A2品种对施

16、氮肥量比较敏感，施氮肥量比不施增产效果极显著，但用量大小之间无显著差异。A3品种对施氮肥量的大小很敏感，三种施氮肥量结果其产量均高于不施氮肥的产量，并且达显著差异，中施氮肥量的苹果产量极显著高于中低、高施氮量的产量。处理间比较 将各处理平均数列入表9-14，并以表9-12的LSR值衡量全部处理平均数间差异的显著性。 表9-14 各处理组合闸的差异显著性测验(LSR法)处理差异显著性5%1%A3B358.67aAA2B347.33 bBA1B345.67 bcBCA2B244.67 bc BCA3B244.00 bcBCA1B443.67 bcBCA1B243.00 bcd BCDA2B442.

17、33 bcd BCDA3B441.33 bcd BCDA1B140.00 cde BCDA2B136.67 de CDA3B134.00 e D 表9-14测验结果表明，A3品种(丙品种)在中施氮量水平条件下产量最高，为最优组合，该处理组合极显著优于其它各处理组合的产量。三个苹果品种不施氮肥产量最低，显著的低于其它施氮肥产量。3、试验结论 本试验品种主效无显著差异；施氮肥量主效有极显著差异，以B3施氮量产量最高，与B1、B2、B4达极显著差异，B2、B4之间无差异，但与B1(不施氮肥)达极显著差异。品种与施氮肥量的互作极显著，以A3品种与B3施氮量结合，才能取得最高产量，A2、A1品种也需B3

18、施氮量，但不如A3品种产量高，3个品种以施氮肥比不施氮肥产量高。第三节 随机区组试验缺区估计在试验过程中，由于某些难以控制的因素影响，如病虫害或其它意外损失造成植物生长不正常，甚至缺株或缺区。在这种情况下，处理和区组的均衡性遭到破坏，方差分析不能按原定系统进行。为解决这一矛盾，可根据一定的统计原理，估算出缺区数据，并使之代替缺值参加分析。一个试验中，在小区面积较大(或株数较多)而缺株不多的情况下，可根据该小的生长结果的植株数，折合成全小区估计数值。如果缺株过多，须将整个小区剔出而另外估计其试验数据。缺区估计是不得已的补救办法，如果缺区较多，则缺区估计不可靠，应作试验失败对待或者除去缺区过多的处

19、理或区组再作分析。在随机区组试验缺区数值，一般用下式估计： (9-7)式中为缺区理论估计值，n为区组数(重复数)k为处理数，为缺区所在的不包括缺区数值在内的处理的总和，为缺区所在，但不包括缺区数值在内的区组总和，为缺区除外的全试验的总和。一、随机区组试验缺失一个数据的估计 例9.3某水稻品比试验，随机区组设计，缺失一区产量，其结果如表9-15，试作估计。表9-15 水稻品比产量(kg/20m2) 品比区组Tt冀粳8号20.317.616.254.1越富17.117.415.950.472719.615.714.449.7涿17918.817.719.255.7681119.115.334.4+

20、中百4号21.221.819.362.3 花764919.318.617.955.8花823221.420.416.658.4156.8129.2+134.8420.8将9-15有关数值代入公式(9-7)，求得缺区估计数值为：将填于表915的的位置，即可进行方差分析。二、试验缺失两个小区数据的估计对于缺失对于缺失两个小区数据的资料，仍可用(97)公式，但需建立二元一次联立方程才能得到估计值。例9.4有一大豆栽培试验，假定缺失两个小区产量，其产量如表9-16，试进行估计。表9-16大豆随机区组试验缺失产量区试验结果(kg/20m2)处理区组TtA8A8925+AB910101140C111113

21、1247D139D1234+DE1091294051394353186根据公式（97）对A 有方程A对D 有方程D整理成二元一次联立方程得：解之得：=6.88 ，=12.44 将=6.88 ，=12.44 填入表916中，即可进行方差分析。三、缺失数据资料的分析与不缺区数据资料分析的区别对于缺区资料的统计分析，与没有缺区资料分析的步骤相同，但在具体分析时其总自由度和误差项自由度应减去缺区数目，如果缺失一个小区数值，误差项和总变异自由度比常规的少1个，如果缺失两个小区数值，其误差项和总变异自由度比常规的少2个。在进行处理间比较时，一般用t试验，对于非缺区处理问的比较，其标准误仍用公式b，毛在缺失一个小区数据的资料分析时，缺区处理与非缺区处理相比较其标准误公式为： （9-8）上式中MSe为误差项均方，n为区组数，k为处理数。在缺失两个小区数据的资料分析时，缺区处理与非缺区处理其标准误公式为： （9-9） 式中为MSe为误差项均方，n1和n2分别表示两个相比较处理的有效重复数。其计算方法是在同一区组内若两个处理都不缺区，则各记为1；若一处理缺区，另一处理不缺区，则缺区处理记0，不缺区处理记(k2)(k1)，其中k为试验处理数目。 例如在表916资料中，A和D比较时 A的有效重复数D的有效重复数故