《spss析因设计》由会员分享,可在线阅读,更多相关《spss析因设计(50页珍藏版)》请在金锄头文库上搜索。
1、spssspss析因设计析因设计方法:方法:方差分析方差分析(F检验)检验)目的:研究多个(包括两个)处理因素对目的:研究多个(包括两个)处理因素对试验对象的试验指标的作用试验对象的试验指标的作用资料:处理因素分几个水平(非定量),资料:处理因素分几个水平(非定量),试验指标定量试验指标定量结果:结果:应用:应用:FactorialdesignANOVA 析因设计的方差分析析因设计的方差分析一、析因设计二、析因设计的方差分析 (一)两因素两水平 (二)两因素三水平 (三)三因素多水平 实例实例实例实例1 1 1 1:甲乙两药治疗高胆固醇血症的疗效(胆固:甲乙两药治疗高胆固醇血症的疗效(胆固:甲
2、乙两药治疗高胆固醇血症的疗效(胆固:甲乙两药治疗高胆固醇血症的疗效(胆固醇降低值醇降低值醇降低值醇降低值mgmgmgmg),问),问),问),问甲乙两药是否有降低胆固甲乙两药是否有降低胆固甲乙两药是否有降低胆固甲乙两药是否有降低胆固醇的作用?醇的作用?醇的作用?醇的作用?两种药间有无交互作用两种药间有无交互作用两种药间有无交互作用两种药间有无交互作用完全随机的两因素完全随机的两因素2222析因设计析因设计甲药甲药 乙药乙药 用用 不用不用用用 64 56 78 44 80 42不用不用 28 16 31 25 23 18 实例实例实例实例2 2 2 2:小鼠种别:小鼠种别:小鼠种别:小鼠种别A
3、 A A A、体重、体重、体重、体重B B B B和性别和性别和性别和性别C C C C对皮内移植对皮内移植对皮内移植对皮内移植SRSSRSSRSSRS瘤瘤瘤瘤细胞生长特征影响的结果(肿瘤体积细胞生长特征影响的结果(肿瘤体积细胞生长特征影响的结果(肿瘤体积细胞生长特征影响的结果(肿瘤体积cmcmcmcm3 3 3 3)问)问)问)问AAAA、B B B B、C C C C各自的主效应如何?各自的主效应如何?各自的主效应如何?各自的主效应如何?三者间有无交互作用?三者间有无交互作用?三者间有无交互作用?三者间有无交互作用?完全随机的三因素完全随机的三因素222222析因设计析因设计实例实例实例实
4、例3 3 3 3:研究小鼠在不同注射剂量和不同注射频次下:研究小鼠在不同注射剂量和不同注射频次下:研究小鼠在不同注射剂量和不同注射频次下:研究小鼠在不同注射剂量和不同注射频次下药剂药剂药剂药剂ACTHACTHACTHACTH对尿总酸度的影响。问对尿总酸度的影响。问对尿总酸度的影响。问对尿总酸度的影响。问AAAA、B B B B各自的主效各自的主效各自的主效各自的主效应如何?应如何?应如何?应如何?二者间有无交互作用?二者间有无交互作用?二者间有无交互作用?二者间有无交互作用?随机配伍的两因素随机配伍的两因素3232析因设计析因设计显著特征显著特征每每个个处处理理是是各各因因素素各各水水平平的的
5、一一种种组组合合,总总处处理理数数为为各各因因素素各各水水平平的的全全面面组组合合数数,即即各各因因素素各各水水平平数数的的乘乘积积。如如两两因因素素析析因因设设计计,设设A因因素素有有I个个水水平平,B因因素素有有J个个水水平平,则则总总处处理理数数G=IJ。在在三三个个因因素素的的析析因因设设计计中中,若若各因素水平为各因素水平为I、J、K,则总处理数,则总处理数G=IJK。 要要求求各各个个处处理理组组内内的的实实验验单单位位数数相相等等(便便于于手手工工计计算算)且且每每组组至至少少有有两两个个实实验验单单位位,否否则则无无法法分分析析因因素素间间的交互作用,故总的实验单位数至少为的交
6、互作用,故总的实验单位数至少为2G。2.实验设计实验设计各因素各水平的全面组合各因素各水平的全面组合设有设有k个因素,每个因素有个因素,每个因素有L1, L2, , Lk个水平,个水平,那么共有那么共有g= L1L2Lk个处理组。个处理组。例如有三个因素,分别是例如有三个因素,分别是A,B,C。 A因素有因素有2水平,水平, B因素有因素有3水平,水平, C因素有因素有2水平,水平, 则处理组则处理组 g=232=12个个确定了处理组数后,将实验对象分配到各组的确定了处理组数后,将实验对象分配到各组的方法可以采用完全随机设计、随机区组设计或方法可以采用完全随机设计、随机区组设计或拉丁方设计拉丁
7、方设计 析因设计析因设计(Factorial Design)是一种多因是一种多因素多水平交叉分组进行全面试验的设计方素多水平交叉分组进行全面试验的设计方法。它可以研究两个或两个以上因素多个法。它可以研究两个或两个以上因素多个水平的效应。水平的效应。 在析因设计中,研究因素的所有可能的水平组合都能被研究到。 例如4个因素同时进行实验,每个因素取两个水平,实验的总组合数为24=16;如果水平为3,则有34=81种组合数。即是这81种组合均进行实验。 所以析因设计可以分析观测指标与研究因素间的复杂关系,包括各因素间的交互作用(Interaction) 。3.析因设计的特点析因设计的特点2个以上(处理
8、)因素(个以上(处理)因素(factor)(分类)(分类变量)变量)2个以上水平(个以上水平(level)2个以上重复(个以上重复(repeat)每次试验涉及全部因素每次试验涉及全部因素,即因素同时施加即因素同时施加观察指标(观测值)为观察指标(观测值)为计量资料计量资料(独立、(独立、正态、等方差)正态、等方差)4.析因设计的有关术语析因设计的有关术语单独效应(单独效应(simple effects):其它因素):其它因素(factor)的水平()的水平(level)固定为某一值时,某)固定为某一值时,某一因素的效应。即其他因素的水平固定时,同一一因素的效应。即其他因素的水平固定时,同一因素
9、不同水平间的差别。因素不同水平间的差别。主效应(主效应(main effects):某因素各单独效应的):某因素各单独效应的平均效应。某一因素各水平间的平均差别。平均效应。某一因素各水平间的平均差别。交互作用(交互作用(Interaction):某一因素的各个单独):某一因素的各个单独效应随着另一因素变化而变化的情况。(如一级效应随着另一因素变化而变化的情况。(如一级交互作用交互作用AB、二级交互作用、二级交互作用ABC)如果不存在交互效应,则只需考虑各因素的如果不存在交互效应,则只需考虑各因素的主效应。主效应。在方差分析中,如果存在交互效应,解释结在方差分析中,如果存在交互效应,解释结果时,
10、要逐一分析各因素的单独效应,找出果时,要逐一分析各因素的单独效应,找出最优搭配。最优搭配。在两因素析因设计时,只需考虑一阶交互效在两因素析因设计时,只需考虑一阶交互效应。三因素以上时,除一阶交互效应外,还应。三因素以上时,除一阶交互效应外,还需考虑二阶、三阶等高阶交互效应,解释将需考虑二阶、三阶等高阶交互效应,解释将更复杂。更复杂。(A+B)=A+B+AB(A+B+C)=A+B+C+AB+AC+BC+ABC5.析因设计的优缺点析因设计的优缺点优点:优点:比单处理因素设计能提供更多的试验信息比单处理因素设计能提供更多的试验信息,可可用来分析全部因素用来分析全部因素主效应主效应,以及因素间各级,以
11、及因素间各级的的交互作用,交互作用,在医学上可用于筛选最佳治疗方案、药物在医学上可用于筛选最佳治疗方案、药物配方、实验条件等研究。配方、实验条件等研究。缺点:当因素增加时,实验组数呈几何倍缺点:当因素增加时,实验组数呈几何倍数增加,数增加,所需试验的次数很多。所需试验的次数很多。不但计算复杂,不但计算复杂,不但计算复杂,不但计算复杂,而且给众多交互作用的解释带来困难。因此,当因素及水平而且给众多交互作用的解释带来困难。因此,当因素及水平而且给众多交互作用的解释带来困难。因此,当因素及水平而且给众多交互作用的解释带来困难。因此,当因素及水平数较多时,一般不采用完全交叉分组的析因设计,而采用正数较
12、多时,一般不采用完全交叉分组的析因设计,而采用正数较多时,一般不采用完全交叉分组的析因设计,而采用正数较多时,一般不采用完全交叉分组的析因设计,而采用正交设计交设计交设计交设计。(一)两因素两水平(一)两因素两水平(22)完全随机析因设计的方差分析完全随机析因设计的方差分析 A缝合方法缝合方法外膜缝合(外膜缝合(a1)束膜缝合(束膜缝合(a2)合计合计B缝合后时间缝合后时间1月月b12月月b21月月b12月月b210301050103020504070307050605060103030302444285212022014026074044001120048001440034800例例例例11
13、-111-111-111-1:研究不同缝合方法及缝合后时间对家兔:研究不同缝合方法及缝合后时间对家兔:研究不同缝合方法及缝合后时间对家兔:研究不同缝合方法及缝合后时间对家兔轴突通过率()的影响,问轴突通过率()的影响,问轴突通过率()的影响,问轴突通过率()的影响,问两种缝合方法间两种缝合方法间两种缝合方法间两种缝合方法间有无差别?缝合后时间长短间有无差别?有无差别?缝合后时间长短间有无差别?有无差别?缝合后时间长短间有无差别?有无差别?缝合后时间长短间有无差别?两者两者两者两者间有无交互作用间有无交互作用间有无交互作用间有无交互作用A A因素(因素(2 2水平)水平)缝合方法缝合方法 B B
14、因素(因素(2 2水平)水平) 缝合后时间缝合后时间缝合后缝合后1 1月月 b1 b1缝合后缝合后2 2月月 b2 b2外膜缝合外膜缝合a1a12424(a1b1a1b1)4444(a1b2a1b2)束膜缝合束膜缝合a2a22828(a2b1a2b1)5252(a2b2a2b2)2因素因素2水平析因实验示意图水平析因实验示意图表表表表 2 2 2 2因素因素因素因素2 2 2 2水平析因水平析因水平析因水平析因试验试验的均数差的均数差的均数差的均数差别别 B因素因素平均平均b2b1b1b2a124443420a228524024平均平均264822a2a1486A因素因素在在a1b1、a1b2
15、、a2b1和和a2b2的四种处理组合中,每个格子均有的四种处理组合中,每个格子均有5个数据,因此它又是重复数相等的析因设计。由于数据按因素个数据,因此它又是重复数相等的析因设计。由于数据按因素A和因素和因素B两个方向交叉分组,故可用双向方差分析。进一步分两个方向交叉分组,故可用双向方差分析。进一步分析处理的单独效应析处理的单独效应(simpleeffect)、主效应、主效应(maineffect)和交互效和交互效应应(interaction)。单独效应:单独效应:其他因素的水平固定时,同一因素不同水其他因素的水平固定时,同一因素不同水平间的差别平间的差别 A A因素因素缝合方法缝合方法 B B
16、因素因素 缝合后时间缝合后时间 平均平均 B B单独效应单独效应1 1月(月(b1b1)2 2月(月(b2b2) (b2b2b1b1)外膜(外膜(a1a1) 24 24 44 44 34 34 20 20束膜(束膜(a2a2) 28 28 52 52 40 40 24 24平均平均 26 26 48 48 22 22A A单独效应单独效应(a2a2a1a1) 4 4 8 8 6 6 表表2-2 222-2 22析因设计均数差别分析析因设计均数差别分析析因试验的均数差别析因试验的均数差别 b因素因素平均平均b2b1b1b2a124443420a228524024平均平均264822a2a1486
17、a因素因素a因素单独效应:比较缝合后两个时间点的轴突通过率因素单独效应:比较缝合后两个时间点的轴突通过率a a(b1(b1固定)固定)=(a2-a1)=28-24=4=(a2-a1)=28-24=4 a(b2固定固定)=(a2-a1)=52-44=8b因素单独效应:因素单独效应:b(a1固定)固定)=(b2-b1)=20b(a2固定固定)=(b2-b1)=24 主效应:主效应:某一因素各水平间的平均差别某一因素各水平间的平均差别A A因素的主效应解释为:束膜缝合与外膜缝合相比(不考虑缝合时间)因素的主效应解释为:束膜缝合与外膜缝合相比(不考虑缝合时间),神经轴突通过率提高了,神经轴突通过率提高
18、了6%6%(40-3440-34)。)。B B因素的主效应解释为:缝合后因素的主效应解释为:缝合后2 2月与月与1 1月相比(不考虑缝合方法),月相比(不考虑缝合方法),神经轴突通过率提高了神经轴突通过率提高了22%22%(48-2648-26)。)。 A A因素因素缝合方法缝合方法 B B因素因素 缝合后时间缝合后时间 平均平均 B B主效应主效应1 1月(月(b1b1)2 2月(月(b2b2)外膜(外膜(a1a1) 24 24 44 44 34 34 束膜(束膜(a2a2) 28 28 52 52 40 40平均平均 26 26 48 48 2222A A主效应主效应 6 62 2(交互作
19、用)(交互作用)交互作用交互作用(interaction) :当某因素的各个单独效应随另一因素水平的变当某因素的各个单独效应随另一因素水平的变化而变化化而变化,且相互间的差别超出随机波动范围时且相互间的差别超出随机波动范围时,则称则称这两个因素间存在交互作用。这两个因素间存在交互作用。若某因素不同水平间的单独效应差因另一因素水若某因素不同水平间的单独效应差因另一因素水平的影响呈较大幅度增加平的影响呈较大幅度增加,并且差别有统计学意义并且差别有统计学意义,可可认为两因素有协同交互作用认为两因素有协同交互作用;若某因素不同水平间的单独效应差因另一因素水若某因素不同水平间的单独效应差因另一因素水平的
20、影响呈较大幅度下降平的影响呈较大幅度下降,并且差别有统计学意义并且差别有统计学意义,可可认为两因素有拮抗交互作用。认为两因素有拮抗交互作用。交互作用交互作用 统计分析时,若存在交互作用,须逐一分析各因统计分析时,若存在交互作用,须逐一分析各因素的单独效应。反之如果不存在交互作用,则两因素的素的单独效应。反之如果不存在交互作用,则两因素的作用相互独立,分析某一因素的作用只需考察该因素的作用相互独立,分析某一因素的作用只需考察该因素的主效应。两因素的交互作用称为一阶交互作用,当因素主效应。两因素的交互作用称为一阶交互作用,当因素个数大于个数大于2 2时,可计算二阶交互作用、三阶交互作用,时,可计算
21、二阶交互作用、三阶交互作用,。A1A210304050B1B2B1B2两因素无交互作用两因素无交互作用反反应应值值20A1A210204050B1B2B2B1两因素有交互作用,为拮抗作用两因素有交互作用,为拮抗作用反反应应值值30析因析因试验的均数差的均数差别 b因素因素平均平均b2b1b1b2a124443420a228524024平均平均264822a2a1486a因素因素本例考察的交互作用为本例考察的交互作用为不同的缝合方法是不同的缝合方法是否影响两个时间点家兔的轴突通过率。否影响两个时间点家兔的轴突通过率。ab交互作用交互作用(2420)4? 当某因素的各个单独当某因素的各个单独效应随
22、另一因素水平的变效应随另一因素水平的变化而变化化而变化, ,且相互间的差且相互间的差别超出随机波动范围时别超出随机波动范围时, ,则称这两个因素间存在交则称这两个因素间存在交互作用。互作用。交互作用交互作用 缝合后2月后束膜缝合与外膜缝合神经轴突通过率的差异,仅比缝合后1月提高了2%, 两条直线相互平行, 表示两因素交互作用很小 ANOVA分析的必要性分析的必要性 A因素(缝合方法)的主效应为6%,B因素(缝合时间)的主效应为22%,AB的交互作用为2%。 以上都是样本均数的比较结果,要推论总体均数是否有同样的特征,需要对试验结果进行方差分析后下结论。 H0:两种缝合方式间轴突通过率相同:两种
23、缝合方式间轴突通过率相同H1:两种缝合方式间轴突通过率不同:两种缝合方式间轴突通过率不同H0:不同时间轴突通过率相同:不同时间轴突通过率相同H1:不同时间轴突通过率不同:不同时间轴突通过率不同H0:缝合方式与时间存在交互作用:缝合方式与时间存在交互作用H1:缝合方式与时间不存在交互作用:缝合方式与时间不存在交互作用0.052222析因设计方差分析时变异分解及计算析因设计方差分析时变异分解及计算 变异分解变异分解A因素(不同方法因素(不同方法+误差)误差)B因素(不同时间因素(不同时间+误差)误差)AB因素(因素(AB交互交互+误差)误差)误差变异(随机因素误差变异(随机因素+未知因素)未知因素
24、)处理组变异处理组变异总总变变异异SS处理处理的析因分解的析因分解 2.2.2.2.处理间离均差平方和处理间离均差平方和处理间离均差平方和处理间离均差平方和3.A3.A3.A3.A因素离均差平方和因素离均差平方和因素离均差平方和因素离均差平方和1.1.1.1.总离均差平方和总离均差平方和总离均差平方和总离均差平方和SS估计量的计算方法:估计量的计算方法:A A4.B4.B4.B4.B因素离均差平方和因素离均差平方和因素离均差平方和因素离均差平方和5.AB5.AB5.AB5.AB交互作用交互作用交互作用交互作用6.6.6.6.误差离均差平方和误差离均差平方和误差离均差平方和误差离均差平方和析因分
25、析结果析因分析结果 两因素方差分析的分析策略小结:两因素方差分析的分析策略小结:1 1 先做两因素方差分析确定是否有交互作用;先做两因素方差分析确定是否有交互作用;2 2 如果没有交互作用,看主效应的差别是否有如果没有交互作用,看主效应的差别是否有统计学意义:若有统计学意义,考察相应统计学意义:若有统计学意义,考察相应的样本均数,确定哪种情况的均数高。的样本均数,确定哪种情况的均数高。3 3 如果有交互作用,则不能分析主效应。而化如果有交互作用,则不能分析主效应。而化为单因素的方差分析为单因素的方差分析( (组数为各个因素的水组数为各个因素的水平数之和平数之和) ),作两两比较。,作两两比较。
26、4 4 在有交互作用的情况下,通过计算样本均数在有交互作用的情况下,通过计算样本均数确认交互作用为协同作用还是拮抗作用。确认交互作用为协同作用还是拮抗作用。如果有交互作用,则:如果有交互作用,则:两个药都用的均数两个药都用的均数AA药的均数药的均数B B药的均数两个药的均数两个药都未用的均数,则称协同作用。药都未用的均数,则称协同作用。两个药都用的均数两个药都用的均数AA药的均数药的均数B B药的均数两个药的均数两个药都未用的均数,则称拮抗作用。药都未用的均数,则称拮抗作用。A A、B B因素因素(二)两因素多水平完全随机析因设计的方差分析例例11-2观察观察A,B两种镇痛药物联合运用在两种镇
27、痛药物联合运用在产妇分娩时的镇痛产妇分娩时的镇痛效果。效果。A药取药取3个剂量:个剂量:1.0m,2.5mg,5.0mg;B药也取药也取3个剂量:个剂量:5g,15g,30g。共。共9个处理组。将个处理组。将27名产妇随机等分为名产妇随机等分为9组,组,每组每组3名产妇,记录每名产妇分娩时的镇痛时间。名产妇,记录每名产妇分娩时的镇痛时间。分析分析A,B两药联合运用的镇痛效果。两药联合运用的镇痛效果。两种药物联合镇痛效果研究两种药物联合镇痛效果研究 例例11-2:观察:观察A、B两药联合应用在产妇分娩时的镇痛时两药联合应用在产妇分娩时的镇痛时间(间(min)完全随机的两因素完全随机的两因素3333析因设计析因设计完全随机的两因素析因设计完全随机的两因素析因设计方差分析表方差分析表A、B两药联合运用的镇痛时间两药联合运用的镇痛时间方差分析结果表方差分析结果表DRUGADRUGB