多因素方差分析(1)

《多因素方差分析(1)》由会员分享，可在线阅读，更多相关《多因素方差分析(1)（45页珍藏版）》请在金锄头文库上搜索。

1、第二章第二章 多因素方差分析多因素方差分析第一节第一节 引言引言 第二节第二节 基本概念基本概念第三节第三节 多因素方差分析多因素方差分析第四节第四节 协方差分析第五节第五节 多因子方差应用实例与计多因子方差应用实例与计算机实现算机实现1概述：概述：单因素方差分析是检验多个样本均数间单因素方差分析是检验多个样本均数间差异有无统计学意义的统计学方法。差异有无统计学意义的统计学方法。在医学领域中，还经常碰到研究多个因在医学领域中，还经常碰到研究多个因素对某个观察指标的作用的问题素对某个观察指标的作用的问题 。多因素方差分析是分析两个及两个以上多因素方差分析是分析两个及两个以上因素对观察指标影响的统

2、计方法。因素对观察指标影响的统计方法。2第一节第一节 引言引言 方差分析中，影响观察指标的因素称为方差分析中，影响观察指标的因素称为因子因子factor);factor);因子所处的状态称为因子的因子所处的状态称为因子的一个水平一个水平level of factor);level of factor);各因子水平各因子水平的组合称为处理的组合称为处理treatment).treatment).3第二节第二节 根本概念根本概念一一 方差分析的基本思想方差分析的基本思想二二 主效应和单独效应主效应和单独效应三三 交互作用交互作用 4一、方差分析根本思想一、方差分析根本思想n 将全部观测值的总变异按

3、影响结果的诸因素分解为相应的假设干局部变异，构造出反映各局部变异作用的统计量，在此根底上，构建假设检验统计量，以实现对总体参数的推断。5主效应和单独效应主效应和单独效应主效应主效应(main effect)：某一因素各个水平间的平均差异：某一因素各个水平间的平均差异 单独效应单独效应(simple effect):指其它因素水平固定在一个水平时，某一因素不同水指其它因素水平固定在一个水平时，某一因素不同水平之间均数的差异。平之间均数的差异。6二、主效应二、主效应n n n 7二因子方差分析二因子方差分析n例：例：A A、B B两药治疗缺铁性贫血两药治疗缺铁性贫血2424例，试验结果如下：例，试

4、验结果如下：四种疗法治疗缺铁性贫血后红细胞增加数四种疗法治疗缺铁性贫血后红细胞增加数1012/L1012/Ln 疗法疗法 红细胞增加数红细胞增加数 总体均数记号总体均数记号 一般疗法一般疗法 0.8 0.9 0.7 0.8 0.9 0.7 一般疗法加一般疗法加A A药药 1.3 1.2 1.1 1.3 1.2 1.1 一般疗法加一般疗法加B B药药 0.9 1.1 1.0 0.9 1.1 1.0一般疗法加一般疗法加A A药加药加B B药药 2.1 2.2 2.0 2.1 2.2 2.0 8研究目的研究目的 本例研究目的之一为本例研究目的之一为A A药的使用是药的使用是否会引起病人的红细胞数变化

5、。否会引起病人的红细胞数变化。研究目的之二为研究目的之二为B B药的使用是否会引药的使用是否会引起病人的红细胞数的变化。起病人的红细胞数的变化。n 9研究目的研究目的 研究目的之三为A药与B药是否有交互作用。所谓有协同作用，是指同时用A、B两药起的作用大于单独用A药和B药的作用之和。所谓有拮抗作用，是指同时用A、B两药起的作用小于单独用A药各B药的作用之和。10研究目的研究目的 不管协同或拮抗作用均意味着A、B药同时使用的作用不等于单独作用之和。两药有无协同作用或拮抗作用，只要检验假设：11两因素有重复数据的方差分析变异分解两因素有重复数据的方差分析变异分解n 方差分析表方差分析表n n 变异

6、来源变异来源 SS df MS F P SS df MS F P处理间模型处理间模型 2.9625 3 0.9875 98.75 0.0000 2.9625 3 0.9875 98.75 0.0000 因子因子A 1.6875 1 1.6875 168.75 0.0000A 1.6875 1 1.6875 168.75 0.0000因子因子B 0.9075 1 0.9075 90.75 0.0000B 0.9075 1 0.9075 90.75 0.0000A A与与B B的交互作用的交互作用 0.3675 1 0.3675 36.75 0.0000 0.3675 1 0.3675 36.75

7、 0.0000 误差误差 0.0800 8 0.0100 0.0800 8 0.0100 总的总的 3.0425 11 0.2766 3.0425 11 0.276612n n n 13注意：注意：当因子当因子A与与B间的交互作用有统计学意义时，对间的交互作用有统计学意义时，对A或或B的单独作用的解释须小心。本例，用的单独作用的解释须小心。本例，用B药时，用药时，用A药病人药病人比不同时用比不同时用A药的病人的红细胞数均数大，不用药的病人的红细胞数均数大，不用B药时，用药时，用A药病人比不同时用药病人比不同时用A药的病人的红细胞数均数也大，故可说药的病人的红细胞数均数也大，故可说明明A药有效。

8、但有时可能出现这种情况，用药有效。但有时可能出现这种情况，用B药时，用药时，用A药病药病人比不同时用人比不同时用A药的病人的红细胞数均数大，不用药的病人的红细胞数均数大，不用B药时，药时，用用A药病人比不同时用药病人比不同时用A药的病人的红细胞数均数小，此时药的病人的红细胞数均数小，此时就不能简单地说就不能简单地说A药有利于病人红细胞数增加，需分别就用药有利于病人红细胞数增加，需分别就用B药和不用药和不用B药两种情况说明药两种情况说明A药的作用。对药的作用。对B作用的作用的作用的作用的解释也是如此。解释也是如此。14三因子方差分析三因子方差分析 例题例题 某研究者以大白鼠作试验，观某研究者以大

9、白鼠作试验，观察指标是肝重与体重之比察指标是肝重与体重之比5%5%，主要想，主要想了解正氟醚对观察指标的作用，同时要考了解正氟醚对观察指标的作用，同时要考察用生理盐水和用戊巴比妥作为诱导药对察用生理盐水和用戊巴比妥作为诱导药对正氟醚毒性作用有无影响，对不同性别大正氟醚毒性作用有无影响，对不同性别大白鼠诱导的作用有何不同，以及对不同性白鼠诱导的作用有何不同，以及对不同性别大白鼠正氟醚的作用是否相同。别大白鼠正氟醚的作用是否相同。15A因子因子 B因子因子 C因子因子 肝重与体重之比肝重与体重之比 总体均数总体均数不用不用 生理盐水生理盐水 雌雌 5.26 5.68 5.83 111不用不用 生理

10、盐水生理盐水 雄雄 5.00 5.52 5.38 112不用不用 戊巴比妥戊巴比妥 雌雌 5.87 5.50 6.20 121不用不用 戊巴比妥戊巴比妥 雄雄 6.13 6.46 5.21 122用用 生理盐水生理盐水 雌雌 5.42 5.60 5.70 211用用 生理盐水生理盐水 雄雄 6.30 7.02 5.90 212用用 戊巴比妥戊巴比妥 雌雌 4.64 4.60 5.44 221用用 戊巴比妥戊巴比妥 雄雄 6.02 5.70 5.48 222161718方差分析的随机效应模型方差分析的随机效应模型一一 随机效应模型随机效应模型 二二 固定效应模型固定效应模型 19一、固定效应和随

11、机效应模型的定义一、固定效应和随机效应模型的定义n 方差分析中的因子有选择型与随机型之分，假方差分析中的因子有选择型与随机型之分，假设数据资料中涉及到因子水平是研究者关心的因子设数据资料中涉及到因子水平是研究者关心的因子水平全体，那么该因子属于选择型因子；相应的模水平全体，那么该因子属于选择型因子；相应的模型称为固定效应模型。型称为固定效应模型。20一、固定效应和随机效应模型的定义一、固定效应和随机效应模型的定义 假设数据资料中涉及到因子水平只是研究者关心的因子水平总体的一个样本，那么该因子属于随机型因子；假设你的研究中有某些因子是随机型因子或全为随机型因子时，方差分析的模型与固定效应模型相同

12、，但关于主效应、和交互效应的假定及F统计量的计算公式有些不同。21n它们的计算公式为：它们的计算公式为：n 22 实例实例 某医院管理者欲了解血压计与量血压某医院管理者欲了解血压计与量血压的医生对血压测定结果是否有影响。他在医的医生对血压测定结果是否有影响。他在医院中随机抽取院中随机抽取3 3台血压计，台血压计，4 4名医生，对名医生，对2424名名体检者测量血压，下面是舒张压的观察结果，体检者测量血压，下面是舒张压的观察结果，请作分析请作分析。23 血压计血压计 医生甲（医生甲（b=1b=1）医生乙医生乙(b=2)(b=2)医生丙医生丙(b=3)(b=3)医生丁医生丁(b=4)(b=4)a=

13、1 60,97 84,63 70,99 74,68a=1 60,97 84,63 70,99 74,68a=2 91,60 85,88 90,74 76,62 a=2 91,60 85,88 90,74 76,62 a=3 85,67 90,71 65,79 75,96 a=3 85,67 90,71 65,79 75,96 24 n 用随机效应模型作为方差分析时，离均差平方和与自由度的计算与固定效应相同，但无效假设与F统计量的计算有所不同。25 n 2627方差分析的混合效应模型方差分析的混合效应模型n例题：设某人研究围产期窒息对新生儿中血中次黄例题：设某人研究围产期窒息对新生儿中血中次黄嘌

14、呤浓度是否有影响，同时还了解新生出生一小时嘌呤浓度是否有影响，同时还了解新生出生一小时内次黄嘌呤浓度是否有变化。他随机抽取围产期窒内次黄嘌呤浓度是否有变化。他随机抽取围产期窒息息9 9名，不窒息的正常新生儿名，不窒息的正常新生儿9 9名作为对照对每名作为对照对每组的组的9 9名新生儿随机安排三个不同时间，测定血中次名新生儿随机安排三个不同时间，测定血中次黄嘌呤浓度如下：黄嘌呤浓度如下：28 因子因子B因子因子A 出生时出生时 出生后出生后20分钟分钟 出生后出生后30分钟分钟对照组对照组 6.20 11.50 14.53 5.80 13.37 11.40 8.25 24.10 12.37围产期

15、围产期 23.06 25.56 10.52窒息组窒息组 21.46 30.40 13.66 11.43 18.19 18.20 29n用混合效应作方差分析时，离均差平方和与自由度的计算与固定效应相同，但无效假设与F统计量不同。它们的计算公式为：n 30几点说明1、每个处理均有观察且有相同观察例数n),这种设计称为完全、平衡设计。假设每个处理均有观察，但观察例数不等，那么属于完全、不平衡设计。2、方差分析要求观察值独立且服从正态分布，还要求各处理组有相同的方差。在这三个条件中，对独立性要求最严，对正态性要求最宽，当各处理组有相同观察例数时。对方差齐性的要求也不严。31n 本章仅介绍二因子、三因子

16、方差分析，二因子方本章仅介绍二因子、三因子方差分析，二因子方差模型中除了各因子的主效应外，还有两因子的一差模型中除了各因子的主效应外，还有两因子的一级交互作用项；三因子模型中除了主效应，每两因级交互作用项；三因子模型中除了主效应，每两因子的交互效应外，还有三因子的二级交互效应；四子的交互效应外，还有三因子的二级交互效应；四因子模型中除主效应，每两因子一级交互效应，三因子模型中除主效应，每两因子一级交互效应，三因子的二级交互效应；还有四因子的三级交互作用因子的二级交互效应；还有四因子的三级交互作用项；五因子及五个以上的因子的模型以此类推。项；五因子及五个以上的因子的模型以此类推。32第四节第四节 协方差分析协方差分析一一 概述概述 二二 协方差的基本相思协方差的基本相思 三三 应用条件应用条件 33概述概述n检验两个或两个以上均数间差异的显著性，可考虑用方差分检验两个或两个以上均数间差异的显著性，可考虑用方差分析。方差分析要求各比较组除了所施加的处理因素不同外，析。方差分析要求各比较组除了所施加的处理因素不同外，其它对观察指标有影响因素的因素齐同或均衡，即要求控制其它对观察指标有影响因素