12_重复测量方差分析报告模板

108 109 112 112 121 120 121 129 115 114 119 128 123 121 124
112 115 119 119 127 118 119 126 111 116 118 121 120 116 118
120 126 124 126 133 131 129 135 123 123 135 148 143 145 142
随机区组设计
处理只能在区组内随 机分配 区组内试验单位彼此 独立
第二节 重复测量数据的两因素两水平分析
高血压患者治疗前后的舒张压（mmHg）
处理组 编号 1 2 前 130 124 后 114 110 编号 11 12 对照组 前 118 132 后 124 122
3
4 5 6
136
128 122 118
5.94 5.49
5.26
5.88 5.43
4.93
5.43 5.32
4.70
5.04 5.04
5
6 7 8
5.71
6.27 5.88 5.32
5.49
6.27 5.77 5.15
5.43
5.66 5.43 5.04
4.93
5.26 4.93 4.48
越是相邻的点，相关系数越大 还是越小？
Correlations T0 T0 Pearson Correlation Sig. (2-tailed) N Pearson Correlation Sig. (2-tailed) N Pearson Correlation Sig. (2-tailed) N Pearson Correlation Sig. (2-tailed) N 1 . 8 .979** .000 8 .937** .001 8 .861** .006 8 T45 .979** .000 8 1 . 8 .882** .004 8 .875** .004 8 T90 .937** .001 8 .882** .004 8 1 . 8 .900** .002 8 T135 .861** .006 8 .875** .004 8 .900** .002 8 1 . 8
第十二章 重复测量设计资料的 方差分析 Repeated measures ANOVA
暨南大学医学院医学统计学教研室 林汉生
重复测量的概念
本章的重复测量实验：
对同一观察对象（病人，受试者、动物、
标本，等）的同一观察指标在不同时间点 上进行的多次测量。
同一受试者的多次测量之间可能存在某种
相关性，用通常的统计方法不能揭示内在 的特点，有时甚至会得出错误的结论。
第一节 重复测量资料的数据特征
前后测量设计 设立对照的前后测量设计 重复测量设计 设立对照的重复测量设计
一、前后测量设计
高血压患者治疗前后的舒张压（mmHg）
编号 1 2 治疗前 130 124 治疗后 114 110 差值 16 14
3
4 5
136
128 122
126
116 102
0.614
0.434
0.314
先将全部受试对象按某种或某些特征分为若干
个区组（block），使每个区组内的观察对象 随机地接受研究因素某一水平的处理。
重复测量数据与随机区组设计的区别
重复测量设计
处理（有对照的重复测 量）是在受试者间随机 分配 区组内试验单位彼此不 独立，同一受试者的测 量结果可能高度相关
T45
T90
T135
**. Correlation is significant at the 0.01 level (2-tailed).
越是相邻的时间点，数据之间的 相关性越大
各列数据处理之间存在不同程度的相关
性，而且越是相邻的时间点，数据之间 的相关性越大。 理论上讲，当重复测量数据之间相关， 用随机区组设计方差分析考察各时间点 之间的差别，会增加I类错误的概率。 如果只重复2次，不涉及多重比较，故不 涉及该问题。
10
12 20
6
7 8 9 10
118Baidu Nhomakorabea
116 138 126 124
100
98 122 108 106
18
18 16 18 18
前后测量设计的统计方法
前提条件：在假定测量时间对观测结果没有
影响时，才能推断处理是否有效。 统计方法：用配对 t 检验或配对秩和检验 高血压患者治疗前后的舒张压平均下降了 16mmHg。 经配对t检验，t =16.18， P＝0.000 能说明治疗有效吗？ 住院休息，环境和情绪的改变？考虑了吗？
表 12-16 不同麻醉诱导时相患者的收缩压（mmHg） 方法 A A A A A B B B B B C C C C C 序号
T0
T1
麻
醉
T2
诱
T3
导
时
相
T4
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
120 118 119 121 127 121 122 128 117 118 131 129 123 123 125
“球对称”检验结果：因只重复两次测量，
不能进行“球对称”检验。
“球对称”假设
球对称假设：各时间点数据的两两相关系
数为常数。 检验。
两个时间点的重复测量，不进行“球对称”
Tests of W ithin-Subjects Effects Measure: MEASURE_1 Source TIME Sphericity Assumed Greenhouse-Geisser Huynh-Feldt Lower-bound TIME * 分组 Sphericity Assumed Greenhouse-Geisser Huynh-Feldt Lower-bound Error(TIME) Sphericity Assumed Greenhouse-Geisser Huynh-Feldt Lower-bound Type III Sum of Squares 1020.100 1020.100 1020.100 1020.100 348.100 348.100 348.100 348.100 333.800 333.800 333.800 333.800 df 1 1.000 1.000 1.000 1 1.000 1.000 1.000 18 18.000 18.000 18.000 Mean Square 1020.100 1020.100 1020.100 1020.100 348.100 348.100 348.100 348.100 18.544 18.544 18.544 18.544 F 55.008 55.008 55.008 55.008 18.771 18.771 18.771 18.771 Sig. .000 .000 .000 .000 .000 .000 .000 .000
治疗前
治疗后
SPSS结果解释
Multivariate Tests b Effect TIME Pillai's Trace Wilks' Lambda Hotelling's Trace Roy' s Largest Root TIME * 分组 Pillai's Trace Wilks' Lambda Hotelling's Trace Roy' s Largest Root a. Exact statistic b. Design: Inte r cept+ 分组 Within Subjects Design: TIME Value .753 .247 3.056 3.056 .510 .490 1.043 1.043 F 55.008a 55.008a 55.008a 55.008a 18.771a 18.771a 18.771a 18.771a Hypothesis df 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 Error df 18.000 18.000 18.000 18.000 18.000 18.000 18.000 18.000 Sig. .000 .000 .000 .000 .000 .000 .000 .000
三、重复测量设计
当前后测量设计的重复测量次数m3时，
称重复测量设计或重复测量数据。 重复测量设计
设立对照的重复测量设计 未设立对照的重复测量设计

1. 设立对照的重复测量设计
将手术要求基本相同的15名患者随机分3
组，在手术过程中分别采用A，B，C三 种麻醉诱导方法，在T0（诱导前）、T1、 T2、T3，T4 五个时相测量患者的收缩压， 数据记录见表。
Tests the null hypothesis that the error covariance matrix of the orthonormalized transformed dependent variables is proportional to an identity matrix. a. May be used to adjust the degrees of freedom for the averaged tests of significance. Corrected tests are displayed in the Tests of Within-Subjects Effects table. b. Design: Inte r cept+分组 Within Subjects Design: TIME
将“治疗前”,“治疗后”选入… 将“分组”选入…
分别 单击Plots 、Options
将分组和time选入，单击Add， 单击Continue
选择“Descriptive statistics” 单击Continue
单击OK
结果解释
Descriptive Statistics 分组 处理组 对照组 Total 处理组 对照组 Total Mean 126.20 124.80 125.50 110.20 120.60 115.40 Std. Deviation 7.084 7.899 7.338 9.307 9.755 10.704 N 10 10 20 10 10 20
117 123 118 120 126 137 133 142 131 133 129 132 136 126 130
2. 未设立对照的重复测量数据
表12-3 受试者血糖浓度（mmol/L） 受试者 编号 1 放置时间（分） 0 5.32 45 5.32 90 4.98 135 4.65
2
3 4
5.32
二、设立对照的前后测量设计
高血压患者治疗前后的舒张压（mmHg）
处理组
编号 1 2 3 4 5 前 130 124 136 128 122 后 114 110 126 116 102 差值 编号 11 12 13 14 15
对照组
前 118 132 134 114 118 后 124 122 132 96 124 差值
6
7 8
118
116 138
100
98 122
16
17 18
128
118 132
118
116 122
9
10
126
124
108
106
19
20
120
134
124
128
确定疗效的前后测量设计必须增 加平行对照
平行对照：处理组和对照组的受试对象
来自同一人群，试验同时进行。 对两组治疗前后的血压差值进行两样本t 检验，信息损失较大。如编号1（130114＝16），编号8（138-122＝16）。
可以分析：
时间主效应：治疗前、后的舒张压均数是
否不同？ 组间主效应：处理组与对照组的舒张压均 数是否不同？ 交互作用：组间的舒张压均数差异是否随 时间改变（或各组舒张压的下降幅度是否 相同）？
AnalyzeGeneral Linear Model Repeated Measures
126
116 102 100
13
14 15 16
134
114 118 128
132
96 124 118
7
8 9 10 均数
116
138 126 124 126.2
98
122 108 106 110.2
17
18 19 20
118
132 120 134 124.8
116
122 124 128 120.6
这是多变量方差分析结果，不是本章讨论内容
Mauchly 's Test of Sphericity Measure: MEASURE_1
b
Epsilona Within Subjects Effect TIME Mauchly's W 1.000 Approx. Chi-Square .000 df 0 Sig. . Greenhous e-Geisser 1.000 Huynh-Feldt 1.000 Lower-bound 1.000
随机区组设计数据
表4-9 区组 1 2 3 4 5 均数 不同药物作用后小白鼠肉瘤重量（g） A药 B药 C药 合计 0.82 0.65 0.51 1.98 0.73 0.54 0.23 1.50 0.43 0.34 0.28 1.05 0.41 0.68 0.21 0.43 0.31 0.24 0.93 1.35