...Repeated Measures ANOVA)——SPSS软件实现
发布网友
发布时间：2024-08-19 04:00
共1个回答
热心网友
时间：2024-08-24 02:02
一、案例介绍
检验科研究血样放置时间对某生化指标浓度检测的影响，采集了10份人体血液标本，分别在放置0分钟(T0)、30分钟(T30)、60分钟(T60)和90分钟(T90)时对该指标的浓度(mmol/L)进行检测，分析放置时间是否对该生化指标检测结果有影响？
创建记录血标本编号的变量“ID”，测量尺度均设为“名义”。创建记录不同时点浓度的变量“T0”、“T30”、“T60”和“T90”，测量尺度设为“标度”。部分数据见图1。本文案例可从“附件下载”处下载。
二、问题分析
本案例的分析目的是比较4个时间点的生化指标浓度是否有差异。由于4个时间点的数据属于重复测量数据，可以使用单因素重复测量方差分析(One-Way Repeated Measures ANOVA)。但需要满足5个条件：
条件1：观察变量唯一，且为连续变量。本研究中观察变量为生化指标浓度，且为连续变量，该条件满足。
条件2：观察变量为重复测量数据，即不满足独立性。本研究中4个时间点测量的生化指标浓度均是针对同一批样本，因此不满足独立性，该条件满足。
条件3：观察变量不存在显著的异常值，该条件需要通过软件分析后判断。
条件4：各水平(时间点)观察变量为正态(或近似正态)分布，该条件需要通过软件分析后判断。
条件5：各重复测量水平观察变量的协方差相等，即满足球形假设，该条件需要通过软件分析后判断。
三、软件操作及结果解读
(一) 适用条件判断
1. 条件3判断(异常值判断)
① 选择“分析”—“描述统计”—“探索”(图2)。
② 在“探索 ”对话框中将表示生化指标浓度的变量“T0”、“T30”、“T60”和“T90”选入右侧“因变量列表”框，点击“确定”(图3)。
2. 条件4判断(正态性检验)
① 选择“分析”—“描述统计”—“探索”(图2)。
② 在“探索”对话框中将表示4个时间结果的变量“T0”、“T30”、“T60”和“T90”选入右侧“因变量列表”框(图3)。
③ 在“图”子对话框中勾选“含检验的正态图”，取消勾选“茎叶图”，其他不变，点击“继续”(图5)后回到“探索”对话框，点击“确定”，则会输出结果。
3. 条件5判断(球形假设检验)
球形假设检验的操作和统计描述及推断在同一对话框中完成，详见下述统计描述及推断部分。
(二) 统计描述及推断
1. 软件操作
①选择“分析”—“一般线性模型”—“重复测量” (图7)。
②出现“重复测量定义因子”对话框，“主体内因子名”处要录入重复测量的因子名称，自定义名称。由于本案例是在4个时间点重复测量，所以此处将“因子”改为“时间”，下方“水平数”是填入测量次数，则填入“4”，然后点击“添加”。最后点击“定义”(图8)。
③出现“重复测量”对话框，将左侧四个变量“T0”、“T30”、“T60”和“T90”选入右侧，则如图9所示。
④在“重复测量”对话框中点击右侧“图”，出现“轮廓图”子对话框，将“时间”选入右侧“水平轴”中，然后点击“添加”，保持选中“折线图”，如图10所示，点击“继续”回到主对话框。
⑤在“重复测量”对话框中点击右侧“EM平均值”，出现“重复测量：估算边际均值”子对话框，将“时间”选入右侧“显示下列各项的均值”列表框中(图11)，勾选下方“比较主效应”，“置信区间调整”下选择“邦弗罗尼”，点击“继续”回到主对话框。
⑥在“重复测量”对话框中点击右侧“选项”，出现“重复测量：选项”子对话框，勾选“描述统计”和“估计效应值”(图12)，点击“继续”回到主对话框后点击“确定”，则得到统计结果。
2. 结果解读
(1) 统计描述
图13列出了T0、T30、T60、T90四个时间点的浓度均值分别为5.66±0.31、5.55±0.34、5.27±0.23、4.87±0.23 mmol/L。
图14提供了各组的“估算边际均值”、“标准误差”及均值的“95% Confidence Interval (95%CI)”。因为不存在协变量，所以估算边际均值和统计描述均值一致。
图15为“MEASURE_1的估算边际均值”图，绘制了各时间点的生化指标浓度变化情况，可见浓度随着时间的增加而降低。
(2) 球形假设检验
图16“莫奇来球形度检验”中给出了球形假设检验结果，W=0.252，P=0.061，>0.05，满足球形假设。因此，本案例可以直接采用“单变量检验/主体内效应检验”非校正法的结果。当违背了球形假设条件时，可采用“多变量检验”结果，或者使用单变量检验校正法的结果，包括格林豪斯-盖斯勒、辛-费德特和下线三种校正方法。单变量检验校正法一般建议采用前两种方法，当epsilon (ε)<0.75时，使用格林豪斯-盖斯勒法，当epsilon (ε)>0.75时，使用辛-费德特法。当违背了球形假设条件时，如果单变量和多变量检验结果不一致，以多变量检验结果为准。
(3) 时间效应
图17“单变量检验/主体内效应检验”结果显示，第一行为“球形假定”结果，即当满足球形性假设时采用的结果，自由度为3，均方为3.709，F=116.103，P<0.001，提示不同时间点该生化指标浓度差异有统计学意义。下面三行提供了如果不满足球形性假设时，单变量检验三种校正法的结果。
图18为“多变量检验”分析结果。
(4) 时间变化趋势和两两比较
以上分析已显示不同时间的生化浓度存在差异，需要进一步探讨浓度变化符合何种变化曲线。由图19“主体内对比检验”可见，“线性”和“二次”关系的P值均<0.05，但“线性”关系的F值大于“二次”关系，且结合轮廓图(图15)可知，线性关系更符合变化趋势。
以上分析已显示不同时间点浓度存在差异，需要进一步两两比较。图20“成对比较”提供了各时间点两两比较的“平均值差值”、“标准误差”、“显著性(校正P值)”和均值的“95% Confidence Interval (95%CI)”。可知，随着时间的延长，各时刻与T0时刻相比，均值差逐渐增大，且均有统计学意义(P<0.05)，均值减少的趋势符合线性变化(F=187.769，P<0.001)。
四、结论
本研究采用单因素重复测量方差分析判断四个时间点的生化指标浓度是否有差异。通过箱式图判断，数据不存在异常值；通过Shapiro-Wilk检验，提示各组数据服从正态分布；球形度检验提示满足球形假设(W=0.252，P=0.061)，使用单变量检验的未校正法进行数据分析。
T0、T30、T60、T90时间点的生化指标浓度均值分别为5.66±0.31、5.55±0.34、5.27±0.23、4.87±0.23 mmol/L。分析提示，不同时间点生化指标浓度差异有统计学意义(F=116.103，P<0.001)。进一步采用“Bonferroni”校正法进行两两比较，可知随着时间的延长，各时刻与T0时刻相比，均值差逐渐增大，且均有统计学意义(P<0.05)，均值减少的趋势符合线性变化(F=187.769，P<0.001)。综上可知，放置时间对该生化指标检测结果具有较大的影响。
五、知识小贴士
(一) 球形度检验
(二) 基线数据处理
