在前面文章中研究的因素可能是多个，但是重复测量因素只有一个，本文介绍存在两个重复测量因素(Two Repeated Measures Factor)但无交互作用的情况及其在SPSS软件中的实现过程。
关键词：SPSS; 重复测量; 重复测量资料; 重复测量方差分析; 两重复测量因素; 球形检验; 交互作用; 主效应; 单独效应
一、案例介绍
评价在某麻醉药物的作用下，不同部位某蛋白酶活性变化情况。在使用该种药物后，检测12只小白鼠用药后10min、20min、30min、40min 4个不同时点A部位和B部位蛋白酶的活性。数据见图1。本案例数据可从“附件下载”处下载。
二、问题分析
本案例资料具有两个重复测量因素，一个是时间因素，有4个水平，即time1~time4，分别代表用药后10min、20min、30min、40min 4个不同时点。另一个是身体不同部位，有2个部位，即A部位和B部位。因此需要使用两重复测量因素方差分析进行数据分析。需要满足以下5个条件：
条件1：观察变量唯一，且为连续变量。本研究中观察变量只有蛋白酶活性，且为连续变量，该条件满足。
条件2：有两个分析因素。本研究有时间因素（time1~time4）和部位因素（A和B）两个因素，该条件满足。
条件3：观察变量为重复测量数据，即不满足独立性。本研究中在4个时点测量的蛋白酶活性均是针对同一批样本，因此不满足独立性；相同个体不同部位也不满足独立性；该条件满足。
条件4：观察变量不存在显著的异常值，该条件需要通过软件分析后判断。
条件5：各组、各水平(时点)观察变量为正态(或近似正态)分布，该条件需要通过软件分析后判断。
关于方差齐性的问题，有无本案例不存在分组因素，因此无须检验方差齐同。
三、软件操作及结果解读
(一) 适用条件判断1
1. 条件4判断(异常值判断)
(1) 软件操作
① 选择“分析”—“描述统计”—“探索”(图2)。
② 在“探索 ”对话框中将8列观察变量全部选入右侧“因变量列表”框，点击“确定”(图3)。
(2) 结果解读
异常值通过箱线图和专业知识进行判断，如图4所示，为A部位time1的箱线图，提示不存在异常值。同样的方法请读者自行操作，依次可判定A、B部位其他时点数据也不存在异常值。
2. 条件5判断(正态性检验)
(1) 软件操作
① 选择“分析”—“描述统计”—“探索”(图2)。
② 在“探索”对话框中将观察变量“time1”“time2”和“time3”选入右侧“因变量列表”框 (图3)。
(步骤①②同条件4判断的软件操作，故条件4和条件5判断可同时在“探索”对话框中进行)。
③ 在“图”子对话框中勾选“含检验的正态图 ”，取消勾选“茎叶图”，其他不变，点击“继续”(图5)，回到“探索”对话框后点击“确定”，则会输出结果。
(2) 结果解读
图6显示了两种正态性检验的结果，Kolmogorov-Smirnov (柯尔莫哥洛夫-斯米诺夫，K-S)检验和Shapiro-Wilk (夏皮罗-威尔克，S-W)检验。K-S检验适用于大样本资料，本案例查看S-W检验结果，可见大多数数据的P值均>0.5，提示所有时间点下两部位的数据服从正态分布。关于正态性检验的注意事项详见推文(医学统计学核心概念及重要假设检验的软件实现(2/4)——正态性假设检验的SPSS实现)。
(二) 软件实现步骤
球形假设检验以及两重复测量因素方差分析时都需要基于重复测量窗口，以下首先介绍调用重复测量窗口的操作步骤。
1. 调用“重复测量窗口”
① 选择“分析”—“一般线性模型”—“重复测量” (图7)。
② 在图8中，“主体内因子名”输入“时间”，“级别数”输入“4”，点击“添加”；再次在“主体内因子名”输入“部位”，“级别数”输入“2”，点击“添加”，出现图9。点击图9中“定义”，出现图10。
③ 出现“重复测量”对话框，将左侧8个变量“A_1”—“B_4”选入右侧“主体内变量”中(图10)。此后的软件操作都基于该窗口。
2. “重复测量窗口”设置
① 在“重复测量”对话框（图10）中点击右侧“模型”，保持默认，即全因子模型，因为需要分析交互作用(图11)。
② 在“重复测量”对话框（图10）中点击右侧“图”，出现“轮廓图”子对话框，将“时间”选入右侧“水平轴”中，将“部位”选入右侧“单独的线条”中，然后点击“添加”，保持选中“折线图”，如图12所示，点击“继续”回到主对话框。
③ 在“重复测量”对话框（图10）中点击右侧“EM 平均值”，出现“重复测量：估算边际均值”子对话框，将“时间”、“部位”和“时间*部位”选入右侧“显示下列各项的均值”列表框中，勾选“比较主效应”，下方“置信区间调整”选则“邦弗罗尼”，如图13所示，点击“继续”回到主对话框。
④ 在“重复测量”对话框（图10）中点击右侧“选项”，出现子对话框，勾选“描述统计”，如图14所示，点击“继续”回到主对话框，点击“确定”，则输出一系列结果。
(三) 球形假设检验
图15“莫奇来球形度检验”中给出了球形假设检验结果，可见“时间”、“时间*部位”的检验结果均<0.05，提示不满足球形假设。当违背了球形假设条件时，可采用“多变量检验”结果，或者使用单变量检验校正法的结果，主要包括格林豪斯-盖斯勒、辛-费德特和下线3种校正方法。单变量检验校正法一般建议采用前两种方法，当epsilon (ε)<0.75时，使用格林豪斯-盖斯勒法，当epsilon (ε)>0.75时，使用辛-费德特法。当违背了球形假设条件时，如果单变量和多变量检验结果不一致，以多变量检验结果为准。
(四) 统计描述与推断
1. 统计描述
图16列出了两部位4个时点的均值。
图17的估算边际均值图绘制了两部位4个时点蛋白酶活性的变化情况，可见两组数据均有下降，但下降的幅度较为相似。
图18提供了两部位4个时点的“估算边际平均值”、“标准误差”及均值的“95%置信区间(95%CI)”。此处均值和图16的统计描述结果一致。
2. 交互作用判断
由于本案例有两个因素(一个为时间因素，另一个为部位)，因此需要首先判断两个因素之间是否存在交互作用。如果交互作用有统计学意义，则需要分析单独效应。
图19为“主体内效应检验”，第一行为“假定球形度”结果，当满足球形性假设时查看该行统计推断结果。此处可见四种方法分析的时间与部位之间的交互作用结果均一致，F时间*部位=1.283，P>0.05，提示时间与部位之间交互作用无统计学意义。
如果使用“多变量检验”分析结果(图20)，与“主体内效应检验”一致，时间与部位之间的交互作用仍然无统计学意义。因此，本案例采用主效应分析结果即可。
(五) 时间效应分析
根据图19“主体内效应检验”结果可知F时间=369.928，P<0.001，提示不同时间点蛋白酶活性差异有统计学意义。如果使用图20“多变量检验”分析结果，也可得到相同结果。
已经发现不同时点的蛋白酶活性存在差异，则需要进一步探讨蛋白酶活性变化符合何种变化曲线。由“主体内对比检验”可知，“线性”关系的F时间=2.344，P<0.001(图21)，结合轮廓图(图17)也可知，线性关系符合变化趋势。
已经发现不同时点蛋白酶活性存在差异，则需要进一步两两比较。图22“成对比较”结果中提供了各时点两两比较的“平均值差值”“标准误差”“显著性(校正P值)”和“95%置信区间(95%CI)”。可知，随着时间的延长，各时点与time1相比，均值差逐渐增大。
(六) 部位效应分析
根据图19“主体内效应检验”结果可知F部位=0.091，P=0.768，提示不同部位间蛋白酶活性差异无统计学意义。如果使用图20“多变量检验”分析结果，也可得到相同结果。
四、结论
本研究采用两重复测量因素方差分析评价在麻醉药物的作用下，A、B机体部位某蛋白酶活性变化情况。数据不存在异常值，整体数据服从正态分布；时间与部位之间交互作用无统计学意义(F时间*部位=1.283，P>0.05)，故进行主效应分析。
时间效应分析显示，不同时点蛋白酶活性差异有统计学意义(F时间=369.928，P<0.001)，其变化符合线性趋势(F时间=2.344，P<0.001)。两两比较发现，随着时间的延长，各时点与time1相比，均值差逐渐增大。部位效应分析显示，不同部位间蛋白酶活性差异无统计学意义(F部位=0.091，P=0.768)。
综上，该麻醉药对A、B两部位的蛋白酶活性影响差异无统计学意义，随着时间的延长，蛋白酶活性逐渐降低。