如何用lmerTest包构建含嵌套随机效应的线性混合模型？

在lmer/lmerTest的语法体系中，嵌套随机效应需通过明确层级归属的写法来体现，核心是表达“某一组别完全包含于另一组别”的关系。以下针对你的需求给出具体修改方案：

假设你的嵌套关系为neighborhood嵌套在school中（即每个学校下包含多个独立社区），可采用两种等价写法：

方式1：斜杠/简化写法 #

model <- lmerTest::lmer(score ~ gender + student + (1|school/neighborhood), data = mydata)

方式2：显式交叉项写法（与方式1完全等价） #

model <- lmerTest::lmer(score ~ gender + student + (1|school) + (1|school:neighborhood), data = mydata)

关键说明： #

• 斜杠写法(1|school/neighborhood)会被自动解析为学校层面的随机截距(1|school)，加上同一学校内不同社区的随机截距差异(1|school:neighborhood)，精准体现嵌套层级。
• 若实际嵌套方向为school嵌套在neighborhood中（每个社区下包含多个学校），只需调换组别顺序：

  # 简化写法
  model <- lmerTest::lmer(score ~ gender + student + (1|neighborhood/school), data = mydata)
  # 等价显式写法
  model <- lmerTest::lmer(score ~ gender + student + (1|neighborhood) + (1|neighborhood:school), data = mydata)
  

• 全程保留lmerTest::lmer调用，即可正常输出固定效应的p值（默认采用Satterthwaite自由度近似，可通过ddf参数调整计算方法）。

内容的提问来源于stack exchange，提问作者Phoebe