R语言线性混合效应模型实战案例

处理分组数据和复杂层次结构的分析师,从嵌入在参与者中的测量,嵌套在州内的县或嵌套在教室内的学生,经常发现他们需要建模工具来反映他们数据的这种结构。

由Kaizong Ye,Coin Ge撰写

在R中,有两种主要的方法来拟合多级模型,这些模型考虑了数据中的这种结构。

介绍

这些教程将向用户展示如何使用lme4R中的包来拟合线性和非线性混合效果模型,以及如何使用rstan以完全适合贝叶斯多级模型。这里的重点是如何使模型适合R而不是模型背后的理论。有关多级建模的背景知识,请参阅参考资料。 

本教程将介绍如何lme4  设置和运行一些基本模型,其中包括:

  • 在R中构造变化的截距,变化的斜率以及变化的斜率和截距模型
  • 从混合效应模型中生成预测和解释参数
  • 广义和非线性多层次模型
  • 完全贝叶斯多级模型适合rstan或其他MCMC方法

设置 环境

在R中开始多级建模很简单。lme4是在R中实现多级模型的规范包,尽管有许多包依赖并增强其功能集,包括贝叶斯扩展。lme4 最近已被重写以提高速度并整合C ++代码库,因此封装的功能有些不断变化。 

要安装lme4,我们只需运行:

读入数据

多级模型适用于特定类型的数据结构,其中单元嵌套在组内(通常为5个以上组),并且我们希望对数据的组结构进行建模。对于我们的介绍性示例,我们将从lme4文档中的一个简单示例开始,并解释模型正在执行的操作。 

模型

让我们首先拟合一个简单的OLS回归 。 

R模型接口非常简单,首先指定因变量,然后是 ~符号,预测变量,每个都被命名。加法符号表明这些被建模为加性效应。最后,我们指定要计算模型的数据。这里我们使用该lm函数执行OLS回归,但R中还有许多其他选项。

如果我们想要提取诸如AIC之类的度量 。

这导致模型拟合较差。现在让我们看一个简单的模型。

拟合不同的 模型

我们的下一步可能是使用分组变量(如学校或班级)来拟合不同的 模型。 

这通常被称为固定效应 。 

学校因子大大提高了我们的模型拟合度。但是,我们如何解释这些影响呢?

在这里,我们可以看到我们有一个完美平衡的设计,在课堂和学校的每个组合中有50个观察结果 。

让我们尝试对这些独特的因素进行建模。 

这非常有用,但如果我们想了解学校的影响和课堂的影响,以及学校和班级的影响,该怎么办? 

探索随机斜率

另一种选择是为每个学校和班级组合建立一个单独的模型。如果我们认为我们的变量之间的关系可能高度依赖于学校和班级组合,我们可以简单地拟合一系列模型并探索它们之间的参数变化:

我们将在未来的教程中更深入地讨论此策略,包括如何在此命令中生成的模型列表中进行性能推断。

建立不同的斜率模型

虽然上述所有技术都是解决这一问题的有效方法,但当我们明确感兴趣的是群体之间的变化时,它们并不一定是最好的方法。这是混合效果建模框架有用的地方。现在我们使用lmer具有熟悉的公式接口的函数, 使用特殊语法指定组级变量:(1|school) ,使lmer拟合具有变量截距组效果的线性模型school

我们可以使用多个组来拟合多个组效果。

最后,我们可以通过以下语法拟合嵌套:

在这里(1|school/class),我们想要为1|学校和学校内的课程设置不同截距的混合效应。

用lmer拟合变化的斜率模型

但是,如果我们想要探索不同学生水平指标的影响,因为它们因教室而异。我们可以拟合不同的斜率模型,而不是按学校(或学校/班级)拟合模型。在这里,我们修改我们的随机效应项,在分组术语之前包含变量:(1 + open|school/class)告诉R拟合变化的斜率和不同的学校和学校类别的截距模型,并允许open变量的斜率因学校而异。

结论

在R语言和生态系统中,拟合混合效应模型和探索组变异非常容易。在以后的教程中,我们将探索模型的比较,使用混合效果模型进行推理,以及创建混合效果模型的图形表示了解它们的效果。


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498