R语言之纵向数据分析:多级线性增长模型2
这篇文章已经是纵向数据分析系列的第三篇了。之前,我们介绍了什么是纵向数据,我们如何把长型的数据集转换成纵向的数据集,并通过建立相应的多级模型进行分析。显然,仅仅介绍基本的多级模型并不足以对虚拟随机对照实验数据集进行分析。这篇文章则延续之前的分析,并介绍一种不错的方法来处理多级模型里的治疗效果问题。
和往常一样,我们还是先创建一个纵向数据集,然后把它转成长型的格式。
这是一个RCT数据集,也就说明,这里潜在存两个以疗效分组的差异。上一次,我们没有考虑异质性,而仅仅对两个组的普通疗效进行了细分。从直觉上看,我们可以添加治疗变量来修复这个模型,从而抓到组间的差异。
这里,我们任然使用lmerTest包是因为它允许修正模型的检验使用自由度。模型的公式和是和上次的一样,除非我们把治疗变量tx作为独立变量。
从上述的summary里的信息可以很明显的看到,time和tx这两个变量都相当显著,而接受B治疗的人群在抑郁症得分上比接受A治疗的人群低2.34分。那么,我们是否可以认为,B治疗方案效果更好?绝对不是。模型m1实际上并不合适对纵向数据集里的组间差异进行比较。m1里的治疗效果,其实也就是一段时间内的平均治疗效果。换句话说,治疗方案A和治疗方案B的疗效所导致的分差在-2.34分,其实只不过是1,2,3时间段内抑郁症均分的差。这是RCT,所以,我们所期待的志愿者参与两组治疗方案里时间段1的差要小于1。我们对两组疗效的均值差异感兴趣,而现在的轨迹在时间段2、3里存在差异。
这个概率可能对初学者来说有点复杂,但是其执行过程很简单,我们只需要添加交互变量time和tx。
在R,我们可以使用:符号来表明变量间的交互关系。因此,我们可以这样设置:time+tx+time:tx。但是,更简短的写法是:timetx。操作符要R涵盖主要变量和治疗疗效。此时,使用操作符就能简单实现各个效应的展示。请记住,操作符也可用于高维变量交互。比如,ABC就要R包含三个主要的效应,三个两两交互的效应,以及一个三元交互效应。
现在,我们从m2模型中获得一个三元交互关系。在这个模型当中,time效应指的就是疗效A(这是特定条件下的疗效和疗效时间的交互),而txB效应就是时间点为0时的效应(数据集里一个虚的时间点)。这个术语解释起来有点复杂:系数值反映了每个时间段内两种治疗方案的疗效上的差异,其条件是取平均时间和治疗效果。口头上说,人们可以理解的是,这表明两组治疗方案是如何随着时间的推移而不同。对于每个单位时间内的增长(即,从时间段1到时间段2),在完成平均治疗效应差的分析以后,接受B治疗方案的参与者理论上在抑郁症得分方面比接受A治疗方案的要低2.88分。所以,在时间段2内,接受B治疗方案的参与者得分差是 -2.88 x 2 + 3.43 = -2.33,而在时间段3内,接受B治疗方案的参与者得分差是-2.88 x 3+ 3.43 = -5.21。我们应当时刻注意,在一个给定的时间段里对比疗效差的时候,主要疗效以及疗效变量间的交互。但是,为什么要这样?你不妨设置分析一下,当你忽略了主要治疗效应的时候,你在时间段1内所得到的疗效差是怎样的情形。
我可以告诉你,模型m2比m1和m0(上一篇文章所采用的模型)都要好,原因在于,这个数据集是我们自己创建的。但在实际生活中,我们会在选择最好的一个模型的时候会遇到很多困难。其中一种选择的方式就是,按照我们之前所讲的,使用summary()函数,以此算出额外变量的统计显著性。其它一个常用的方法就是使用方差分析(ANOVA)表(仅仅针对嵌套模型),而这种方法适合用于这个模型。
这里,我们使用anova()对m0、m1和m2模型进行方差分析,原因在于我们把m0嵌套于m1,然后把m1嵌套于m2。如果还有另外的子集,我们就把模型嵌套到那个子集上。换句话说,含有独立变量A、B的模型则嵌套于含有变量A、B、C的模型中。
如果你做过方差分析,那这个方差分析表对于你来说应该是很熟悉了。Df列表明了模型里的自由度,它简单的反映了这个案例的参数估计。例如,自由度为4的模型m0表明了这个模型预测了4个参数(1.修正截距效应系数,2.时间效应,3.随机截距方差,4.残差)。
AIC和BIC是两个常用的拟合指数。AIC和BIC都考虑到两个因素:模型的数据拟合效果有多好,和这个模型复不复杂。AIC和BIC的不同之处在于模型复杂度的罚分上(BIC对模型复杂度的罚分影响更大)。我发现,AIC理解起来更容易,因为它渐进的使用留一交叉检验法(LOOCV)来预测运行效果。我们应当理解交叉检验和预测性能的概念,随后,如果你对此不太理解,你可以先把它跳过。
LogLik列模型参数的似然对数。从根本上说,它是模型里数据拟合效果的指标。其偏差是−2×logLik−2×logLik。为什么我们需要它?因为两个嵌套模型的偏差在原假设的条件下遵循卡方检验(假设两个嵌套模型的偏差一样)。此时,我们可以把它用作模型拟合差的测试。偏差的差则在Chisq列显示。卡方检验的统计相关自由度则在Chi Df(简单来说,就是额外参数的个数)。最后一行的结果显而易见:它算出了测试两个嵌套模型的差的p值。虽然是这样,我们还是要警惕ANOVA的结果,因为在众多的候选模型中,它会导致错误类型1上升极快。
由于模型m2比其它两个候选模型AIC和BIC都要好,而对于方差分析,我们可以认为这个模型是目前为止我们拟合的最好的了。当然,这个模型还有许多要改进的地方。后续,我们可能还可以在模型的性能上进行提升,并作图进行预测。
数据分析咨询请扫描二维码
CDA持证人Louis CDA持证人基本情况 我大学是在一个二线城市的一所普通二本院校读的,专业是旅游管理,非计算机非统计学。毕业之 ...
2024-12-18最近,知乎上有个很火的话题:“一个人为何会陷入社会底层”? 有人说,这个世界上只有一个分水岭,就是“羊水”;还有人说,一 ...
2024-12-18在这个数据驱动的时代,数据分析师的技能需求快速增长。掌握适当的编程语言不仅能增强分析能力,还能帮助分析师从海量数据中提取 ...
2024-12-17在当今信息爆炸的时代,数据分析已经成为许多行业中不可或缺的一部分。想要在这个领域脱颖而出,除了热情和毅力外,你还需要掌握 ...
2024-12-17数据分析,是一项通过科学方法处理数据以获取洞察并支持决策的艺术。无论是在商业环境中提升业绩,还是在科研领域推动创新,数据 ...
2024-12-17在数据分析领域,图表是我们表达数据故事的重要工具。它们不仅让数据变得更加直观,也帮助我们更好地理解数据中的趋势和模式。相 ...
2024-12-16在当今社会,我们身处着一个飞速发展、变化迅猛的时代。不同行业在科技进步、市场需求和政策支持的推动下蓬勃发展,呈现出令人瞩 ...
2024-12-16在现代商业世界中,数据分析师扮演着至关重要的角色。他们通过解析海量数据,为企业战略决策提供有力支持。要有效完成这项任务, ...
2024-12-16在当今数据爆炸的时代,数据分析师是组织中不可或缺的导航者。他们通过从大量数据中提取可操作的洞察力,帮助企业在竞争激烈的市 ...
2024-12-16在现代企业中,数据分析师扮演着至关重要的角色。他们不仅负责处理和分析大量的数据,还需要将这些分析结果转化为切实可行的商业 ...
2024-12-16在当今的大数据时代,数据分析已经成为推动企业战略的重要组成部分。无论是金融、医疗、零售,还是制造业,各个行业对数据分析的 ...
2024-12-16在当今这个以数据为驱动力的时代,数据分析领域正在迅速扩展与发展。随着大数据、人工智能和机器学习技术的不断进步,数据分析已 ...
2024-12-16在信息爆炸和数据驱动的时代,数据分析专业是否值得一选成为许多人思考的议题。无论是刚刚迈入大学校门的新生,还是考虑职业转型 ...
2024-12-16适合数据分析专业学生的实习岗位有很多,以下是一些推荐: 阿里巴巴数据分析岗位实习:适合经济、统计学、数学及计算机专业的 ...
2024-12-16在数据科学领域,探索实习机会是一个理想的学习和成长方式。实习不仅可以提供宝贵的实践经验,还能帮助学生发展关键的数据分析技 ...
2024-12-16在当今信息驱动的时代,数据分析不仅成为了企业决策的重要一环,还催生了各种职业机会。从技术到业务,数据分析专业的就业岗位种 ...
2024-12-16在现代企业中,数据分析师被誉为“数据探险家”,他们通过揭示隐藏在数据背后的故事,帮助公司优化业务策略和做出明智的决策。然 ...
2024-12-16在大数据崛起的时代,数据分析师被誉为企业的“幕后英雄”。他们通过解读数据,揭示隐藏的真相,为企业战略提供重要的指导。这份 ...
2024-12-16在这个信息大爆炸的时代,数据分析师成为了企业中的“福尔摩斯”,他们能够从庞杂的数据中提取关键洞察,为业务发展提供坚实支持 ...
2024-12-16在这个数据为王的现代社会,数据分析师如同企业的导航员,洞悉数据背后所隐藏的商业机会和战略优势。然而,成为一名优秀的数据分 ...
2024-12-16