616-第六章 方差分析.ppt
《616-第六章 方差分析.ppt》由会员分享,可在线阅读,更多相关《616-第六章 方差分析.ppt(53页珍藏版)》请在三一文库上搜索。
1、第六章 方差分析,第一节 方差分析的基本原理 第二节 多重比较 第三节 方差分析的基本假定和数据转换,第一节 方差分析的基本原理,所谓方差分析(analysis of variance) ,是关于k(k3)个样本平均数的假设测验方法,是将总变异剖分为各个变异来源的相应部分,从而发现各变异原因在总变异中相对重要程度的一种统计分析方法。,假设测验的依据是:扣除了各种试验原因所引起的变异后的剩余变异提供了试验误差的无偏估计 。,方差是平方和除以自由度的商。,一、自由度和平方和的分解,设有k组数据,每组皆具n个观察值,则该资料共有nk个观察值,其数据分组如表6.1。,表6.1 每组具n个观察值的k 组
2、数据的符号表,在表6.1中,总变异是nk个观察值的变异,故其自由度 v = nk1,而其平方和SST则为:,(61),其中的C称为矫正数:,(62),对于第 i 组的变异,有,从而总变异(61)可以剖分为:,(63),即 总平方和=组内(误差)平方和+处理平方和,组间变异由k个 的变异引起,故其自由度 v =k1 , 组间平方和 SSt 为:,组内变异为各组内观察值与组平均数的变异,故每组具有自由度 v =n1和平方和 ;而资料共有k 组,故组内自由度 v = k (n1) ,组内平方和 SSe 为:,(65),(64),因此,得到表6.1类型资料的自由度分解式为:,(66),总自由度DFT
3、=组间自由度DFt +组内自由度DFe,求得各变异来源的自由度和平方和后,进而可得:,(67),例6.1 以A、B、C、D 4种药剂处理水稻种子,其中A为对照,每处理各得4个苗高观察值(cm),其结果如表6.2,试分解其自由度和平方和。,表6.2 水稻不同药剂处理的苗高(cm),根据(66)进行总自由度的剖分: 总变异自由度 DFT=(nk1)=(44)1=15 药剂间自由度 DFt=(k1)=41=3 药剂内自由度 DFe=k(n1)=4(41)=12,根据(63)进行总平方和的剖分:,或,或 药剂A内: 药剂B内: 药剂C内: 药剂D内:,所以,进而可得均方:,二、F分布与F测验,在一个平
4、均数为 、方差为 的正态总体中,随机抽取两个独立样本,分别求得其均方 s12 和 s22,将 s12 和 s22 的比值定义为F:,(68),此F值具有s12 的自由度 v1 和 s22 的自由度 v2。,所谓F分布,就是在给定的 v1 和 v2 下按上述方法从正态总体中进行一系列抽样,就可得到一系列的F 值而作成一个分布。,F分布下一定区间的概率可从已制成的统计表查出。,F分布曲线特征: (1)具有平均数 =1 (2)取值区间为0,; (3)某一特定曲线的形 状则仅决定于参数 v1和 v2 。 在 v1=1或 v1=2时,F分布 曲线是严重倾斜成反向J型;,当 v13时,曲线转为偏态(图6.
5、1)。,图6.1 F分布曲线 (随v1和v2的不同而不同),F测验需具备条件: (1)变数y遵循正态分布N( , ), (2) s12 和 s22 彼此独立 。,在F 测验中,如果作分子的均方小于作分母的均方,则F0.05,应接受H0。,例6.2 测定东方红3号小麦的蛋白质含量10次,得均方 s12 =1.621;测定农大139小麦的蛋白质含量5次,得均方 s22 =0.135。试测验东方红3号小麦蛋白质含量的变异是否比农大139为大。,假设H0:东方红小麦总体蛋白质含量的变异和农大139一样,即 ,对 。,显著水平 =0.05,v1=9,v2 =4时,F0.05 =6.00。,测验计算: F
6、 =1.621/0.135=12.01,此FF0.05,即P0.05。,推断:否定H0,接受HA,即东方红3号小麦蛋白质含量的变异大于农大139。,例6.3 在例6.1算得药剂间均方st2=168.00,药剂内均方se2=8.17,具自由度 v1=3,v2=12。试测验药剂间变异是否显著大于药剂内变异?,假设 对,显著水平 =0.05, F0.05=3.49。,测验计算:F =168.00/8.17=20.56,查附表5 v1 =3,v2=12时 F0.05 =3.49,F0.01=5.95,实得FF0.01F0.05 。,推断:否定 ,接受 ;即药剂间变异显著地大于药剂内变异,不同药剂对水稻
7、苗高是具有不同效应的。,例6.1和例6.3的分析结果可以归纳在一起,列出方差分析表,如表6.3所示。,表6.3 水稻药剂处理苗高方差分析表,第二节 多重比较,所谓多重比较(multiple comparisons)是指一个试验中k个处理平均数间可能有k(k1)/2个比较,亦称为复式比较。,多重比较有多种方法,本节将介绍常用的三种: 最小显著差数法 复极差法( q法) Duncan氏新复极差法,一、最小显著差数法,最小显著差数法(least significant difference,简称LSD法), 法实质上是第五章的t 测验。 其程序是: (1)在处理间的F测验为显著的前提下,计算出显著水
8、平为 的最小显著差数 ; (2)任何两个平均数的差数( ),如其绝对值 ,即为在水平上差异显著;反之,则为在水平上差异不显著。,已知:,若|t| , 即为在 水平上显著。,因此,最小显著差数为:,(69),当两样本的容量n相等时,,在方差分析中,上式的se2有了更精确的数值 MSe(因为此自由度增大),因此(69)中 的为:,(610),例6.4 试以LSD法测验表6.2资料各种药剂处理的苗高平均数间的差异显著性。,由(例6.3)计算得F=20.56为显著,MSe=8.17,DFe=12,,故,由附表4,v =12时,t0.05 =2.179,t0.01=3.055,故 LSD0.05 =2.
9、1792.02=4.40(cm) LSD0.01=3.0552.02=6.17(cm),然后将各种药剂处理的苗高与对照苗高相比,差数大于4.40cm为差异显著;大于6.17cm为差异极显著。,二、q法,q测验是Student-Newman-Keul基于极差的抽样分布理论提出来的,或称复极差测验,有时又称SNK测验或NK测验。,q法是将一组k个平均数由大到小排列后,根据所比较的两个处理平均数的差数是几个平均数间的极差分别确定最小显著极差值 的。 q测验因是根据极差抽样分布原理的,其各个比较都可保证同一个 显著水平。,q测验尺度值构成为:,(611),(612),式中2pk,p是所有比较的平均数按
10、大到小顺序排列所计算出的两极差范围内所包含的平均数个数(称为秩次距)。 SE为平均数的标准误,可见在每一显著水平下该法有 k1个尺度值。,平均数比较时,尺度值随秩次距的不同而异。,例6.5 试对表6.2资料的各平均数作q测验。,由6.1资料得:,查附表7 q值表,当DF=12时,p=2,3,4的 值,并由(611)计算出尺度值 ,列于表6.4。,表6.4 表6.2资料 值的计算(q测验),由表6.2可知, =29cm, =23cm, =18cm, =14cm。:,由此可得到,-,-,-,-,三、新复极差法,新复极差法是D.B. Duncan(1955)基于不同秩次距p下的最小显著极差变幅比较大
11、而提出的,又称最短显著极差法( shortest significant ranges,SSR )。,查得 后,有,(613),此时,在不同秩次距p下,平均数间比较的显著水平按两两比较是 ,但按p个秩次距则为保护水平,例6.6 试对表6.2资料的各平均数作新复极差测验。,已知 =29cm, =23cm, =18cm, =14cm, MSe=8.17,,查附表8,得值,由(613)算得在p=2,3,4时的值(表6.5),即为测验不同p时的平均数间极差显著性的尺度值。,表6.5 表6.2资料LSR值的计算(新复极差测验),当p=2时, =6(cm) 5水平显著; =5(cm) 5水平显著; =4(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 616-第六章 方差分析 616 第六
链接地址:https://www.31doc.com/p-3022276.html