第十章一年多点试验资料的方差分析.ppt

第十章 一年多点试验资料的方差分析,一、品种多点试验资料的方差分析,设有v个品种，在L个地点做比较试验每个地点皆设r个重复，按随机区组设计进行试验，则第i个品种(i=1, 2, …, v)在第j个地点(j=1, 2, …, u)，第k区组(k=1, 2, …, r)的观测值为xijk，如下表： 表10－1 一年多点试验数据表,它的数学模型为： xijk＝μ＋ti＋Lj＋(tv)ij＋rjk＋eijk 式中μ为群体的平均值，ti为品种i的效应，Lj为地点j的效应，(tv)ij为品种×地点互作效应，rjk为地点内的区组效应，eijk为随机误差由此，可以得到一年多点区域试验的方差分析表(表10-2)表10-2 一年多点试验资料的方差分析,二、品种多点试验结果统计分析示例,设有一个早稻品种多点试验，供试品种四个（V＝4），以V1、V2, V3, V4表示，其中V4品种为对照，三次重复（r＝3），以I、II、III表示，随机区组试验设计，分别在五个试验点（L＝5）同时进行，以L1, L2, L3, L4, L5表示，小区面积为100m2，试验小区产量结果（kg）列于表10-31．各试验点品种比较试验的方差分析,表10-3 早稻5点试验各试验点小区产量（kg）,(1) L1试验点品种比较试验的方差分析,分别对表10-3早稻多点试验各试验点小区产量结果（kg）进行方差分析，计算出各试验点相应的平方和、自由度和均方。

矫正数C＝T2/Vr=2042/4×3=3468 总平方和SST=x2－C=（172＋102＋……＋122）－C＝438 区组平方和SSr=Tr2－C＝(602+682+762)/4-C=32 品种平方和SSv =Tv2－C ＝(692+332+572+452)/3－C＝240,误差平方和SSe＝438－32－240＝166 总dfT＝Vr－1＝4×3－1＝11 区组 dfr＝r－1=3－1=2 品种dfv=V－1=4－1=3 误差dfe＝(v－1)(r－1)＝(4－1)(3－1)＝6 区组均方MSr＝32/2=16 品种均方MSv＝240/3=80 误差均方MSe＝166/6=27.67,(2)其他试验点品种比较试验的方差分析,同理分别对L2、L3、L4、L5试验点进行方差分析（具体计算方法同L1）各试验点的自由度都相同，均方值只需将相应的平方和除以自由度，其各试验点方差分析结果列于表10－4表10-4 各试验点方差分析结果,2．各试验点误差均方同质性测验,对品种多点试验结果进行联合分析时，通常要对各试验点误差均方进行同质性（齐性）测验，只有当各试验点误差均方差异不显著时，才能将各试验点的试验结果合并分析，否则，不宜合并。

对各试验点的误差均方行同质性测验1）Bartlett 2测验方法,设有L个独立误差均方估计值S12，S22，…，SL2，其相应自由度分别为V1，V2，……VL，那么合并方差S2为：,*,Bartlett 2值为：,如果求出的2值小于查2表的临界02 值， 则说明各误差均方同质，即可将多点试验结果合并进行联合分析；若求得的2值大于查2表的临界02 若值， 则说明各误差均方不同质，需要对试验数据进行适当的数据处理，通常可剔除个别“特殊”的试点，或将原始数据作平方根或对数转换，获得一个同质的方差，再合并进行联合分析当然在对试验结果统计分析要求不太严格的情况下，也可以不进行各试验点误差均方同质性测验，直接将各试验点的结果合并进行联合分析2）Bartlett 2测验的计算,可由表10-4列成表10-5形式进行计算,本例Bartlett 2测验计算（L=5）查2表的自由度＝L－1＝5－1＝4，得0.05显著水平临界值20.05(4)＝9.49 Bartlett 2测验结果，实得2 ＝6.54小于临界值20.05(4) ＝9.49 2测验不显著，可以认为各试验点误差均方同质，可以将各试验点的结果合并进行联合分析。

3. 品种多点试验结果的联合分析,（1）联合分析的平方和与自由度的计算 将5个试验点的试验结果，合并成为表10－3形式，根据表10－2计算各变异来源的平方和与自由度 矫正数C＝6782(4×5×3)＝7661.4 总平方和SST＝1944.60，总自由度dfT＝4×5×3－1＝59 区组平方和SSr＝(602＋682＋……＋332)/4－C= 108.5，区 组dfr＝5×（3－1）＝10 品种与试验点处理组合平方和SSvl＝(692＋332＋……＋332)/3－C＝1422.6,试验点平方和SSL＝(2042＋1592＋……＋1172)/4×3－C＝689.1，dfL＝5－1＝4 品种平方和SSv＝(1682＋1322＋2312＋1472)/5×3－C＝379.8，dfv＝4－1＝3 品种×试验点互作SSvl＝1422.6―689.1―379.8＝353.7，互作df＝12 误差平方和SSe＝1944.6―108.5―1422.6＝413.5，dfe ＝5(4－1)(3－1)＝30,(2）列方差分析表，进行F测验表10－6 早稻品种多点试验方差分析,F测验,品种多点试验的主要目的在于鉴定参试品种的优劣及其适应区域，而对试验点间的产量差异和试验点内区组间的差异不感兴趣，所以在品种多点试验资料联合分析时，只作品种以及品种×试验点互作的F测验。

一般品种多点试验，品种为固定模型，而试验点和区组往往是随机模型，故品种多点试验为混合模型本例按表10-6所列的均方进行F测验方差分析中处理效应的分类： 固定效应：在单因素试验的方差分析中，把k个处理看作k个明晰的总体如果研究的对象只限于这k个总体的结果，而不需推广到其它总体；研究目的在于推断这k个总体平均数是否相同，即在于检验k个总体平均数相等的假设H0：μ1=μ2=…=μk；H0被否定，下步工作在于作多重比较；重复试验时的处理仍为原k个处理这样，则k个处理的效应(如=μi-μ)固定于所试验的处理的范围内，处理效应是固定的随机效应：在单因素试验中，k个处理并非特别指定，而是从更大的总体中随机抽取的k个处理而已，即研究的对象不局限于这k个处理所对应的总体的结果，而是着眼于这k个处理所在的更大的总体；研究的目的不在于推断当前k个处理所属总体平均数是否相同，而是从这k个处理所得结论推断所在更大总体的变异情况，检验的假设一般为处理效应方差等于零，即H0： =0；如果H0被否定，进一步的工作是估计 ；重复试验时，从更大的总体随机抽取新的处理这样，处理效应是随机的按处理效应的类别来划分方差分析的模型，在单因素试验时，有2种，即固定模型和随机模型；在多因素试验时，则有3种，即固定模型、随机模型和混合模型。

若各试验因素水平的效应均属固定，则称之为固定模型一般品种比较试验、肥料试验等均属固定模型若各试验因素水平的效应均属随机，则称之为随机模型随机模型在遗传、育种和生态试验研究方面有广泛的应用例如，为研究中国早稻产量变异情况，从大量早稻品种中随机抽取部分品种为代表进行试验，从试验结果推断中国早稻产量变异情况，这就属于随机模型在多因素试验时，若各试验因素水平的效应既有固定的、也有随机的，则称之为混合模型混合模型在试验研究中是经常采用的例如，进行多年、多点品种区域试验，品种效应、地点效应是固定的，而年份效应是随机的由于模型不同，方差分析中各项期望均方的计算也有所不同，因而F检验时分母项均方的选择也有所不同 就试验资料的具体统计分析过程而言，这三种模型的差别并不太大，但从解释和理论基础而言，它们之间是有很重要的区别的品种间F＝品种均方（品种×试验点互作均方）＝126.629.475=4.295* 查F表，可得F0.05,(3.12) =3.49, F0.01(3,12)=5.95 品种间实得F值为4.295，大于F0.05,(3.12) =3.49，则表明供试品种产量间存在着显著的差异，表明总的说来各品种平均产量间存在着真正的差异。

但是F测验不能具体指出究竟哪几个品种产量间有真正差异，还必须进一步作品种平均产量间的多重比较品种×试验点互作F＝（品种×试验点均方）误差均方＝29.4513.783=2.139* 查F表可得F0.05(12.30=2.09，F0.05(12.30=2.84， 品种×试验点互作实得F值 2.139，大于F0.05＝2.09，则表明品种×试验点互作亦达到显著，说明不同品种在不同的试验条件下的表现存在差异，这种差异主要是品种基因型与环境互作造成的所以当品种×试验点互作显著或极显著时，有必要测定品种的稳定性这里还需说明一点，有些教材或参考资料对品种多点试验（包括品种区域试验）的联合分析，是按固定模型的期望均方进行F测验 按固定模型则品种间F＝品种均方误差均方＝26.613.783=9.185 品种间实得F值达到极显著水平，表明供试品种产量间存在着极显著差异；品种×试验点互作F值的计算采用混合模型3）品种平均产量间的多重比较,品种平均产量间的多重比较采用新复极差测验 计算标准误Sx＝[(品种×试验点均方)(试验点数目×重复次数)]－2 ＝(29.75（5×3）)＝1.40 查新复极差测验5％和1％SSR值表，并计算各个LSR值，列于表10－7。

查SSR值表的自由度为品种×试验点的自由度，本例dfvL＝12，K＝2，3，4表10－7 LSR值的计算,③品种平均数间的多重比较,按品种小区平均产量的高低依次排列，然后逐个侧验各品种平均产量间的差异显著性，测验结果用标记字母法表示于表10－8表明，V3, V1品种间产量差异不显著，V1, V4, V2品种间产量差异亦不显著；本试验只有V3品种的产量显著高于V4对照品种，极显著高于V2品种，其他品种与V4对照品种均无显著差异根据试验结果，建议V3品种进一步扩大试种、示范、推广4. 品种稳定性测定,当品种多点试验联合分析品种×试验点互作达显著或极显著时，就要测验品种的稳定性在品种多点试验（包括品种区域试验）中，品种稳定性主要指产量稳定性，所谓品种产量稳定性，就是指品种在不同的环境条件下，能够保持产量的稳定状态 (1）由表10-1资料计算出每一品种在每一试验点的小区平均产量，列于表10-9表10-9 品种与试验点相应的平均产量和环境指数,环境指数：品种多点试验中将每一个试验点所有参试品种的平均产量，作为该试验点的环境指数2）计算每一品种的平均产量（Y），回归于环境指数（x）的回归系数（b）、回归截距（a）及其直线回归方程式。

现以计算V1品种的b、a和直线回归方程为例，说明计算步骤和方法∑x＝17.0＋13.25＋8.0＋8.5＋ 9.75＝56.5 x＝56.5/5=11.3 ∑x2＝17.02＋13.252＋8.02＋ 8.52＋ 9.752＝695.875 ∑y＝23.0＋13.0＋5.0＋7.0＋8.0＝56.0 y＝56.0/5＝11.2 ∑y2＝23.02＋13.02＋5.02＋7.02＋8.02＝ ∑xy＝(17.0×23.0)＋(13.25×13.0)＋(8.0×5.0)＋(8.5×7.0)＋(9.75×8.0)＝740.75,,,,,(8) a=y-bx=11.2-1.88×11.3＝－10.04 (9) V1品种的直线回归方程式为 ŷ＝－10.04＋1.88x V2品种为：ŷ＝2.84＋0.53x V3品种为： ŷ＝4.99＋0.92x V4品种为： ŷ＝2.21＋0.67x,,,(3) 品种产量稳定性统计数b的解释,可将上述4个参试品种的平均产量和回归系数b列成表10－10形式， 当用回归系数b来测定品种产量的区域稳定性时，若某品种b＝1，表示该品种具有平均稳定性；b1，表示低于平均稳定性，其实际含义是该品种在有利的环境条件下，具有较大的增产潜力。

表10-10 4个早稻品种的小区平均产量和回归系。