问三级部分嵌套模型

EN

我正在用R和lme4模拟随时间变化的群体心理治疗对象。我的数据有以下结构：

• 主题(id)
• 时间(等距重复测量的代码1-10 )
• 结果(对每一次重复测量)
• 治疗(心理治疗/等候名单控制为0/1)

我的第一个带有随机斜坡和拦截的两级模型工作得很好，而且很简单：

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

现在我想知道我是否应该使用一个三层部分嵌套模型，因为小组心理治疗对象是嵌套在治疗师(在那里有几个治疗师提供治疗)，但对照是非嵌套的。我猜我至少应该解释治疗师的主要影响，如der，Moerbeek & Van der Leeden (2010)所述。

De Jong，K.，Moerbeek，M. & 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)

他提供的数据实际上是相同的，但他在模型中添加了“组”变量。我不明白为什么这是因为治疗/控制分组和治疗/非治疗分组是一样的。如果被试接受治疗比他在实验组，如果不是他在对照组。您将如何编写这个三层部分嵌套模型？我知道这是一个交叉验证的论坛更多的问题，我张贴在那里没有回应，我真的感到困惑。谢谢。