《方差分析基本原理.ppt》由会员分享,可在线阅读,更多相关《方差分析基本原理.ppt(119页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于方差分析的基本原理第一张,PPT共一百一十九页,创作于2022年6月1、u或或t 检验过程烦琐检验过程烦琐 例如,一试验包含例如,一试验包含5个处理,采用个处理,采用t检验法要检验法要进行进行 =10次两两平均数的差异显著性检验;次两两平均数的差异显著性检验;若有若有k个处理,则要作个处理,则要作 =k(k-1)/2 次类似的次类似的检验。检验。第二张,PPT共一百一十九页,创作于2022年6月2、无统一的试验误差,误差估计的准确性和、无统一的试验误差,误差估计的准确性和检验的灵敏性低检验的灵敏性低(1)(1)t 检验要进行两两比较,每次仅用检验要进行两两比较,每次仅用2 2个样本信息个样
2、本信息估计总体方差,误差估计的估计总体方差,误差估计的准确性低准确性低(2)k个处理平均值的自由度为个处理平均值的自由度为k(n-1),而,而t 检验查检验查t值的自由度为值的自由度为2(n-1),从而,从而降低了检验的灵敏性降低了检验的灵敏性两两比较合并均方两两比较合并均方:k 个样本合并均方个样本合并均方:第三张,PPT共一百一十九页,创作于2022年6月3、t检验增大犯检验增大犯错误的概率错误的概率 t 检验时对具有不同秩次的平均数采用同一个检验时对具有不同秩次的平均数采用同一个t,会增大犯会增大犯 错误的概率,降低推断的可靠性。错误的概率,降低推断的可靠性。2个平均数比较:个平均数比较
3、:=0.055个平均数比较:个平均数比较:=1-(1-0.05)10=0.401310个平均数比较:个平均数比较:”=1-(1-0.05)45=0.9006第四张,PPT共一百一十九页,创作于2022年6月因此,多个平均数的差异显著性检验不宜用因此,多个平均数的差异显著性检验不宜用t(或或u)检验,须采用检验,须采用方差分析法方差分析法。方差分析方差分析(analysis of variance)由英国统计学家由英国统计学家R.A.Fisher于于1923年提出的。年提出的。将多个样本将多个样本(处理处理)的观测值作为一个总体,的观测值作为一个总体,用方差来表示变异,把引起事件总的变异分解为用
4、方差来表示变异,把引起事件总的变异分解为各种因素的变异,并对每个因素引起的变异作数各种因素的变异,并对每个因素引起的变异作数量估计,从而说明各因素的变异幅度及其在总变量估计,从而说明各因素的变异幅度及其在总变异中的重要程度;并用剩余变异无偏估计随机误异中的重要程度;并用剩余变异无偏估计随机误差,进而比较处理均值间的差异。差,进而比较处理均值间的差异。第五张,PPT共一百一十九页,创作于2022年6月有关术语:1、试验指标(experimental index)为衡量试验结果的好坏或处理效应的高低,在试验中具体测定的性状或观测的项目称为试验指标。由于试验目的不同,选择的试验指标也不相同。在畜禽、
5、水产试验中常用的试验指标有:日增重、产仔数、产奶量、产蛋率、瘦肉率、某些生理生化和体型指标(如血糖含量、体高、体重)等。2、试验因素(experimental factor)试验中所研究的影响试验指标的因素叫试验因素。如研究如何提高猪的日增重时,饲料的配方、猪的品种、饲养方式、环境温湿度等都对日增重有影响,均可作为试验因素来考虑。当试验中考察的因素只有一个时,称为单因素试验;若同时研究两个或两个以上的因素对试验指标的影响时,则称为两因素或多因素试验。试验因素常用大写字母A、B、C、等表示。第六张,PPT共一百一十九页,创作于2022年6月3、因素水平(level of factor)试验因素所
6、处的某种特定状态或数量等级称为因素水平,简称水平。研究某种饲料中4种不同能量水平对肥育牛瘦肉率的影响,这4种特定的能量水平就是饲料能量这一试验因素的4个水平。因素水平用代表该因素的字母加添足标1,2,来表示。如A1、A2、,B1、B2、,等。4、试验处理(treatment)事先设计好的实施在试验单位上的具体项目叫试验处理,简称处理。在单因素试验中,实施在试验单位上的具体项目就是试验因素的某一水平。例如进行饲料的比较试验时,实施在试验单位(某种畜禽)上的具体项目就是喂饲某一种饲料。所以进行单因素试验时,试验因素的一个水平就是一个处理。在多因素试验中,实施在试验单位上的具体项目是各因素的某一水平
7、组合。例如进行3种饲料和3个品种对猪日增重影响的两因素试验,整个试验共有33=9个水平组合,在多因素试验时,试验因素的一个水平组合就是一个处理。第七张,PPT共一百一十九页,创作于2022年6月5、试验单位(experimental unit)在试验中能接受不同试验处理的独立的试验载体叫试验单位。在畜禽、水产试验中,一只家禽、一头家畜、一只小白鼠、一尾鱼,即一个动物;或几只家禽、几头家畜、几只小白鼠、几尾鱼,即一组动物都可作为试验单位。试验单位往往也是观测数据的单位。6、重复(repetition)在试验中,将一个处理实施在两个或两个以上的试验单位上,称为处理有重复;一处理实施的试验单位数称为
8、处理的重复数。例如,用某种饲料喂4头猪,就说这个处理(饲料)有4次重复。第八张,PPT共一百一十九页,创作于2022年6月二、方差分析的基本原理二、方差分析的基本原理1、线性模型与基本假定线性模型与基本假定假设某单因素试验有假设某单因素试验有k个处理,每个处理有个处理,每个处理有n次次重复,共有重复,共有nk个观测值。这类试验资料的数据模个观测值。这类试验资料的数据模式如表式如表7.1所示。所示。第九张,PPT共一百一十九页,创作于2022年6月表表7.1 k个处理每处理有个处理每处理有n个观测值的数据模式个观测值的数据模式x1.x2.xi.xk.x注:注:xij指第指第i个处理第个处理第j个
9、观察值(个观察值(i=1k;j=1n)第十张,PPT共一百一十九页,创作于2022年6月xij 可以分解为可以分解为:其中其中:表示全试验观测值总体的平均数;表示全试验观测值总体的平均数;i 是是 第第 i 个个 处理的效应处理的效应(treatment effects)表示处理)表示处理i对试验结对试验结果产生的影响。果产生的影响。显然有显然有ij 是试验误差,相互独立,且服从是试验误差,相互独立,且服从 正态分布正态分布N(0,2)。)。(7-1)式叫做)式叫做 单因素试验的单因素试验的线性模型线性模型,亦称,亦称数学模型数学模型。在这个模型中在这个模型中xij表示为总平均数表示为总平均数
10、、处理效应、处理效应 i、试验误、试验误差差ij之和。之和。(7-1)(7-2)第十一张,PPT共一百一十九页,创作于2022年6月 由由ij 相相 互独立且服从正态分布互独立且服从正态分布 N(0,2),可),可知各处理知各处理i(i=1,2,k)所属总体亦应具正态性,所属总体亦应具正态性,即服从正态分布即服从正态分布N(i,2)。尽管各总体的均数。尽管各总体的均数 i 可可以不等或相等,以不等或相等,2则必须是相等的。所以,单因则必须是相等的。所以,单因素试验的数学模型可归纳为:素试验的数学模型可归纳为:效应的可加性效应的可加性(additivity)、分布的正态性分布的正态性(norma
11、lity)、方差的同质性方差的同质性(homogeneity)。这也。这也是进行其它类型方差分析的前提或基本假定。是进行其它类型方差分析的前提或基本假定。第十二张,PPT共一百一十九页,创作于2022年6月若将表若将表7.1中的观测值中的观测值 xij(i=1,2,k;j=1,2,n)的数)的数据结构(模型)用样本符号来表示,则据结构(模型)用样本符号来表示,则(7-3)第十三张,PPT共一百一十九页,创作于2022年6月 (7-1)、()、(7-3)两式告诉我们:)两式告诉我们:每每 个个 观观 测测 值值 都包含处理效应都包含处理效应(i 或或 ),与误差与误差(ij 或或 ),故,故kn
12、个观测值的总变异个观测值的总变异可分解为可分解为处理处理间间的变异的变异和和处理处理内内的变异的变异两部分。两部分。第十四张,PPT共一百一十九页,创作于2022年6月药剂药剂 A(xA(x1 1)B(xB(x2 2)C(xC(x3 3)D(xD(x4 4)19 19212120202222 23 23242418182525 21 21272719192727 13 13202015152222总和总和T Ti.76 76929272729696平均平均 19 19232318182424观察值观察值【例7.1】-2 2 -3 3 第十五张,PPT共一百一十九页,创作于2022年6月药剂药剂
13、 A(xA(x1 1)B(xB(x2 2)C(xC(x3 3)D(xD(x4 4)1919212120202222232324241818252521212727191927271313202015152222总和总和T Ti.7676929272729696平均平均1919232318182424观察值观察值第十六张,PPT共一百一十九页,创作于2022年6月2、自由度与平方和的剖分自由度与平方和的剖分 在方差分析中是用样本方差即在方差分析中是用样本方差即均方均方(MS)来度量资)来度量资料的变异程度的。料的变异程度的。表表7.1中全部观测值的总变异可以用总均方来度量。中全部观测值的总变异可
14、以用总均方来度量。将总变异分解为处理间变异和处理内变异,就是要将总变异分解为处理间变异和处理内变异,就是要将将 总总 均方均方 分解为处理间均方和处理内均方。但这种分分解为处理间均方和处理内均方。但这种分解是通过将总均方的分子解是通过将总均方的分子称为总离均差平方和,简称为总离均差平方和,简称为总平方和,剖分成处理间平方和与处理内平方和称为总平方和,剖分成处理间平方和与处理内平方和两部分;将总均方的分母两部分;将总均方的分母称为总自由度,剖分成处称为总自由度,剖分成处理间自由度与处理内自由度两部分来实现的。理间自由度与处理内自由度两部分来实现的。第十七张,PPT共一百一十九页,创作于2022年
15、6月(1)平方和的剖分平方和的剖分 在表在表7.1中,反映中,反映 全部观测值总变异的总平方和是各全部观测值总变异的总平方和是各观测值观测值xij与总平均数的离均差平方和,记为与总平均数的离均差平方和,记为SST。即即(7-4)第十八张,PPT共一百一十九页,创作于2022年6月因为因为第十九张,PPT共一百一十九页,创作于2022年6月其中其中所以所以(7-5)(7-5)式中,)式中,为各处理平均数与总平均为各处理平均数与总平均数的离均差平方和与重复数数的离均差平方和与重复数n的乘积的乘积,反映了重复,反映了重复 n 次的处理间变异次的处理间变异,称为处理间平方和,记为,称为处理间平方和,记
16、为SSt,即,即:(7-6)第二十张,PPT共一百一十九页,创作于2022年6月(7-5)式中,式中,为各处为各处 理内离均差平方和之理内离均差平方和之和,反映了各处理内的变异即误差,称为处理内平和,反映了各处理内的变异即误差,称为处理内平方和或误差平方和,记为方和或误差平方和,记为SSe,即,即:(7-7)第二十一张,PPT共一百一十九页,创作于2022年6月于是有于是有 SST=SSt+SSe (7-8)这个关系式中三种平方和的简便计算公式如下:这个关系式中三种平方和的简便计算公式如下:其中,其中,C=T2/(kn)称为矫正数。称为矫正数。(7-9)CTnSSikit-=S=2.11CxS
17、SijnjkiT-=SS=211第二十二张,PPT共一百一十九页,创作于2022年6月(2)自由度的剖分自由度的剖分 在计算总平方和时,资料中的各个观测值要受在计算总平方和时,资料中的各个观测值要受 这一条件的约束,故总自由度等于资料中观测值的总这一条件的约束,故总自由度等于资料中观测值的总个数减个数减1,即,即kn-1。总自由度记为总自由度记为dfT,即,即:dfT=kn 1 (7-10)第二十三张,PPT共一百一十九页,创作于2022年6月在计算处理间平方和时,各处理均数要受在计算处理间平方和时,各处理均数要受这一条件的约束,故处理间自由度为处这一条件的约束,故处理间自由度为处理数减理数减
18、1,即,即k-1。处理间自由度记为处理间自由度记为dft,即:,即:dft=k-1 (7-11)在计算处理内平方和时,要受在计算处理内平方和时,要受k个条件的约束,即个条件的约束,即,(i=1k)。故处理内自由度为资料中观测值的。故处理内自由度为资料中观测值的总个数减总个数减k,即,即kn-k。处理内自由度记为。处理内自由度记为dfe,即:,即:dfe=kn k=k(n-1)(7-12)第二十四张,PPT共一百一十九页,创作于2022年6月因为因为所以所以 综合以上各式得:综合以上各式得:(7-13)(7-14)第二十五张,PPT共一百一十九页,创作于2022年6月即:即:总均方一般不等于处理
19、间均方加处理内均方。总均方一般不等于处理间均方加处理内均方。各部分平方和除以各自的自由度便得到总均各部分平方和除以各自的自由度便得到总均方、处理间均方和处理内均方,方、处理间均方和处理内均方,分别记为分别记为 MST(或(或 )、)、MSt(或(或 )和)和MSe(或(或 )。)。(7-15)(3)均方的计算均方的计算 第二十六张,PPT共一百一十九页,创作于2022年6月计算均方后,通过比较处理间均方相对误差均方的计算均方后,通过比较处理间均方相对误差均方的大小即可判断处理效应所引起的变异所占比重,从大小即可判断处理效应所引起的变异所占比重,从而可以判断试验是否存在明显处理效应。而可以判断试
20、验是否存在明显处理效应。第二十七张,PPT共一百一十九页,创作于2022年6月三、三、方差分析中的方差分析中的F 测验测验方差分析的方差分析的F 检验用于测验某项变异因素的效应检验用于测验某项变异因素的效应是否真实存在是否真实存在将要测验的那一项变异因素的均方作分子,另将要测验的那一项变异因素的均方作分子,另一项变异因素(例如一项变异因素(例如误差项误差项)的均方作分母()的均方作分母(具具体情况与所用试验设计和模型有关体情况与所用试验设计和模型有关)如果如果F0.05,应该,应该接受接受H0第二十八张,PPT共一百一十九页,创作于2022年6月F测验需具备的条件:测验需具备的条件:(1)被抽
21、样总体的变数)被抽样总体的变数x服从正态分布,即服从正态分布,即xN(,2 2)(2)和和 彼此相互独立。彼此相互独立。注:当试验资料不符合这些条件时,需要作适注:当试验资料不符合这些条件时,需要作适当转换。当转换。第二十九张,PPT共一百一十九页,创作于2022年6月第二节 方差分析的一般步骤一、一、平方和平方和与与自由度自由度的的分解分解二、二、列列出出方差分析表方差分析表,作,作F测验测验三、若三、若F检验显著,则进行检验显著,则进行多重比较多重比较四、四、结果结果的的解释解释/说明说明第三十张,PPT共一百一十九页,创作于2022年6月【例例】以A(CK)、B、C、D共4种药剂处理水稻
22、种子,测得苗高结果如下表:第三十一张,PPT共一百一十九页,创作于2022年6月一、一、平方和平方和与与自由度自由度的的分解分解1、平方和的计算方法一:方法一:方法二:方法二:,C=T2/(kn)CTnSSikit-=S=2.11CxSSijnjkiT-=SS=211第三十二张,PPT共一百一十九页,创作于2022年6月方法二:C=T2/(kn)=3362/(44)=7056CxSSijnjkiT-=SS=211CTnSSikit-=S=2.11=(182+212+322)C=602=(722+922+562+1162)/4C=504SSe=SST SSt=602504=98第三十三张,PPT
23、共一百一十九页,创作于2022年6月2、自由度的分解tTedfdfdf-=-=tkdf1Tkndf-=1=441=15=41=3=153=12第三十四张,PPT共一百一十九页,创作于2022年6月二、二、列列出出方差分析表方差分析表,作,作F测验测验表表7.5 药剂处理后水稻苗高的方差分析表药剂处理后水稻苗高的方差分析表MS168.008.17F20.56显著F值F0.05=3.49F0.01=5.95*第三十五张,PPT共一百一十九页,创作于2022年6月三、若三、若F检验显著,则进行检验显著,则进行多重比较多重比较1、最小显著差数最小显著差数法(法(LSD法法)到底哪些处理间存在真实差异?
24、如何判断?到底哪些处理间存在真实差异?如何判断?多重比较多重比较第三十六张,PPT共一百一十九页,创作于2022年6月将这一判断标准记作:将这一判断标准记作:第三十七张,PPT共一百一十九页,创作于2022年6月LSD0.05=2.1792.02=4.40(cm)LSD0.01=3.0552.02=6.17(cm)第三十八张,PPT共一百一十九页,创作于2022年6月2、最小显著极差最小显著极差法(法(LSR法)法)Least Significant Ranges根据根据极差极差抽样分布原理,将一组抽样分布原理,将一组k个平均数由大个平均数由大到小排列后,依所比较的两个处理平均数的差数是几到小
25、排列后,依所比较的两个处理平均数的差数是几个平均数间的极差分别确定最小显著极差个平均数间的极差分别确定最小显著极差LSR值的值的多重比较方法。多重比较方法。用于比较任意两个处理间的平均数的差异显著性,两平用于比较任意两个处理间的平均数的差异显著性,两平均数的差数的绝对值大于或等于均数的差数的绝对值大于或等于LSR则差异显著差异显著与与LSD法相比,不同处理间的差数进行比较时需要根据差数所法相比,不同处理间的差数进行比较时需要根据差数所包含的平均数个数包含的平均数个数(m)不同确定不同的比较标准不同确定不同的比较标准第三十九张,PPT共一百一十九页,创作于2022年6月LSR法根据其比较标准此方
26、法包括两种:法根据其比较标准此方法包括两种:复极差法复极差法(q法)法)新复极差法新复极差法(SSR法,法,shortest significant ranges)*这里的q、SSR由、df、m三因素确定。*SE指标准误,随检验的对象不同(平均数或总指标准误,随检验的对象不同(平均数或总和数)而不同。和数)而不同。第四十张,PPT共一百一十九页,创作于2022年6月第四十一张,PPT共一百一十九页,创作于2022年6月 计算计算SE(比较四种药剂处理的平均数)(比较四种药剂处理的平均数)据据ve=12查附表查附表6得得SSR,计算,计算LSR,并列表并列表第四十二张,PPT共一百一十九页,创作
27、于2022年6月 按从大到小的顺序排列各处理平均数,并按从大到小的顺序排列各处理平均数,并 用应定方法标识其差异显著性用应定方法标识其差异显著性P=4P=3P=2P=2P=3P=2第四十三张,PPT共一百一十九页,创作于2022年6月字母标记法字母标记法 原则:原则:凡是两个平均数无共同字母则表示差异显著,凡是两个平均数无共同字母则表示差异显著,只要有一个字母相同就说明其间差异不显著只要有一个字母相同就说明其间差异不显著字母标记法标记步骤字母标记法标记步骤(0.050.05标记小写字母,标记小写字母,0.010.01标记大写字母标记大写字母):a)首先将全部平均数从大到小依次排列首先将全部平均
28、数从大到小依次排列b)在最大的平均数标上字母在最大的平均数标上字母ac)将该平均数与以下的平均数依次比较,凡将该平均数与以下的平均数依次比较,凡不显著的都标上相同字母不显著的都标上相同字母a,直到显著,标,直到显著,标上下一字母上下一字母b3、多重比较结果的表示方法多重比较结果的表示方法第四十四张,PPT共一百一十九页,创作于2022年6月d)以标上字母以标上字母b的平均数作为标准,从下向上与的平均数作为标准,从下向上与上方各自比其大的平均数比较,凡是不显著上方各自比其大的平均数比较,凡是不显著的都标以的都标以b,直到显著为止,直到显著为止注意:向上比较显著时不标新字母,向下比较注意:向上比较
29、显著时不标新字母,向下比较显显著时著时要要标以新字母标以新字母e)再以标有再以标有b的最大平均数为标准,向下依次的最大平均数为标准,向下依次比较,直到显著,标以字母比较,直到显著,标以字母c,再以标,再以标c平均平均数向上比较直到显著为止数向上比较直到显著为止f)重复上述步骤,直到所有平均数都比较完重复上述步骤,直到所有平均数都比较完第四十五张,PPT共一百一十九页,创作于2022年6月用字母标记法标记用字母标记法标记四种药剂试验的平四种药剂试验的平均数差异关系均数差异关系abcAcABBCC第四十六张,PPT共一百一十九页,创作于2022年6月四种药剂处理效应间存在显著差异四种药剂处理效应间
30、存在显著差异在在0.050.05显著性水平上有显著性水平上有D D与其它药剂之间与其它药剂之间差异显著,差异显著,B B药剂显著高于药剂显著高于A A和和C C药剂,而药剂,而A A和和C C药剂间差异不显著药剂间差异不显著而在而在0.010.01显著性水平上,显著性水平上,D D与与A A、C C药剂处理药剂处理间差异极显著,间差异极显著,B B药剂则极显著高于药剂则极显著高于C C药剂药剂的处理效果。的处理效果。四、四、结果结果的的解释解释/说明说明第四十七张,PPT共一百一十九页,创作于2022年6月另另一一例例aabbbccccccdddd第四十八张,PPT共一百一十九页,创作于202
31、2年6月a)凡是与对照比较或与预定的比较对象的比较,凡是与对照比较或与预定的比较对象的比较,一般可选用一般可选用LSD法法关于多重比较方法的选择关于多重比较方法的选择b)选用时考虑到否定一个正确的选用时考虑到否定一个正确的H0和接受一和接受一个不正确的个不正确的H0的相对重要性的相对重要性第四十九张,PPT共一百一十九页,创作于2022年6月显著尺度:显著尺度:LSD法最低,法最低,SSR法次之,法次之,q法最高法最高犯第一类错误概率:犯第一类错误概率:LSD法最大,法最大,SSR法次之,法次之,q法法最小最小*所以,当试验结论事关重大或要求较严格时,所以,当试验结论事关重大或要求较严格时,宜
32、采用宜采用q测验,一般试验可采用测验,一般试验可采用LSD法法第五十张,PPT共一百一十九页,创作于2022年6月第三节 数学模型一、一、方差分析的模型方差分析的模型据其中据其中 i 的不同可将模型分为两类:的不同可将模型分为两类:(1)固定模型固定模型(2)随机模型随机模型第五十一张,PPT共一百一十九页,创作于2022年6月1、固定模型固定模型 指试验的各处理均抽自特定的处理总体,这些总体遵循指试验的各处理均抽自特定的处理总体,这些总体遵循N(m mi,s s )的分布,处理效应的分布,处理效应 i=m mi m m 是固定的是固定的,试验的目的在于试验的目的在于研究研究 i 的大小,如果
33、重复做试验,所用处理不变,试验的假设的大小,如果重复做试验,所用处理不变,试验的假设为:为:m mm m=iiH或或:00e2*在在 F 测验测验 H0 被否定后,进一步的工作是多重比较,用于比被否定后,进一步的工作是多重比较,用于比较每个处理的效应大小较每个处理的效应大小第五十二张,PPT共一百一十九页,创作于2022年6月*每次重复试验时所用材料和试验处理都是固定的,即每次重复试验时所用材料和试验处理都是固定的,即 i =m mi-m m是固定的是固定的*固定模型所做出的推断仅限于试验所用处理范围之内固定模型所做出的推断仅限于试验所用处理范围之内,不可推而广之,不可推而广之*用于效应比较一
34、类的试验,包括多数栽培试验(如肥料、用于效应比较一类的试验,包括多数栽培试验(如肥料、密度、农药密度、农药.)以及很多室内效应比较试验等)以及很多室内效应比较试验等第五十三张,PPT共一百一十九页,创作于2022年6月 指试验的各处理是从同一个总体指试验的各处理是从同一个总体N(0,)抽得的一组随抽得的一组随机样本,机样本,i =m mi m m 是随机变量,随试验的不同而不同。试验的是随机变量,随试验的不同而不同。试验的目的不在于研究目的不在于研究 i的大小,而在于研究的大小,而在于研究 i的变异度的变异度 ;如果如果重复做试验,所用处理会改变。重复做试验,所用处理会改变。试验的假设为:试验
35、的假设为:2、随机模型随机模型2 s s00220=s ss s:AHH*在在 F 测验测验 H0 被否定后,进一步的工作是被否定后,进一步的工作是是计算是计算 的大小,然后估计遗传力等遗传参数的大小,然后估计遗传力等遗传参数2 s s2 s s第五十四张,PPT共一百一十九页,创作于2022年6月*每次重复试验时都要再次随机抽取样本,因此每次重复试验时都要再次随机抽取样本,因此 i =m mi m m 会发生变化会发生变化*随机模型所做出的推断不仅限于试验所用处理,其目的在于随机模型所做出的推断不仅限于试验所用处理,其目的在于推断抽出这些处理的整个总体,因而其推断可以在一定条推断抽出这些处理
36、的整个总体,因而其推断可以在一定条件下推而广之件下推而广之*在育种和生态一类的试验中应用较为广泛,目的在于总在育种和生态一类的试验中应用较为广泛,目的在于总体特征的研究体特征的研究第五十五张,PPT共一百一十九页,创作于2022年6月两类模型的比较两类模型的比较第五十六张,PPT共一百一十九页,创作于2022年6月判断两类模型的标准:判断两类模型的标准:(1)(1)看处理因素的水平是否可以完全控制:看处理因素的水平是否可以完全控制:随机模型的处理因素水平不可以控制随机模型的处理因素水平不可以控制固定模型的处理因素水平可以完全控制固定模型的处理因素水平可以完全控制温度菌落生长的影响化肥植物产量的
37、提高农家肥植物产量的提高第五十七张,PPT共一百一十九页,创作于2022年6月判断两类模型的标准判断两类模型的标准:*欲通过试验看一个地区辣椒产量的整体水平欲通过试验看一个地区辣椒产量的整体水平(2)看研究目的:看研究目的:固定模型:特定个体的比较固定模型:特定个体的比较随机模型:总体变异情况的研究随机模型:总体变异情况的研究*欲通过试验从一批辣椒品种中选出几个产量水平较欲通过试验从一批辣椒品种中选出几个产量水平较高的品种在某地区推广高的品种在某地区推广如前面四种药剂试验的例子如前面四种药剂试验的例子第五十八张,PPT共一百一十九页,创作于2022年6月 第三节 单因素方差分析适用于适用于单因
38、素完全随机试验设计单因素完全随机试验设计第五十九张,PPT共一百一十九页,创作于2022年6月一、一、每个处理每个处理观察值观察值数目数目相等相等的单的单向分组资料向分组资料第六十张,PPT共一百一十九页,创作于2022年6月【例例】表 水稻盆栽施肥试验产量结果(g/盆)第六十一张,PPT共一百一十九页,创作于2022年6月表表 不同肥料处理后水稻产量的方差分析表不同肥料处理后水稻产量的方差分析表LSR:LSD:第六十二张,PPT共一百一十九页,创作于2022年6月资料的LSR值(n=15)多重比较结果aabbbAAABBBCCc第六十三张,PPT共一百一十九页,创作于2022年6月资料的LS
39、R值(n=15)多重比较结果aaabbbcAAABBBCC第六十四张,PPT共一百一十九页,创作于2022年6月LSD:多重比较结果v=15时,时,t0.05=2.131;t0.01=2.947LSD0.05=3.91;LSD0.01=5.40第六十五张,PPT共一百一十九页,创作于2022年6月二、每个处理每个处理观察值观察值数目数目不不相等相等的单的单向分组资料向分组资料第六十六张,PPT共一百一十九页,创作于2022年6月每个处理观察值数目不相等时每个处理观察值数目不相等时平方和平方和与与自由度自由度的分解的分解第六十七张,PPT共一百一十九页,创作于2022年6月每个处理观察值数目不相
40、等时的每个处理观察值数目不相等时的多重比较多重比较第六十八张,PPT共一百一十九页,创作于2022年6月 5 5个不同品种羊的育肥试验,后期个不同品种羊的育肥试验,后期3030天增重天增重(kg)(kg)如表如表7.127.12所示。试比较品种间增重有无差异。所示。试比较品种间增重有无差异。【例例】第六十九张,PPT共一百一十九页,创作于2022年6月1、计算各项平方和与自由度、计算各项平方和与自由度第七十张,PPT共一百一十九页,创作于2022年6月2、列方差分析表,进行、列方差分析表,进行F测验测验方差分析显示方差分析显示F=5.99F0.01(4,20),P0.01,表明品种间,表明品种
41、间差异极显著,此试验为固定模型,下一步需要进行品种差异极显著,此试验为固定模型,下一步需要进行品种间的多重比较间的多重比较第七十一张,PPT共一百一十九页,创作于2022年6月3、多重比较、多重比较假设用假设用SSR法(可用其它多重比较方法)法(可用其它多重比较方法)因为各处理重复数不等,应先计算出平均重复次数因为各处理重复数不等,应先计算出平均重复次数n0来代来代替标准误替标准误SE中的中的n 第七十二张,PPT共一百一十九页,创作于2022年6月羊品种平均增重多重比较梯形表羊品种平均增重多重比较梯形表 第七十三张,PPT共一百一十九页,创作于2022年6月三、系统分组资料三、系统分组资料第
42、七十四张,PPT共一百一十九页,创作于2022年6月第七十五张,PPT共一百一十九页,创作于2022年6月系统分组资料系统分组资料平方和与自由度平方和与自由度的的分解分解第七十六张,PPT共一百一十九页,创作于2022年6月二级分组资料的二级分组资料的方差分析方差分析及其及其期望均方期望均方*若组间差异显著,则进行多重比较若组间差异显著,则进行多重比较*若亚组间差异不显著,可将亚组间和亚组内的平方若亚组间差异不显著,可将亚组间和亚组内的平方和进行合并,求其合并均方:和进行合并,求其合并均方:*首先对亚组间进行首先对亚组间进行F 测验测验*若亚组间差异显著,然后对组间进行若亚组间差异显著,然后对
43、组间进行F 测验测验第七十七张,PPT共一百一十九页,创作于2022年6月系统分组资料的多重比较系统分组资料的多重比较F测验中的备测验中的备比量比量被比较平均数所包含的观被比较平均数所包含的观察值个数察值个数第七十八张,PPT共一百一十九页,创作于2022年6月【例例】表 4种培养液下株高增长量(mm)第七十九张,PPT共一百一十九页,创作于2022年6月【例例】资料的方差分析表资料的方差分析表第八十张,PPT共一百一十九页,创作于2022年6月第四节 两向分组资料的方差分析两向分组两向分组(交叉分组交叉分组)资料资料:试验包括两个因素试验包括两个因素两个因素的每个水平均衡相遇两个因素的每个水
44、平均衡相遇其水平组合总数即为处理数其水平组合总数即为处理数按组合观察值是否有重复其方差分析分两种情况按组合观察值是否有重复其方差分析分两种情况组合内有重复观察值组合内有重复观察值组合内无重复观察值组合内无重复观察值第八十一张,PPT共一百一十九页,创作于2022年6月A因素水平数为因素水平数为i=1,2,aB因素水平数为因素水平数为j=1,2,b共有处理组合数共有处理组合数ab个个,即观察值总数为即观察值总数为ab个个一、组合内无重复观测值的两向资料一、组合内无重复观测值的两向资料x.1xa.x.bx.jx.2x1.xi.x2.x.第八十二张,PPT共一百一十九页,创作于2022年6月1、模型
45、模型第八十三张,PPT共一百一十九页,创作于2022年6月2、平方和与自由度的分解平方和与自由度的分解第八十四张,PPT共一百一十九页,创作于2022年6月3、方差分析及其期望均方方差分析及其期望均方第八十五张,PPT共一百一十九页,创作于2022年6月4、多重比较多重比较第八十六张,PPT共一百一十九页,创作于2022年6月【例例7.9】(P119)采用采用5种生长素处理豌豆,未处理为对照,待种子发种生长素处理豌豆,未处理为对照,待种子发芽后,分别每盆中移植芽后,分别每盆中移植4 组,每组组,每组6盆,试验共有盆,试验共有4组组24盆,盆,并按组排于温室中,使同组各盆条件一致。当各盆见第并按
46、组排于温室中,使同组各盆条件一致。当各盆见第一朵花时记录总结间数,见下表:一朵花时记录总结间数,见下表:第八十七张,PPT共一百一十九页,创作于2022年6月【例7.9】资料的方差分析表第八十八张,PPT共一百一十九页,创作于2022年6月二、组合内有重复观测值的两向资料二、组合内有重复观测值的两向资料A因素水平数为因素水平数为i=1,2,aB因素水平数为因素水平数为j=1,2,b每个处理重复观察值数为每个处理重复观察值数为k=1,2,n共有处理组合数共有处理组合数ab个个,则观察值总数为则观察值总数为abn个个第八十九张,PPT共一百一十九页,创作于2022年6月1、模型模型第九十张,PPT
47、共一百一十九页,创作于2022年6月2、平方和与自由度的分解平方和与自由度的分解第九十一张,PPT共一百一十九页,创作于2022年6月3、方差分析及其期望均方方差分析及其期望均方因不同模型因不同模型而异而异第九十二张,PPT共一百一十九页,创作于2022年6月4、多重比较多重比较第九十三张,PPT共一百一十九页,创作于2022年6月【例例】表7.15 3种肥料施于3种土壤的小麦产量第九十四张,PPT共一百一十九页,创作于2022年6月表7.16 表7.15资料的方差分析表肥类平均数的比较组合平均数的比较第九十五张,PPT共一百一十九页,创作于2022年6月表7.17 表7.15资料的多重比较(
48、肥类均值)第九十六张,PPT共一百一十九页,创作于2022年6月表7.18 表7.15资料的多重比较(组合均值)第九十七张,PPT共一百一十九页,创作于2022年6月如互作不存在,两个因素各最好水平的组合为如互作不存在,两个因素各最好水平的组合为最优处理最优处理如存在互作,两个因素各最好水平的组合不一如存在互作,两个因素各最好水平的组合不一定为最优处理定为最优处理第九十八张,PPT共一百一十九页,创作于2022年6月第五节 方差分析的基本假定和数据转换一、方差分析的三个基本假定方差分析的三个基本假定第九十九张,PPT共一百一十九页,创作于2022年6月1、处理效应和环境处理效应和环境(误差误差
49、)效应具有效应具有“可加性可加性”第一百张,PPT共一百一十九页,创作于2022年6月SST=SSt+SSe第一百零一张,PPT共一百一十九页,创作于2022年6月效应的效应的可加性可加性与与倍加性倍加性处理与环境彼此独立,处理与环境彼此独立,其效应是线性可加的其效应是线性可加的处理与环境不是彼此独立的,处理与环境不是彼此独立的,其效应非线性可加其效应非线性可加第一百零二张,PPT共一百一十九页,创作于2022年6月2、试验误差是独立的随机变量,且符合、试验误差是独立的随机变量,且符合N(0,2)“正态性正态性”*误差的随机性是误差的随机性是F测验的基本假定测验的基本假定*ij 相互之间是独立
50、的相互之间是独立的 严格遵循试验设计的随机性严格遵循试验设计的随机性 要杜绝系统误差,严格试验管理与环要杜绝系统误差,严格试验管理与环境控制境控制*ij 与与 i 之间也是独立的,如存在关联,则会表之间也是独立的,如存在关联,则会表现为现为 ij 不符合正态分布不符合正态分布第一百零三张,PPT共一百一十九页,创作于2022年6月3、所有试验处理必须具有共同的误差方差、所有试验处理必须具有共同的误差方差“同质性同质性”即假定各处理的即假定各处理的 ij N(0,s s 2 2)如果计算所得处理内方差相差比较悬殊,可以通过如果计算所得处理内方差相差比较悬殊,可以通过F测测验或验或 2 2检验其是