方差分析

《方差分析》由会员分享，可在线阅读，更多相关《方差分析（73页珍藏版）》请在金锄头文库上搜索。

1、第六章 方差分析,在实际研究过程中，对三个及三个以上多个样本平均数进行比较，用t或u检验会出现的问题：,第六章 方差分析,在实际研究过程中，对三个及三个以上多个样本平均数进行比较，用t或u检验会出现的问题：,所出现的问题,检验过程繁锁 无统一的试验误差，误差估计的精确性和检验的灵敏性低 推断的可靠性降低，犯I类错误的概率增加,方差分析(analysis of variance, ANOVA)又叫变量分析，是英国著名统计学家R.A. Fisher于1923年提出的。方差分析是将所有处理的观测值作为一个整体，一次比较就对所有各组间样本平均数是否有差异作出判断。如果差异不显著，则认为他们都是相同的，

2、如果差异显著，再进一步比较是哪组数据与其他数据不同。,第六章 方差分析,方差分析 是对两个或多个样本平均数差异学著性检验的方法，是将测量数据的总变异按照变异来源分解为处理效应和试验误差，并做出其数量估计。,分类,单因素方差分析 二因素方差分析 多因素方差分析,基本步骤,分解 F检验 多重比较,第六章 方差分析,方差分析 是对两个或多个样本平均数差异性作显著性检验的方法，是将测量数据的总变异按照变异来源分解为处理效应和试验误差，并做出其数量估计。,多重比较,最小显著差数法 最小显著极差法,第六章 方差分析,方差分析 是对两个或多个样本平均数差异性作显著性检验的方法，是将测量数据的总变异按照变异来

3、源分解为处理效应和试验误差，并做出其数量估计。,多重比较,最小显著差数法 最小显著极差法,基本假定条件,正态性 可加性 方差同质性,数据缺失,可利用误差平方和最小的原则对其进行弥补,第六章 方差分析,1. 方差分析的基本原理,1. 相关术语,1) 试验因素,试验中所研究的影响试验指标的原因或原因组合称为试验因素(experimental factor)或处理因素(treatment factor)，简称因素或因子。按性质不同，可分为可控因素和非控因素。,可控因素(controllable factor),非控因素(uncontrollable factor),试验因素常用大写字母A、B、C等来

4、表示。,第六章 方差分析,1. 方差分析的基本原理,1. 相关术语,2) 因素水平,每个试验因素的不同状态(处理的某种特定状态或数量上的差别)称为因素水平(level of factor)，简称水平。,水平常用大写字母添加下标1、2、3等来表示，如A1、A2、A3、 ，B1、 B2、 B3 ， C1、 C2、C3 等来表示。,1) 试验因素,第六章 方差分析,1. 方差分析的基本原理,1. 相关术语,2) 因素水平,试验处理(experimental treatment)通常也称为处理(treatment)，指对受试对象给予的某种外部干预，是试验实施因子水平的一个组合。,1) 试验因素,3)

5、试验处理,单因素处理(single factor treatment),多因素处理(multiple factor treatment),4) 试验单位,5) 重复,二因素处理(single factor treatment),2. 方差分析的基本原理,观测值的差异,处理效应,误差效应,方差分析的基本思想就是将测量数据的总变异按照变异原因不同进行分解，分为处理效应和试验误差，并对其进行数量估计。,判定标准？,方差，即均方(mean squares),3. 数学模型,固定模型(fixed model),随机模型(random model),混合模型(mixed model),不同的模型在平方和及

6、自由度的计算上是相同的，但在进行假设检验时，F值的计算公式是不同的，具有不同的侧重点,4. 平方和与自由度的分解,总平方和(total sum of squares,总自由度(total degree of freedom),1) 平方和的分解,2) 自由度的分解,3) 算计方差,5. 统计假设的显著性检验-F检验,处理内方差SSe,也就是饲料内方差可以估计误差方差,处理间方差SSt,饲料间方差，可以估计不同饲料喂养增重的差异。,F检验(F-test),被检验因素的均方,误差均方,6. 多重比较,1) 最小显著差数法(LSD法),2)最小显著极差法(LSR法),为了比较不同处理平均数两两间差异

7、的显著性，每个处理的平均数都要与其他处理的平均数进行比较，统计上把多个平均数两两间的相互比较称为多重比较(multiple comparisons),是Fisher最早用于检验所有总体均数间两两相等假设的方法，其实质是两个平均数相比较的t检验法。,是在一定的显著水平上，根据极差范围内所包含的处理数据(也称为秩次距)M的不同而采用不同的显著差数标准进行比较。其又可分为新复极差检验和q检验。,最小显著差数法(LSD),首先计算出达到差异显著的最小差数，记为LSD； 然后与两个处理平均数的差值进行比较； 如果处理的平均数差值大于LSD，则认为在给定的显著水平上有显著性差异，反之，则认为没有显著性差异

8、。,多重比较结果表示方法,标记字母法 梯形法,标记字母法，首先将全部平均数从大到小依次排列，然后在最大的平均数上标记字母a，将该平均数以下各平均数相比，凡相差不显著的(LSD0.05)都标上字母a，直至某个与之相差显著的则标以字母b。再以标有b的平均数为标准与各个平均数比较，凡差数差异不显著的在字母a后再续标字母b，直至差异显著的平均数标以c,然后重复上述工作，直到最小的平均数有标记为止。凡标有一个相同字母的即为差异不显著。,多重比较结果表示方法,标记字母法 梯形法,梯形法，首先将各处理的平均数差数按梯形列于表中，并将这些差数的LSD值比较。若差数LSD0.05说明处理平均数间的差异达到显著水

9、平，在右上角标上“*”；差数 LSD0.01，说明处理平均数间差异达到极显著水平，在差数的右上角标“*”，反之如差数小于LSD说明差异不显著。,利用LSD法进行多重比较的步骤， 1) 计算最小显著差数LSD0.05和LSD0.01； 2) 列出平均数的多重比较表，表中各处理按其平均数从大到小依次进行排列。 3) 将两两平均数的差数与 LSD0.05和LSD0.01 进行比较，作出统计推断。,不足，在差数比较时没有考虑到大小排列上的次序，仍有推断可靠性低,犯I类错误概率增加的问题。,最小显著极差法(LSR法),相对于LSD法，该类方法同时考虑到平均数差数的秩序问题，然后进行不同的显著性标准进行比

10、较，因此比LSD法更为可靠。 LSR检验又可分为新复极差检验和q检验,新复极差检验是由邓肯提出的，因此又称Duncan法或SSR法；,q检验法也称为Student-Newman-Keuls法，SNK,新复极差检验，新复极差检验是由邓肯于1955年提出的，又称Duncan法，或SSR法。其比较步骤如下：,1） 按相比较的样本容量计算平均数标准误，当n1=n2=n时，有：,2） 查s2e所具有自由度dfe和比较所含平均数个数M时的SSR值，然后计算最小显著极差：,3） 将各平均数按大小顺序排列，用各个M值的LSRa值即可检验各平均数间差异的显著性。：,对上面的各组平均值作新复极差检验：,查附表获得

11、df=16,M=2时，SSR的值：,查附表获得df=16,M=2时，SSR的值：SSR0.05=3.0, SSR0.01=4.13，则：,当A1与A4比较时，M=2，差值(32.0)31.02且小于42.70，所以两者差异达到显著水平； 当A1与A2比较时，M=3，差值(49.0)44.57，所以差异达到极显著水平； 同理可将本个平均数一一用此方法进行比较和显著性水平判断。,LSR法,LSD法,新复极差检验(new multiple range test)是由邓肯于1955年提出，因此又称Duncan法或SSR(shortest significant ranges)法。,q检验法也称为Stu

12、dent-Newman-Keuls (SNK)检验，是以统计量q的概率分布为基础的，方法与新复极差检验相似，其区别仅在于计算最小显著极差时不是查SSR，而是相q值。,对三种检验方法的比较中发现，当样本数(处理数)k=2时，LSD法、SSR法和q检验的显著尺度是相同的，当k3时三种显著尺度便不相同：LSD法SSR法q法 因此 对于精度要求高的研究-q检验 一般试验-SSR检验,方差分析的基本步骤,将样本数据的总平方和与总自由度分解为各变异因素的平方和与自由度； 列方差分析表进行F检验，分析各变异因素在总变异中的重要程度； 若F检验显著，对各处理平均数进行多重比较,第六章 方差分析,2. 单因素方

13、差分析,按重复数分类,单因素试验(single factor experiment)资料的方差分析是比较简单的一种，目的在于正确判断试验因素各水平的相对效果。在单因素方差分析中，根据组内观测数目(重复数)是否相同，可分为组内观测次数相等的方差分析和组内观测资料不等的方差分析,组内观测次数相等的方差分析,组内观测次数不相等的方差分析,一、 组内观测次数相等的方差分析,K组资料中，每一处理组皆含有n个观测值，其方差分析方法前面已做过介绍，此处只给出相关的计算公式：,例：测定东北、内蒙古、河北、安徽、贵州5个地区黄鼬冬季针毛的长度，每个地区随机抽取4个样本，测定的结果列于表中，试比较各地区其针毛长度

14、的差异显著性。,组间有极显著性差异！,哪一组？,计算结果：,多重检验1.: (最小显著差数法),当自由度为15时，可以得到t0.05、t0.01水平的值分别为：2.131， 2.947,多重检验2.: (最小显著极差法 LSR法),1) 按相比较的样本容量计算平均数标准误,2) 根据整体方差的自由度和各个M值，查出所对应的SSR值，然后计算出最小显著极差：,3) 将各平均数按大小顺序排列，用各个M值的LSR值来检查各平均数间差异的显著性。,查表，当df=15,M=2时，SSR0.05和0.01水平的值分别为：3.01，4.17。,多重检验3.: (q法),二、 组内观测次数不相等的方差分析,当

15、每组资料中，每一处理组的观测值个数n不同时，上面所述的方差分析方法仍然可用，只是总观测数要重新计算，不是nk。相对应的，计算平方和的公式也稍有改变。,多重比较-平均数的标准误-平均样本量n0,第六章 方差分析,3. 二因素方差分析,按重复数分类,二因素试验(two factor experiment)方差分析中，需要对因素的主效和因素间的互作进行分析。因素间的交互作用显著与否关系到主效的利用价值，有时互作效应相当大，甚至可以忽略主效应。二因素间是否存在互作可根据专门的统计方法或专业知识进行判断。,无重复观测值的二因素方差分析,具有重复 观测值的二因素方差分析,一、无重复观测值的二因素方差分析,

16、即在试验中，依据经验或专业知识判断二因素无互作时，每个处理可不设重复。假定A因素有a个水平，B因素有b个水平，每个处理组合只有一个观测值。其分组资料模式如下：,在无重复观测值的二因素试验资料中，A因素的每一个水平可看作有b个重复，B因素的每一个水平可看作有a个重复。因此每个观测值既受A因素影响，又受到B因素的影响。,若因素间不存在互作时，则二因素方差分析观测值的线性模型为：,其中， 、 分别是A因素和B因素的效应，可以是固定的，也可以是随机的，且 。 是随机误差，彼此独立且服从 。,(1) 平方和的分解,(2) 自由度的分解,(2) 各项方差的计算,例：将一种生长激素配成M1、M2、M3、M4、M5五种浓度，并用H1、H2、H3三种时间浸渍某大豆种子，出苗45天后得各处理每一植株的平均干物重(g)，结果列于下表，试作方差分析并进行多重比较。,(1) 平方和的分解,(2) 自由度的分解,