详细方差分析PPT课件

.,1,第三章方差分析Chapter3ANOVA(AnalysisofVariance),.,2,方差分析是判断多组数据（K3）之间平均数差异是否显著的一种假设测验方法。2个样本平均数可用t或U测验的方法来评定其差数的显著性。如果有K个平均数，且K3，若仍然用两两比较的方法来测验，则需要作K(K-1)/2次测验，如果K10，则需要45次测验，不但测验程序繁琐，而且在理论上，其显著水平已经扩大了。因此，对于多样本平均数的假设测验，需采用一种更为合适的统计方法，即方差分析法（Fisher,1923）。,第三章方差分析,方差是平方和除以自由度的商。,.,3,方差分析是将总变异分裂为各个因素的相应变异，作出其数量估计，从而发现各个因素在变异中所占的重要程度，而且除了可控制因素所引起的变异后，其剩余变异又可提供试验误差的准确而无偏的估计，作为统计假设测验的依据。,第三章方差分析,例如，若有5组数据要比较，则共需要比较(54)/2=10次。若H0正确，每次接受的概率为10.95，10次都接受的概率为0.95100.60，因此，=10.600.40，即犯第一类错误的概率为0.40，这显然是不能接受的。,.,4,本章主要内容：第一节方差分析的基本原理和方法。第二节单向分组资料的方差分析。第三节两向分组资料的方差分析。,第三章方差分析,.,5,第一节方差分析的基本原理和方法,1.自由度和平方和的分解2.F分布(FDistribution)3.多重比较(multiplecomparisons)4.方差分析的基本假定5.数据转换,第三章方差分析,.,6,1、自由度和平方和的分解设有K组样本，每样本均具有n个观察值，则该资料共有nk个观察值，数据如下表。,表每组具n个观察值的k组样本的符号表,第一节方差分析的基本原理和方法,Xij，i=1,2,k，j=1,2,n。,.,7,总变异是nk个观察值的变异，故其自由度为nk1，平方和SST为：,式中，C称为矫正数。,总平方和(SST),第一节方差分析的基本原理和方法,.,8,总平方和SST组内平方和SSe处理平方和SSt,总平方和SST的计算：,.,9,组内的变异为各组内观察值与组平均数的相差，故每组具有n1个自由度，平方和为，而总共有k组资料，故组内自由度为k（n1），而组内平方和SSe为：,第一节方差分析的基本原理和方法,上述总变异的自由度和平方和可分解为组间和组内两个部分。组间变异即k个平均数的变异，故其自由度为k1，平方和SSt为：,.,10,因此，上述资料的自由度和平方和的分解式为：总自由度组间自由度组内自由度（nk-1）（k1）+k(n-1)总平方和组间平方和组内平方和,第一节方差分析的基本原理和方法,.,11,均方的计算：,第一节方差分析的基本原理和方法,.,12,方差分析表,第一节方差分析的基本原理和方法,.,13,例1：测定东小麦品种东方红3号的蛋白质含量（）10次，得其平均数为14.3，方差为1.621；测定农大139号的蛋白质含量5次，得其平均数为11.7，方差为0.135。试测验东方红3号小麦蛋白质含量的变异是否比农大139为大。假设：H0：1222；HA：1222。显著水平：0.05,DF1=9,DF2=4时,F0.05,(9,4)6.00。,推断：此FF0.05，所以，P0.05接受HA，即东方红3号小麦蛋白质含量的变异大于农大139。,第一节方差分析的基本原理和方法,.,14,分析：两样本分别来自于两个不同的总体，总体方差均为未知，不能假设1222。可采用近似t分布两尾测验的方法。假设：H0：12；HA：12。显著水平：0.05。,回顾t测验法：,东方红3：均数：14.3，方差：1.621，n1=10农大139：均数：11.7，方差：0.135，n2=5,计算；两个样本的样本容量不同，需转换自由度。,.,15,推断：接受HA，否定H0，即两品种蛋白质含量有极显著差异。在12时的t测验，如果两个样本的样本容量相同n1=n2=n，则在t测验时，可不必进行自由度的转换，可直接取自由度为n1。,查表，t0.05,112.301。计算值|t|=5.98t0.05,11，故P22。显著水平：0.05,DF1=3,DF2=12时,F0.05,(3,12)3.49。,第一节方差分析的基本原理和方法,.,17,自由度分解：总变异自由度44115药剂间自由度413药剂内自由度4（41）12平方和分解：SST=222SSt=104SSe=SST-SSt=222-104=118均方：ST2=222/15=14.80St2=104/3=34.67Se2=118/12=9.83,其中，Se2为4种药剂内变异的合并均方，是试验误差的估计值；药剂均方St2则为试验误差加上不同药剂对苗高的效应。,第一节方差分析的基本原理和方法,.,18,推断：接受HA，即测验药剂间变异显著地大于药剂内变异，不同药剂对水稻苗高具有不同效应。,查表5（F值表）：自由度（3；12）,第一节方差分析的基本原理和方法,F.05=3.49；F.01=5.95,.,19,方差分析表,平方和自由度均方FF0.05SSt=1043St2=104/3=34.67St2/Se2=3.53*3.49SSe=SST-SSt=11812Se2=118/12=9.83SST=22215ST2=222/15=14.80,第一节方差分析的基本原理和方法,.,20,2.F分布FDistribution,第一节方差分析的基本原理和方法,.,21,此F值具有S12的自由度1和S22的自由度2。如果我们在给定的1和2下进行一系列抽样，就可得到一系列的F值，这一系列的F值呈F分布。理论统计研究证明，F分布具有平均数F1和取值区间为【0，】的一组曲线，而某一特定的曲线的形状则仅决定于参数1和2。11或12时，F分布曲线呈反向“J”型；当13时，曲线呈偏态。,定义：在一个平均数为，方差为的正态总体中，随机抽取两个独立样本，并求得其均方S12和S22，我们将这两个均方的比值定义为F。,FDistribution,.,22,因自由度不同的F分布曲线,FDistribution,当11或12时，F分布曲线呈反向“J”型；当13时，曲线呈偏态。,f(F),F,.,23,F分布下一定区间的概率可从已制成的统计表查出。附表5系各种v1和v2下右尾概率=0.05和=0.01时的临界F值（一尾概率表）。如查附表5，v1=3，v2=12时，F0.05=3.49，F0.01=5.95，即表示如以v1=3（n1=4）、v2=12（n2=13）在一正态总体中进行连续抽样，则所得F值大于3.49的仅有5%，而大于5.95的仅有1%。所以附表5的数值实际是专供测验S12的总体方差12是否显著大于S22的总体方差22而用的。（H0：1222；HA：1222）。在作F则验时，应以取大值的均方（S12）作分子、取小值的均方（S22）作分母计算F值。若所得FF0.05或F0.01。则该F值即为在=0.05或=0.01水平上显著，应否定H0，接受HA；若所得FF0.05，则接受H0。,FDistribution,.,24,在方差分析的体系中，F测验某项变异因素的效应或方差是否真实存在。所以在计算F值时，总是将要测验的那一项变异因素的均方作分子，而以另一项变异因素（如试验误差项）的均方作分母。这个问题与方差分析的模型和各项变异来源的期望均方有关。在此测验中，如果作分子的均方小于作分母的均方，则F1;此时不必查F表即可确定P0.05,应接受H0。,F测验需具备：（1）变数x遵循正态分布N(,2)（2）S12和S22彼此独立两个条件。当资料不符合这些条件时，需作适合转换。,FDistribution,.,25,3.多重比较(multiplecomparisons),在上例中，接受了HA，仅是指出了东方红3号小麦蛋白质含量的变异大于农大139的。但是，是否各个平均数彼此间都有显著差异呢？还是仅有一部分平均数间有显著差异而另一部分平均数间没有显著差异？仅根据上述分析结果是无法确定的。要明确各个平均数彼此间的差异显著性，还必须对各平均数进行多重比较。,3.1、最小显著差数测验法3.2、最小显著极差测验法(1)Duncans新复极差测验法（Duncan,1955）(2)q测验3.3、比较方法的选择,.,26,3.1最小显著差数测验法leastsignificantdifference，简称LSD法。用此法测验多个平均数时，首先算得平均数差数的标准误：,式中，Se2为方差分析时的误差均方值，n为样本容量。然后查t表得Se2所具有自由度下两尾概率值为的临界t值t，计算得最小显著差数：,若两个平均数的差数LSD，即为在水平上显著。,multiplecomparisons,.,27,例3：以A、B、C、D4种药剂处理水稻种子，各药剂处理后的苗高平均数依次为19、23、18、24cm，作多重比较。,假设：H0:B=A,C=A,D=A;HA:BA,CA，DA。显著水平：0.05,multiplecomparisons,已经算得Se2=9.83，A为对照。,SSDFMSFF0.05104334.673.53*3.49118129.832221514.8,.,28,DF12时，显著水平：0.05查t分布表，t0.052.179，,可以看出，只有XD与对照在0.05上有显著差异，其余两个药剂和对照无显著差异。,multiplecomparisons,SSDFMSFF0.05104334.673.53*3.49118129.832221514.8,注：用LSD法测验多个样本的所有平均数间的差异显著性是不合理的，因为LSD实质是t测验。,.,29,3.2最小显著极差测验法leastsignificantranges，简称LSR法。此法的特点是不同平均数间的比较采用不同的显著差数标准，克服了LSD法的局限性，可用于多样本平均数间的差异显著性比较。这里主要介绍两种类型：（1）Duncans新复极差测验法（Duncan,1955）又称最短显著极差（Shortestsignificantranges,SSR）。,式中，SE为平均数的标准误；Se2为误差均方，n为样本容量。,multiplecomparisons,.,30,查SSR表，查得Se2所具有自由度下，P2，3，k时的SSR值，其中P为两极差间所包含的平均数个数。根据上述公式利用SSR值计算最小显著极差LSR值。具体做法是：将各平均数按大小顺序排列，用各个P的LSR值测验平均数极差的显著性，凡两极差LSR者为接受H0；凡两极差LSR者为接受HA。,multiplecomparisons,.,31,例3：以A、B、C、D4种药剂处理水稻种子，各药剂处理后的苗高平均数依次为19、23、18、24cm，作多重比较。由于已经算得Se2=9.83，且A为对照。假设：H0:B=A,C=A,D=A;HA:BA,CA，DA。显著水平：0.05,查SSR表，DF12，P2时，SSR0.053.08，LSR0.051.573.084.84同理可得，DF12，P3时，SSR0.053.23，LSR0.051.573.235.07同理可得，DF12，P4时，SSR0.053.33，LSR0.051.573.335.23,multiplecomparisons,SSDFMSFF0.05104334.673.53*3.49118129.832221514.8,.,32,平均数从大到小排序：D24B23A19C18,D与B比：242315.23;显著B与A比：231944.84;不显著B与C比：231855.07;不显著A与C比：191814.84;不显著,比较依据:（显著水平：0.05）P2时，LSR0.051.573.084.84P3时，LSR0.051.573.235.07P4时，LSR0.051.573.335.23,multiplecomparisons,.,33,平均数从大到小排序：D24B23A19C18,D与B比：242317.35;不显著B与A比：231946.78;不显著B与C比：231857.14;不显著A与C比：191816.78;不显著,比较依据:（显著水平：0.01）P2时，LSR0.051.574.326.78P3时，LSR0.051.574.557.14P4时，LSR0.051.574.687.35,结论：4个药剂处理水稻的苗高的显著差异来源于处理D和C在0.05水平的显著差异，其余皆差异不显著。,multiplecomparisons,.,34,(2)q测验q测验与新复极差测验相似，其区别仅在于计算最小显著极差LSR值时不是查SSR表，而是查q表，采用下式计算：,multiplecomparisons,.,35,例3：以A、B、C、D4种药剂处理水稻种子，各药剂处理后的苗高平均数依次为19、23、18、24cm，作多重比较。由于已经算得Se2=9.83，且A为对照。假设：H0:B=A,C=A,D=A;HA:BA,CA，DA。显著水平：0.05,multiplecomparisons,.,36,DF=12,P=2、2、4，q0.05。查q表DF=12,P=2时,q0.05=3.08，LSR0.05=1.573.08=4.84同理可得，DF=12,P=3时,q0.05=3.77，LSR0.05=1.573.77=5.92同理可得，DF=12,P=4时,q0.05=4.20，LSR0.05=1.574.20=6.59按同样的方法，可计算出LSR0.01值。详细结果见下表：,表LSR值,multiplecomparisons,.,37,平均数从大到小排序：D24B23A19C18,D与B比：242314.84，不显著D与A比：241955.92，不显著D与C比：241866.59，不显著B与A比：231944.84，不显著B与C比：231855.92，不显著A与C比：191814.84，不显著,结果表明，各平均数之间差异均不显著。,multiplecomparisons,.,38,多重比较结果的表示方法(1)、标记字母法：,表新复极差测验结果的字母标记,D与B比：242315.23，显著B与A比：231944.84，不显著B与C比：231851，则接受HA:0。比较处理效应的试验都应该用固定模型。,4、方差分析的基本假定,.,54,随机模型(randommodel),例：研究水稻杂交F5代系间单株干草重量的遗传变异，随机抽取76个系进行测验，每系取2个样品测定干草重(g/株)。测定结果的方差分析表如下：,变异SSDFMS期望均方系间5459.257572.79系内1350.527617.772,随机模型中i是从总体中随机抽出，服从N(0,2)。这里为随机效应的方差。,4、方差分析的基本假定,.,55,随机模型的F测验：,查表：当n1=75,n2=76时，F.05=1.48；F.01=1.74该例F=4.09，说明系间差异大于系内变异。,若处理效应=0（H0：1=2=k），则F的期望值1。该例中F1，则接受HA:0。,4、方差分析的基本假定,.,56,变异SSDFMS期望均方系间5459.257572.79系内1350.527617.772,该例F1，说明存在，即系间差异存在。,进一步分析系间差异。,这里27.51表示系间差异，即系间遗传变异。2代表环境条件所导致的变异，记作。,代表系间表型变异。,数量遗传学中的遗传率(h2)为：,即F5代家系的表型变异中有60是归属于遗传变异的原因。,4、方差分析的基本假定,.,57,固定模型与随机模型的区别,.,58,方差分析的基本假定,方差分析是建立在一定的线性模型的基础上的。它具有三类原因或效应：（1）处理原因或效应，（2）环境原因或效应，（3）试验误差（这是处理内和环境内的其他非可控因素的变异），故其线性模型为x=+i+j+ij建立这一模型，有如下3个基本假定:,4、方差分析的基本假定,.,59,方差分析的基本假定处理效应与环境效应应该是“可加性”的。对于非可加性资料，一般需作对数转换或其他转换，使其效应变为可加性，才能符合方差分析的线性模型。试验误差应该是随机的、彼此独立的，而且作正态分布，具有平均数为零。N(0,2)所有试验处理必须具有共同的误差方差，即误差同质性假定。,4、方差分析的基本假定,.,60,(1)处理效应与环境效应等应该是“可加性”,依据（xij-)=（i+j+ij）上式两边各取平方求其总和，则得平方和为：(x-)2=bi2aj2ij2因为三类原因均各自独立，所以右边有三个乘积和，即、和，皆为零值。因而得到总平方和等于处理效应平方和加环境效应平方和再加上试验误差平方和。,4、方差分析的基本假定,.,61,可加性特性是方差分析的主要特性，是根据线性模型而产生的必然结果。当从样本估计时，则为(x-x0)2=b(xi.-x0)2+a(x.j-x0)2+(x-xi.-x.j+x0)2或SST=SSA+SSB+SSe由于方差分析具有效应必须可加的假定，故必然导致试验中的=0，=0和=0。一般言之，即各种效应总和与试验误差总和皆等于零。,4、方差分析的基本假定,.,62,（2）试验误差应该是随机的、彼此独立的，而且作正态分布，具有平均数为零。以为多样本的F测验是假定k个样本从k个正态总体中抽取的，所以一定是随机性的。在田间试验中，处理安排在每一区组中均用独立的随机步骤决定而不用顺序排列；这些措施都是为了保证各个误差的彼此独立性和随机性。顺序排列设计的主要缺点是不能获得无偏的试验误差估计，以致方差分