1、.第五章第五章 专业试验统计专业试验统计.一、生物测定n生物测定，又称为半数效量，是生物、医学研究中最普遍生物测定，又称为半数效量，是生物、医学研究中最普遍的一类实验，通常包括供试因子的一类实验，通常包括供试因子(如化学或生物性杀虫剂、杀菌剂、如化学或生物性杀虫剂、杀菌剂、除草剂、植物生长调节剂、新医药、保健品或可能具有某种活性的物质或产除草剂、植物生长调节剂、新医药、保健品或可能具有某种活性的物质或产品品)和受试靶体和受试靶体(如一种昆虫、植物或小白鼠之类的实验动物如一种昆虫、植物或小白鼠之类的实验动物)。这类实。这类实验有定性与定量水平之分。定性测定较为简单，如仅测定验有定性与定量水平之分
2、。定性测定较为简单，如仅测定某种微生物是否对某种动物有致病性，其结果仅仅回答是某种微生物是否对某种动物有致病性，其结果仅仅回答是或非两种可能；定量测定则较为复杂，通常包括供试因子或非两种可能；定量测定则较为复杂，通常包括供试因子的一系列剂量，每剂量要处理合理数量的一系列剂量，每剂量要处理合理数量(越多实验误差越小越多实验误差越小)的的受试靶体并且整个实验要持续观察一段时间。受试靶体并且整个实验要持续观察一段时间。.一、生物测定n定性数据概率分析定性数据概率分析q例：根据甘蔗工业科学研例：根据甘蔗工业科学研究所用农药液剂、烟草粉究所用农药液剂、烟草粉肥皂液和烟脉肥皂液三种肥皂液和烟脉肥皂液三种药
3、液对棉蚜进行毒力测定药液对棉蚜进行毒力测定的数据的数据(王鉴明，王鉴明，1985)，进行机值分析并检验三种进行机值分析并检验三种毒力回归直线的平行性。毒力回归直线的平行性。q数据编辑格式：第数据编辑格式：第1列为农列为农药药剂药药剂(依次从依次从1，2， )编号、第编号、第2列为农药处理浓列为农药处理浓度，第度，第3、4列分别为供试列分别为供试虫数和死亡虫数。定义数虫数和死亡虫数。定义数据块。据块。.一、生物测定定义数据块定义数据块选择选择“专业统计专业统计生生物测定物测定计数型机值分析计数型机值分析”.一、生物测定n定性数据概率分析定性数据概率分析q系统提示确认浓度是否经过对数转换。如果已经
4、转换则按系统提示确认浓度是否经过对数转换。如果已经转换则按“Y”再再回车，若没有转换则直接回车。然后，根据系统提示输入待求的回车，若没有转换则直接回车。然后，根据系统提示输入待求的致死剂量效应指标，致死剂量可同时输入致死剂量效应指标，致死剂量可同时输入10、50、90 等，各值等，各值间用空格隔开，再按回车。间用空格隔开，再按回车。q如果直接回车，系统将取默认值如果直接回车，系统将取默认值LD50、LD95 。回车后可立即得到。回车后可立即得到回归独立曲线图和分析结果。回归独立曲线图和分析结果。.一、生物测定结果分析解释结果分析解释 拟合程度卡方值及其显著水平拟合程度卡方值及其显著水平。如显著
5、水。如显著水平大于平大于0.05，则表明所建立的毒力回归曲线，则表明所建立的毒力回归曲线模型是合适的，即否定所求模型异质性的假模型是合适的，即否定所求模型异质性的假设；反之，则认为所建立的模型不合适。设；反之，则认为所建立的模型不合适。本例中三种农药各自的本例中三种农药各自的概率分析模型的拟合度概率分析模型的拟合度以第一、二个较差，第以第一、二个较差，第三个较好，卡方检验三个较好，卡方检验x2值分别为值分别为36.63，140.35和和5.82。.一、生物测定 对多种药剂进行机值分析并对多种药剂进行机值分析并检验毒力回归直线平行性检验毒力回归直线平行性。可以。可以其方差分析结果推断：当方差分其
6、方差分析结果推断：当方差分析表中的析表中的F 值的显著水平小于等值的显著水平小于等于于0.05 时，表明毒力回归直线时，表明毒力回归直线互相不平行。互相不平行。本例中，害虫对三种农药反应的本例中，害虫对三种农药反应的差异极显著差异极显著(F =10.522，大于，大于Fa=0.01=8.05)，所以这三条回归，所以这三条回归线的斜率不同，即不存在共同斜线的斜率不同，即不存在共同斜率。率。 上述方差分析只能揭示毒力上述方差分析只能揭示毒力回归直线的总的趋势，而两种药回归直线的总的趋势，而两种药剂其致死剂量之间的差异，剂其致死剂量之间的差异，DPS系统则采用致死剂量比率测定方系统则采用致死剂量比率
7、测定方法来检验：当致死剂量比率的法来检验：当致死剂量比率的95%置信区间包含置信区间包含1，则两种药，则两种药剂的某致死剂量之间的差异不显剂的某致死剂量之间的差异不显著。著。 注意：有时因异质因子较大，注意：有时因异质因子较大，计算的置信限会溢出，以致于不计算的置信限会溢出，以致于不能给出其置信区间。能给出其置信区间。.一、生物测定n定量数据概率分析定量数据概率分析q数据编辑格式是每一个浓度数据编辑格式是每一个浓度在一行中按浓度和死亡在一行中按浓度和死亡(抑制抑制)率顺序填入。率顺序填入。q对多种药剂，且每种药剂有对多种药剂，且每种药剂有几个剂量水平处理的试验结几个剂量水平处理的试验结果进行机
8、值分析时，数据编果进行机值分析时，数据编辑格式为：第辑格式为：第1 列为药剂列为药剂(依依次从次从1，2， )编号，第编号，第2列列为药剂浓度，第为药剂浓度，第3 列为死亡列为死亡(抑制抑制)率。率。.一、生物测定n定量数据概率分析定量数据概率分析q定义数据块后，选择定义数据块后，选择“专业统计专业统计生物测定生物测定数量型数量型机值分析机值分析”，系统提示确认浓度是否已经过对数转换。，系统提示确认浓度是否已经过对数转换。如果已经转换则按如果已经转换则按“Y”再回车，若没有转换则直接回再回车，若没有转换则直接回车。此后，系统再提示用户输入待求的致死剂量效应车。此后，系统再提示用户输入待求的致死
9、剂量效应指标。致死剂量可同时输入指标。致死剂量可同时输入10、50、90 等，各值间等，各值间用空格隔开，再按回车。如果直接回车，系统将取默用空格隔开，再按回车。如果直接回车，系统将取默认值认值LD50。回车后瞬即得到分析结果。回车后瞬即得到分析结果。.一、生物测定分析结果解释分析结果解释 毒力回归方程方差分析毒力回归方程方差分析F 值及其显著水值及其显著水平：数量型数据毒力回归方程诊断采用方平：数量型数据毒力回归方程诊断采用方差分析方法进行，这是数量型资料和计数差分析方法进行，这是数量型资料和计数资料在生物测定分析进行模型诊断方面的资料在生物测定分析进行模型诊断方面的差异。当差异。当F 值的
10、显著水平小于等于值的显著水平小于等于0.05 时，时，表明了所求的毒力回归曲线是合适的；反表明了所求的毒力回归曲线是合适的；反之则表明所建立的模型不合适。之则表明所建立的模型不合适。 对多种药剂进行对多种药剂进行机值分析，机值分析，DPS系系统采用致死剂量比统采用致死剂量比率测定方法来检验率测定方法来检验两种药剂在致死剂两种药剂在致死剂量方面的差异：当量方面的差异：当致死剂量比率的致死剂量比率的95%置信区间包含置信区间包含1，则两种药剂的某致则两种药剂的某致死剂量之间的差异死剂量之间的差异不显著。不显著。.一、生物测定n时间时间-剂量剂量-死亡率模型分析死亡率模型分析q对于包含供试品剂量和受
11、试生物在试验期间死亡对于包含供试品剂量和受试生物在试验期间死亡(反应反应)率的生测率的生测数据分析，一般用概率分析的方法，即将累计死亡率进行概率值数据分析，一般用概率分析的方法，即将累计死亡率进行概率值转换后分别对剂量或时间作线性回归分析，通过回归参数的估计转换后分别对剂量或时间作线性回归分析，通过回归参数的估计而建立直线回归模型，从而估计剂量效应而建立直线回归模型，从而估计剂量效应(如如LD50)或时间效应或时间效应(如如LT50)。但这样的分析方法导致时间与剂量的效应相互分开，无法。但这样的分析方法导致时间与剂量的效应相互分开，无法使所建模型充分体现实验数据的完整性。因此，有必要将时间和使
12、所建模型充分体现实验数据的完整性。因此，有必要将时间和剂量的效应统一到同一个模型中来。时间剂量的效应统一到同一个模型中来。时间-剂量剂量-死亡率死亡率(time-dose-mortality，缩为，缩为TDM)模型就是这样的模型，按其数学结模型就是这样的模型，按其数学结构又称互补重对数模型构又称互补重对数模型(complementary log-log model)，简称，简称CLL 模型。模型。.一、生物测定n时间时间-剂量剂量-死亡率模型分析死亡率模型分析q如：用乙酰甲胺磷的如：用乙酰甲胺磷的5 个浓度处理西红柿上某害虫。个浓度处理西红柿上某害虫。 处理后试虫处理后试虫死亡过程缓慢，故适合
13、用死亡过程缓慢，故适合用TDM 模型分析数据。模型分析数据。q分析前，先将数据编辑定义。其数据编辑和存放格式为：数据块分析前，先将数据编辑定义。其数据编辑和存放格式为：数据块第一列为供试药剂的浓度，以后各列为用药后不同时间的试虫存第一列为供试药剂的浓度，以后各列为用药后不同时间的试虫存活数。活数。q然后执行然后执行“专业统计专业统计生物测定生物测定时间剂量死亡率模型时间剂量死亡率模型”功能功能项项.一、生物测定输出结果包括：输出结果包括：模型系数的协方差模型系数的协方差矩阵；矩阵； 所求系数即各参数所求系数即各参数的估计值、参数标准的估计值、参数标准误、误、t 测验值、显著水测验值、显著水平；
14、平；模型拟合值、残差、模型拟合值、残差、标准残差。标准残差。结果表列中有些行显示结果表列中有些行显示“数数据不合理，未进行处理据不合理，未进行处理”的的字样，这主要因为在最高剂字样，这主要因为在最高剂量下试虫在整个实验观察完量下试虫在整个实验观察完成前的数个时间区间已经全成前的数个时间区间已经全部死亡，即累计死亡率在前部死亡，即累计死亡率在前一个时间区间已达一个时间区间已达100%，此，此后各时间区间的条件死亡率后各时间区间的条件死亡率在数学上无法计算在数学上无法计算(0/0)。每。每当这种情况出现时，当这种情况出现时，DPS系系统自动将该数据视为无效，统自动将该数据视为无效，但截至前一时间区
15、间的拟合但截至前一时间区间的拟合仍然有效，并且不影响其余仍然有效，并且不影响其余各剂量数据的继续拟合。各剂量数据的继续拟合。结果中，如果皮尔逊卡方检结果中，如果皮尔逊卡方检验能通过，即其显著水平大验能通过，即其显著水平大于等于于等于0.05，则模型的拟合，则模型的拟合效果很好。而实际上，一般效果很好。而实际上，一般皮尔逊卡方检验不能通过皮尔逊卡方检验不能通过(显显著水平小于著水平小于0.05)，这表明模，这表明模型异质性的假设成立。但型异质性的假设成立。但Nowierski 等提出用等提出用Hosmer 和和Lemeshow 改进的方法进改进的方法进行拟合度测验，如行拟合度测验，如Hosmer
16、-Lemeshow 检验的卡方值的检验的卡方值的显著水平大于等于显著水平大于等于0.05，则，则可认为数据拟合是成功的。可认为数据拟合是成功的。.n新引入、新育成的品种或品系，在大面积推广之前，都需要经过区域新引入、新育成的品种或品系，在大面积推广之前，都需要经过区域试验，来鉴定其产量水平和适应性，为其推广和合理利用提供依据。试验，来鉴定其产量水平和适应性，为其推广和合理利用提供依据。n区域试验中需要研究的主要效应有：区域试验中需要研究的主要效应有：q 品种效应品种效应，即品种的产量或品质效应，因区域试验中供试品种是一定，即品种的产量或品质效应，因区域试验中供试品种是一定的，故品种效应是固定型
17、效应；的，故品种效应是固定型效应；q 地区效应地区效应，即地区之间土壤类型、耕作制度、生产水平和管理方法等，即地区之间土壤类型、耕作制度、生产水平和管理方法等可以预见的环境差异对品种的影响效应，一般亦属固定型效应；可以预见的环境差异对品种的影响效应，一般亦属固定型效应；q 年际效应年际效应，即不同年份的温度、雨量、偶然性灾害等难以预知的随机，即不同年份的温度、雨量、偶然性灾害等难以预知的随机性环境差异对品种的影响效应，一般属随机型效应；性环境差异对品种的影响效应，一般属随机型效应；q 品种品种地区的互作效应地区的互作效应，即研究品种对于可预知的环境差异是否具有，即研究品种对于可预知的环境差异是
18、否具有特殊的适应性，若品种与地区互作效应显著，说明品种对地区有特殊适应特殊的适应性，若品种与地区互作效应显著，说明品种对地区有特殊适应性，反之则说明适应性广泛，故该互作效应一般亦属固定型效应；性，反之则说明适应性广泛，故该互作效应一般亦属固定型效应；q 品种品种年份的互作效应年份的互作效应，反映品种对难以预测的环境差异是否有特殊，反映品种对难以预测的环境差异是否有特殊适应性，一般属随机型效应；适应性，一般属随机型效应；q 品种品种地点地点年份的互作效应年份的互作效应，反映品种与地点的互作效应是否随难，反映品种与地点的互作效应是否随难以预知的环境变化而改变，因此该效应属于随机型效应。以预知的环境
19、变化而改变，因此该效应属于随机型效应。二、品种区域试验.二、品种区域试验n评定一个品种的应用价值，主要考虑以下三个效应值：评定一个品种的应用价值，主要考虑以下三个效应值：q 品种效应；品种效应；q 品种品种地区点的互作效应及方差；地区点的互作效应及方差；q 品种品种年份的互作效应及方差。年份的互作效应及方差。q品种效应显著而互作效应小的品种是具有广泛适应性的丰产型品品种效应显著而互作效应小的品种是具有广泛适应性的丰产型品种，适于大面积推广；而互作效应显著的品种具有特殊适应性种，适于大面积推广；而互作效应显著的品种具有特殊适应性(如如对环境条件有特殊要求对环境条件有特殊要求)，只能在特定地区推广
20、才能发挥增产作用。，只能在特定地区推广才能发挥增产作用。.（一）一年多点试验稳定性分析n例如，对例如，对5个品种个品种5个地点的品种区个地点的品种区域试验结果进行稳定性分析，得到域试验结果进行稳定性分析，得到的结果如图。的结果如图。n分析时按图中的阴影部分编辑并用分析时按图中的阴影部分编辑并用鼠标选中待分析的数据。鼠标选中待分析的数据。n然后执行然后执行“专业统计专业统计品种区域试品种区域试验验 品种稳定性分析品种稳定性分析”功能，这功能，这时系统会要求你输入地点数、品种时系统会要求你输入地点数、品种数数(这里分别是这里分别是5、5)。输入后系统。输入后系统立即给出立即给出George C.
21、C. Tai 模型的模型的品种稳定性分析图。品种稳定性分析图。.（一）一年多点试验稳定性分析ni用以量测第用以量测第i 个品种对环境效应个品种对环境效应的线性响应，而的线性响应，而i是表示以误差均是表示以误差均方为单位的线性回归离差。方为单位的线性回归离差。n此图表示此图表示与与的分布图。图中两条的分布图。图中两条正交轴分别以正交轴分别以、值表示，如果值表示，如果i =0为真，那么将有为真，那么将有95的的 落在落在双曲线内。第一、二条垂直线是双曲线内。第一、二条垂直线是0=1时时90置信区间的上、下限，置信区间的上、下限，若途中还有其余垂直线则是若途中还有其余垂直线则是0 1时的置信上限。这
22、些直线把时的置信上限。这些直线把、区区分成不同的稳定性区域。分成不同的稳定性区域。n从图中各个品种分布位置可以看出，从图中各个品种分布位置可以看出，供试品种供试品种3 表现最稳定，其次是品表现最稳定，其次是品种种1 和品种和品种4；而品种；而品种2 的稳定性最的稳定性最差。差。a .（一）一年多点试验稳定性分析结束图形分析，返回到编辑结束图形分析，返回到编辑状态，则可得到几种方法的状态，则可得到几种方法的稳定性分析结果。稳定性分析结果。分析结果的第一部分是各个分析结果的第一部分是各个试验点误差均方的同质性检试验点误差均方的同质性检验。一般要求是卡方值较小，验。一般要求是卡方值较小，概率概率p
23、值大于值大于0.05。本例中。本例中的的p=0.7717，表明各个试验，表明各个试验点的误差均方是同质的，可点的误差均方是同质的，可以进行进一步分析。以进行进一步分析。随后，系统给出了各个试点、随后，系统给出了各个试点、品种的均值，并给出了各个品种的均值，并给出了各个试点参试品种的方差分析结试点参试品种的方差分析结果和几个试点的联合方差分果和几个试点的联合方差分析结果。从这里可以初步分析结果。从这里可以初步分析各个试点各个品种的表现析各个试点各个品种的表现情况，以及他们的综合表现情况，以及他们的综合表现情况。情况。.（二）一年多点区域试验的统计分析n例如，在例如，在4 个地点进行个地点进行6
24、个品种的区域品种比较试验，每个品种的区域品种比较试验，每个地点个地点2 个区组，获得一组试验资料。现对该试验资料进个区组，获得一组试验资料。现对该试验资料进行一年多点试验的方差分析，分析前先按下图方式输入、行一年多点试验的方差分析，分析前先按下图方式输入、编辑数据：编辑数据：.（二）一年多点区域试验的统计分析n或者，按如下格式编辑或者，按如下格式编辑定义数据块。定义数据块。n定义数据块后，选择定义数据块后，选择“专业统计专业统计品种区域品种区域试验试验一年多点一年多点”分析分析项，即可得到分析结果：项，即可得到分析结果：.（二）一年多点区域试验的统计分析.（三）多年多点品种区域试验的统计分析n
25、品种区域试验常常需要在多个地点连续进行数年，品种区域试验常常需要在多个地点连续进行数年，以便为品种的评定和应用提供更广泛的信息，尤以便为品种的评定和应用提供更广泛的信息，尤其是品种对不同年份随机变化的气候条件适应性其是品种对不同年份随机变化的气候条件适应性的信息。的信息。n按地点、品种、区组顺序输入试验数据并按规定按地点、品种、区组顺序输入试验数据并按规定格式编辑和定义数据块，系统执行运算后，将输格式编辑和定义数据块，系统执行运算后，将输出计算结果。输出结果包括双向列表、方差分析出计算结果。输出结果包括双向列表、方差分析表、表、f 统计量、多重比较和稳定性分析结果以及统计量、多重比较和稳定性分
26、析结果以及最后的综合评价。最后的综合评价。.（三）多年多点品种区域试验的统计分析n如有如有5 个玉米品种在个玉米品种在4 个地个地点连续进行了点连续进行了3 年试验，年试验，每次试验设每次试验设2 个区组。个区组。n试验数据按图格式输入并试验数据按图格式输入并编辑定义成数据块，注意编辑定义成数据块，注意地点、品种、年份、重复地点、品种、年份、重复的排列格式。的排列格式。n进入专业统计菜单，选择进入专业统计菜单，选择“品比试验品比试验多年多点多年多点”项，按提示输入各参数后项，按提示输入各参数后回车执行运算，分析结束回车执行运算，分析结束时输出如下结果：时输出如下结果：.（三）多年多点品种区域试
27、验的统计分析.（三）多年多点品种区域试验的统计分析.（四）品种区域试验AMMI 模型分析n作物品种区域试验旨在鉴定品种的丰产性、稳定性和适应性。参加区作物品种区域试验旨在鉴定品种的丰产性、稳定性和适应性。参加区试的品种在不同地点的产量往往表现不一致。这表明品种的基因型和试的品种在不同地点的产量往往表现不一致。这表明品种的基因型和环境互作环境互作(ge)效应的存在。以往对这种互作效应大多采用线性模型进效应的存在。以往对这种互作效应大多采用线性模型进行分析，但线性模型一般仅能解释很少一部分交互作用的变化。近年行分析，但线性模型一般仅能解释很少一部分交互作用的变化。近年来一种更为有效的来一种更为有效
28、的加性主效应乘积交互作用加性主效应乘积交互作用(AMMI)模型模型已开始被用于已开始被用于多年多点的区域试验资料的分析。该模型与方差分析模型、线性回归多年多点的区域试验资料的分析。该模型与方差分析模型、线性回归模型相比，应用范围更广且更有效。模型相比，应用范围更广且更有效。AMMI 模型的分析结果可以用直模型的分析结果可以用直观简洁的图形表达和解释。观简洁的图形表达和解释。AMMI 分析不仅可以对品种进行更可靠的分析不仅可以对品种进行更可靠的稳定性分析，而且可以鉴别一些具有特殊品种基因型和环境稳定性分析，而且可以鉴别一些具有特殊品种基因型和环境(ge)互作互作效应的基因型，为针对某一特殊环境的
29、特殊适应性的品种的育种提供效应的基因型，为针对某一特殊环境的特殊适应性的品种的育种提供有价值的信息。该模型的主要特点是将方差分析和主成分分析有机地有价值的信息。该模型的主要特点是将方差分析和主成分分析有机地结合在一起。结合在一起。.（四）品种区域试验AMMI 模型分析nDPS平台的数据格式及其编辑平台的数据格式及其编辑q数据存放格式要求以行代表环境数据存放格式要求以行代表环境(地点地点)，以列代表基因，以列代表基因型型(品种品种)，按试验地点的顺序，输完一个地点的数据后，按试验地点的顺序，输完一个地点的数据后再输下一个。再输下一个。q若一个品比试验有若一个品比试验有e 个地点、个地点、g 个品
30、种、个品种、r 个区组个区组(重复重复)，其资料整理排列顺序如下表。其资料整理排列顺序如下表。.（四）品种区域试验AMMI 模型分析.（四）品种区域试验AMMI 模型分析n应用示例应用示例q假设对假设对3 个品种在个品种在3 个地点、个地点、4 次重复的品种比较试验中次重复的品种比较试验中获得的资料，进行获得的资料，进行AMMI 模模型分析。型分析。q按右图格式输入数据并编辑按右图格式输入数据并编辑定义成数据块。定义成数据块。q选择专业统计选择专业统计品比试验品比试验有重复有重复AMMI 模型分析项，模型分析项，系统给出系统给出AMMI 模型分析图模型分析图形界面，其中包括双标图分形界面，其中
31、包括双标图分析析(图图1)、品种适应性分析、品种适应性分析(图图2)和各品种和各品种AMMI 排序分排序分析析(图图3)。.（四）品种区域试验AMMI 模型分析图图1解释解释AMMI1双标图的原则：双标图的原则：在水平方向，品种和地点的分布范围表明相应的主效应。在水平方向，品种和地点的分布范围表明相应的主效应。在垂直方向，如果以在垂直方向，如果以 IPCA1=0 作一条水平线，则品种与同在此水平线一作一条水平线，则品种与同在此水平线一侧的地点有正的互作，与另一侧的地点的互作则为负。侧的地点有正的互作，与另一侧的地点的互作则为负。靠近水平线的品种为较为稳定的品种。所以，高产稳产的品种对应最右边靠
32、近水平线的品种为较为稳定的品种。所以，高产稳产的品种对应最右边且靠近水平线的品种图标。且靠近水平线的品种图标。图图1 AMMI 双标图双标图 其纵、横轴都是平均产量、离差、其纵、横轴都是平均产量、离差、IPCA1（主成分分析）（主成分分析）、IPCA2 和和IPCA3，用户可以指定任意两个变量，作出双标图。双标图有助于揭示品种和环境之用户可以指定任意两个变量，作出双标图。双标图有助于揭示品种和环境之间的相互关系。间的相互关系。AMMI 双标图是以品种和地点的主效应值或平均产量为双标图是以品种和地点的主效应值或平均产量为x 坐标，坐标，y 坐标为品种坐标为品种和地点的和地点的IPCA1 值所作的
33、图形。这一双标图同时表达了品种的丰产性和值所作的图形。这一双标图同时表达了品种的丰产性和GE 交互作用的信息。交互作用的信息。.（四）品种区域试验AMMI 模型分析图图2n图图2 AMMI 品种适应图品种适应图q忽略地点主效应，用品忽略地点主效应，用品种主效应和种主效应和IPCA1 值估值估计品种的产量计品种的产量。q这里不需要考虑地点效这里不需要考虑地点效应是因为它只是将产量应是因为它只是将产量直线沿着产量坐标一块直线沿着产量坐标一块移上或移下，而我们只移上或移下，而我们只对每一地点所处最高的对每一地点所处最高的直线感兴趣。这样，每直线感兴趣。这样，每一品种对应一条直线，一品种对应一条直线，
34、在同一图上画出所有品在同一图上画出所有品种对应的直线，并在种对应的直线，并在x 轴轴内侧上标出的是参试地内侧上标出的是参试地点代码，则处于最上端点代码，则处于最上端的折线对应的品种即为的折线对应的品种即为有最佳适应性的品种。有最佳适应性的品种。q本例中，从产量水平角本例中，从产量水平角度来看，在参试地点环度来看，在参试地点环境范围内，境范围内，V1 和和V3 两两个品种表现较佳。个品种表现较佳。.（四）品种区域试验AMMI 模型分析图图3n图图3 AMMI 排序排序图图q以品种的平均产量为以品种的平均产量为x 轴，以轴，以IPCA1为为y 轴。轴。q当指定某一品种后，作出通过该点的两条曲线，这
35、时通过它们的空间位置可当指定某一品种后，作出通过该点的两条曲线，这时通过它们的空间位置可以分析其品种特性：即指定品种特性优于那些落在它左边的品种，如该图表以分析其品种特性：即指定品种特性优于那些落在它左边的品种，如该图表示品种示品种V1 好于好于V2。.（四）品种区域试验AMMI 模型分析n应用示例应用示例q图形分析完成后，获得如下输出结果：图形分析完成后，获得如下输出结果：n 平均产量两向列表平均产量两向列表(列表示品种，行表示地点列表示品种，行表示地点)；n 离差两向列表离差两向列表(列表示品种，行表示地点列表示品种，行表示地点)；n 各种方差分析表：品种试验方差分析表，品比试验结果各种方
36、差分析表：品种试验方差分析表，品比试验结果线性回归分析模型方差分析表，品比试验结果线性回归分析模型方差分析表，品比试验结果AMMI分析；分析；n 品种和环境各个变量的平均值、品种和环境各个变量的平均值、 离差、主成分值及特离差、主成分值及特征向量；征向量；n AMMI模型拟合两向列表模型拟合两向列表( 列表示品种，行表示地点环列表示品种，行表示地点环境境)。.（四）品种区域试验AMMI 模型分析n无区组品比试验的无区组品比试验的AMMI 模型分析模型分析q当参与区试的品种没有设置区组时，品比试验的当参与区试的品种没有设置区组时，品比试验的AMMI 模型模型分析步骤如下：分析步骤如下：n数据编辑
37、格式：以行代表环境数据编辑格式：以行代表环境(地点地点)，以列代表基因型，以列代表基因型(品种品种)，按试验地点的顺序，输完一个地点的数据后再输下一个地点的按试验地点的顺序，输完一个地点的数据后再输下一个地点的数据，在一个地点中依次输入各品种的数据。数据，在一个地点中依次输入各品种的数据。n定义数据块定义数据块.（四）品种区域试验AMMI 模型分析n无区组品比试验的无区组品比试验的AMMI 模型分析模型分析q当参与区试的品种没有设置区组时，品比试验的当参与区试的品种没有设置区组时，品比试验的AMMI 模型分析步骤如下：模型分析步骤如下：n进入专业统计菜单，选择进入专业统计菜单，选择“品种比较试
38、验品种比较试验无重复无重复AMMI 模型分析模型分析”功能项，按系统提示输入试验点数功能项，按系统提示输入试验点数(E)、试验、试验品种数品种数(G)，然后运行并输出分析结果，然后运行并输出分析结果n分析结果包括各方差分析表、品种和地点的两向列表、分析结果包括各方差分析表、品种和地点的两向列表、各个阶次各个阶次AMMI 模型预测估计值等。模型预测估计值等。.n王宏富（1962-），男，山西晋城人，教授，硕士生导师，作物栽培与耕作学系主任。中国耕作制度研究会理事，山西省质量技术监督局农业标准专家组成员，山西省作物学会理事。1984年毕业于山西农业大学，后留校任教，主要从事作物化控与逆境生理、农田杂草与防除、计算机在农业中的应用、农业结构调整等方面的教学与研究工作，为博士、硕士和本专科生主讲旱地农业专题、高级耕作学、农业传播技术与应用、农学概论、计算机在农业中的应用、杂草学、水土保持农学、无公害农产品生产技术等课程。 主持和参加国家“948”项目、教育部项目、省攻关项目、省自然基金项目、省开发项目与教学项目等多项，获省级研究成果一等奖1项、二等奖4项、三等奖1项；国家级出版社出版专著1部、21世纪教材1部（副主编）、国家统编教材1部（副主编）、十一五教材1部（副主编）、参编十一五教材1部，在国内外专业刊物上发表了学术论文30余篇。
