重复测量方差分析在抑郁症研究中的应用

重复测量设计(repeated measurement design)在医学、生物学研究中较为常见,即在给予种或多种处理后,在多个时间点上从同一个受试对象重复获得指标的观察。
重复测量研究的目的是探讨同一研究对象在不同时间点某指标的变化情况,例如患者在治疗后(或手术后)1天、2天、1周、2周等,各时间点上某指标的变化。今天,我们将通过一个案例介绍重复测量方差的分析步骤及结果解读。

一、重复测量方差案例分析

抑郁症是一种严重影响全球数百万人的心理障碍,在心理学和医学研究中,了解治疗方法的效果以及病情变化如何随时间推移而变化,是至关重要的。在这个过程中常用的统计方法即为重复测量方差分析。接下来以其在抑郁症研究中的应用为例,帮助读者理解其原理、步骤及结果解读。

1、案例背景

本研究共有20名抑郁症患者,其中男性10名,女性10名,服用抗抑郁药物进行抑郁症的治疗。在第0周时对患者采用SF-36量表进行生命质量的测量,并在治疗的第4周、第8周时进行随访,同时用同样的量表进行测量。

本案例对20名患者基期(第0周)、第4周、第8周时的生命质量进行重复测量方差分析,收集部分数据如下:

(数据杜撰无实际意义)

2、SPSSAU软件操作

将整理好的数据上传至SPSSAU系统,在【实验/医学研究】模块选择分析方法【重复测量方差】,将变量拖拽到右侧相应分析框中,选择事后多重比较方法,本案例使用Bonferroni校正法;点击【开始分析】,操作如下图:

3、重复测量方差分析结果解读

重复测量方差分析包括三个步骤,接下来分别从球形度检验、组内效应分析、组间效应分析三个方面进行分析。

1)球形度检验

球形度检验通过检验数据不同时间点构成的协方差矩阵是否符合球形特征,来检验组内各水平之间之间的自相关程度。如果p


SPSSAU输出球形度检验结果如下:

从上表可以看到,球形度检验p值大于0.05,所以数据通过球形度检验。如果没有通过,应该选择相应的校正结果进行分析。接下来进行组内效应分析。
【提示】如果组内项水平数量为2(比如两个时间点)则不能进行球形度检验;如果有缺失数据,可能导致无法进行球形度检验。
2)组内效应分析

组内效应分析是对同一患者不同时间点测得生命质量进行差异性分析,SPSSAU输出组内效应分析结果如下

球形度检验通过,所以查看满足球形检验时的分析结果。从上表查看时间点的主效应和性别:时间点二者交互效应的显著性情况。分析上表可知,时间点主效应对应p值小于0.05,呈现出显著性;性别:时间点二者交互效应对应p值大于0.05,没有呈现显著性,说明不存在交互效应。

  • 事后多重比较

不同时间点之间呈现出显著性差异,可以使用事后多重比较进行时间之间两两比较,查看具体哪两组时间的换着生命质量得分之间具有显著差异。具体使用的方法为分析前选择的Bonferroni校正法进行事后多重比较,SPSSAU输出时间点的事后多重比较分析结果如下:

分析上表可知,三次时间点测量的患者生命质量得分两两之间均呈现出显著差异性。具体差异可查看SPSSAU输出的可视化图形及表格如下图:

从上图及表格数据可清晰看出,20名患者从初期生命质量平均得分为9分,到第4周生命质量平均得分提高到44.45分,第8周时生命质量平均得分已经提高到70.65分。说明随着治疗时间的推移,患者的生命质量得分显著提升,即抑郁情况逐渐转好。

3)组间效应分析

组间效应分析是指对不同性别的患者生命质量的差异性分析,SPSSAU输出组间效应分析结果如下:

分析上表可知,不同性别患者生命质量得分并未呈现出显著差异性,对应p值大于0.05,说明男性患者和女性患者的生命质量不存在差异性。

总结上述分析结果可知:

  • 组内效应分析时,不同时间点之间患者生命质量呈现出显著性差异,进一步通过事后多重比较进行对比,两两时间点之间均存在显著差异,且随着治疗时间的推移,患者的生命质量得分显著提升,抑郁情况逐渐转好。

  • 组间效应分析时,不同性别患者生命质量并未呈现差异性。

  • 时间点与性别之间不存在交互效应。SPSSAU输出汇总效应分析结果如下:

二、重复测量方差简单效应分析

当交互效应存在时,可继续进行简单效应分析。简单效应指A因素在某水平时,B因素不同水平之间差异比较。

比如本案例中,控制性别为男性或女性,比较不同时间点两两之间患者生命质量的差异;或者控制时间点为第0周、第4周、第8周,比较男性、女性之间生命质量的差异。


进行简单效应分析的前提条件是需要交互效应显著,如果交互效应不显著一般不需要进行简单效应分析。本案例交互效应不显著,所以不进行简单效应分析。如果需要进行简单效应分析,在SPSSAU系统重复测量方差分析中勾选【简单效应】即可输出简单效应分析结果。

假如本案例继续进行简单效应分析,SPSSAU输出简单效应分析结果如下:

分析上表可知,控制时间点为第0周、第4周、第8周,对比男性和女性之间生命质量均为呈现出差异性。控制性别为女性、男性,对比两两时间点之间生命质量均呈现出显著差异性。

若交互效应存在,也可通过下图及平均值表格进行具体差异性比较:

由于本案例并不存在交互效应,故不进行深入分析。

三、重复测量设计理论知识

1、重复测量设计

重复测量设计是指对同一受试对象的同一指标在不同时间点上进行多次测量,分析其不同时段的变化规律。重复测量设计通常使用重复测量方差分析进行研究。

由于同一受试者各个时间点的观察值间具有相关性,因此其数据分析与常规的统计方法存在较大区别,比如方差分析、t检验等都要求各样本数据之间相互独立,所以不能使用这类分析方法进行分析。

2、重复测量数据格式

1)数据格式

重复测量数据是指同一批样本(病例)在不同的时间点测量了多次数据,因此重复测量数据的特殊之处在于一定会有ID号(即样本或者病例号),以及时间点数据,如下图:


同一个ID会有多个时间点的数据,比如上图中有12个样本(12个ID号),并且测量5个时间点。那么就一定会有12*5=60行数据。同一个ID号会重复5次,同一个时间点会重复12次。

2)组内项&组间项

重复测量方差分析时涉及两个重要的术语名词,分别是组内项和组间项。

  • 组内项:表示对同一对象测试多次的标识项,比如上图中对同一样本在不同时间点的记录,称为组内项。

  • 组间项:表示不同样本组别的标识项,比如上图中对不同药物类型的记录,称为组间项。

在进行重复测量方差分析时,一定要有组内项(即不同时间点的数据,时间点要大于等于3次);但并不一定存在组间项(可有可无,组间项可同时存在多个)。

四、其他常见问题

1如果有缺失数据如何办?

重复测量方差要求为平衡数据,即不能有缺失数据,比如有12名被试,每名被试重复3次,则一定需要有12*3=36行数据;如果有某个被试数据缺失,那么可将该被试筛选出去;也可以先对数据进行填补(通常使用平均值)后再进行分析。
2)事后多重比较的类型选择说明?

如果选中事后多重比较,SPSSAU默认会对组内项和组间项均地事后多重比较。具体多重比较类型上,通常建议使用Bonferroni校正法较优。如果各组别样本不同时可使用scheffe,如果各组别样本完全相同可使用tukey法等。

3)为什么没有呈现球形度检验值?

如果组内项仅2个水平(比如仅测量2个时间点的数据),此时不需要进行球形度检验。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,772评论 6 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,458评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,610评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,640评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,657评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,590评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,962评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,631评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,870评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,611评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,704评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,386评论 4 319
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,969评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,944评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,179评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 44,742评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,440评论 2 342

推荐阅读更多精彩内容

  • 重复测量方差分析流程: 一、案例背景 当前有一项关于两种药物剂型在服用后,血液中浓度有无不同的研究。研究共有18名...
    spssau阅读 358评论 0 1
  • 一、前期准备 1.研究目的 重复测量方差:相关领域(比如医学研究)时,常常需要对同一观察单位重复进行多次测量,比如...
    spssau阅读 947评论 0 0
  • 在某些实验研究中,常常需要考虑时间因素对实验的影响,当需要对同一观察单位在不同时间重复进行多次测量,每个样本的测量...
    spssau阅读 548评论 0 2
  • 重复测量方差分析 MathJax.Hub.Config({tex2jax:{inlineMath:[['']]}}...
    毛线东东a阅读 41,465评论 0 20
  • 一、重复测量方差说明 在数据分析中,尤其是在医学研究领域,往往存在大量重复数据,比如想要观测某药物疗效或者胎儿生长...
    spssau阅读 257评论 0 2