Cla*_*eri 5 nested r lme4 multi-level
I am modeling change over time in group psychotherapy subjects using R and lme4. My data have the following structure:
My first two-level model with random slopes and intercepts works well and is simple:
lmer(outcome ~ time * treatment + (time | subject), data=data, REML=FALSE)
Now I was wondering if I should use a three-level partially nested model because the group psychotherapy subjects are nested within therapist (there where several therapist providing the treatment), but controls were non-nested. I am guessing I should account at least for main effects of therapists as argued by De Jong, Moerbeek & Van der Leeden (2010).
De Jong, K., Moerbeek, M., & Van der Leeden, R. (2010). A priori power analysis in longitudinal three-level multilevel models: an example with therapist effects. Psychotherapy Research, 20(3), 273-284.
I have found a very useful resource in the following link under "Partially nested models": http://rpsychologist.com/r-guide-longitudinal-lme-lmer
The author gives the following code for a model virtually identical to the one I would like to test:
lmer(outcome ~ time * treatment + (1 | group:subject) + (0 + time | therapist:subject) + (0 + time:treatment | group) + (0 + treatment | group), data=data)
The data he presents is virtually identical but he adds in the model the "group" variable. I don't understand why this is because treatment/control grouping is the same as treatment/nontreatment grouping. If a subject received treatment than he is in the experimental group, if not he in the control group. How would you write this three-level partially nested model? I know this is more a question for Crossvalidated forum, I posted it there with no response and I really feel confused. Thank you.
把事情简单化。none只需在没有治疗师的受试者中添加一个虚拟治疗师即可。然后拟合下面的模型。
lmer(outcome ~ time * treatment + (time | therapist/subject), data=data)
治疗师none对治疗感到困惑waiting list。治疗师是一种随机效应，因此受到惩罚。治疗具有固定效果，不会受到惩罚。因此，所有信息都将用于治疗效果，waiting list而治疗师效果none将为零。