多因素实验设计的方差分析

《多因素实验设计的方差分析》由会员分享，可在线阅读，更多相关《多因素实验设计的方差分析（37页珍藏版）》请在装配图网上搜索。

1、,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,多因素实验资料的方差分析,SPSS,实现,邹莉玲,PH.D,同济大学医学院,主要内容,（一）析因设计的方差分析,（二）正交设计的方差分析,（三）嵌套设计的方差

2、分析,（四）裂区设计的方差分析,析因设计,（,factorial design,）,：,是将两个或多个实验因素的各水平进行组合，对各种可能的组合都进行实验，从而探讨各实验因素的,主效应,（,main effect,），以及各因素间的,交互作用,（,interaction,）的研究设计类型。,例,1,为最简单的析因设计，即两因素两水平，记作,2,2,或,22,的析因设计,。,（一）析因设计的方差分析,实例分析,例,1,：某研究人员为了解,升白细胞药物（,A,）,和,纯苯（,B,）,对大鼠吞噬指数的影响，以及两者同时使用的作用。将,20,只性别相同、体重相近的大鼠，按,A,、,B,两因素有无分为,

3、a,1,b,1,、,a,1,b,2,、,a,2,b,1,、,a,2,b,2,四组（,1,表示用药，,2,表示不用）。测得吞噬指数结果见表,9-1,。,若,i,：表示因素,A,的水平（,i=1,，,2,，,，,a,），,j,：表示因素,B,的水平（,j=1,，,2,，,，,b,），,k,：表示因素,A,和因素,B,各水平组合下的观察单位数（,k,=1,，,2,，,，,n,）。,单独效应、主效应和交互效应,单独效应：,指其他因素水平固定时，同一因素不同水平的差异。如,A,因素固定在,1,水平时，,B,因素的单独效应为,-1.8100,。,主效应：,某一因素各水平间的差异。如,A,的主效应为,0.0

4、53,。,交互效应：,当某一因素的各单独效应随另一因素变化而变化时，称这两个因素间存在交互效应。如,AB,的交互效应：,AB=(a1b1-a2b1)-(a1b2-a2b2)/2=(0.0960-0.0100)/2=0.0430,。,研究目的,当研究的因素不止一个时，这种研究设计就称为多因素的实验设计。其方法有很多种，析因设计就是其中的一种。,研究目的,：,不仅分析单个因素不同水平效应之间的差异，还要知道两个因素各水平间效应的相互影响。,分析方法：,采用多因素方差分析。,变异分解：,自由度的分解：,方差分析的基本思想,方差分析的基本步骤,建立检验假设，确定检验水准,对于因素,A,（升白细胞药物）

5、：,H,0,：给药与不给药的大鼠吞噬指数的总体均数相等,H,1,：给药与不给药的大鼠吞噬指数的总体均数不等,对于因素,B,（纯苯染毒）：,H,0,：染毒与不染毒的大鼠吞噬指数的总体均数相等,H,1,：染毒与不染毒的大鼠吞噬指数的总体均数不等,对于交互作用,AB,：,H,0,：因素,A,和因素,B,无交互作用,H,1,：因素,A,和因素,B,有交互作用,均取,计算检验统计量,表,9-4,例,1,的两因素析因设计方差分析表,确定,P,值，作出推断结论,（,1,）,AB,交互效应的,P0.05,，提示按,0.05,的检验水准，接受,H,0,假设，即还不能认为,AB,两因素间存在交互作用。,（,2,）

6、,A,因素主效应的,P0.05,，提示不能认为给予升白细胞药物对大鼠吞噬细胞指数有影响。,（,3,）,B,因素主效应的,PGeneral Linear Model-Univariate,单变量,：,-,因变量名称：,y,-,固定因子（处理因素）：,A,、,B,2,）,点击,“,模型,”,按钮，弹出重复度量模型对话框。,-,模型：,全因子,，即分析所有主效应及交互效应（系统默认）。,-,平方和：类型,III,（系统默认）。,3,）,点击,“,图,/plot,”,按钮，弹出交互作用轮廓图对话框。,交互作用轮廓图是将各因素不同水平组合的均值在二维图形上标出，以直观描述交互效应。,-,水平轴：因素,A

7、,-,单图（线段，,separate lines,）：因素,B,-,多图（分图，,separate plots,）：无,4,）,Post Hoc,（对比）按钮：,用于某处理因素多个水平间的多重比较。本例的研究因素均为两水平，所以无需此步骤。,5,）,Options,（选项）按钮,-,显示均值：输出所选因素的均数、标准误、可信区间。,-,输出：包括描述统计、参数估计、方差齐性检验等供选择项。本例不选择任何选项。,主要输出结果,1,）,均值估计,：,2,）,方差分析表,：包括处理因素主效应和交互效应比较。,3,）,交互效应轮廓图,：,结论：总的模型拟合效果理想，,R,2,=0.990,。,（,1,

8、）,AB,交互效应的,P0.05,，提示按,0.05,的检验水准，接受,H,0,假设，即还不能认为,AB,两因素间存在交互作用。,（,2,）,A,因素主效应的,P0.05,，提示不能认为给予升白细胞药物对大鼠吞噬细胞指数有影响。,（,3,）,B,因素主效应的,P0.01,，提示染毒对吞噬指数有影响，可以降低大鼠吞噬指数。,交互效应轮廓图中，两条直线几乎平行，提示,A,、,B,两因素的交互效应不显著。反之，若两条直线交叉，则提示可能存在交互效应。,小 结,析因设计的,优点,：,全面高效性，可以对各因素的不同水平进行组合，对各因素不同水平主效应进行分析的同时，还可以对交互效应进行分析；通过比较各实

9、验组合，还可以寻求最佳组合。,析因设计的,缺点,：,工作量大，含有较多因素和水平的实验一般不用完全交叉分组的析因设计，而采用非全面试验的正交设计，可以大幅度减少实验次数。,两因素析因设计与随机区组设计方差分析的区别：,后者每个组合下的数据无重复，不能分析交互效应。,析因设计资料分析：应先分析交互效应。,若交互效应有统计学意义，要逐一分析各因素的单独效应，即固定一个因素对其他因素进行分析；反之，若交互效应无统计学意义，则因素间的作用相互独立，直接分析各因素的主效应。,（二）正交设计的方差分析,正交设计：,是利用一套规范化的,正交表,，使每次试验的各因素及其水平得到合理安排的,一种,高效,、,多因

10、素,实验设计。,常用于：寻找疗效好的药物配方、医疗仪器多个参数的优化组合、生物体的培养条件等最优搭配方案的研究。,方差分析基本思想,将总的离均差平方和（,SS,）分解成各因素及各交互作用的离均差平方和，构造,F,统计量，对各因素是否对实验效应指标具有显著影响作,F,检验。,例,11-4,（,P202,）：研究雌螺产卵的最优条件，在,20cm,2,的泥盒里饲养同龄雌螺,10,只，试验条件有,4,个因素（表,11-15,），每个因素,2,个水平。试在考虑温度与含氧量对雌螺产卵有交互作用的情况下安排正交试验。,表,11-15,影响雌螺产卵数的试验因素及其水平,水平,试 验 因 素,温度,(,),含氧

11、量,(,),含水量,(,),pH,值,A,B,C,D,1,5,0.5,10,6.0,2,25,5.0,30,8.0,SPSS,操作过程：,（,1,）建立数据文件（例,11-04.sav,）,（,2,）过程：,分析,analyze,常规线性模型,General Linear Model,单变量,Univariate,因变量：,X,固定因子：,a,b,ab,c,d,模型：定制,a,b,ab,c,d,（主效应）,平方和：类型,III,，,在模型中包含截踞,SPSS,结果输出：,（三）嵌套设计的方差分析,嵌套设计：一个因素的不同水平分别与另一个因素的不同水平发生组合。或者说,B,因素的不同水平是嵌套在

12、,A,因子内的。,实例分析,例,11-6,（,P208,）：实验甲乙丙三种催化剂在不同温度下对某化合物的转化作用。将催化剂作为一级实验因素（,I=3,），温度作为二级实验因素（,J=3,），采用嵌套设计，每个处理重复,2,次（,n=2,）,实验结果见表,11-25,。试做方差分析。,表,11-25,某化合物的转化率（,%,）,催化剂,温度（）,A,B,C,70,80,90,55,65,75,90,95,100,实验,结果（,X,）,82,91,85,65,62,56,71,75,85,84,88,83,61,59,60,67,78,89,实例,11-6,数据的,方差分析表,表,11-26,某化

13、合物转化率的方差分析表,变异来源,自由度,SS,MS,F,P,总变异,17,2406.5,催化剂,2,1956.0,978.0,177.82,0.01,温度,/,催化剂,6,401.0,66.8,12.15,0.01,误差,9,49.5,5.5,SPSS,操作过程：,（,1,）建立数据文件（例,11-06.sav,）,（,2,）过程：,分析,analyze,常规线性模型,General Linear Model,单变量,Univariate,因变量：,X,（转化率,%,）,固定因子：,a,（催化剂）,b,（温度）,模型：定制,a,b,（主效应）,平方和：,类型,I,，,在模型中包含截踞,SPS

14、S,结果输出：,（四）裂区设计的方差分析,裂区设计：,析因设计的一种特殊形式。是将全区试验分解成多个裂区组安排析因处理。,实验设计方法：,1.,完全随机,IJ,裂区设计,（,1,）将一级实验单位随机等分成,I,组，每组例数为,n,（,2,），分别接受,a,1,a,2,a,3,a,i,各水平的处理。,（,2,）分别将各一级实验单位内的二级实验单位随机分配接受,b,1,b,2,b,3,b,j,的处理。,2.,随机区组,IJ,裂区设计：（略）,实例分析,例,11-7,（,P211,）：一种全身注射抗毒素对皮肤损伤的保护作用实验，将,10,只,家兔,随机等分为两组，一组注射抗毒素，一组注射生理盐水作对

15、照。分组后，每只家兔取,甲乙两部位,分配注射低浓度和高浓度毒素，观察指标为皮肤受损直径（,mm,），实验结果见表,11-31,。试作方差分析。,表,11-31,家兔皮肤损伤直径（,mm,）,注射药物（,A,）,随机化后家兔编号,毒素浓度（,B,）,家兔小计,低浓度（,b1,）,高浓度（,b2,）,抗毒素（,a1,）,1,15.75,19.00,34.75,生理盐水（,a2,）,6,合计,10,179.00,214.50,393.50,表,11-32,家兔皮肤损伤直径的方差分析表,变异来源,自由度,SS,MS,F,P,两级单位总计,（,19,）,（,146.1735),家兔间,（,9,）,(81

16、.0125),注射药物（,A,）,1,63.0125,63.0125,28.01,0.01,一级单位误差,8,18.0000,2.2500,部位间,（,10,）,(65.1250),毒素浓度（,B,）,1,63.0125,63.0125,252.05,0.05,二级单位误差,8,2.0000,0.2500,SPSS,操作过程：,（,1,）建立数据文件（例,11-07.sav,）,（,2,）过程：,分析,analyze,常规线性模型,General Linear Model,单变量,Univariate,因变量：,diameter,（损伤直径）,固定因子：,drug,（药物）,cons,（浓度）,随机因子：,No,（家兔编号）,模型：定制,drug,cons,drug*cons,平方和：类型,III,，,在模型中包含截踞,点击粘贴,Paste:,将弹出的程序中的最后一行改为：,“,/DESIGN=drug No(drug)cons cons*drug.,”,SPSS,结果输出：,