多因素试验结果的统计分析.ppt

《多因素试验结果的统计分析.ppt》由会员分享，可在线阅读，更多相关《多因素试验结果的统计分析.ppt（103页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第八章多因素试验结果的统计分析,二因素随机区组设计的结果分析三因素随机区组设计的结果分析二因素裂区设计的结果分析,8.1多因素随机区组试验的统计分析,一、二因素随机区组试验结果的方差分析1、二因素随机区组试验的线性模型和期望均方设有A、B两个试验因素，A因素有a个水平,B因素有b个水平,采用随机区组设计，重复r次。该二因素试验共有ab个水平组合，每一个水平组合有r个观察值。则该试验共有rab个观察值。,在二因素试验中，由于有两个试验因素，其处理效应由三部分构成，即：,故二因素随机区组试验中每一观察值的线性模型为：,式中，h=1,2，r；i=1,2,a；j=1,2，b;,因此，在可加性的假设下，

2、二因素随机区组试验结果的总变异可分解为区组间、处理间和试验误差三部分，而处理又可分解为A因素、B因素和AB互作三个部分。,二因素随机区组试验设计可参照单因素随机区组试验进行，唯一不同点是二因素随机区组试验把各因素不同水平组合当作单因素试验中的处理看待，并按随机的原则排列在各区组。,假定有一个A、B二因素试验，a=3,b=4,随机区组设计，重复两次r=2，该试验共有12个水平组合.,因重复2次，故应先划分为两个区组；又因有12个水平组合，故每区组划分为12个试验小区。,肥瘦,2、二因素随机区组试验的结果分析,单因素随机区组试验,二因素随机区组试验,二因素随机区组设计的平方和与均方,变异期望均方来

3、源固定模型随机模型区组2e+ab22e+ab2A2e+rb2A2e+r2(AB)+rb2AB2e+ra2B2e+r2(AB)+ra2BAB2e+r2(AB)2e+r2AB误差2e2e,二因素随机区组的期望均方,对于多因素试验而言，效应模型的不同将导致F测验的方法不同。对固定模型来说，各变异项的均方除误差均方即构成相应的F测验。但对随机模型来说，区组变异和互作变异用误差均方进行F测验；而A、B的变异则应用互作项的均方进行F测验。,当选用固定模型：,其F值都是以误差项的均方为分母的。,当选用随机模型：,测验,应以误差项均方为分母；,而测验,需以互作项的均方为分母。,二因素随机区组与单因素随机区组的

4、差别：二因素试验的处理项可以再分解为A因素水平间、B因素水平间和AB互作三部分，因此二因素处理项的平方和与自由度亦可作相应的剖分：SSt=SSA+SSB+SSAB(ab-1)=(a-1)+(b-1)+(a-1)(b-1)其中，SSt为处理平方和；SSAB为互作项平方和;SSA为A因素平方和；SSB为B因素平方和。,二、二因素随机区组试验结果的分析实例(固定模型）【例8.1】玉米品种与施肥二因素随机区组试验，A因素有A1，A2，A3(a=3)三个品种，B因素有B1，B2，B3(b=3)三个施肥水平，重复3次(r=3)，小区计产面积20m2，田间排列和小区产量(kg)如图8.1，试作分析。,玉米品

5、种与施肥随机区组试验田间排列和小区产量,(1)结果整理将试验所得结果按处理和区组两向分组整理成表；,资料处理与区组两向表,资料品种(A)与施肥(B)两向表B1B2B3TAA1453828111A2435251146A3532824105TB141118103362,再按品种(A)和施肥(B)作两向分组整理成表。,(2)自由度和平方和的分解在表中，Tr为区组总和，Tt为处理总和，TA为各品种总和，TB为各施肥水平总和，T为全试验总和。,自由度的分解,矫正数C=T2/rab=3622/(333)=4853.48,平方和的分解,误差SSe=SST-SSr-SSt=436.52-35.19-338.5

6、2=62.81,对处理SSt进行再分解可得：,AB互作,B因素,A因素,玉米品种与施肥二因素试验的方差分析（区组随机，处理固定）,（3）方差分析和F测验,这里对A和B两因素都取固定模型，区组则取随机模型，根据表8.1中所列各变异来源EMS的组成情况，各F值的计算都应用误差项MS作分母，取显著水平=0.05。,表8.4的F测验结果表明，区组间、品种间、施肥水平间以及品种与施肥互作间的差异均达5%显著水平，所以应进一步进行差异显著性测验。,(4)差异显著性测验(SSR)计算过程已多次重复叙述，本章只列出各自的标准误、LSR和显著性测验结果，查表计算过程从略。,k=2时，LSR0.05，16=1.9

7、8(kg)，LSR0.01，16=2.73(kg)。k=3时，LSR0.05，16=2.08(kg)，LSR0.01，16=2.88(kg)。,以各品种的小区平均产量为单位进行新复极差测验，求得：,品种间比较,品种小区平均产量(kg)差异显著性5%1%A216.22aAA112.33bBA311.67bB,表8.5三个品种小区平均产量新复极差测验,推断：品种A2与A1、A3的差异达=0.01水平，A1与A3间差异不显著。因此，A2品种平均产量最高，极显著地优于A1、A3品种。,施肥水平间比较仍以施肥水平的小区平均产量进行比较，求得：因施肥水平间与品种间的标准误相同，故LSR也一样。显著性测验结

8、果如表8.6。,施肥量小区平均产量(kg)差异显著性5%1%B213.11bABB311.44bB,表8.6三种施肥量小区平均产量间差异显著性测验,推断：以上分析表明施B1水平的肥量，各品种的平均产量最高，显著高于B2、B3，并与B3差异达极显著。,水平组合间的比较AB的F测验显著、说明不同品种要求的施肥量是不相同的。因此，还需要比较两因素水平组合之间的差异显著性。简便而常用的方法是对A各水平下B间(或B各水平下A间)作多重比较。,各品种在不同施肥水平下的小区平均产量比较对各Ai，算得AB的标准误为,k=2,LSR0.05,16=3.44(kg)，LSR0.01,16=4.73(kg)。k=3

9、,LSR0.05,16=3.61(kg)，LSR0.01,16=4.97(kg)。,施肥平均(kg)差异显著性水平产量5%1%B115.00aAB212.67abABB39.3bB,表8.7各品种在不同施肥水平下的差异显著性,施肥平均(kg)差异显著性水平产量5%1%B217.33aAB317.00aAB114.33aA,A2品种,A1品种,比较结果列于表8.7，,施肥平均(kg)差异显著性水平产量5%1%B117.67aAB29.33bBB38.00bB,A3品种,从表8.7可看出,A1品种以B1施肥量产量最高，它与B2无显著差异，但与B3差异极显著；A2品种以B2施肥量产量最高，但与B3、

10、B1差异不显著；A3品种以B1施肥量最优，并与B2、B3有极显著差异。,另一种方法是直接比较全部九个处理的差异显著性新复极差法的显著标准,9个处理间的差异显著性（SSR）,(5)试验结论参试品种水平间有显著差异:以A2平均产量最高，与A1，A3均有极显著差异。施肥量水平间有显著差异:以B1产量为最高，与B2、B3有显著差异，并与B3达极显著差异。品种与施肥量互作显著:A1，A3品种应取B1施肥量为优，A2品种在3种施肥量下产量差异不显著。,二因素随机区组试验结果的分析实例(随机模型）,例：将水稻的3个不同细胞质源的不育系（A1、A2、A3）和五个恢复系杂交，配成15个F1。采用随机区组设计，重

11、复2次，小区计产面积6平方米。其田间排列和小区稻谷产量列于下表。试作分析。,处理和区组两向表,因素A和因素B的两向表,变异期望均方来源固定模型随机模型区组2e+ab22e+ab2A2e+rb2A2e+r2(AB)+rb2AB2e+ra2B2e+r2(AB)+ra2BAB2e+r2(AB)2e+r2AB误差2e2e,二因素随机区组的期望均方,方差分析表,解得：,故产量性状的广义、狭义遗传力为：,二、三因素随机区组试验结果的方差分析1.三因素随机区组试验的线性模型与期望均方设试验有A，B，C三个因素，各具有a，b，c个水平，且将abc个处理安排在r个随机完全区组内，则该试验有rabc个观察值。由于

12、有A，B，C三个因素参加试验，处理效应则为三个因素的总效应，它可以分解为Ai，Bj，Ck，(AB)ij，(AC)ik，(BC)jk和(ABC)ijk七个部分。,各观察值的线性模型为：,式中：h=1，2，r；i=1，2，a；j=1，2，b；k=1,2，c。,根据上述分析，三因素随机区组试验的总变异可分解为区组、处理和误差，而处理间又分解为A因素，B因素，C因素，AB，AC，BC及ABC互作七项。各项变异来源的名称、符号、计算公式、期望均方均分别列在表8.8。,表8.8的期望均方提供了三因素随机区组试验正确进行F测验的依据，测验时仍按分子均方EMS比分母均方EMS只多一个分量的原则。当有现成适合的

13、被比量时，直接计算F值进行F测验。当无现成适合的被比量时，可通过有关项均方相加来近似测验。,例如，在随机模型中，测验H0：2A=0对HA：2A0，,于是，由F=s21/s22可测验H0：2A=0对HA：2A0。,再将AB和AC的均方相加得s22=s2AB+s2AC，估计22+2r2ABC+rc2AB+rb2AC。,先将A和ABC的均方相加得s21=s2A+s2ABC，估计22+2r2ABC+rc2AB+rb2AC+rbc2A，,此F值的近似自由度为：,以后遇到随机模型以及混合模型的某些变异来源没有现成适合的被比量时，F测验均可仿此进行近似测验。但当期望均方中某些分量的H0被接受，使该变异来源可

14、以直接测验时，则不应使用上述近似方法。,2、三因素随机区组试验结果的分析示例【例8.2】有一玉米三因素随机区组试验，A因素为品种，有A1，A2，A3三个水平(a=3)；B因素为密度，有B1，B2二个水平(b=2)；C因素为施肥量，有C1，C2二个水平(c=2)；重复三次，小区计产面积44.4m2。试验处理与代号、田间排列和小区产量(kg)见图8.2，试作分析。,图8.2玉米三因素随机区组试验的田间排列和小区产量,(1)结果整理先将试验结果按区组与处理两向分组整理成表8.9。,8.9图8.2资料的区组与处理两向表区组A1A2A3B1B2B1B2B1B2C1C2C1C2C1C2C1C2C1C2C1

15、C2Tr333230293028272628262524338343330292928292727252625342323229282927282828262524336Tt999789868883848183777673T=1016,A，C两向表A，B两向表B，C两向表A1A2A3TCA1A2A3TBC1C2TBC1188172159517B1196171160527B1270257527C2183164150497B2175165149489B2249240489TA3713363091016TA3713363091016TC5194971016,表8.10图8.2资料二因素两项表,再按任

16、二因素两向分组整理成表8.10的，。,(2)自由度与平方和的分解r=3,a=3,b=2,c=2总变异：dfT=rabc-1=3322-1=35,区组变异：dfr=r-1=2,处理间变异：dft=abc-1=11,A品种间：dfA=a-1=2,B密度间：dfB=b-1=1,C施肥间：dfC=C-1=1,AB:df(AB)=(a-1)(b-1)=2,AC:df(AC)=(a-1)(c-1)=2,BC:df(BC)=(c-1)(b-1)=1,ABC:dfABC=(a-1)(b-1)(c-1),误差：dfe=(r-1)(abc-1),(3)方差分析表和F测验这三个试验因素都取固定模型，所以各项均方都可

17、与误差均方相比而得出F值。,表8.11玉米三因素试验的方差分析表变异来源DFSSMSFF0.05F0.01区组21.550.781.63处理11226.2220.57A(品种)2161.0580.53167.76*3.445.72B(密度)140.1140.1183.56*4.307.94C(施肥)113.4413.4428.00*4.307.94AB29.734.8710.14*3.445.72AC20.730.371BC10.450.451ABC20.710.361误差2210.450.48总变异35238.22,F测验表明，A，B，C三个因素主效以及一级互作AB均达到极显著，其余皆不显著

18、。所以应进一步对A，B，C以及AB水平组合作多重比较。,(4)各因素及水平组合的显著性测验(SSR)品种间比较,A3品种的亩产量=3091.25=386.3(千克/亩),A2品种的亩产量=3361.25=420.0(千克/亩),A1品种的亩产量=3711.25=463.8(千克/亩),每个TA是rbc=322=12个小区的产量总和，故将小区总产量折算为亩产量的改算系数为，cf=666.67/(1244.4)=1.25,以亩产量为单位比较,k=2时，LSR0.05,22=8.79；LSR0.01,22=11.97k=3时，LSR0.05,22=9.24；LSR0.01,22=12.51以上述尺度

19、去测验三个品种亩产量的差数，其测验结果列于表8.12。,表8.12三个品种亩产量的新复极差测验品种亩产量(kg)差异显著性5%1%A1468.8aAA2420.0bBA3386.3cC,推断：A1，A2，A3三个品种间亩产量的差异均达到极显著，A1品种为最优。,密度间和施肥量间比较,当变异来源的自由度为1时，t测验、SSR测验和F测验的结果完全相同。所以，对密度间、施肥水平间可以根据F测验直接作出判断：B1极显著地高于B2，C1极显著地高于C2，而不必再作测验。,因为AB交互作用极显著，说明不同品种适宜不同的密度，所以还需测验AB水平组合之间的差异显著性。将各AB水平组合的亩产量列于表8.13

20、，根据如下计算做显著性测验。,AB水平组合间的比较,cf=666.67/(644.4)=2.50,k=2时，LSR0.05,22=12.4，LSR0.01,22=16.9,k=5时，LSR0.05,22=13.7，LSR0.01,22=18.5,k=4时，LSR0.05,22=13.4，LSR0.01,22=18.1,k=3时，LSR0.05,22=13.1，LSR0.01,22=17.7,A1B1的亩产量=小区总产（99+97）cf,A3B2的亩产量=小区总产（76+73）cf,A3B1的亩产量=小区总产（83+77）cf,A2B2的亩产量=小区总产（84+81）cf,A2B1的亩产量=小区

21、总产（88+83）cf,A1B2的亩产量=小区总产（89+86）cf,品种密度亩产量(kg)差异显著性5%1%A1B1490.0aAA1B2437.5bBA2B1427.5bBCA2B2412.5cCDA3B1400.0dDA3B2372.5eE,表8.13不同品种不同密度的新复极差测验,比较结果表明，在5%显著水平下，除A1B2与A2B1之间差异不显著外，其余品种密度间均存在显著差异，在1%显著水平下，A1B2与A2A1，A2B1与A2B2，A2B2与A3B1之间的差异不显著，其余品种密度间存在极显著差异。,(5)试验结论本试验品种、密度、施肥量间均有极显著差异，品种应选A1，密度应选B1，

22、施肥量应选C1。但AB互作显著，根据AB水平组合间显著性测验结果，A1B1较其他品种密度极显著增产，与C1组合为A1B1C1是本试验的最佳组合。,第二节裂区试验的方差分析一、二因素裂区试验结果的方差分析1二因素裂区试验的线性模型和期望均方设有A，B二因素，A因素为主处理，具a个水平，B因素为副处理，具b个水平，有r个完全区组，则全试验共有rab个观察值。,样本中任意观察值的线性可加模型为：,式中：h=1,2，r;i=1,2，a;j=1,2，b；,h,估计区组效应,Ai,估计主区因素A的主效，,(e1)hi,估计主区误差，,Bj,估计副区因素B的主效，,（AB）ij,估计二因素互作效应,(e2)

23、hij,估计副区误差。,由此不难看出，二因素裂区试验与二因素随机区组试验在分析上大致相同，不同之处仅在于将误差分成了主区误差和副区误差两部分，分别用来测验主处理、副处理以及主与副处理互作的的显著性。,二因素裂区设计的方差分析,当选用固定模型时，Ea可以测验区组间和主处理（A）水平间均方的差异显著性；Eb可用以测验副处理（B）水平间和AB互作间均方的差异显著性。,多重比较时的标准误,2、二因素裂区试验结果的分析示例【例8.3】设有一水稻施N量(A)、品种(B)试验，主处理为A，有A1，A2，A3(a=3)三个水平，副处理为B，有B1，B2，B3，B4(b=4)四个水平，裂区设计，重复三次(r=3

24、)，副区计产面积13.34m2，其田间排列和小区产量(单位：kg)列于图8.3，试作分析。P129图8.3水稻施N量与品种裂区试验田间排列和小区产量,(1)试验结果的整理先将图8.3的试验结果按区组和处理作两向分组整理成表8.16.再按A和B因素作两向分组整理成表8.17。,表8.16图8.3资料处理与区组两向表主处理A副处理B区组TABTAA1B114141341B213121237B311101132B415181952Tm535455162A2B119212262B220222163B319232466B425262374Tm839290265A3B116181953B217162053

25、B321252773B420192160Tm747887239Tr210224232T=666,表8.17图8.3资料的A和B两向表B1B2B3B4TA平均数A14137325216213.50A26263667426520.08A35353736023919.92TB156153171186T=666平均数17.317.019.020.618.5,(2)自由度与平方和的分解表8.16和表8.17中，Tr为各区组总和，TAB为各处理总和。TA为A因素各水平总和，TB为B因素各水平总和，Tm为各主区总和，T为全试验总和。,自由度的分解,总变异dfT=rab-1=334-1=35,主区部分：,主区

26、dfm=ra-1=33-1=8,区组dfr=r-1=3-1=2,A因素dfA=a-1=3-1=2,误差dfEa=(r-1)(a-1)=22=4,副区部分：,B因素dfB=b-1=4-1=3,A与B互作df(AB)=(a-1)(b-1)=23=6,误差dfEb=a(r-1)(b-1)=323=18,区组SSr=T2r/ab-C=(2102+2242+2322)/(34)-12321=20.67,平方和的分解,矫正数C=T2/rab=6662/(334)=12321,总变异SST=x2-C=142+142+212-12321=739,主区总变异SSm=T2m/b-C=(532+542+872)/4

27、-12321=512,处理SSAB=T2AB/r-C=(412+372+602)/3-12321=669,A因素SSA=T2A/rb-C=(1622+2652+2392)/(34)-12321=478.17,主区误差SSEa=SSm-SSr-SSA=512-20.67-478.17=13.16,平方和的分解,副区误差SSEb=SST-SSm-SSB-SSAB=739-512-77-113.83=36.17,由表8.17可求得：,AB互作SSAB=SSAB-SSA-SSB=669-478.17-77=113.83,B因素SSB=T2B/ra-C=(1562+1532+1712+1862)/(33

28、)-12321=77,平方和的分解,(3)方差分析表和F测验将上述计算结果录入表8.18，因本试验研究各指定处理的效应，故按固定模型作F测验。结果表明，施N量(A)间，品种(B)间，以及AB互作均达极显著。,表8.18图8.3资料的方差分析表变异来源DFSSMSFF0.01主区部分区组220.6710.34A2487.17239.0972.67*.18.00Ea413.163.29总变异8512.00副区部分B377.0025.6712.77*5.09AB6113.8318.979.44*4.01Eb1836.172.01总变异35739.00,(4)主效的显著性测验在此以小区平均产量用新复极

29、差法进行测验。主处理(施N量)间比较k=2时，LSR0.05，4=2.06，LSR0.01，4=3.41。k=3时，LSR0.05，4=2.10,LSR0.01，4=3.56。,以上述尺度测验施N量(A)的各水平间差数，其测验结果列于表8.19，表明A2，A3的小区平均产量都极显著高于A1；A2与A3之间无显著差异。,副处理(品种)间比较见表8.20。以上述尺度测验品种(B)的各水平小区平均产量间的差异，结果列于表8.21。表明B4品种最优，它显著高于B3品种，极显著高于B1和B2品种；B3的小区平均产量显著高于B1和B2品种，B1和B2品种间差异不显著。,表8.19三种施N量小区平均产量的S

30、SR测验施氮量小区平均产量差异显著性5%1%A320.08aAA219.92aAA113.50bB,表8.20副处理比较的LSRk234LSR0.05，181.401.481.52LSR0.01，181.932.022.07表8.21品种(B)小区平均产量间SSR测验品种小区平均数差异显著性5%1%B420.67aAB319.00bABB117.33cBB217.00cB,同一主处理(施N量A)内不同副处理(品种)间比较采用SSR测验，见表8.22。根据表8.15的有并以此为尺度测量不同施N水平下，各品种小区平均产量的差异显著性于表8.23中。结果表明，在A1和A2施N量下以B4品种为最优，在

31、A3施N量下以B3品种为最优。,表8.22互作比较的LSR值k234LSR0.05，182.4322.9952.629LSR0.01，183.3333.4973.587,表8.23各施N量下不同品种的小区平均产量及其差异显著A1施N量A2施N量A3施N量品种产量差异显著性品种产量差异显著性品种产量差异显著性5%1%5%1%5%1%B417.33aAB424.67aAB324.33aAB113.67bBB322.00bABB420.00bBB212.33bcBB221.00bBB117.67bBB310.67cBB120.67bBB217.67bB,表8.24各处理组合平均数比较的LSR值计算表

32、k23456789SSR0.053.313.433.503.543.573.603.613.62SSR0.014.935.115.205.265.305.315.365.39LSR0.052.923.023.083.123.153.173.183.19LSR0.014.344.504.504.644.674.714.724.75,表8.25各处理组合平均数比较(SSR测验)水平组合平均产量(kg)差异显著性5%1%A2B424.67aAA3B324.33aAA2B322.00abABA2B221.00bABA2B120.67bcABA3B420.00bcABA3B117.67cBA3B217.67cBA1B417.33cBA1B113.67dBCA1B212.33dCA1B310.67dC,