《方差分析⑵》PPT课件.ppt

试验统计学试验统计学第四章第四章概率论与数理统计的基础知识概率论与数理统计的基础知识本课程使用区靖祥编本课程使用区靖祥编著的著的试验统计学试验统计学一书作为课本。全程一书作为课本。全程为为50学时，占学分。学时，占学分。第二章第二章常用的试验设计常用的试验设计第三章第三章试验数据的整理试验数据的整理第五章第五章参数区间估计参数区间估计第八章第八章常用试验设计的资料分析常用试验设计的资料分析第六章第六章统计假设测验统计假设测验第七章第七章方差分析方差分析第九章第九章直线相关与回归直线相关与回归第一章第一章绪论绪论第十章第十章 协方差分析协方差分析第二节第二节 处理平均数间的多重比较处理平均数间的多重比较 第一节第一节 方差分析的基本原理方差分析的基本原理第三节第三节 方差分量的估计方差分量的估计 第七章第七章方差分析方差分析第四节第四节 单向分类资料的方差分析单向分类资料的方差分析第五节第五节 两向分类资料的方差分析两向分类资料的方差分析第六节第六节 系统分组资料的方差分析系统分组资料的方差分析第七节第七节 方差分析的基本假设和数据转换方差分析的基本假设和数据转换二、二、组内观察值数目不等的单向分类资料组内观察值数目不等的单向分类资料l单向分类资料单向分类资料(one-wayclassifications)是指那些只含有一个可控因素的资料。通常这个可控因素也就是考察因素。l根据具体的数据结构又分为两种情况：第四节第四节 单向分类资料的方差分析单向分类资料的方差分析一、一、组内观察值数目相等的单向分类资料组内观察值数目相等的单向分类资料二、二、每处理组合内有多于一个观察值的两向分类资料每处理组合内有多于一个观察值的两向分类资料l两向分类资料两向分类资料(two-wayclassifications)是指那些含有两个可控因素的资料。根据具体的数据结构又分为两种情况：第五节第五节 两向分类资料的方差分析两向分类资料的方差分析一、一、每处理组合内只有一个观察值的两向分类资料每处理组合内只有一个观察值的两向分类资料l在动物学的试验中，由于每一个母本不可能同时与若干个父本交配，所以不能采用交叉式设计的杂交方式，转而使用另一种被称为巢式设计巢式设计（nesteddesign）的杂交方案。巢式设计的杂交方案以及相类似的试验设计所取得到数据资料称为系统分组的系统分组的资料资料(hierarchalclassification)。第六节第六节 系统分类资料的方差分析系统分类资料的方差分析本节介绍对最简单的系统分组的资料的统计分析方法，这类资料为组内又分亚组的单向分类资料组内又分亚组的单向分类资料组内又分亚组的单向分类资料组内又分亚组的单向分类资料l方差分析的三个基本假定：方差分析的三个基本假定：第七节第七节 基本假设和数据转换基本假设和数据转换在结束本章之前，再稍微详细地讨论一下这些假定。数据中的各种效应应该具有“可加性可加性”；所有处理应该具有相同的误差方差误差方差，即具有“同质同质 性性”。误差应该是“随机、独立随机、独立”的，并且具有“平均数为平均数为0、方差为方差为 的正态分布的正态分布”；l方差分析的三个基本假定：第七节第七节 基本假设和数据转换基本假设和数据转换那么，怎么样的数据是可加性的呢？数据中的各种效应应该具有“可加性可加性”；线性可加模型是方差分析的基础，只有当数据具有可只有当数据具有可加性时，总平方和才能分解为各项平方和之和加性时，总平方和才能分解为各项平方和之和；并且，由由于数据具有可加性，又必然导致各项效应之和为于数据具有可加性，又必然导致各项效应之和为0 0，所有，所有误差之和为误差之和为0 0。以单向分类资料为例，因为数学模型为：，因此才有：和，即SST SSt SSe。也因此有。此式左边为0，因此右边自然也应该等于0。于是有和。表表7.63可加性资料与非可加性资料的比较可加性资料与非可加性资料的比较(a)可加性资料可加性资料(b)倍加性资料倍加性资料(c)对倍加性资对倍加性资料取对数后料取对数后处处理理分组分组处处理理分组分组处处理理分组分组121212A1020A1020A1.00 1.30B3040B3060B1.48 1.78l方差分析的三个基本假定：方差分析的三个基本假定：第七节第七节 基本假设和数据转换基本假设和数据转换误差应该是“随机、独立随机、独立”的，并且具有“平均数为平均数为0、方差为 的正态分布正态分布”；首先首先，在数学模型中的误差效应必须是随机的在数学模型中的误差效应必须是随机的，因为数据中的k 个处理仅仅是从所研究的k 个总体中随机抽取出来的k 个样本，而F 测验正是通过样本统计量对总体参数进行判断的手段。其次其次，在观察这个个体时的误差误差误差误差与观察另一个个体时的误差应该是无关误差应该是无关的，即误差彼此之间是相互独立的。前面已讨论过，计算F值的两个方差，所来自的（亚）总体应该是正态正态正态正态分布的。l方差分析的三个基本假定：方差分析的三个基本假定：第七节第七节 基本假设和数据转换基本假设和数据转换 因为在方差分布中将k个样本的“组内平方和”和“组内自由度”合并为整个试验的“组内平方和”和“组内自由度”，并利用它们算出的“组内均方”来估计试验误差，其前提必须是各处理的方差是相等的，不相等不相等不相等不相等怎么能合并呢？怎么能合并呢？怎么能合并呢？怎么能合并呢？资料中各组的方差是否相等可以通过BartlettBartlett卡方测验卡方测验卡方测验卡方测验来检验。所有处理应该具有相同的误差方差，即具有“同质性同质性”。l当试验资料不符合上述假定时要先对数据进行一些适当的处理，然后用经过处理的数据进行方差分析。第七节第七节 基本假设和数据转换基本假设和数据转换常用的数据变换方法有四种，分别适用下述的四种情况：剔除一些表现“特殊”的观察值、处理或重复。用同一个体或小区的重复观察值的平均数进行方差 分析。对需要分析的资料进行研究，了解它们不符合哪个 基本假定，然后针对性地采用下述数据转换方法中 的一种，先对数据进行某种尺度变换，用经变换的 数据进行方差分析(及多重比较)，而在对分析结果 进行解释时，再反代换为原来的尺度。第七节第七节 基本假设和数据转换基本假设和数据转换当数据的 时，各观察值的方差近似与其平 均数成比例关系：即平均数越大，方差越大。这时 宜采用平方根转换，即当有部分观察值小于10时，采用公式：或服从泊松分布泊松分布泊松分布泊松分布的计数资料宜采用这种转换。通常认为计数资料，如每一个显微镜视野中的细菌数、每土方中的昆虫幼蛹数等，都服从泊松分布。第七节第七节 基本假设和数据转换基本假设和数据转换服从二项分布的百分数二项分布的百分数二项分布的百分数二项分布的百分数资料宜采用这种变换。例如昆虫死亡率、产品合格率、种子发芽率等都属于这种资料。请注意上式方根号内xijk是百分数，例如当xijk75.5%时，xijk60.33。当数据的时，宜采用反正弦转换，即。有些百分率不是用计数数据换算得到的，而是用度量数据换算得到的，如玉米种子中的蛋白质含量、花生油中的不饱和脂肪酸含量等，不应当作离散性随机变量看待。不必要采用反正弦代换。第七节第七节 基本假设和数据转换基本假设和数据转换当资料不近似服从正态分布或资料不接近于加性模型时，宜采用这种转换。当数据的 时，各观察值的方差近似与其 平均数的平方成比例关系。这种资料宜采用对数转 换。即当有部分数据小于10时，采用公式：第七节第七节 基本假设和数据转换基本假设和数据转换 当数据的 时，各观察值的标准差近似与其 平均数的平方成比例关系。这种资料宜采用倒数转 换，即 ，当有部分数据小于10时，采用公式。第七节第七节 基本假设和数据转换基本假设和数据转换l 如果能事先把您的离散型数据的各处理平均数和方差计算出来，看看它们之间的关系，就可以很容易地决定采用何种尺度变换方法。如果无法决定那种变换更加合适时，可以只变换少数几种处理的数据（这种处理的原数据最好能包括有大、中、小的数字），然后再看看经变换的数据的平均数与方差间有没有上述的关联性，这些关联性最小的那一种变换往往就是最好的变换方法。五个病人在三种不同室温下淋巴细胞玫瑰花簇玫瑰花簇形成率（按200个淋巴细胞计算）如所示。试测验不同室温下这种形成率之间是否有显著差异。黄成达编制黄成达编制31.255.235.415.055.534.5516.065.534.5440.049.034.0336.058.034.0249.048.040.013037oC2025oC46oC室温病人编号表表7.64不同室温下淋巴细胞玫瑰花簇形成率（不同室温下淋巴细胞玫瑰花簇形成率（%）黄成达编制黄成达编制黄成达编制黄成达编制46oC54321病人编号表表7.65经反正弦代换后的数据经反正弦代换后的数据()33.378648.014336.501935.638222.786548.157735.970437.859423.578254.029635.970439.775739.231544.427035.668540.714036.869949.603435.668542.504144.427043.853839.23153037oC2025oC室温黄成达编制黄成达编制黄成达编制黄成达编制1049.223814总变异总变异46.3924371.13918误差误差0.452320.982183.92834病号间病号间8.64914.45906.4036*297.0782594.15642室温间室温间F0.01F0.05FMSSSdf变异变异表表7.66例例7.14的方差分析表的方差分析表黄成达编制15.230314.43835.004.74310.32619.93023.393.262LSR0.01LSR0.05SSR0.01SSR0.05g表表7.67查得的查得的SSR值及各种值及各种LSR值的计算值的计算表表7.687.68 对不同室温的进行比较的梯形表对不同室温的进行比较的梯形表 室温 33.37863037oC3.123336.501946oC11.5124*14.6357*48.01432025oC平均数黄成达编制人类人类T淋巴细胞表面具有绵羊红细胞淋巴细胞表面具有绵羊红细胞(SRBC)的的受体，故能与受体，故能与SRBC结合形成花环，目前常用结合形成花环，目前常用E花花环；目前常用环；目前常用E花环形成试验来检测花环形成试验来检测T淋巴细胞数淋巴细胞数量，作为了解人体细胞免疫功能状态的指标之一，量，作为了解人体细胞免疫功能状态的指标之一，总总E花环（花环（Et）试验应用最广，反映）试验应用最广，反映T细胞活性的细胞活性的是活性花环试验是活性花环试验(Ea)。E-花环试验作为体外测定细胞免疫方法之一，花环试验作为体外测定细胞免疫方法之一，其形成百分率下降是细胞免疫功能降低的反映，临其形成百分率下降是细胞免疫功能降低的反映，临床上已用于恶性肿瘤、白血病、自身免疫性疾病、床上已用于恶性肿瘤、白血病、自身免疫性疾病、免疫缺陷病、器官移植排斥反应等方面的研究，作免疫缺陷病、器官移植排斥反应等方面的研究，作为了解这些疾病患者机体细胞免疫状态、预后及疗为了解这些疾病患者机体细胞免疫状态、预后及疗效观察的指标。效观察的指标。T淋巴细胞玫瑰花环：淋巴细胞玫瑰花环：左左.（甲紫染色，（甲紫染色，800)：可见：可见2个被羊红细胞个被羊红细胞包围而形成玫瑰环的包围而形成玫瑰环的T淋巴细胞；淋巴细胞；右右.（吖啶橙染色）：图中可见（吖啶橙染色）：图中可见3个染成橙黄色荧光的个染成橙黄色荧光的T淋巴细胞被羊红细胞所包围。淋巴细胞被羊红细胞所包围。淋淋巴巴细细胞胞玫玫瑰瑰花花环环电电镜镜照照片片一、组内观察值数目相等的单向分类资料一、组内观察值数目相等的单向分类资料如果资料中含有如果资料中含有k组数据，每组含组数据，每组含n个观察值，个观察值，全部共有全部共有nk个观察值。那么，此类资料观察值的数个观察值。那么，此类资料观察值的数学模型为：学模型为：（i1，2，k，j1，2，n）因此总变异分解为两大部分：因此总变异分解为两大部分：组间变异组间变异和和组内变异组内变异。方差分析表如表所示。方差分析表如表所示。黄成达编制黄成达编制dfTnk1总变总变异异MSedfek(n1)误误差差FMSt/MSeMStdftk1组间组间随机模型随机模型固定模型固定模型期望均方期望均方(EMS)F值值均方均方平方和平方和自由度自由度变变异异来源来源表表7.14组组内内观观察察值值数目相等的数目相等的单单向分向分类资类资料的方差分析表料的方差分析表事实上，在前几节中所举的例子都是这种资料。事实上，在前几节中所举的例子都是这种资料。这里，我们另举一个例子，并且把整个计算过程系这里，我们另举一个例子，并且把整个计算过程系统地陈述一遍。以后在各种不同试验设计的分析方统地陈述一遍。以后在各种不同试验设计的分析方法中，也都将按这样的方式进行陈述。法中，也都将按这样的方式进行陈述。例例在一个塑料大棚内进行番茄无土栽培试验，全部采用同一品种，5种不同的培养液，每种培养液观察4株。试验指标为单株产量，结果如表所示，如果培养液的效应为固定模型，试对5种培养液的效应进行显著性测验。数据整理：数据整理：黄成达编制黄成达编制黄成达编制265520013898520总总和和20640016228021162221E3115376386812427333331D2812544315411230252730C24921623189626212425B2711664293610826283024A观察值观察值培养液培养液表表7.15番茄无土栽培试验的产量数据番茄无土栽培试验的产量数据平方和及自由度的分解：平方和及自由度的分解：总自由度总自由度dfT观察值总数观察值总数120119处理间自由度处理间自由度dft处理数处理数1514误差自由度误差自由度dfedfTdft19415矫正数矫正数C.T.总和平方观察值总数总和平方观察值总数5202/2013520总平方和总平方和=13898-13520=378总平方和总平方和=各观察值的平方之和各观察值的平方之和C.T.处理间平方和处理间平方和=55200/4-13520=280误差平方和误差平方和37828098误差平方和误差平方和=总平方和总平方和各项已知因素的平方和各项已知因素的平方和将各项自由度和平方和填入方差分析表。将各项自由度和平方和填入方差分析表。黄成达编制37819总变异总变异6.5339815误差误差4.893.0610.72*702804处理间处理间F0.01F0.05F均方均方平方和平方和自由度自由度变异来源变异来源表表7.16例例7.7的方差分析表的方差分析表多重比较（以多重比较（以Duncan法为例）：法为例）：计算标准误：计算标准误：其中其中MSe为方差分析表中的误差均方；为方差分析表中的误差均方；n 是计算所比较的平均数时用到的观察值数目。是计算所比较的平均数时用到的观察值数目。查表并计算各种LSR值（误差自由度dfe=15）：黄成达编制5.854.234.583.3155.754.154.503.2545.584.044.373.1635.323.854.173.012LSR0.01LSR0.05SSR0.01SSR0.05g表表7.17查得的查得的SSR值及各种值及各种LSR值的计算值的计算3.列梯形表进行比较：列梯形表进行比较：黄成达编制5.854.234.583.3155.754.154.503.2545.584.044.373.1635.323.854.173.012LSR0.01LSR0.05SSR0.01SSR0.05g表表7.17查得的查得的SSR值及计算得到的值及计算得到的LSR值值黄成达编制黄成达编制20E4*24B37*27A148*28C347*11*31D平均数平均数处理处理表表7.18多重比较的梯形表多重比较的梯形表作统计推断：作统计推断：黄成达编制黄成达编制20E4*24B37*27A148*28C347*11*31D平均数平均数处理处理表表7.18多重比较的梯形表多重比较的梯形表培养液培养液D与培养液与培养液E、B之间有显著差异之间有显著差异(=0.01)；培养液培养液C、A与培养液与培养液E之间有显著差异之间有显著差异(=0.01)；培养液培养液B与培养液与培养液E之间有显著差异之间有显著差异(=0.05)；。其余各处理平均数之间没有显著差异。其余各处理平均数之间没有显著差异。二、组内观察值数目不等的单向分类资料二、组内观察值数目不等的单向分类资料如果资料中含有如果资料中含有k组数据，其中第组数据，其中第i组含组含ni个观个观察值，全部共有个观察值。此类资料观察值的数学察值，全部共有个观察值。此类资料观察值的数学模型为：模型为：（i1，2，k，j 1，2，ni）因此总变异分解为两大部分：因此总变异分解为两大部分：组间变异组间变异和和组内变异组内变异。方差分析表如表所示方差分析表如表所示dftk-1黄成达编制黄成达编制黄成达编制SSTdfT总变总变异异MSeSSedfe误误差差FMSt/MSeMStSStdft组间组间随机模型随机模型固定模型固定模型期望均方期望均方(EMS)F值值MSSSdf变变异异来源来源表表7.19组组内内观观察察值值数目不等的数目不等的单单向分向分类资类资料的方差分析表料的方差分析表其中其中为各组样本含量的调和平均数，即为各组样本含量的调和平均数，即但是如果利用随机模型对但是如果利用随机模型对进行估计，即进行估计，即时，可能会出现方差估计值为负的奇怪现象。这时时，可能会出现方差估计值为负的奇怪现象。这时可以使用可以使用代替代替进行计算。进行计算。各组内观察值数目相等的试验资料称为各组内观察值数目相等的试验资料称为平衡资平衡资料料；各组内观察值数目不等的试验资料称为；各组内观察值数目不等的试验资料称为不平不平衡资料衡资料。利用平衡资料所得到的参数估计值是无。利用平衡资料所得到的参数估计值是无偏的，利用不平衡资料所得到的参数估计值不是偏的，利用不平衡资料所得到的参数估计值不是无偏的。一般地说，如果是经过精心设计来安排无偏的。一般地说，如果是经过精心设计来安排试验的话，应尽可能安排平衡试验，避免采用不试验的话，应尽可能安排平衡试验，避免采用不平衡试验设计。如果是进行一些野外调查，由于平衡试验设计。如果是进行一些野外调查，由于特殊条件的要求或限制，不得不采用不平衡设计特殊条件的要求或限制，不得不采用不平衡设计时，仍应特别注意估计结果的有偏性。时，仍应特别注意估计结果的有偏性。例例7.8调查四个水库中的氯离子浓度，按水库的调查四个水库中的氯离子浓度，按水库的水面面积大小取不同数目的样本。数据如表所示，水面面积大小取不同数目的样本。数据如表所示，如果各水库的效应为固定模型，试测验各水库间氯如果各水库的效应为固定模型，试测验各水库间氯离子的浓度之间是否有显著差异。离子的浓度之间是否有显著差异。表表7.20四个水库氯离子浓度的数据资料四个水库氯离子浓度的数据资料水水库库观观察察值值1213141514161491088910992101112131211141011121312黄成达编制黄成达编制3603371131328总和128648747261080088480895675716371413721382987ni水库水库表表7.20四个水库氯离子浓度的数据资料的整理四个水库氯离子浓度的数据资料的整理数据处理：数据处理：平方和及自由度的分解：平方和及自由度的分解：总自由度总自由度dfT=观察值总数观察值总数1=281=27处理间自由度处理间自由度dft=处理数处理数1=41=3误差自由度误差自由度dfe=dfTdft=273=24矫正数矫正数C.T.=总和平方总和平方/观察值总数观察值总数=3132总平方和总平方和=37113498.8929=212.1071处理间平方和处理间平方和36033498.8929104.1071误差平方和误差平方和=212.1071104.1071=108列方差分析表：列方差分析表：表表7.21例例7.8的方差分析表的方差分析表变因变因自由度自由度平方和平方和均方均方FF0.05F0.01处理间处理间3104.1071 34.70247.712*3.014.72误差误差24108.00004.5000总变异总变异27212.1071多重比较（以多重比较（以Duncan法为例）：法为例）：先计算样本容量的调和平均数先计算样本容量的调和平均数或或：计算标准误：计算标准误：表表7.22查得的查得的SSR值及各种值及各种LSR值的计算值的计算gSSR0.05SSR0.01LSR0.05LSR0.0122.923.962.3533.19233.074.112.4743.33743.154.242.5393.417查表并计算各种查表并计算各种LSR值（误差自由度值（误差自由度dfe=24）如表所示：如表所示：a.第第号水库与第号水库与第、号水库的氯离子浓度有极显著差异；号水库的氯离子浓度有极显著差异；b.第第号水库与第号水库与第号水库的氯离子浓度有显著差异；号水库的氯离子浓度有显著差异；c.其余各水库的氯离子浓度之间没有显著差异。其余各水库的氯离子浓度之间没有显著差异。列梯形表进行比较：列梯形表进行比较：黄成达编制黄成达编制黄成达编制911023*1224*5*14平均数平均数处理处理表表7.23例例7.8的多重比较梯形表的多重比较梯形表作统计推断：作统计推断：一、每处理组合内只有一个观察值的两向分一、每处理组合内只有一个观察值的两向分类资料类资料如果资料中有两个可控因素如果资料中有两个可控因素A和和B，其中，其中A有有a个水平，个水平，B有有b个水平，于是共有个水平，于是共有ab个处理组合。个处理组合。每个处理组合含每个处理组合含1个观察值，全部共有个观察值，全部共有ab个观察值。个观察值。为了陈述上的方便，将数据结构列出如表所示。为了陈述上的方便，将数据结构列出如表所示。其中处理组合其中处理组合AiBj的观察值记为的观察值记为xij；第；第i个个A水水平的观察值之和及平均数分别记为平的观察值之和及平均数分别记为，第，第j个个B水平的观察值之和及平均数分别记为水平的观察值之和及平均数分别记为，全部观察值的总和及平均数分别记为全部观察值的总和及平均数分别记为黄成达编制黄成达编制黄成达编制黄成达编制黄成达编制xabxajxa2xa1Aaxibxijxi2xi1Aix2bx2jx22x21A2x1bx1jx12x11A1BbBjB2B1B因素因素(j)A因素因素(i)表表7.24每处理组合只含一个观察值的两向分类资料每处理组合只含一个观察值的两向分类资料此类资料观察值的数学模型为：此类资料观察值的数学模型为：(i1，2，a；j1，2，b)因此总变异分解为三大部分：因此总变异分解为三大部分：A因素各水平间的变异因素各水平间的变异B因素各水平间的变异因素各水平间的变异剩余的变异剩余的变异(即误差变异即误差变异)。方差分析表如表方差分析表如表所示。所示。表表7.25每每处处理理组组合内只有一个合内只有一个观观察察值值的的两向分两向分类资类资料的方差分析表料的方差分析表变变异异来源来源自由度自由度平方和平方和均方均方期望均方期望均方(EMS)固定模型固定模型随机模型随机模型混合模型混合模型A间间dfASSAMSAB间间dfBSSBMSB误误差差DfeSSeMSe总变总变异异dfT SSTdfAa1dfBb1dfe(a1)(b1)dfT ab1表中列出了固定模型和随机模型的期望均方，如果两个处表中列出了固定模型和随机模型的期望均方，如果两个处理效应中有一个理效应中有一个(如如A)为固定，另一个为固定，另一个(如如B)为随机，则其为随机，则其模型称为模型称为(A固定固定B随机的随机的)混合模型混合模型(mixedmodel)。对各效。对各效应进行应进行F测验时采用的测验时采用的F值计算公式应视各项均方的期望值计算公式应视各项均方的期望值而定，例如在固定模型中，测验值而定，例如在固定模型中，测验A间差异时，间差异时，FMSA/MSe，测验，测验A间差异时，间差异时，FMSB/MSe。例例六个水稻品种六个水稻品种(A1、A2、A3、A4、A5和和A6)栽植在四种不同的土壤类型栽植在四种不同的土壤类型(B1、B2、B3和和B4)中，中，产量数据如表所示，如果品种和土壤类型都是固产量数据如表所示，如果品种和土壤类型都是固定效应，试对资料进行适当的分析。定效应，试对资料进行适当的分析。黄成达编制黄成达编制黄成达编制6361.563.065.062.5151236937839037562.5063.2563.5061.7566.2560.75250253254247265243表表7.26例例7.9的产量资料及数据整理的产量资料及数据整理59626564A664626562A561626764A460636460A365686765A260616260A1B4B3B2B1土壤类型（土壤类型（B）品种品种（A）数据处理：如表所示数据处理：如表所示平方和及自由度的分解：平方和及自由度的分解：总自由度总自由度dfT=观察值总数观察值总数1=241=23品种间自由度品种间自由度dfA=品种数品种数1=61=5土壤类型间自由度土壤类型间自由度dfB=土壤类型数土壤类型数1=41=3误差自由度误差自由度dfe=dfTdfAdfB=2353=15误差平方和误差平方和=1427139=32列方差分析表：列方差分析表：变异来源变异来源自由度自由度平方和平方和均方均方FF0.05F0.01品种间品种间57114.200 6.656*2.901 4.556土壤类型间土壤类型间33913.000 6.094*3.287 5.417误差误差15322.133总变异总变异23142表表例的方差分析表例的方差分析表统计推断：统计推断：品种间差异极显著品种间差异极显著；土壤类型间差异；土壤类型间差异也极显著也极显著。多重比较多重比较在本例中，品种和土壤类型都是固定效应；在本例中，品种和土壤类型都是固定效应；F测验表明，品种间差异极显著，因此要对不同品种测验表明，品种间差异极显著，因此要对不同品种的平均数进行多重比较；同时土壤类型间差异也极的平均数进行多重比较；同时土壤类型间差异也极显著，如果研究目的要求对土壤差异进行分析，也显著，如果研究目的要求对土壤差异进行分析，也应对它进行多重比较。现以应对它进行多重比较。现以Duncan法为例说明之。法为例说明之。对对品种间差异品种间差异的多重比较：的多重比较：计算标准误：计算标准误：其中b为每品种的观察值数目，在本题为土壤类型种类数。查查Duncan表得表得SSR值并计算各种值并计算各种LSR值值(误差自由度dfe=15)如表。表表7.28查得的查得的SSR值及各种值及各种LSR值的计算值的计算gSSR0.05SSR0.01LSR0.05LSR0.0123.014.172.203.0533.164.372.313.1943.254.502.373.2953.314.582.423.3463.364.642.453.39列梯形表进行比较，如表所示。列梯形表进行比较，如表所示。60.75A11.0061.75A30.751.7562.50A60.751.502.50*63.25A50.251.001.752.75*63.50A42.75*3.00*3.75*4.50*5.50*66.25A2平均数平均数品种品种表表7.29对品种平均数的进行比较的梯形表对品种平均数的进行比较的梯形表根据多重比较做出统计推断：根据多重比较做出统计推断：a.品种品种A2的产量最高，并且与其它所有品的产量最高，并且与其它所有品b.种之间都有显著或极显著差异；种之间都有显著或极显著差异；b.品种品种A4与品种与品种A2、品种、品种A1间有显著差异；间有显著差异；c.品种品种A5与品种与品种A1间有显著差异；间有显著差异；d.其余各品种之间没有显著差异。其余各品种之间没有显著差异。对对土壤类型间差异土壤类型间差异的多重比较：的多重比较：计算标准误：计算标准误：其中a为每土壤类型的观察值数目，在本题为品种的数目。查查Duncan表得表得SSR值并计算各种值并计算各种LSR值值(误差自由度dfe=15)如表。表表7.30查得的查得的SSR值及各种值及各种LSR值的计算值的计算gSSR0.05SSR0.01LSR0.05LSR0.0123.014.171.792.4933.164.371.882.6143.254.501.942.68列梯形表进行比较，如表所示。列梯形表进行比较，如表所示。61.5B4162.5B10.51.563.0B32*2.5*3.5*65.0B2平均数平均数土壤类型土壤类型表表7.31对土壤类型平均数进行比较的梯形表对土壤类型平均数进行比较的梯形表根据多重比较做出统计推断：根据多重比较做出统计推断：a.土壤类型土壤类型B2与其它所有类型之间都有显著与其它所有类型之间都有显著b.或极显著差异；或极显著差异；b.其余土壤类型之间没有显著差异。其余土壤类型之间没有显著差异。如果土壤类型虽然是可控因素，但不是考察因如果土壤类型虽然是可控因素，但不是考察因素，那么，尽管土壤类型之间的素，那么，尽管土壤类型之间的F测验显著也可以测验显著也可以不进行多重比较（因为你对它的差异与否不感兴趣）不进行多重比较（因为你对它的差异与否不感兴趣）。二、每处理组合内有多于一个观察值的两向二、每处理组合内有多于一个观察值的两向分类资料分类资料如果资料中有两个可控因素如果资料中有两个可控因素A和和B，其中，其中A有有a个水平，个水平，B有有b个水平，于是共有个水平，于是共有ab个个处理组合处理组合。每个处理组合含每个处理组合含n个观察值，全部共有个观察值，全部共有abn个观察个观察值。值。当每组合含有多于一个观察值时，可以通过方当每组合含有多于一个观察值时，可以通过方差分析考察因素间的交互作用，我们记因素差分析考察因素间的交互作用，我们记因素A和因和因素素B之间的交互作用为之间的交互作用为(AB)。xabnxab2xab1xajnxaj2xaj1xa2nxa22xa11Aaxibnxib2xib1xijnxij2xij1xi2nxi22xi11Aix2bnx2b2x2b1x2jnx2j2x2j1x22nx222x211A2x1bnx1b2x1b1x1jnx1j2x1j1x11nx112x111A1BbBjB1B因素因素(j)A因素因素(i)表表7.32每每处处理理组组合内含有合内含有n个个观观察察值值的两向分的两向分类资类资料料为了陈述上的方便，将数据结构列于表。为了陈述上的方便，将数据结构列于表。xijk处理组合处理组合AiBj的第的第k个观察值记；个观察值记；上表其中：上表其中：A因素第因素第i个水平的观察值总和；个水平的观察值总和；B因素第因素第j个水平的观察值的平均数；个水平的观察值的平均数；B因素第因素第j个水平的观察值总和；个水平的观察值总和；A因素第因素第i个水平的观察值的平均数；个水平的观察值的平均数；全部观察值的部和数；全部观察值的部和数；全部观察值的平均数。全部观察值的平均数。此类资料观察值的数学模型为：此类资料观察值的数学模型为：（i1，2，a；j1，2，b；k1，2，n）从以上数学模型看出，对于这类资料的分析可从以上数学模型看出，对于这类资料的分析可以分两步进行：以分两步进行：先按模型先按模型将总变异分解为处将总变异分解为处理理(组合组合)间变异和误差变异；这是一个单向分类的间变异和误差变异；这是一个单向分类的方差分析。方差分析表如表方差分析。方差分析表如表7.33所示。如果这一步所示。如果这一步的的F测验不显著，就不必进行第二步。分析就此结测验不显著，就不必进行第二步。分析就此结束。束。表表7.33第一步的方差分析表：第一步的方差分析表：将将总变总变异分解异分解为处为处理理(组组合合)间变间变异和异和误误差差变变异异变变异来源异来源自由度自由度平方和平方和均方均方F值值处处理理间间dftab1SStMStMSt/MSe误误差差dfeab(n1)SSeMSe总变总变异异dfTabn1SST 如果第一步的如果第一步的F测验显著，应进一步按模型测验显著，应进一步按模型将处理变异分解为将处理变异分解为A因素各水因素各水平间的变异、平间的变异、B因素各水平间的变异和交互作用因素各水平间的变异和交互作用(AB)引起的变异。引起的变异。该步骤的方差分析表如表所示，在第二步中的该步骤的方差分析表如表所示，在第二步中的F测验值如何计算需要视研究目的所确定的模型而测验值如何计算需要视研究目的所确定的模型而定。定。表表7.34第二步的方差分析表：第二步的方差分析表：将将处处理理间变间变异分解异分解为为A间变间变异、异、B间变间变异和交互作用异和交互作用变变异来异来源源自由度自由度平方和平方和均方均方F值值A间间dfAa1SSAMSA视视模模型而型而确定确定B间间dfBb1SSBMSBAB互作互作 dfAB(a1)(b1)SSABMSAB处处理理间间dft ab1SStMSt将以上两步合并，得到总的方差分析表，如表将以上两步合并，得到总的方差分析表，如表所示。表中列出了固定模型和随机模型的期望均方。所示。表中列出了固定模型和随机模型的期望均方。如果两个处理效应中有一个如果两个处理效应中有一个(例如例如A)为为固定模型固定模型，另一个另一个(例如例如B)为为随机模型随机模型，这时的模型称为，这时的模型称为(A固固定定B随机的随机的)混合模型混合模型(mixedmodel)。如表也列出。如表也列出了了(A固定固定B随机的随机的)混合模型的期望均方。混合模型的期望均方。SSTdfT 总变总变异异MSeSSedfe误误差差MSABSSABdfABAB互作互作MSBSSBdfBB间间MSASSAdfAA间间视视模模型型而而定定MStSStdft 处处理理间间A固定固定B随机随机随机模型随机模型固定模型固定模型期望均方期望均方(EMS)F值值均方均方平平方方和和自自由由度度变变异异来源来源表表7.35每每处处理理组组合内含有多于一个合内含有多于一个观观察察值值的的两向分两向分类资类资料的料的总总方差分析表方差分析表21340例例三个水稻品种三个水稻品种(A1、A2、A3)种在四种不同的种在四种不同的土壤类型土壤类型(B1、B2、B3和和B4)中，每个组合种了两个中，每个组合种了两个小区。产量数据已经过简化，列于表。试对资料进小区。产量数据已经过简化，列于表。试对资料进行适当的分析。行适当的分析。数据处理：数据处理：5.258.55.58.04.06.07.03.53.55.04.03.54.5146228912125664144196494910064498175014565130347298252952342541126171116812147710879847567354534979367426345观观察察值值B4B3B2B1B4B3B2B1B4B3B2B1合合计计A3A2A1处处理理组组合合表表7.36例例7.10的的产产量量资资料及数据整理料及数据整理为了计算因素间的交互作用，将各处理组合之为了计算因素间的交互作用，将各处理组合之和列成和列成AB二向表如表所示。二向表如表所示。6.505.004.25表表7.37AB二向表二向表6.55.55.04.0408615211089900576146253338135419454601462126393330242704730521711168A3160043840121477A211562943410879A1B4B3B2B1土壤土壤类类型（型（B）品种品种（A）第一步：将总变异分解为处理间变异和误差变第一步：将总变异分解为处理间变异和误差变异两部分异两部分全试验中，有全试验