R 三层部分嵌套模型

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

• 受试者（身份证）
• 时间（等距重复测量代码1-10）
• 结果（每次重复测量）
• 治疗（心理治疗/等待名单控制为0/1）

我的第一个具有随机坡度和截距的两级模型运行良好，非常简单：

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

现在我想知道我是否应该使用一个三级部分嵌套模型，因为团体心理治疗受试者嵌套在治疗师中（在那里有几个治疗师提供治疗），但对照组没有嵌套。我猜我至少应该解释一下德容、莫尔贝克和范德莱登（2010）提出的治疗师的主要影响

De Jong，K.，Moerbeek，M.，和Van der Leeden，R.（2010）。纵向三水平多水平模型中的先验功率分析：以治疗师效应为例。心理治疗研究，20（3），273-284

我在“部分嵌套模型”下的以下链接中找到了一个非常有用的资源：

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

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

---

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

保持简单。只需在没有治疗师的受试者中添加一个虚拟治疗师

none

。然后拟合下面的模型

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

---

治疗师

无

与治疗

等待名单

混淆。治疗师是一种随机效应，因此会受到惩罚。治疗是一种固定的效果，不会受到惩罚。因此，

等待名单

的所有信息将进入治疗效果，

无

的治疗师效果将为零。

为什么不写信给相应的作者，特别是关于组和治疗之间的差异的问题？这似乎是对的，而且非常简单。我是传销新手，这就是为什么我想不起这个。我还是觉得自己很笨。这很有道理。非常感谢。