r - 3 级逻辑回归

Question

我正在尝试使用 3 级逻辑回归分析具有二元响应的纵向数据。我的数据的层次结构看起来像这个例子 - 学生(1 级)嵌套在一个类(class)(2 级)中，而这些类(class)嵌套在一个学校(3 级)内。

任何人都可以在 SAS 上给我建议一些合适的读物吗？或 R执行 3 级混合效应逻辑回归的代码？我不是在寻求有关 3 级逻辑回归背后的理论的帮助。任何带有 SAS 的示例或 R代码会有很大帮助。

Answer 1

在 R 中，您可以使用 glmer()来自 的函数lme4 混合效应模型包。您在 formula= 中指定效果选项。对于嵌套的随机效果，您可能需要以下内容:

model <- glmer(formula = response ~ (1|student|class/block), family=binomial)

R-sessions上有一些多层次模型的例子。

并且您可以使用 family=binomial 拟合非线性响应，例如二项式.