常用试验设计及其统计分析-OK课件.ppt

1、12022-9-29生物学试验研究中，一项工作要取得客观、理想的结果，必须做到试验方案设计合理，精心组织操作，采用相应的统计方法对试验结果进行分析。本章主要讨论试验设计的基本原理和常用试验设计方法及其统计分析。22022-9-293制定制定试验试验方案的要点方案的要点2022-9-29 广义的试验设计广义的试验设计是指整个研究课题的设计，包括试验方案的拟订，试验单位的选择，分组的排列，试验过程中生物性状和试验指标的观察记载，试验资料的整理、分析等内容；狭义的试验设计狭义的试验设计则仅是指试验单位的选择、分组与排列方法。42022-9-29合理的试验设计对科学试验是非常重要的。它不仅能够节省人力

2、、物力、财力和时间，更重要的是它能够减少试验误差，无偏估计误差，提高试验的精确度，取得真实可靠的试验资料，为统计分析得出正确的判断和结论打下基础。52022-9-291 1、试验目的要明确、试验目的要明确明确选题，制定合理的实验方案。一是要抓住当时生产实践和科学实验中急需解决的问题，二是要照顾到长远和不久的将来可能突出的问题。2 2、试验条件要有代表性、试验条件要有代表性试验条件应能代表将来准备推广试验结果的地区的自然条件、经济和社会条件。62022-9-293 3、试验结果要可靠、试验结果要可靠 试验结果的可靠程度主要用准确度与精确度进行描述。准确度指观察值与真值的接近程度，由于真值是未知数

3、，准确度不容易确定，故常设置对照处理，通过与对照相比以了解结果的相对准确程度。精确度是指试验中同一性状的重复观察值彼此接近的程度，即试验误差的大小，它是可以计算的。试验误差越小，处理间的比较越精确。72022-9-294 4、试验结果要能够重演、试验结果要能够重演 指在相同条件下，再次进行试验或实验，应能获得与原试验相同的结果。注意保持试验条件的一致性。82022-9-29 试验设计包括三个基本组成部分，即：处理因素、受试对象和处理效应。1.1.处理因素处理因素 一般是指对受试对象给予的某种外部干预(或措施)，称为处理因素，简称处理处理。试验因素：试验因素：在科学试验中，被变动的并设有待比较的

4、一组处理的因子称为试验因素，简称因素或因子（factor)。92022-9-29 水平：水平：试验因素的量的不同级别或质的不同状态称为水平（level）。试验水平可以是定性的，如供试的不同品种，具有质的区别，称为质量水平质量水平；也可以是定量的，如N肥的施用量，具有量的差异，称为数量水平数量水平。单因素试验单因素试验(single-factor experiment)：整个试验中只变更、比较一个试验因素的不同水平，其它作为试验条件的因素均严格控制一致的试验。102022-9-29 多因素试验多因素试验(multiple-factor or factorial experiment)：在同一试验

5、方案中包含两个或两个以上的试验因素，各个因素都分为不同水平，其它试验条件严格控制一致的试验。2.2.受试对象受试对象 受试对象是处理因素的客体，实际上就是根据研究目的而确定的观测总体。试验指标：试验指标：用于衡量试验效果的指标性状称为试验指标试验指标(experimental indicator)。112022-9-293.3.处理效应处理效应 处理效应是处理因素作用于受试对象的反应，是研究结果的最终体现。简单效应简单效应(simple effect)：同一因素内两种水平间试验指标的差数。主要效应主要效应(main effect)：一个因素内各简单效应的平均数，又称平均效应，简称主效。交互作用

6、交互作用(interaction effect)：两个因素简单效应间的平均差异称为交互作用效应，简称互作122022-9-2913PNN1N2平均N2N1P11016136P21824216平均14206P2P18880，0/2=02022-9-2914PNN1N2平均N2N1P11016136P218282310平均14228P2P1812104/2=22022-9-2915PNN1N2平均N2N1P11016136P21814164平均14151P2P182310/2=5两个因素间的互作称为一级互作一级互作，三个因素间的互作称为二级互作二级互作，余类推。一级互作易于理解，实际意义明确。二级

7、以上的互作较难理解，实际意义不大。2022-9-297.1.4 7.1.4 制定试验方案的要点制定试验方案的要点 试验方案：是根据试验目的和要求所拟进行比较的一组试验处理(treatment)的总称。1 1、明确试验目的、明确试验目的 通过回顾以往的研究进展、调查研究、文献探索等明确试验的目的，形成对所研究主题及外延的设想，使待拟订的试验方案能针对主题确切而有效地解决问题。162022-9-29 2 2、根据试验目的确定恰当的供试因素及水平、根据试验目的确定恰当的供试因素及水平 供试因素不宜过多，应该抓住1-2个或少数几个主要因素解决关键性问题。每因素的水平数目也不宜多，且各水平间距要适当，使

8、各水平能明确区分，并把最佳水平范围包括在内。例如通过喷施矮壮素控制玉米株高，其浓度试验设置为50、100、150、200、250mg/L等5个水平，其间距为50mg/L。如果将间距缩小为10mg/L，水平数猛增到20个。172022-9-29 这会导致两方面的问题：一是实验无法进行；二是受误差影响不容易发现试验效应的规律。3、试验方案中应包括对照水平或处理试验方案中应包括对照水平或处理(check,CK)对照是试验中比较处理效应的基准。品种比较试验中常统一规定同生态区内使用的对照品种。4、注意比较间的唯一差异性原则，才能正确注意比较间的唯一差异性原则，才能正确解析出试验因素的效应解析出试验因素

9、的效应。182022-9-29 例如，在对小麦进行叶面喷施P肥的试验中，如果只设施P（A）与不施P（B）两个处理，因为P肥是兑在水中然后喷到小麦叶面上的，两者的差异可能有P的作用，也可能有水的作用，无法将它们区分开。如果再加入一个喷施等量清水的处理（C），则P和水的作用可分别从A与C及B与C的比较中解析出来。5 5、正确处理试验因素与试验条件间的关系、正确处理试验因素与试验条件间的关系（1）试验因素的表现受试验条件的制约（2）注意试验条件的代表性与典型性192022-9-29 6 6、尽量用多因素试验、尽量用多因素试验（1）在同一试验中提供了比单因素试验更多的效应估计；（2）误差自由度多，试验

10、精确度提高。202022-9-29 1.1.试验误差的来源试验误差的来源(1)试验材料固有的差异 如基因型不一致、种子生活力有差异、秧苗素质有差异等(2)试验条件不一致 如各试验单位所处的外部环境不一致。田间试验中农事操作和管理技术的不一致。（3)操作技术不一致(4)偶然因素的影响。212022-9-29 2.2.控制试验误差的途径控制试验误差的途径(1)选择同质一致的试验材料。(2)改进操作和管理技术，使之标准化。(3)精心选择试验单位。各试验单位的性质和组成要求均匀一致。(4)采用合理的试验设计。222022-9-29进行试验设计的目的，在于降低试验误差，无偏估计误差，提高试验的准确度与精

11、确度，使试验结果正确可靠。为了有效地控制和降低误差，试验设计必须遵循下面三条基本原则。232022-9-291.1.重复重复 定义：重复重复(replication)在试验中同一处理设置的试验单位数。作用：作用：(1)估计误差 i=yi 式中为总体平均数，是一个无法得到的理论值。在实际工作中，通常用样本的平均数来估计。而24_y2022-9-29 单一小区所得数值易受特别高或低的肥力的影响，多次重复所估计的处理效应（平均数）可以抵消部分误差的影响，使处理间的比较更加可靠。25nyyniii1_2022-9-29 2.随机 随机随机 (random)：指一个重复中每个处理都有同等的机会设置在任何

12、一个试验单位上，避免任何主观成见。作用：使估计的误差无偏。方法：(1)抽签法(2)随机数字表（P351）处理在9个以内，直接用随机数字表。中任意指定页中的任意一行的数字次序即可。例如：有8个处理，分别用1、2、3、4、5、6、7、8代表。在随机数字表中得到一行随机数字为：262022-9-29 5264862339，9718302620 去掉序列中的0、9和重复数字，得到：52648371 这就是8个处理在区组内的排列顺序，即第一小区安排5号处理，第二小区安排2号处理，第三号小区安排6号处理，余类推。多于9个的处理，从随机数字表中任意行开始，每次取两位数。如12个处理，可查任何一页的一行，去掉

13、00、97、98、99后，凡大于12的数均被12除后得余数，将重复数字划去，即得到12个处理排列的次序。272022-9-29 去掉00、97、98、99这几个值是为了保证每个处理都有相同的次数被取到，12个处理，从01到96这些数字中，每个处理都可能取8次。例如：从随机表中取得97、39、24、89、90、89、86、49、15、18、25、43、80、74、30、41、67、36、43、58、42、07、04、25、17、54、60、88、49、34、42等随机数，除去97，大于12的数用12除后取余数，将重复数字划去，所得随机排列为：3 3、1212、5 5、6 6、2 2、1 1、7

14、7、8 8、1010、4 4、9 9、1111282022-9-29 3.3.局部控制局部控制 将整个试验环境分解成若干个相对一致的小环境(称为区组、窝组或重复)，再在小环境内分别配置一套完整的处理，在局部对非处理因素进行控制。作用：降低试验误差。方法：在田间试验中将试验田划分成等于重复数的区组，区组内的肥力水平尽可能保持一致；在温室试验中，将区组安排在同一光照水平上；在微生物接种试验中，将接种时间安排为区组。292022-9-29三个基本原则的关系和作用重复30无偏的试验误差估计随机局部控制降低试验误差2022-9-29312022-9-29completely random design)

15、1、特点：使用了试验设计三个原则中的两个（重复、随机），能够得到无偏的误差估计值，但控制试验环境误差的能力不强。2、常用于试验环境因素相当均匀的场合，如实验室培养试验、网室温室的盆钵试验。322022-9-29 3 3、设计示例、设计示例 有三种生长激素，分别用A、B、C代替，测定其对小麦株高的影响，包括对照（用等量的清水）在内，共4个处理，进行盆栽试验，每盆小麦为一个单元，每处理用4盆（重复4次）共16盆。第一步：用数字代表处理 A：14，B：58，C：912，CK：1316 第二步：抽签或查随机数字表，得到随机数字 14、9、7、1、5、12、16、3、11、8、4、2、6、13、10、1

16、5332022-9-29 第三步：将随机数字对应的处理安排到相应的盆内。34ckCBABCckACBAABckCck2022-9-29 第五章方差分析的例子用的就是完全随机设计，请参见教材第六章第四节等有关内容，这里从略。352022-9-29362022-9-29 1、特点：特点：使用了田间试验设计三个原则，并根据“局部控制”的原则，将试验地按肥力程度划分为等于重复数的区组，一区组安排一重复，区组内各处理独立地随机排列。是田间试验最常用的设计。2、优缺点：、优缺点：优点：（1）设计简单，容易掌握；（2）富于伸缩性，单因素、多因素以及综合性试验都能用；（3）能提供无偏的误差估计，并有效减小单向

17、的肥力差异，降低误差；372022-9-29（4）对试验地要求不严，必要时，不同的区组可以分散设置在不同地段上。缺点：（1）设计不允许处理数太多，一般不超过20个；（2）只能在一个方向上控制土壤差异。382022-9-29 3、设计示例设计示例(1)8个处理，4次重复，共32个小区。39肥力梯度IIIIIIIV251483765142867364537218452413782022-9-29（2）16个处理，3次重复，小区布置成两排40肥力梯度IIIIII13810715 149613416 1121252022-9-29(3)区组布置在不同的地块上41IIIIII2022-9-29一、单因素

18、随机区组试验结果的方差分析一、单因素随机区组试验结果的方差分析 可将处理看作A因素，区组看作B因素，其余部分则为试验误差。设试验有k个处理，n个区组，则自由度与平方和的分解为：nk-1=(n-1)+(k-1)+(n-1)(k-1)总自由度=区组自由度+处理自由度+误差自由度422022-9-2943kntrktnrknyyyyyynyykyy1122121112)()()()(总平方和=区组平方和+处理平方和+误差平方和例12.3P228有一小麦品比试验，其有A、B、C、D、E、F、G、H8个品种(k=8)，其中A是标准品种(ck)，采用随机区组设计，重复3次(n=3)，小区计产面积25m2，

19、其产量如下，试作分析。2022-9-2944品 种区 组TtIIIIIIA10.99.112.232.210.7B10.812.314.037.112.4C11.112.510.534.111.4D9.110.710.129.910.0E11.813.916.842.514.2F10.110.611.832.510.8G10.011.514.135.611.9H9.310.414.434.111.4Tr83.191.0103.9T=278.010.411.413.0tyry2022-9-29 1.1.自由度与平方和分解自由度与平方和分解(1)自由度的分解 总DFT=nk-1=(38)-1=23

20、 区组DFR=n-1=3-1=2 品种DFt=k-1=8-1=7 误差DFe=(n-1)(k-1)=(3-1)(8-1)=14(2)平方和的分解452022-9-294617.32208327822nkTC61.844.141.99.102222CCySST56.2789.1030.911.83)(222212CCkTyykSSrnrR2022-9-29 2.F 测验4708.3431.341.372.32)(222221CCnTyynSStktt97.2208.3456.2761.84tRTeSSSSSSSS变异来源DFSSMSFF0.05区组间227.5613.788.403.74品种间7

21、34.084.872.972.77误 差1422.971.64总变异2384.612022-9-293.3.品种平均数的比较品种平均数的比较本例目的是测验各供试品种是否与标准品种A有显著差异，宜应用LSD法。48)(05.1364.12221kgnMSseyy由于=14时，t0.05=2.145,t0.01=2.977,故LSD0.05=1.052.145=2.25(kg)LSD0.01=1.052.977=3.132022-9-29 各品种产量与对照相比的差异显著性49品种差异E14.2 3.5B12.41.7G11.91.2H11.40.7C11.40.7F10.80.1A(CK)10.7

22、D10.0-0.7ty2022-9-29二、二因素随机区组试验结果的方差分析 设有A和B两个试验结果，各具a和b个水平，那么共有ab个处理组合，作随机区组设计，有r次重复，则该试验有rab个观察值。其自由度与平方和分解为：abr-1=(r-1)+(ab-1)+(r-1)(ab-1)总自由度=区组自由度+处理自由度+误差自由度502022-9-295121212121)()()()(yyyyyyryyabyykljabrjklabklrjabrjkl总平方和=区组平方和+处理平方和+误差平方和接下来，对处理项进行再分解 ab-1=(a-1)+(b-1)+(a-1)(b-1)处理自由度=A的自由度

23、+B的自由度+AB自由度2022-9-295221211221)()()()(yyyyryyrayyrbyyrlkabklblakabkl SSt =SSA +SSB +SSAB2022-9-29 二因素随机区组试验自由度的分解53变异来源DFSS区 组r-1处理组合ab-1Aa-1Bb-1AB(a-1)(b-1)误 差(r-1)(ab-1)总变异rab-1CabTSSrR/2CrTSSABt/2CrbTSSAA/2CraTSSBB/2BAtABSSSSSSSStRTeSSSSSSSSCySST22022-9-29例13.1P249有一早稻二因素试验，A因素为品种，分A1(早熟)、A2(中熟)

24、、A3(晚熟)三个水平(a=3),B因素为密度，分B1(16.56.6cm)、B2(16.59.9cm)、B3(16.513.2cm)三个水平(b=3)，共ab=33=9个处理，重复3次(r=3),小区计产面积20m2。其田间排列和小区产量(kg)如下图，试作分析。54A1B18A2B27A3B310A2B38A3B28A1B36A3B17A1B27A2B19A2B37A3B27A1B27A3B17A1B35A2B19A2B29A3B39A1B18A3B16A1B36A2B18A1B26A2B26A3B39A1B18A2B36A3B28IIIIII2022-9-291.资料整理55处理处理区组

25、区组I区组区组II区组区组IIITABA1B188824A1B277620A1B365617A2B199826A2B279622A2B387621A3B377621A3B287823A3B3109928Tr706863T=201(1)区组与处理的两向表2022-9-29(2)品种(A)和密度(B)的两向表56 BAB1B2B3TAA124201761A226222169A320232871TB706566T=2012.自由度与平方和分解33.149633320122rabTC67.409882222CCySST2022-9-29 在A、B因素两向表的基础上对处理平方和进行再分解5789.233

26、6368702222CCabTSSrR00.3032820242222CCrTSSABt78.700.3089.267.40tRTeSSSSSSSS2022-9-295821.2256.123.600.3056.13366657023.63371696122222222BAtABBBAASSSSSSSSCCraTSSCCrbTSS2022-9-293.方差分析表和F 测验59变异来源DFSSMSFF0.05区组间22.891.452.963.63处理(组合)间830.003.757.652.59 品种26.233.126.373.63 密度21.560.781.593.63 品种密度422.2

27、15.5511.333.01误差167.780.49总变异2640.674.差异显著性测验(1)品种间比较)(233.0949.0kgrbMSSEe2022-9-29 三个品种小区产量的新复极差测验60pSSR0.05SSR0.01LSR0.05LSR0.0123.004.130.700.9633.154.340.731.01品种产量差异显著性5%1%A37.9aAA27.7aABA16.8bB2022-9-29(2)品种密度互作61pSSR0.05SSR0.01LSR0.05LSR0.0123.004.131.211.6733.154.341.271.75)(404.0349.0kgrMSS

28、Ee各品种在不同密度下的小区平均产量及差异显著性品种产量差异显著性5%1%B18.0aAB26.7bABB35.7bB品种产量差异显著性5%1%B18.7aAB27.3bABB37.0bB品种产量差异显著性5%1%B39.3aAB27.7bABB16.7bBA1品种A2品种A3品种2022-9-295.试验结论 本试验品种主效有显著差异，以A3产量最高，与A1有显著差异，而与A2差异不显著。密度主效无显著差异。但品种与密度的互作极显著，A3品种需要用B3密度，A2品种需用B1密度，才能取得最高产量。622022-9-29632022-9-29 将k个不同符号排成k列，使每个符号在每一行、每一列

29、都仅出现一次的方阵，叫kk 拉丁方拉丁方 1、特点：将处理从纵横两个方向排列成区组，具有双向局部控制的能力，因而有较高的精确度。2、优缺点：精度高，但缺乏伸缩性，因为在设计中，重复数必须等于处理数，两者相互制约。642022-9-293、使用范围：只限于48个处理，不能象随机区组那样区组可以分开，故在田间试验时要求有整块平坦的土地。在动物实验中，如要控制来自两个方向的系统误差，且在动物头数较少情况下，常采用这种设计方法。652022-9-294.设计示例研究5种不同饲料(分别用1,2，3,4，5号代表)对乳牛产乳量影响试验，选择5头乳牛，每头乳牛的泌乳期分为5个阶段，随机分配5个饲料的5个水平

30、。由于乳牛个体及牛的泌乳期不同对产乳量都会有影响，故可以分别把其作为区组设置(牛号用I至V表示，为横向区组；泌乳期用一月至五月表示，为纵向区组)，采用一个55的拉丁方设计。662022-9-2967（1）选择标准方标准方：第一直行和第一横行均为顺序排列的拉丁方。教材P25表2.1中列出了（44）（88）的选择标准方。A B C D EB A E C DC D A E BD E B A CE C D B A2022-9-29(2)按随机数字1、4、5、3、2调整直行(对泌乳期区组进行随机)68A B C D EB A E C DC D A E BD E B A CE C D B A选择标准方AB

31、CDE第1行DCEAB第4行EDBCA第5行CEABD第3行BADEC第2行2022-9-29(3)按随机数字5、1、2、4、3调整横行(对牛号区组进行随机)69A D E C BB C D E A C E B A DD A C B EE B A D C调整直行后的拉丁方E B A D C第5行A D E C B第1行B C D E A第2行D A C B E第4行C E B A D第3行2022-9-29(4)按随机数字2A、5B、4C、1D、3E,安排饲料。70E B A D CA D E C BB C D E AD A C B EC E B A D调整横行后的拉丁方35214213455

32、41321245343521一I二三四五IIIIIIVV纵向区组横向区组2022-9-29拉丁方试验中行、列均为区组，因此在试验结果统计分析中比随机区组多一项区组间变异，即总变异可分解为处理间、行区组间、列区组间和试验误差四个部分。其自由度与平方和的分解为：71)2)(1()1()1()1(12kkkkkk总自由度=行自由度+列自由度+处理自由度+误差自由度2022-9-29722112212112)2()()()()(22yyyyyykyykyykyyckkrtkckkr SST=SS行 +SS列 +SSt +SSe1.结果整理将乳牛试验资料按横行、纵行，并计算总和，整理成表1，饲料处理的总

33、和与平均数列于表2。2022-9-2973月份一二三四五TrIIIIIIIVVE 300D 420B 350A 280C 400A 320C 390E 360D 400B 380B 390E 280D 400C 390A 350C 390B 370A 260E 280D 430D 380A 270C 400B 370E 32017801730177017201880Tc175018501810173017408880(T)表1 饲料类型对乳牛产乳量影响的试验结果(kg)2022-9-29 2.自由度和平方和的分解：74饲料饲料ABCDE总和总和Tt14801860197020301540888

34、0296372394406308ty表2 饲料的总和(Tt)与平均数()ty124-4-4-2441-5)(41-5)(41-5)(241-55etrcTdfdfdfdfdf误差自由度自由度饲料处理自由度乳牛横行自由度月份纵行总自由度2022-9-297563224-320420300-31541765588802222T2CcySSC总平方和矫正数214451740.185017502222CCkTSScc纵行(月份)平方和2022-9-2976横行(乳牛)平方和322451880.173017802222CCkTSSrr处理(饲料)平方和5050451540.186014802222CCk

35、TSStt2022-9-293.列方差分析表，计算F 值：77误差平方和7352505043224214463224trcTeSSSSSSSSSS变异来源变异来源DFSSMSFF0.05F0.01纵行纵行(月份月份)间间横行横行(乳牛乳牛)间间处理处理(饲料饲料)间间误误 差差4441221443224505047352536.00806.0012626.00612.6720.61*3.265.41总总 变变 异异24632242022-9-29784.比较各处理间的差异，采用q 测验。)(07.11567.612kgkMSSEep2345q0.05q0.01LSR0.05LSR0.013.0

36、84.3234.09647.8223.775.0441.73455.7934.205.5046.49449.9264.515.8449.92664.649表3 表1资料q测验的LSR值2022-9-2979饲料名称饲料名称平均产乳量平均产乳量差异显著性差异显著性0.050.01DCBEA406394372308296aaabbAAABB表4 不同饲料的乳牛产乳量比较(q测验)2022-9-29802022-9-29 正交设计是一种研究多因素试验的设计方法。在多因素试验中，随着试验因素和水平数的增加，处理组合数将急剧增加。例如，3因素3水平，就有33=27个处理组合，4因素4水平，就有44=25

37、6个处理组合。正交试验是利用到现在一套规格化的表格正交表，科学合理地安排试验。这种设计的特点是在试验的全部处理组合中，仅挑选部份有代表性的水平组合(处理组合)进行试验。通过部份实施了解全面试验情况，从中找出较优组合。812022-9-29 例如，要进行一个4因素3水平的多因素试验，如果全面实施需要34=81个处理。但是采用一张L9(34)的正交表安排试验，则只要9个处理组合就够了。822022-9-29 正交表是正交设计的基本工具。在正交设计中，安排试验、分析结果，均在正交表上进行。教材附表13(P380)给出了常用的正交表。现以L9(34)正交表为例，说明正交表的概念与特点。L表示一张正交表

38、，括号内下面的3表示因素的水平数，3的右上方为指数4，表示最多可以安排因素(包括互作)的个数。L右下角的数字9表示试验次数(水平组合数)832022-9-2984列列 号号ABCD水平组合水平组合1234试试验验号号123456789111222333123123123123231312123312231A1B1C1D1A1B2C2D2A1B3C3D3A2B1C2D3A2B2C3D1A2B3C1D2A3B1C3D2A3B2C1D3A3B3C2D1L9(34)正交表2022-9-29 1.正交表的两个性质：(1)每一列中不同数字出现的次数相等。(2)在任三列中，将同一横行的两个数字看成有序数对时

39、，每一数出现的次数相等。上表中有序数对共有9种：(1,1)、(1,2)、(1,3)、(2,1)、(2,2)、(2,3)、(3,1)、(3,2)、(3,3)，它们各出现一次，也就是说每个因素的每个水平与另一个因素的各个水平各碰到一次，也仅碰到一次，表明任何两因素的搭配是均衡的。由于正交表的这两个特点，所以用正交表安排的试验具有均衡分散均衡分散和整齐可比整齐可比的特性：852022-9-29(1)均衡分散：是说明正交表挑出来的这部分水平组合，在全部可能的水平组合中分布均匀，因此代表性强，能较好地反映全面情况。例如，对L9(34)正交表而言，如有三个因素，则全面试验为33=27次，它们的水平组合为：

40、8632111CCCBA32121CCCBA32131CCCBA32112CCCBA32122CCCBA32132CCCBA32113CCCBA32123CCCBA32133CCCBA2022-9-29（2）整齐可比：由于正交表中各因素的水平是两两正交的，因此，任一因素任一水平下都必须均衡地包含其它因素的各水平。例如，A1、A2、A3条件下各有三种B 水平，三种C 水平，即：873322111CBCBCBA3322112CBCBCBA3322113CBCBCBA2022-9-29 1.确定试验因素和水平数确定试验因素和水平数。例：为了解决花菜留种问题，科技人员考察了浇水、施肥、病害防治和移入温

41、室时间对花果留种的影响，进行了4 个因素各2水平的正交试验。各因素及水平见下表：88因 子水 平 1水 平 2A：浇水次数B：喷药次数C：施肥次数D：进室时间不干死为原则，整个生长期只浇12次水发现病害即喷药开花期施用硫酸铵11月初根据生长需要量和自然条件浇水，但不过湿。每半月喷一次发根期、抽苔期、开花和结实期各施一次11月15日2022-9-29 2.选用合适的正交表选用合适的正交表 其原则是既要能安排下全部试验因素，又要使部分试验的水平组合数尽可能的少。各试验因素的水平数减1之和加1，即为需要做的最少试验次数，若用互作，需要再加上互作的自由度。892022-9-29 本例试验最少需要做的试

42、验次数=(2-1)4+1=5，然后从2n因素正交表中选用处理组合数稍多于5的正交表安排试验，据此选用L8(27)正交表。3.进行表头设计，列出试验方案 所谓表头设计，就是把试验中挑选的各因素填到正交表的表头各列。902022-9-29表头设计的原则是：不要让主效间、主效与互作间有混杂现象。由于正交表中一般都有交互列，因此当因素少于列数时，尽量不在交互列中安排试验因素，以防发生混杂；当存在交互作用时，需查交互作用表，将交互作用安排在合适的列上，如本例若只考虑A、B因素间的互作，其表头设计如下：91列号1 2 3 4 5 6 7因子A B AB C AC D2022-9-29 表头设计好后，把该正

43、交表L8(27)中各列水平号换成各因素的具体水平就成为试验方案。92试验号试验号1列列(浇水次数浇水次数)2列列(喷药次数喷药次数)4列列(施肥方法施肥方法)7列列(进室时间进室时间)123456781 浇水浇水12次次1 浇水浇水12次次1 浇水浇水12次次1 浇水浇水12次次2 需要就浇需要就浇2 需要就浇需要就浇2 需要就浇需要就浇2 需要就浇需要就浇1 发病喷药发病喷药1 发病喷药发病喷药2 半月喷药一次半月喷药一次2 半月喷药一次半月喷药一次1 发病喷药发病喷药1 发病喷药发病喷药2 半月喷药一次半月喷药一次2 半月喷药一次半月喷药一次1 开花施开花施2 施施4次次1 开花施开花施2

44、 施施4次次1 开花施开花施2 施施4次次1 开花施开花施2 施施4次次1 11月初月初2 11月月15日日2 11月月15日日1 11月初月初2 11月月15日日1 11月初月初1 11月初月初2 11月月15日日2022-9-29 1.1.正交试验结果的直观分析正交试验结果的直观分析93A1B2AB3C4A C5D7种子产量12345678111122221122112211222211121212121212212112212112350325425425200250275375T1T2152511001125150013251300125013751400122513001325T=2

45、625381.25275.00281.25375.00331.25325.00312.50343.75350.00306.25325.00331.25R106.25-93.756.25-31.2543.75-6.25水平列号试验号21yy2022-9-29 1.1.逐列计算各因素同一水平之和：逐列计算各因素同一水平之和：第1 列A因素各水平之和：941100375275250200152542542532535021TT第2 列B因素各水平之和：1500375275425425112525020032535021TT2.2.逐列计算各水平的平均数：逐列计算各水平的平均数：第1 列A因素各水平的

46、平均数分别为：00.275281100225.38128152522211nTynTy2022-9-29 3.3.逐列计算各水平均数的极差：逐列计算各水平均数的极差：第一列A因子各水平平均数的极差为：9525.10600.27525.38121yyR4.4.比较极差，确定各因子或交互作用结果的影比较极差，确定各因子或交互作用结果的影响：响：浇水次数(A)和喷药次数(B)的极差|R|分居第一、第二位，是影响花菜种子产量的关键性因子。2022-9-295.5.水平选优与组合选优：水平选优与组合选优：根据各试验因子的总和或平均数可看出A取A1，B取B2，C取C2，D取D2为好。在没有互作的情况下，花

47、菜留种最好的栽培管理方式为：A1 B2 C2 D2。96从表中看出AC 对产量的影响较大，因此A和C 选那个水平应根据A与C 的最好组合来决定。2022-9-29975.3122/)375250(:5.2372/)275200(:0.3752/)425325(:5.3872/)425350(:22122111CACACACA在考虑 AC 互作的情况下，花菜留种的最适条件应为：A1 B2 C1 D2。2022-9-29 2.2.正交试验结果的方差分析正交试验结果的方差分析(1)(1)平方和与自由度的分解：平方和与自由度的分解：983.7031.38288/)12251400(1.788/)130

48、01325(1.788/)13251300(1.19538/)13751250(1.175788/)15001125(1.225788/)11001525(9.46796375.3253501.86132882625222222222222ACABDCBATeACABDCBATSSSSSSSSSSSSSSSSSSSSSSSSSSSSCCySSnTC2022-9-292.2.列方差分析表进行列方差分析表进行F F测验测验9911111117112718ACABDCBATeACABDCBATdfdfdfdfdfdfdfdfdfdfdfdfdfdfdf变异来源DFSSMSFF 0.05F 0.01浇

49、水次数喷药次数施肥方法进室时间浇水次数喷药次数浇水次数施肥方法试验误差111111122578.117578.11953.178.178.13828.1703.322578.117578.11953.178.178.13828.1703.332.1024.992.78115.44161161161161161161405405405405405405总变异746796.32022-9-29F测验各项变异来源均不显著，这是由于误差自由度太小所致。解决这个问题的根本办法是试验设置重复。折中的办法可以将F值小于1的变异项的平方和和自由度与误差项的平方和和自由度合并，作为试验误差平方和的估计值。100

50、5.8591.781.783.703ABDeeSSSSSSSS2022-9-29101变异来源变异来源DFSSMSFF 0.05F 0.01浇水次数浇水次数喷药次数喷药次数施肥方法施肥方法浇水次数浇水次数施肥方法施肥方法试验误差试验误差1111322578.117578.11953.13828.1859.522578.117578.11953.13828.1286.578.18*61.35*6.8213.36*10.1310.1310.13 10.1334.1234.1234.1234.12总变异总变异746796.3花菜留种正交试验的方差分析(去掉F1因子后)2022-9-293.3.互作分