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数据分析师在中国航信高科技产业园进行了面向测试度量的数据分析培训课程,培训人数近2 ...
2024-05-01CDA数据分析师走进深圳迈瑞生物医疗电子股份有限公司,在迈瑞总部展开了为期两天的培训,本次课程参训人员线上及线下近百人, ...
2024-05-01CDA数据分析师在合肥市对合肥阳光新能源科技有限公司开展了为期8天的企业内训。 合肥阳光新能源科技 ...
2024-05-01CDA数据分析师走进海尔大学,进行了《数据治理与数据中台建设的道与术》专题培训,培训现场爆满,近百人参加了此次培训。 ...
2024-05-01在中国银行苏州分行培训中心开始数据分析师培训,此次培训课程共10天内容,包括Excel、MySQL、概率论与数理统计、SPSS等内容, ...
2024-05-01从实际的业务需求出发,结合行业的典型应用特点,围绕实际的商业问题,探讨数据挖掘、机器学习模型在金融领域的应用,包括获客、信用评分、细分画像、交叉销售、反欺诈、违规识别、时序预测、运筹优化、流程挖掘九个方面,形成 ...
2024-05-01本次培训课程为线上+线下的模式,由于学员编程能力不一、部分学员没有编程基础,故提供统计学、python基 ...
2024-05-01华夏银行信用卡中心-机器学习培训 1、课程亮点 取材于业界一流企业和顶级咨询公司的行业实践;已经被证明是人人 ...
2024-05-01主 题:数据中台建设及数据分析应用主题分享 1. 数据中台市场洞察 2. 主流数据中台产品比较 3. 某企业数据中 ...
2024-05-01围绕“数据驱动”战略,全力打造我行 300 人数字化人才梯队,着力培养数字化管理人才、大数据专业团队 ...
2024-05-01在当今数据驱动的商业环境中,数据分析成为了企业决策的重要依据。通过对大量数据的收集、处理和分析,企业能够更好地理解市场 ...
2024-04-29在人工智能(AI)的世界里,提示词(Prompt)是一种强大的工具,它能够引导AI按照用户的需求产生特定的输出。本文将深入探讨AI ...
2024-04-29CDA立足未来职场,拓展前沿视野——对外经贸大学保险学院举办“三全育人大讲堂”分享行业最新动态。 ...
2024-04-294月2日,CDA数据分析师创始发起人兼协会理事长赵坚毅博士受邀在浙江万里学院举办了一场以“数字化能力在职场中的作用” ...
2024-04-29随机森林(Random Forests)现在机器学习中比较火的一个算法,是一种基于Bagging的集成学习方法,能够很好地处理分类和回归的问 ...
2022-12-23方差分析是数据分析中常用的一种统计分析方法,接下来让我们简单了解一下方差分析的基本思想和原理吧。 方差分析(Analysis ...
2022-12-23来源:关于数据分析与可视化 关于streamlit-aggrid 数据排序 表格样式的调整 数据 ...
2022-08-03作者:麦叔 定义 「把上面晦涩的概念汇成一句话就是:」 ❝ 回调函数就是一个被作为参 ...
2022-08-03现今,高学历人群日益增多,物以稀为贵的高学历光环淡去。无论本科生还是研究生,甚至博士生,求职竞争力都大不如前,就业压力越来越大。
2022-06-01某家企业10个人面试,有9个本科生……如何脱颖而出,除得体的举止和良好的沟通力外,证书成重要筹码,这也是很多人考证的关键所在。
2022-04-14