第4章 方差分析

.,第四章多个样本均数比较的方差分析,.,2020/4/27,.,3,主要内容,第一节方差分析的基本概念第二节完全随机设计的方差分析第三节随机区组设计的方差分析第四节拉丁方设计资料的方差分析第五节两阶段交叉设计资料的方差分析第六节多个样本均数间的多重比较第七节多样本方差齐性检验,.,第一节方差分析的基本概念,什么是方差？离均差离均差平方和SS方差（2S2）均方（MS）标准差：S自由度：关系：MS=SS/,一、方差分析的几个名词,.,原始数据对于总均数的离差平方和可分解成几个部分，每一部分代表一种变异的来源，自由度也可得到相应的分解。通过对不同来源变差的均方的比较，以判断某种变异原因的存在与否。,二、方差分析的基本思想,.,三、统计量F的计算及其意义,F=MS组间/MS组内自由度：组间=组数-1组内=N-组数通过这个公式计算出统计量F，查表求出对应的P值，与进行比较，以确定是否为小概率事件。,.,第二节完全随机设计的方差分析(单因素方差分析),.,何谓完全随机设计,是一种将研究对象随机地分配到处理因素各水平组的单因素设计方法,其研究目的是推断处理因素不同水平下的试验结果的差异是否有统计学意义,即该处理因素是否对试验结果有本质影响.,.,例4-1某医师为了研究一种降血脂新药的临床疗效，按统一纳入标准选择120名患者，采用完全随机设计方法将患者等分为4组进行双盲试验。问如何进行分组？,分组方法：先将120名高血脂患者从1开始到120编号，见表4-2第一行；从随机数字表中的任一行任一列开始，如第5行第7列开始，依次读取三位数作为一个随机数录于编号下见表4-2第二行；然后将全部随机数从小到大编序号（数据相同的按先后顺序编序号），将每个随机数对应的序号记在表4-2第三行；规定序号1-30为甲组，序号31-60为乙组，序号61-90为丙组，序号91-120为丁组，见表4-2第四行。,.,表4-2120名患者等分为4组的完全随机设计分组结果,.,例4-2某医生为了研究一种降血脂新药的临床疗效，按统一纳入标准选择120名高血脂患者，采用完全随机设计方法将患者等分为4组（具体分组方法见例4-1），进行双盲试验。6周后测得低密度脂蛋白作为试验结果，见表4-3。问4个处理组患者的低密度脂蛋白含量总体均数有无差别？,.,表4-34个处理组低密度脂蛋白测量值（mmol/L）,.,Xij第i个组的第j个观察值I=1,2,gJ=1,2,nini第i个处理组的例数ni=N,各种符号的意义,.,数据结构,.,.,，即4个试验组的整体均数相等H1：4个试验组的整体均数不全相等按表4-4中的公式计算各离均差平方和SS、自由度、均方MS和F值。=102.91+81.46+80.94+58.99=324.30=367.85+233.00+225.54+132.13=958.52C=(324.30)2/120=876.42SS总=958.52-876.42=82.10,=120-1=119,.,方差分析表见表4-5。按水准,=3,=116查附表3的F界值表,得结论:按水准,拒绝H0,接受H1,认为4个处理组患者的低密度脂蛋白总体均数不全相等,即不同剂量药物对血脂中低密度脂蛋白降低有影响.,.,表4-5例4-2的方差分析表表4-4完全随机设计资料的方差分析表,.,在满足下列条件下：1）各处理组所代表的总体服从正态分布；2）各处理组所代表的总体相互独立；3）各处理组所代表的总体方差齐性；4）各处理组所代表的总体均数相等。,均方（MS)之比的意义,.,F分布及其上侧F界值,.,第三节随机区组设计资料的方差分析,也叫配伍组设计的多个样本均数的比较（双因素方差分析）,应用分层的思想，事先将全部受试对象按某种或某些特性分为若干个区组，区组内的试验对象均衡，区组之间试验对象具有较大的差异，这样利用区组控制非处理因素的影响，并在方差分析时将区组间的变异从组内变异中分解出来。当区组间差别有统计学意义时，这种设计的误差比完全随机设计小，试验效率得以提高。是配对（两个组）资料的扩充（两个以上组）。,.,例4-3如何按随机区组设计,分配5个区组的15只小白鼠接受甲、乙、丙三种抗癌药物？,方法为先将小白鼠的体重从轻到重编号，体重相近的3只小白鼠配成一个区组，见表4-6。然后在随机数字表（附表15）中任选一行一列开始的2位数作为1个随机数，如从第8行第3列开始纪录，见表4-6；在每个区组内将随机数按大小排序，各区组中内序号为1的接受甲药，序号为2的接受乙药，序号为3的接受丙药，分配结果见表4-6。,.,表4-615只小白鼠分5个区组的随机区组设计分配结果,.,例4-4某研究者采用随机区组设计进行实验，比较三种抗癌药物对小白鼠肉瘤的抑瘤效果，先将15只染有肉瘤小白鼠按体重大小配成5个区组，每个区组内3只小白鼠随机接受三种抗癌药物（具体分配方法见例4-3），以肉瘤的重量为指标，试验结果见表4-9。问三种不同药物的抑瘤效果有无差别？,.,表4-9三种不同药物作用后小白鼠肉瘤重量(g),.,分析变异,总变异:组间变异（因素一的变异）:误差（组内）变异（随机误差）:配伍间变异（因素二的变异）:,.,统计量F的计算,F1=MS组间/MS组内F2=MS区间/MS组内自由度：组间=组数-1=3-1=2区间=区数-1=8-1=7组内=（组数-1）（区数-1）=14都是和随机误差进行比较。,.,H0：即三种不同药物作用后小白鼠肉瘤重量的总体均数相等H1：三种不同药物作用后小白鼠肉瘤重量的总体均数不全相等按表4-8中的公式计算各离均差平方和SS、自由度、均方MS和F值。C=6.812/15=3.0917SS总=3.6245-3.0917=0.5328,=15-1=14,.,表4-8随机区组设计资料的方差分析表,.,按水准,查附表3的F界值表,得按水准,拒绝H0,接受H1,认为三种不同药物作用后小白鼠肉瘤重量的总体均数不全相等,即不同药物的抑瘤效果有差别.,.,F分布及其上侧F界值,方差分析结果提供了各组均数间差别的总的信息，但尚未提供各组间差别的具体信息，即尚未指出哪几个组均数间的差别具有或不具有统计学意义。为了得到这方面的信息，可进行多个样本间的两两比较。,注意,.,第四节拉丁方设计资料的方差分析,.,完全随机设计只涉及到一个处理因素。随机区组设计涉及一个处理因素、一个区组因素（配伍因素）。倘若实验研究涉及一个处理因素和两个控制因素，且每个因素的水平数相等时，此时可采用拉丁方设计来安排试验。,.,例4-5某研究者为了比较甲、乙、丙、丁、戊、己6种药物给家兔注射后产生的皮肤疱疹大小（mm2），采用拉丁方设计，选用6只家兔，并在每只家兔的6个不同部位进行注射。试验结果见表4-11，试做方差分析。,.,表4-11例4-5的拉丁方设计与试验结果（皮肤疱疹大小，mm2）,.,，即6种药物注射后家兔产生皮肤疱疹大小的总体均数相等：6种药物注射后家兔产生皮肤疱疹大小的总体均数不全相等，即6只家兔皮肤疱疹大小的总体均数相等6只家兔皮肤疱疹大小的总体均数不全相等，即6个注射部位皮肤疱疹大小的总体均数相等6个注射部位皮肤疱疹大小的总体均数不全相等据例4-5的结果表4-11，按表4-12中的公式计算各离均差平方和SS、自由度、均方MS和F值，计算过程从略，将计算结果按表4-12汇总，得方差分析表，见表4-13。,.,拉丁方设计资料采用三向分类的方差分析，总变异可分解为处理组变异、行区组变异、列区组变异和误差4部分，其方差分析表见表4-12。表4-12拉丁方设计资料的方差分析表,.,表4-13例4-5结果(表4-11)的方差分析表,.,按水准，查附表3的F界值表，按水准，拒绝，接受，即6种药物注射后家兔产生皮肤疱疹大小的总体均数不全相等。按水准，不拒绝，还不能认为6只家兔皮肤疱疹大小的总体均数不全相等；按水准，不拒绝，还不能认为6个注射部位皮肤疱疹大小的总体均数不全相等。,.,第五节两阶段交叉设计资料的方差分析,.,两阶段交叉设计,.,例4-6表4-14是A、B两种闪烁液测定血浆中的交叉试验结果。第阶段1、3、4、7、9号用A测定，2、5、6、8、10号用B测定；第阶段1、3、4、7、9号用B测定，2、5、6、8、10号用A测定。试对交叉试验结果进行方差分析。表4-14两种闪烁液测定血浆中的交叉试验,.,.,将以上结果列成方差分析表，见表4-15。表4-15例4-6的方差分析表,结论：还不能认为A和B两种闪烁液的测定结果有差别；可认为测定阶段对测定结果有影响；可认为各受试者的值不同。,.,第六节多个样本均数间的多重比较,.,能否用t检验或u检验？4个样本均数两两比较,若用t检验,共需做6次比较,若每次不犯第一类错误的概0.95，则10次都不犯一类错误的概率不是0.95，而是：0.956=0.74,这样总的检验水准变为0.402，犯一类错误的概率较0.05有了较大的提高。,.,最小有意义t（LSD-t）检验：检验K组中某一对或几对在专业上有特殊意义的均数（DAB=XA-XB）的总体水平是否为0。,.,LSD-t检验(Leastsignificantdifference),.,例4-7对例4-2资料，问高血脂患者的降血脂新药2.4克组、4.8克组、7.2克组与安慰剂组的低密度脂蛋白含量总体均数有无差别?,降血脂新药2.4克组与安慰剂组的比较:,有统计学意义。可认为降血脂新药2.4g组的低密度脂蛋白含量总体均数低于安慰剂组。同样，对降血脂新药4.8g组、7.2g组,可算出LSD-t为-4.29、-8.59,都得P0.001。可认为降血脂新药4.8g组、7.2g组的低密度脂蛋白含量总体均数低于安慰剂组。,.,Dunnett-t检验:K-1个实验组与一个对照组均数差别的多重比较,.,例4-8对例4-2资料，问高血脂患者的三个不同剂量降血脂新药组与安慰剂的低密度脂蛋白含量总体均数是否有差别？,，,.,（双侧），得有统计学意义。可认为降血脂新药各组与安慰剂组的低密度脂蛋白含量总体均数有差别，即三个不同剂量降血脂新药都可降低血脂中低密度脂蛋白含量。,.,.,例4-9对例4-4资料，问三种不同药物的抑瘤效果两两之间是否有差别？,，将三个样本均数由小到大排列，并编组次：均数0.3140.4340.614组别C药B药A药组次123MS误差=0.0096,。,.,按水准，三个对比组中的1与3组、2与3组拒绝H0，接受H1，有统计学意义，可认为C药与A药间、B药与A药间抑瘤效果有差别。但1与2组不拒绝H0，还不能认为C药与B药间抑瘤效果有差别。,.,第七节多样本方差比较的Bartlett检验和Levene检验,.,方差分析的应用条件：各样本相互独立，均服从正态分布并且方差具有齐性(h