三级部分嵌套模型

Question

我使用R和lme4模拟团体心理治疗受试者随时间的变化。 我的数据具有以下结构：

• subject（id）
• 时间（等间隔重复测量的代码1-10）
• 结果（针对每个重复的措施）
• 治疗（0/1用于心理治疗/等候名单控制）

我的第一个带有随机斜率和截距的两级模型效果很好而且很简单：

lmer(outcome ~ time * treatment + (time | subject), data=data, REML=FALSE)

现在我想知道我是否应该使用三级部分嵌套模型，因为小组心理治疗科目嵌套在治疗师（有几位治疗师提供治疗的地方），但控制是非嵌套的。我猜我应该至少考虑到治疗师的主要影响，如De Jong，Moerbeek＆amp; Van der Leeden（2010）。

De Jong，K.，Moerbeek，M。，＆amp; Van der Leeden，R。（2010）。纵向三级多层次模型中的先验功效分析：具有治疗效果的示例。心理治疗研究，20（3），273-284。

我在“部分嵌套模型”下的以下链接中找到了一个非常有用的资源： http://rpsychologist.com/r-guide-longitudinal-lme-lmer

作者给出了以下代码，该代码几乎与我想要测试的模型完全相同：

lmer(outcome ~ time * treatment + (1 | group:subject) + (0 + time | therapist:subject) + (0 + time:treatment | group) + (0 + treatment | group), data=data)

他提供的数据实际上是相同的，但他在模型中添加了“组”变量。我不明白为什么这是因为治疗/对照组与治疗/非治疗组相同。如果受试者接受治疗而不是他在实验组，如果不是他在对照组。你会如何编写这个三级部分嵌套模型？我知道这对于Crossvalidated论坛来说更像一个问题，我在那里发布了没有回复的内容，我真的感到困惑。谢谢。

Answer 1

保持简单。只需在没有治疗师的情况下为受试者添加假治疗师none。然后适合下面的模型。

lmer(outcome ~ time * treatment + (time | therapist/subject), data=data)

治疗师none与治疗waiting list混淆。治疗师是一种随机效应，因而受到惩罚。治疗是固定的，不会受到惩罚。因此，所有信息都将转至waiting list的治疗效果，none的治疗效果将为零。