《方差分析教学》PPT课件.ppt
《《方差分析教学》PPT课件.ppt》由会员分享,可在线阅读,更多相关《《方差分析教学》PPT课件.ppt(87页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第五章 方方 差差 分分 析析1西南科技大学生命科学与工程学院周海廷制作*5.1 方差分析的基本原理方差分析的基本原理5.1.1 5.1.1 自由度和平方和的分解自由度和平方和的分解自由度和平方和的分解自由度和平方和的分解5.1.2 5.1.2 F F分布与分布与分布与分布与F F测验测验测验测验2西南科技大学生命科学与工程学院周海廷制作*上章介绍了一个或两个样本平均数的假设测验方法。本章将介绍k(k3)个样本平均数的假设测验方法,即方差分析方差分析(analysis of variance)。这种方法的基本特点是:将所有将所有k k个个样本的观察值和平均数作为一个整体加以考虑,样本的观察值和
2、平均数作为一个整体加以考虑,把观察值总变异的自由度和平方和分解为不同变把观察值总变异的自由度和平方和分解为不同变异来源的自由度和平方和,进而获得不同变异来异来源的自由度和平方和,进而获得不同变异来源的总体方差估计值。源的总体方差估计值。3西南科技大学生命科学与工程学院周海廷制作*其中,扣除了各种试验原因所引起的变异扣除了各种试验原因所引起的变异后的剩余变异提供了试验误差的无偏估后的剩余变异提供了试验误差的无偏估计,作为假设测验的依据计,作为假设测验的依据。4西南科技大学生命科学与工程学院周海廷制作*5.1.1 自由度和平方和的分解自由度和平方和的分解 方差方差是平方和除以自由度的商。要将一个试
3、验资料的总变异分解为各个变异来源的相应变异,首先必须将总自由度和总平方和分解为各个必须将总自由度和总平方和分解为各个变异来源的相应部分变异来源的相应部分。因此,自由度和平方和的分解是方差分析的第一步。下面我们首先用一个例子来说明这一问题。5西南科技大学生命科学与工程学院周海廷制作*例5.1以A、B、C、D4种药剂处理水稻种子,其中A为对照,每处理各得4个苗高观察值(cm),试分解其自由度和平方和。药剂苗高观察值总和Ti平均数 A182120137218B202426229223C101517145614D2827293211629T=336=216西南科技大学生命科学与工程学院周海廷制作*1
4、1、总变异、总变异 把表中的全部观察值作为一个组看待即把4个处理(4组、每组有4个观察值)合并成一组,共有24个观察值,根据前面讲过的计算平方和的公式,可以计算出总变异的平方和和自由度7西南科技大学生命科学与工程学院周海廷制作*自由度DFT=nk-1=44-1=15。表中的每一个观察值,即包括有处理的效应(不同药剂对苗高的影响)又受到误差的影响。其中:称为矫正数矫正数,用C表示。8西南科技大学生命科学与工程学院周海廷制作*2 2、误差效应、误差效应表中处理内(组内)各观察值之间,若不存在误差,则各观察值应该相等,由于误差是客观存在的,因而处理内(组内)各观察值之间必然是有差异的,因此,可以用组
5、内(处理内)的差异度量误差效应:9西南科技大学生命科学与工程学院周海廷制作*药剂A内:药剂B内:药剂C内:药剂D内:10西南科技大学生命科学与工程学院周海廷制作*从理论上讲,这4个误差平方和除以相应的自由度得的误差均方都可以作为总体误差方差的无偏估计值。但是,用它们的加权平均值来估计总体误差方差,则效果更佳。所以:每个组内(处理内)的自由度为:n-1=4-1=3,11西南科技大学生命科学与工程学院周海廷制作*所以误差的自由度为:DFe=k(n-1)=4(4-1)=123 3、处理效应、处理效应如果没有处理效应,表中各个处理(组)平均数来度量处理效应。从理论上讲均应该相等,因此可以用12西南科技
6、大学生命科学与工程学院周海廷制作*需要注意的是,系样本平均数的方差,为了进行正确的F 测验,必须使它们都是估的估值,而则是是的估值。因而,处理(组间)平方和计同一参数应为:13西南科技大学生命科学与工程学院周海廷制作*14西南科技大学生命科学与工程学院周海廷制作*本例中 平方和:602=504+98自由度:15=3+12因此误差平方和可以采用简单的办法计算SSe=SST-SSt=602-504=98。进而可得均方:15西南科技大学生命科学与工程学院周海廷制作*平方和与自由度的分解归纳为下表变异来源DFSSMS处理间(组间)k-1MSt误差(组内)k(n-1)MSe总变异kn-1将上述例子推广到
7、一般,设有k组数据,每组皆具n个观察值,则资料共有nk个观察值,其数据分组如表6.1(P99)。16西南科技大学生命科学与工程学院周海廷制作*5.1.2 F分布与分布与F测验测验一、一、F F 分布分布 在一个平均数为、方差为2的正态总体中随机抽取两个独立样本,分别求得和,将和的比值定义为F:其均方17西南科技大学生命科学与工程学院周海廷制作*按上述方法从正态总体中进行一系列抽样,就可得到一系列的F值而作成一个F分布。它是具平均数F=1和取值区间为0,的一组曲线;而某一特定曲线的形状仅决定于参数1和2。F分布下一定区间的概率可从已制成的统计表中查出。附表5给出了各种1和2下右尾概率=0.05和
8、=0.0118西南科技大学生命科学与工程学院周海廷制作*二、二、F F 测验测验 在方差分析的体系中,F测验可用于检测某项变异因素的效应或方差是否存在。所以在计算F值时,总是将要测验的那一项变异因素的均方作分子,而以另一项变异(如误差项)作分母。时的临界F 值。其值是专供测验的总体方是否显著大于的总体方差差而设计的(H0:对HA:)。19西南科技大学生命科学与工程学院周海廷制作*F测验需具备的条件:(1)变数y遵循N(,2);例6.3 在例6.1中算得药剂间均方=168.00,药剂内均方=8.17,具有自由度1=3,2=12。试测验药剂间变异是否显著大于药剂内变异?假设H0:对HA:=0.05
9、和彼此独立。(2)20西南科技大学生命科学与工程学院周海廷制作*查附表5在1=3,2=12时 F0.05=3.49,F0.01=5.95实得F F0.01 P0.01测验计算:21西南科技大学生命科学与工程学院周海廷制作*将例6.1和例6.3的分析结果归纳在一起,列出方差分析表如下:变异来源DFSSMSF显著F值药剂处理间3504168.0020.56F0.05=3.49药剂处理内(误差)12988.17F0.01=5.95总变异15602水稻药剂处理苗高方差分析表推断:否定H0:,接受HA:22西南科技大学生命科学与工程学院周海廷制作*5.2 多重比较多重比较5.2.1 5.2.1 最小显著
10、差数法最小显著差数法最小显著差数法最小显著差数法5.2.2 5.2.2 q q 法法法法5.2.3 5.2.3 新复极差法新复极差法新复极差法新复极差法5.2.4 5.2.4 多重比较方法的选择多重比较方法的选择多重比较方法的选择多重比较方法的选择23西南科技大学生命科学与工程学院周海廷制作*5.2.1 最小显著差数法最小显著差数法最小显著差数法最小显著差数法(least significant differrence,简称LSD法)24西南科技大学生命科学与工程学院周海廷制作*例6.4 试以LSD法测验各种药剂处理的苗高平均数之间的差异显著性。由附表4,=12时,t0.05=2.179,t
11、0.01=3.055故 LSD0.05=2.1792.02=4.40(cm)LSD0.01=3.0552.02=6.17(cm)25西南科技大学生命科学与工程学院周海廷制作*处理苗高平均数差异显著性0.050.01D29B23A18C14不同药剂处理水稻苗高平均数比较(LSD法)abccAABBCC26西南科技大学生命科学与工程学院周海廷制作*5.2.2 q 法法 q q 测验方法测验方法是将k个平均数由大到小排列后,根据所比较的两个处理平均数的差数是几个平均数间的极差分别确定最小显著极差LSR值的。27西南科技大学生命科学与工程学院周海廷制作*例6.5 试以q法测验各种药剂处理的苗高平均数之
12、间的差异显著性。查附表7,得到当DF=12时,p=2,3,4的q值28西南科技大学生命科学与工程学院周海廷制作*LSR值Pq 0.05q 0.01LSR0.05LSR0.0123.084.324.406.1833.775.045.397.2144.205.506.017.87处理苗高平均数差异显著性0.050.01D29aAB23 b ABA18 c BCC14 c C不同药剂处理水稻苗高平均数比较(q法)29西南科技大学生命科学与工程学院周海廷制作*5.2.3 新复极差法新复极差法 新复极差法新复极差法,又称最短显著极差法(shortest significant range),与q法相似。
13、计算LSR值查的是SSR值(附表8)而不是q表。LSRLSR 值值PSSR 0.05SSR 0.01LSR0.05LSR0.0123.084.324.406.1833.234.554.626.5143.334.684.766.6930西南科技大学生命科学与工程学院周海廷制作*5.2.4 多重比较方法的选择多重比较方法的选择1、试验事先确定比较的标准,凡是与对照相比较,或与预定要比较的对象比较,一般可选用最小显著差数法;2、根据否定一个正确的H0和接受一个不正确的H0的相对重要性来决定。31西南科技大学生命科学与工程学院周海廷制作*方差分析的基本步骤:方差分析的基本步骤:(1)分解平方和与自由度
14、;(2)F测验;(3)平均数的多重比较。32西南科技大学生命科学与工程学院周海廷制作*5.3 方差分析的线性模型与期望均方方差分析的线性模型与期望均方5.3.1 5.3.1 方差分析的线性数学模型方差分析的线性数学模型方差分析的线性数学模型方差分析的线性数学模型5.3.2 5.3.2 期望均方期望均方期望均方期望均方33西南科技大学生命科学与工程学院周海廷制作*5.3.1 方差分析的线性数学模型方差分析的线性数学模型 方差分析是建立在一定的线性可加模型的基础上的。所谓线性可加模型是指总体每一个变量可按其变异的原因分解成若干个线性组成部分,它是方差分析的基础。34西南科技大学生命科学与工程学院周
15、海廷制作*表6.1数据的线性模型可表示为:式中,为总体平均数,i为试验处理效应,ij为随机误差具有N(0,2)。35西南科技大学生命科学与工程学院周海廷制作*在以样本符号表示时,样本的线性组成为:是的无偏估计值,36西南科技大学生命科学与工程学院周海廷制作*37西南科技大学生命科学与工程学院周海廷制作*5.3.2 期望均方期望均方 在线性可加模型中,由于对i有不同解释产生了固定模型(I)和随机模型(II)。一、固定模型一、固定模型(fixed model(fixed model)指试验的各处理都抽自特定的处理总体,其处理效应i=(i-)是一个固定的常量,我们的目的就在于研究i,所测验的假设是H
16、0:i=0或H0:i=。38西南科技大学生命科学与工程学院周海廷制作*一般的栽培和饲养试验,如肥料试验、药效试验、密度试验、饲料试验、品种试验等均属于固定模型。例6.8以5个水稻品种作大区比较试验,每品种作3次取样,测定其产量,所得数据为单向分组资料。本试验需明确各品种的效应,故为固定模型,方差分析和期望均方的参数列入下表:39西南科技大学生命科学与工程学院周海廷制作*5个水稻品种产量的方差分析与期望均方表变异来源DFSSMS期望均方(EMS):固定模型品种间品种内41087.624.021.92.40固定模型的处理效应(本例为品种效应)i属于固定效应,固定效应的方差用表示。固定模型的F测验4
17、0西南科技大学生命科学与工程学院周海廷制作*二、随机模型二、随机模型(random model)(random model)指试验中的各处理皆是抽自N(0,)的一组随机样本,因而处理效应i是随机的,它会因试验的不同而不同;故我们的目的不在于研究i而在于研究i的变异度。随机模型在遗传、育种和生态的研究试验方面有较广泛的用处。41西南科技大学生命科学与工程学院周海廷制作*例6.9研究籼粳杂交F5代系间单株干草重的遗传变异,随机抽取76个系进行试验,每系随机取2个样品测定干草重(g/株)。因这76个系是随机抽取的样本,要从这些样本来估计F5代系间单株干草重的遗传变异,故这是随机模型。其方差分析的结果
18、如下:变异来源DFMS期望均方(EMS):固定模型系统间系统内757672.7917.7742西南科技大学生命科学与工程学院周海廷制作*随机模型的F 测验本例中系统内 MS 估计了2,因而;系统间MS估计了因而43西南科技大学生命科学与工程学院周海廷制作*这是测验处理效应的变异度,而不是测验处理效应本身。本例F72.79/17.77=4.09F0.05,说明单株干草重存在遗传变异。44西南科技大学生命科学与工程学院周海廷制作*5.4 单向分组资料的方差分析单向分组资料的方差分析5.4.1 5.4.1 组内观察值数目相等的单向分组内观察值数目相等的单向分组内观察值数目相等的单向分组内观察值数目相
19、等的单向分组资料的方差分析组资料的方差分析组资料的方差分析组资料的方差分析5.4.25.4.2组内又分亚组的单向分组资组内又分亚组的单向分组资组内又分亚组的单向分组资组内又分亚组的单向分组资料的方差分析料的方差分析料的方差分析料的方差分析45西南科技大学生命科学与工程学院周海廷制作*5.4.1 组内观察值数目相等的单向分组资料的方组内观察值数目相等的单向分组资料的方差分析差分析每组具n个观察值的k组数据的符号表组别观察值(yij,i=1,2,k;j=1,2,n)总和 平均 均方1y11y12y1jy1nT12y21y22y2jy2nT2:iyi1yi2yijyinTi:kyk1yk2ykjyk
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 方差分析教学 方差分析 教学 PPT 课件
限制150内