第11章多因素实验的方差分析方法:方差分析(F检验)目的:研究多个(包括两个)处理因素对试验对象的试验指标的作用资料:处理因素分几个水平(非定量),试验指标定量结果:应用:FactorialdesignANOVA析因设计的方差分析一、析因设计二、析因设计的方差分析(一)两因素两水平(二)两因素三水平(三)三因素多水平实例1:甲乙两药治疗高胆固醇血症的疗效(胆固醇降低值mg%),问①甲乙两药是否有降低胆固醇的作用?②两种药间有无交互作用完全随机的两因素2×2析因设计甲药乙药用不用用645678448042不用281631252318实例2:小鼠种别A、体重B和性别C对皮内移植SRS瘤细胞生长特征影响的结果(肿瘤体积cm3)问①A、B、C各自的主效应如何?②三者间有无交互作用?种别A体重(g)雄性雌性昆明种24~250.70690.18850.78540.34030.35810.250313~151.08380.95500.94250.92150.33350.8514泸白种24~250.06280.47120.09420.08800.04710.175913~150.01260.25130.00940.36760.01250.1327性别完全随机的三因素2×2×2析因设计实例3:研究小鼠在不同注射剂量和不同注射频次下药剂ACTH对尿总酸度的影响。问①A、B各自的主效应如何?②二者间有无交互作用?配伍组编号日注射量AB1(少)B2(多)1A133.633.0237.130.5334.133.3434.634.41A233.028.5229.531.8329.229.9430.728.31A331.430.7228.328.2328.928.4428.630.6注射次数B随机配伍的两因素3×2析因设计显著特征每个处理是各因素各水平的一种组合,总处理数为各因素各水平的全面组合数,即各因素各水平数的乘积。如两因素析因设计,设A因素有I个水平,B因素有J个水平,则总处理数G=I×J。在三个因素的析因设计中,若各因素水平为I、J、K,则总处理数G=I×J×K。要求各个处理组内的实验单位数相等(便于手工计算)且每组至少有两个实验单位,否则无法分析因素间的交互作用,故总的实验单位数至少为2G。2.实验设计各因素各水平的全面组合设有k个因素,每个因素有L1,L2,…,Lk个水平,那么共有g=L1×L2×…×Lk个处理组。例如有三个因素,分别是A,B,C。A因素有2水平,B因素有3水平,C因素有2水平,则处理组g=2×3×2=12个确定了处理组数后,将实验对象分配到各组的方法可以采用完全随机设计、随机区组设计或拉丁方设计析因设计(FactorialDesign)是一种多因素多水平交叉分组进行全面试验的设计方法。它可以研究两个或两个以上因素多个水平的效应。在析因设计中,研究因素的所有可能的水平组合都能被研究到。例如4个因素同时进行实验,每个因素取两个水平,实验的总组合数为24=16;如果水平为3,则有34=81种组合数。即是这81种组合均进行实验。所以析因设计可以分析观测指标与研究因素间的复杂关系,包括各因素间的交互作用(Interaction)。3.析因设计的特点2个以上(处理)因素(factor)(分类变量)2个以上水平(level)2个以上重复(repeat)每次试验涉及全部因素,即因素同时施加观察指标(观测值)为计量资料(独立、正态、等方差)4.析因设计的有关术语单独效应(simpleeffects):其它因素(factor)的水平(level)固定为某一值时,某一因素的效应。即其他因素的水平固定时,同一因素不同水平间的差别。主效应(maineffects):某因素各单独效应的平均效应。某一因素各水平间的平均差别。交互作用(Interaction):某一因素的各个单独效应随着另一因素变化而变化的情况。(如一级交互作用AB、二级交互作用ABC…)如果不存在交互效应,则只需考虑各因素的主效应。在方差分析中,如果存在交互效应,解释结果时,要逐一分析各因素的单独效应,找出最优搭配。在两因素析因设计时,只需考虑一阶交互效应。三因素以上时,除一阶交互效应外,还需考虑二阶、三阶等高阶交互效应,解释将更复杂。(A+B)=A+B+AB(A+B+C)=A+B+C+AB+AC+BC+ABC5.析因设计的优缺点优点:比单处理因素设计能提供更多的试验信息,可用来分析全部因素主效应,以及因素间各级的交互作用,在医学上可用于筛选最佳治疗方案、药物配方、实验条件等研究。缺点:当因素增加时,实验组数呈几何倍数增加,所需试验的次数很多。不但计算复杂,而且给众多交互作用的解释带来困难。因此,当因素及水平数较多时,一般不采用完全交叉分组的析因设计,而采用正交设计。(一)两因素两水平(2×2)完全随机析因设计的方差分析A缝合方法外膜缝合(a1)束膜缝合(a2)合计B缝合后时间1月b12月b21月b12月b210301050103020504070307050605060103030302444285212022014026074044001120048001440034800例11-1:研究不同缝合方法及缝合后时间对家兔轴突通过率(%)的影响,问①两种缝合方法间有无差别?缝合后时间长短间有无差别?②两者间有无交互作用iTix2ixA因素(2水平)缝合方法B因素(2水平)缝合后时间缝合后1月b1缝合后2月b2外膜缝合a124(a1b1)44(a1b2)束膜缝合a228(a2b1)52(a2b2)2因素2水平析因实验示意图表2因素2水平析因试验的均数差别B因素平均b2-b1b1b2a124443420a228524024平均264822a2-a1486A因素在a1b1、a1b2、a2b1和a2b2的四种处理组合中,每个格子均有5个数据,因此它又是重复数相等的析因设计。由于数据按因素A和因素B两个方向交叉分组,故可用双向方差分析。进一步分析处理的单独效应(simpleeffect)、主效应(maineffect)和交互效应(interaction)。单独效应:其他因素的水平固定时,同一因素不同水平间的差别A因素缝合方法B因素缝合后时间平均B单独效应1月(b1)2月(b2)(b2-b1)外膜(a1)24443420束膜(a2)28524024平均264822A单独效应(a2-a1)486表2-22×2析因设计均数差别分析析因试验的均数差别b因素平均b2-b1b1b2a124443420a228524024平均264822a2-a1486a因素a因素单独效应:比较缝合后两个时间点的轴突通过率a(b1固定)=(a2-a1)=28-24=4a(b2固定)=(a2-a1)=52-44=8b因素单独效应:b(a1固定)=(b2-b1)=20b(a2固定)=(b2-b1)=24主效应:某一因素各水平间的平均差别A因素的主效应解释为:束膜缝合与外膜缝合相比(不考虑缝合时间),神经轴突通过率提高了6%(40-34)。B因素的主效应解释为:缝合后2月与1月相比(不考虑缝合方法),神经轴突通过率提高了22%(48-26)。A因素缝合方法B因素缝合后时间平均B主效应1月(b1)2月(b2)外膜(a1)244434束膜(a2)285240平均264822A主效应62(交互作用)交互作用(interaction):当某因素的各个单独效应随另一因素水平的变化而变化,且相互间的差别超出随机波动范围时,则称这两个因素间存在交互作用。若某因素不同水平间的单独效应差因另一因素水平的影响呈较大幅度增加,并且差别有统计学意义,可认为两因素有协同交互作用;若某因素不同水平间的单独效应差因另一因素水平的影响呈较大幅度下降,并且差别有统计学意义,可认为两因素有拮抗交互作用。交互作用统计分析时,若存在交互作用,须逐一分析各因素的单独效应。反之如果不存在交互作用,则两因素的作用相互独立,分析某一因素的作用只需考察该因素的主效应。两因素的交互作用称为一阶交互作用,当因素个数大于2时,可计算二阶交互作用、三阶交互作用,……。A1A210304050B1B2B1B2两因素无交互作用反应值20A1A210204050B1B2B2B1两因素有交互作用,为拮抗作用反应值30析因试验的均数差别b因素平均b2-b1b1b2a124443420a228524024平均264822a2-a1486a因素本例考察的交互作用为不同的缝合方法是否影响两个时间点家兔的轴突通过率。ab交互作用(24-20)=4?当某因素的各个单独效应随另一因素水平的变化而变化,且相互间的差别超出随机波动范围时,则称这两个因素间存在交互作用。交互作用22122111()()2(84)22ABabababab缝合后2月后束膜缝合与外膜缝合神经轴突通过率的差异,仅比缝合后1月提高了2%,两条直线相互平行,表示两因素交互作用很小ANOVA分析的必要性A因素(缝合方法)的主效应为6%,B因素(缝合时间)的主效应为22%,AB的交互作用为2%。以上都是样本均数的比较结果,要推论总体均数是否有同样的特征,需要对试验结果进行方差分析后下结论。•H0:两种缝合方式间轴突通过率相同•H1:两种缝合方式间轴突通过率不同H0:不同时间轴突通过率相同H1:不同时间轴突通过率不同H0:缝合方式与时间存在交互作用H1:缝合方式与时间不存在交互作用α=0.052×2析因设计方差分析时变异分解及计算•变异分解•A因素(不同方法+误差)•B因素(不同时间+误差)•AB因素(AB交互+误差)•误差变异(随机因素+未知因素)处理组变异EAMSMSF1EBMSMSF2EABMSMSF3总变异()()ABABEABABESSSSSSSSSSSSSS处理总误差处理总误差SS处理的析因分解2.处理间离均差平方和3.A因素离均差平方和1.总离均差平方和SS估计量的计算方法:742020/7403480022CXSST2620526014022012022222CCnTSSiiB18052)260140()220120()(21222221CCAArSSNXC2)(A480026207420BTESSSSSS4.B因素离均差平方和242052)260220()140120()(21)(222221CCBBrBSS5.AB交互作用2024201802620)()()(BSSASSSSABSS处理6.误差离均差平方和析因分析结果两因素方差分析的分析策略小结:1先做两因素方差分析确定是否有交互作用;2如果没有交互作用,看主效应的差别是否有统计学意义:若有统计学意义,考察相应的样本均数,确定哪种情况的均数高。3如果有交互作用,则不能分析主效应。而化为单因素的方差分析(组数为各个因素的水平数之和),作两两比较。4在有交互作用的情况下,通过计算样本均数确认交互作用为协同作用还是拮抗作用。如果有交互作用,则:两个药都用的均数A药的均数+B药的均数-两个药都未用的均数,则称协同作用。两个药都用的均数A药的均数+B药的均数-两个药都未用的均数,则称拮抗作用。A、B因素(二)两因素多水平完全随机析因设计的方差分析例11-2观察A,B两种镇痛药物联合运用在产妇分娩时的镇痛效果。A药取3个剂量:1.0m,2.5mg,5.0mg;B药也取3个剂量:5μg,15μg,30μg。共9个处理组。将27名产妇随机等分为9组,每组3名产妇,记录每名产妇分娩时的镇痛时间。分析A,B两药联合运用的镇痛效果。两种药物联合镇痛效果研究例11-2:观察A、B两药联合应用在产妇分娩时的镇痛时间(min)完全随机的两因素3×3析因设计完全随机的两因素析因设计方差分析表A、B两药联合运用的镇痛时间方差分析结果表MultipleCompar