重复测量方差分析结果解读一、方法概述重复测量方差分析Repeated Measures ANOVA是一种用于分析同一受试对象在多个时间点或条件下重复观测数据的统计方法。与独立样本方差分析不同重复测量设计考虑了同一受试对象多次测量之间的相关性将总变异分解为组间效应不同处理组之间的差异、组内效应不同时间点之间的差异以及交互效应组间因素与组内因素的联合作用三个部分。该方法需要满足球形度假定即各时间点之间差异的方差相等。当球形度假定不满足时需采用Greenhouse-GeisserGG或Huynh-FeldtHF校正方法调整自由度。本分析使用SPSSAU软件完成重复测量方差分析的全部计算与检验过程。在SPSSAU【实验/医学研究】模块选择【重复测量方差】将变量拖拽至右侧对应分析框操作如下图二、效应分析结果总览表1展示了重复测量方差分析的效应分析结果包括截距、组间因素group、组内因素time以及交互效应group:time四个效应的平方和、自由度、均方、F值、p值及效应量指标。从表1可知截距效应高度显著F8861.487p0.001偏eta^20.994说明因变量的总体均值显著不为零。组间因素group的主效应不显著F0.122p0.7290.05偏eta^20.002效应量极小表明两组之间的总体差异无统计学意义。组内因素time的主效应高度显著F163.752p0.001偏eta^20.759说明不同时间点之间因变量存在极其显著的差异且时间因素能够解释因变量75.9%的组内变异效应量属于大效应。交互效应time x group同样高度显著F14.879p0.001偏eta^20.222表明时间因素对因变量的影响模式在两组之间存在显著差异即两组随时间变化的趋势不同。三、球形度检验在进行组内效应分析之前需要检验球形度假定。球形度检验的原假设为各时间点间差异的方差矩阵满足球形性。表2报告了球形度检验结果。球形度检验结果显示time效应的Mauchly W值为0.869对应的p值为0.2130.05不能拒绝球形度假设表明数据满足球形度假定。group:time交互效应的球形度检验结果相同W0.869p0.213同样满足球形度假定。因此后续组内效应分析可直接使用满足球形检验的结果无需进行GG或HF校正。GG校正系数为0.915HF校正系数为0.971均接近1进一步佐证了球形度假定的合理性。四、组内效应分析表3展示了组内效应的详细分析结果包括满足球形检验、GG校正和HF校正三种情况下的F值、自由度和p值。由于球形度检验通过p0.2130.05应优先参考满足球形检验一栏的结果。time主效应在三种校正方法下均高度显著F163.752p0.001偏eta^20.759表明不同时间点之间因变量的差异极其显著时间因素解释了因变量75.9%的组内变异。group:time交互效应在三种校正方法下同样高度显著F14.879p0.001偏eta^20.222说明两组随时间变化的趋势存在显著差异交互效应解释了22.2%的组内变异。三种校正方法下的p值均为0.000结果高度一致进一步增强了结论的稳健性。五、事后多重比较由于time主效应显著且包含4个水平时间点需进行事后多重比较以明确各时间点之间的具体差异。表4报告了time因素6对两两比较的结果。事后多重比较结果显示time因素的4个时间点之间两两差异均达到极其显著的水平均p0.001。从均值差值来看时间1到时间4的均值呈现持续递增趋势时间1与时间2之间的均值差为-3.951时间2与时间3之间为-3.619时间3与时间4之间为-4.823。时间1与时间4之间的均值差最大-12.392t-20.415表明起始时间点与最终时间点之间的差异最为突出。所有比较的t值绝对值均大于5p值均小于0.001说明因变量在各时间点之间呈持续且显著的递增趋势每两个相邻时间点之间均存在统计学意义上的显著差异。group因素仅有2个水平事后比较结果显示两组之间的差异不显著均值差-0.573t-0.349p0.729与主效应分析中group效应不显著的结论一致。六、简单效应分析由于交互效应time x group显著需进一步进行简单效应分析以揭示交互作用的具体模式。表5报告了在time的每个水平上两组之间的差异比较。从表5可知在时间1时两组之间的差异具有统计学意义均值差3.889t2.165p0.0350.05group1的均值显著高于group2。而在时间2p0.789、时间3p0.138和时间4p0.104时两组之间的差异均不具有统计学意义。这表明两组在起始时间点存在基线差异但随着时间推移两组之间的差异逐渐缩小并趋于一致。表6报告了在group的每个水平上各时间点之间的两两差异比较。从表6可以看出两组在各时间点之间的变化模式存在明显差异。对于group1时间1与时间2之间的差异不显著p0.275但从时间2开始各时间点之间的差异均达到显著水平p0.01说明group1在时间2之后才开始出现显著变化。对于group2所有时间点之间的两两差异均极其显著均p0.01且均值差普遍大于group1如时间1 vs 时间4group2为-15.830group1为-8.954说明group2在各时间点上的变化幅度更大、变化更为剧烈。这一结果与交互效应显著的分析结论一致两组的时间变化轨迹存在本质差异。七、效应图分析图1和图2分别展示了time和group主效应的均值变化趋势。图1 time主效应均值图从图1可以看出因变量的均值随时间呈持续递增趋势从时间1的约71.4逐步上升至时间4的约83.8。四个时间点之间的递增趋势较为均匀相邻时间点之间的增幅约为3.5~4.8个单位与事后多重比较中各时间点两两差异均显著的结果相一致。图2 group主效应均值图从图2可以看出两组的均值非常接近group1约77.1group2约77.6两条线几乎重合直观反映了group主效应不显著p0.729的结果。两组之间的差异仅为0.573个单位在实际意义上几乎可以忽略。图3 time x group交互效应图图3展示了两组在各时间点上的均值变化趋势。group1从时间1的73.3缓慢上升至时间4的82.3整体呈平稳递增趋势group2从时间1的69.4快速上升至时间4的85.3增速明显快于group1。两条折线呈现非平行趋势group2的斜率更大说明group2随时间的增长幅度更大。这种非平行的变化模式正是交互效应显著的直观体现——时间因素对因变量的影响因组别不同而存在差异。八、描述统计与均值汇总表7至表9分别报告了time各水平、group各水平以及time与group交叉分组的描述统计结果。表7显示因变量均值从时间1的71.381持续上升至时间4的83.773各时间点标准差在6.4~6.8之间离散程度较为稳定。表8显示两组的总体均值极为接近77.073 vs 77.646但group2的标准差9.442明显大于group16.412说明group2组内个体间的变异程度更大。表9的交叉均值进一步揭示了两组的变化差异。在时间1时group1的均值73.326高于group269.436存在约3.89个单位的基线差异。到时间2时两组均值趋于接近75.094 vs 75.569差异几乎消失。到时间3和时间4时group2的均值反超group180.310 vs 77.59085.267 vs 82.280且差距逐步扩大。这一变化模式清晰地展示了两组不同的发展轨迹group1起点较高但增速平缓group2起点较低但增速更快最终实现了反超。九、结论本研究采用重复测量方差分析探讨了因变量在不同时间点和不同组别之间的差异。主要结论如下第一球形度检验通过W0.869p0.213满足重复测量方差分析的前提假定。第二time主效应极其显著F163.752p0.001偏eta^20.759因变量随时间呈持续递增趋势各时间点之间两两差异均显著。第三group主效应不显著F0.122p0.729两组的总体均值无显著差异。第四time x group交互效应显著F14.879p0.001偏eta^20.222简单效应分析表明group2随时间的增长幅度显著大于group1两组的发展轨迹存在本质差异。建议在实际干预中关注不同组别的动态变化特征针对group2增速更快的特点制定差异化的策略。