实习五-ANOVA及双因素方差分析

1、实习五-ANOVA及双因素方差分析 实习五实习五 实习五-ANOVA及双因素方差分析 实习目的：实习目的： 实习五-ANOVA及双因素方差分析 方差分析（方差分析（ANOVA）用）用 于进行两组及多组间样本于进行两组及多组间样本 均数的比较。均数的比较。 实习五-ANOVA及双因素方差分析 总变异总变异= 随机变异随机变异 + 处理因素导致的变异处理因素导致的变异 总变异总变异 = 组内变异组内变异 + 组间变异组间变异 方差分析的基本思想是基于变异分方差分析的基本思想是基于变异分 解的原理：解的原理： 需要证明：组间变异远远大于组内变异，则说明处理因需要证明：组间变异远远大于组内变异，则说明

2、处理因 素的存在，如果两者相差无几，则说明该影响不存在。素的存在，如果两者相差无几，则说明该影响不存在。 实习五-ANOVA及双因素方差分析 单因素方差分析单因素方差分析 双因素方差分析双因素方差分析 实习五-ANOVA及双因素方差分析 适用条件：适用条件： 各样本是相互独立的随机样本；各样本是相互独立的随机样本； 各样本来自正态分布；各样本来自正态分布； 各样本方差相等（方差齐性检验）各样本方差相等（方差齐性检验） 。 实习五-ANOVA及双因素方差分析 统计量统计量F 的计算及其意义的计算及其意义 F=MS组间 组间/MS组内组内 自由度：自由度： 组间 组间=组数 组数-1 组内 组内=

3、N-组数 组数 通过这个公式计算出统计量通过这个公式计算出统计量F，求出对应的，求出对应的P值，值， 与与 进行比较，进行比较，做出推论做出推论。 如果如果p ，说明，说明得到大于现有统计量得到大于现有统计量F值的可能性值的可能性p 大于大于 ，不属于小概率事件，则不拒绝，不属于小概率事件，则不拒绝H0，差别无统计学，差别无统计学 意义，结论是各总体均数相等。意义，结论是各总体均数相等。 如果如果p ，说明得到等于或大于现有统计量的可能性，说明得到等于或大于现有统计量的可能性 p小于小于 ，可判断为小概率事件，则拒绝，可判断为小概率事件，则拒绝H0，接受，接受H1，差，差 别有统计学意义，结论

4、是各总体均数不相等。别有统计学意义，结论是各总体均数不相等。 实习五-ANOVA及双因素方差分析 SPSS操作步骤：操作步骤： 实习五-ANOVA及双因素方差分析 选入需要分析选入需要分析 的变量的变量 选入需要比较的分选入需要比较的分 组因素组因素(only one) u对均数的变动趋势进行检验；对均数的变动趋势进行检验； u定义根据研究目的需要进行的某些精确两两比较定义根据研究目的需要进行的某些精确两两比较 实习五-ANOVA及双因素方差分析 练习练习5.1因变量和因素变量分别选入：因变量和因素变量分别选入： 实习五-ANOVA及双因素方差分析 定义是否在方差分析中进行趋势检验定义是否在方

5、差分析中进行趋势检验 实习五-ANOVA及双因素方差分析 用于选择各组间进行两用于选择各组间进行两 两比较的方法两比较的方法 实习五-ANOVA及双因素方差分析 提供了当各组方差齐时可用的两两比较的方法：提供了当各组方差齐时可用的两两比较的方法： 运用最广泛运用最广泛最敏最敏 感感,每每 组之组之 间都间都 比较比较 实习五-ANOVA及双因素方差分析 实习五-ANOVA及双因素方差分析 输出常用统计描输出常用统计描 述指标述指标 方差齐性检验方差齐性检验 用各组均数做图，同时可辅助对均数间趋势做出判断。用各组均数做图，同时可辅助对均数间趋势做出判断。 实习五-ANOVA及双因素方差分析 结果

6、解释： 单因素方差分析：必须方差齐性。下述结果表明确实可单因素方差分析：必须方差齐性。下述结果表明确实可 以继续进行以继续进行ANOVE分析。分析。 方差齐性检验结果：方差齐性检验结果： 实习五-ANOVA及双因素方差分析 A AN NO OV VA A dna 47.785223.89311.648.000 43.074212.051 90.86023 Between Groups Within Groups Total Sum of SquaresdfMean SquareFSig. 单因素方差分析结果可见，单因素方差分析结果可见，F=11.648,P0.001。因此可认为。因此可认为 三

7、组鼠脾中三组鼠脾中DNA含量不全相同含量不全相同(差别有统计学意义，但具体差别有统计学意义，但具体 谁与谁不同，需要看下面的结果。谁与谁不同，需要看下面的结果。)。 实习五-ANOVA及双因素方差分析 Multiple ComparisonsMultiple Comparisons Dependent Variable: dna 2.2321*.7412.007.6913.774 3.3194*.6959.0001.8724.767 -2.2321*.7412.007-3.774-.691 1.0873.7218.147-.4142.588 -3.3194*.6959.000-4.767-1.

8、872 -1.0873.7218.147-2.588.414 (J) group 2 3 1 3 1 2 (I) group 1 2 3 LSD Mean Difference (I-J)Std. ErrorSig.Lower BoundUpper Bound 95% Confidence Interval The mean difference is significant at the .05 level. *. LSD法两两比较结果法两两比较结果 上表可见：显示了每两组之间的具体的差别。与后面的 SNK方法完全不同，例如在这里： 2与3之间比较是0.147，没有统计学差异。 实习五-AN

9、OVA及双因素方差分析 S-N-K法两两比较结果法两两比较结果 dnadna 911.656 712.743 814.975 .1461.000 group 3 2 1 Sig. Student-Newman-Keulsa,b N12 Subset for alpha = .05 Means for groups in homogeneous subsets are displayed. Uses Harmonic Mean Sample Size = 7.916. a. The group sizes are unequal. The harmonic mean of the group s

10、izes is used. Type I error levels are not guaranteed. b. SPSS在用在用S-N-K法进行两两比较时，如果有差异，则只法进行两两比较时，如果有差异，则只 会告诉你会告诉你P值小于预定的界值（默认为值小于预定的界值（默认为0.05），而不会给），而不会给 出具体的出具体的P值。值。这个跟前面的给出具体的组与组之这个跟前面的给出具体的组与组之 间差异值的间差异值的LSD方法不一样。方法不一样。 从表中可见，自发性（从表中可见，自发性（group2）和移植）和移植 性白血病鼠性白血病鼠（group3）被分在一组，被分在一组， 与正常鼠处于不同亚

11、组，因此，正常鼠和与正常鼠处于不同亚组，因此，正常鼠和 其他两种鼠脾其他两种鼠脾DNA含量间有差异。含量间有差异。 亚组内亚组内DNA含量差含量差 异无统计学意义异无统计学意义 实习五-ANOVA及双因素方差分析 l将结果整理成方差分析的步骤和方法将结果整理成方差分析的步骤和方法 方差分析：方差分析： 1.建立检验假设建立检验假设 2.确定检验水平确定检验水平 3.计算检验统计量计算检验统计量F值，并将方差分析的结果整理成表格形式（方差值，并将方差分析的结果整理成表格形式（方差 分析表）：分析表）： 变异来源（组间、组内、总）、离均差平方和、自由度、均方、变异来源（组间、组内、总）、离均差平方

12、和、自由度、均方、F值、值、P值值 4.确定确定P值值 5.做出统计推断做出统计推断 方差齐性检验方差齐性检验 各组均数间的两两比较各组均数间的两两比较 实习五-ANOVA及双因素方差分析 (2)方差分析：方差分析： 建立假设：建立假设： H0 ：？总体均数相等：？总体均数相等 1 = 2 = 3 H1 ：？总体均数：？总体均数不全不全相等相等 确定显著性水平，用确定显著性水平，用 表示表示，常取，常取0.05。 计算统计量计算统计量F：F=MS组间 组间/MS组内组内=？ ？ 计算概率值计算概率值P：p=？ 做出推论：做出推论： (1)方差齐性检验：方差齐性检验：F=?，p=? 与与 =0.

13、05的关系，各组方差齐？。的关系，各组方差齐？。 (3)各组均数间的两两比较：各组均数间的两两比较： 采用？两两比较的方法，？组与？组均数有差别采用？两两比较的方法，？组与？组均数有差别 (p0.05)；再写出相应的专业结论再写出相应的专业结论。 实习五-ANOVA及双因素方差分析 双因素方差分析双因素方差分析(two-way ANOVA) 配伍组设计配伍组设计(或双因素无重复试验设计或双因素无重复试验设计):是配对设计的扩是配对设计的扩 展展 在实验研究中将实验动物按窝别等特征配伍在实验研究中将实验动物按窝别等特征配伍,再随机分配到各再随机分配到各 处理组中处理组中; 对同一受试对象不同时间

14、点上（即多余两个时间点）的观察对同一受试对象不同时间点上（即多余两个时间点）的观察 研究研究; 给同一样本不同的处理给同一样本不同的处理; 各处理组样本数相同各处理组样本数相同如果是重复测量数据，就是说如果是重复测量数据，就是说 在不同时间点不能有缺省值。在不同时间点不能有缺省值。 无法考察交互作用和方差齐性，所以不用管方差齐不齐。无法考察交互作用和方差齐性，所以不用管方差齐不齐。 实习五-ANOVA及双因素方差分析 统计量统计量F 的计算及其意义的计算及其意义 F1=MS组间 组间/MS组内组内 F2=MS配伍间 配伍间/MS组内组内 都是和随机误差进行比较。都是和随机误差进行比较。 自由度

15、：自由度： 组间 组间 = 组数 组数-1 配伍间 配伍间 = 区组数 区组数-1 组内 组内 = (组数 组数-1)(区组数区组数-1) 实习五-ANOVA及双因素方差分析 练习练习5.2分析分析: 处理因素：即研究的主要因素，为不同处理因素：即研究的主要因素，为不同 时期（变量名：时期（变量名：time）；）； 配伍因素：不同患者治疗前、中、后期配伍因素：不同患者治疗前、中、后期 各时点血液中粘蛋白含量不同，即研究各时点血液中粘蛋白含量不同，即研究 对象的不同。对象的不同。 实习五-ANOVA及双因素方差分析 操作步骤：从操作步骤：从analyze-普通线性模型（普通线性模型（GLM）-进

16、入，见下图。进入，见下图。 实习五-ANOVA及双因素方差分析 因变量因变量 固定变量固定变量 自定义方差分自定义方差分 析模型析模型 协变量协变量:与自变量和因变量可能有关系的连续性变量与自变量和因变量可能有关系的连续性变量 随机因素变量随机因素变量 加权最小二乘法的权重系数加权最小二乘法的权重系数 实习五-ANOVA及双因素方差分析 将分析主效应将分析主效应 的变量选入的变量选入 分析所有分分析所有分 类变量的主类变量的主 效应和交互效应和交互 作用。而作用。而 Custom只是只是 比较主效应，比较主效应， 不看交互作不看交互作 用。用。 方差分析进行变方差分析进行变 异分解的方法异分解

17、的方法 实习五-ANOVA及双因素方差分析 用于对精细趋势用于对精细趋势 检验和精确两两检验和精确两两 比较的选项进行比较的选项进行 定义，比较专业，定义，比较专业， 少用。少用。 实习五-ANOVA及双因素方差分析 各组均数间的多重比较：各组均数间的多重比较： 实习五-ANOVA及双因素方差分析 选入需要估计均数的因素选入需要估计均数的因素 和交互作用和交互作用 进行均数间的多重比较进行均数间的多重比较 输出非常重要指标输出非常重要指标 实习五-ANOVA及双因素方差分析 第一行（校正的模型）：所用方差分析模型的检验第一行（校正的模型）：所用方差分析模型的检验上上 述结果显示有差异，方差不齐。述结果显示有差异，方差不齐。 第二行：截距，分析中无实际意义，可忽略。第二行：截距，分析中无实际意义，可忽略。 第三行：配伍因素第三行：配伍因素不同对象的粘蛋白也有差异不同对象的粘蛋白也有差异 第四行：变量（治疗不同时期粘蛋白含量有差别，具体是第四行：变量（治疗不同时期粘蛋白含量有差别，具体是 哪两个或几个间存在差别，需要后面两两比较来回答。）哪两个或几个间存在差别，需要后面两两比较来回答。） 第五行：误差；第六行：合计；第七行：校正的合计第五行：误差；第六行：合计；第七行：校正的合计 实习五-AN