R 混合效应模型

我正在尝试用lmer创建一个混合效果模型。子地块应与地块嵌套。处理应与子地块嵌套。所以在子地块中有3到7个处理，在一个地块中总是有3个子地块。我创建了以下模型：

model <- lmer(Depth ~ Mass + (1|Plot:SubPlot:Treatment), data=mydata)

---

有什么想法可以继续吗？

如果您包含一个简单的示例输入，可以用来测试和验证可能的解决方案，那么就更容易帮助您了。我添加了一些数据。有帮助吗？

治疗

在您的数据中缺失。您的错误消息表示每组中没有足够的观察结果。事实上，

Plot:Subplot:Treatment

结构在每组中只进行了1次观察。如果只使用2级分组变量而不是3级分组变量，则模型将运行（但不一定收敛）。处理列的名称不正确，现在它是正确的。

Plot SubPlot Treatment Depth Mass
1   1   1   0   21
1   1   2   4   50
1   1   3   4.5 78
1   1   4   5.5 103
1   1   5   6   128
1   1   6   6   147
1   1   7   6   172
1   2   1   3   21
1   2   2   4.5 49
1   2   3   6.5 77
1   2   4   7   102
1   2   5   8   127
1   2   6   9   146
1   2   7   10.5    171
1   3   1   3   21
1   3   2   1.5 49
1   3   3   1.5 77
1   3   4   1.5 102
1   3   5   1.5 127
1   3   6   1.5 146
1   3   7   1.5 171
2   1   1   3   21
2   1   2   5   50
2   1   3   5   78
2   1   4   7   103
2   1   5   9   128
2   1   6   9.5 146
2   1   7   10  171
2   2   1   1.5 21
2   2   2   4   50
2   2   3   5   78
2   2   4   9   103
2   2   5   10  128
2   2   6   10.5    146
2   2   7   10.5    171
2   3   1   0   21
2   3   2   0   50
2   3   3   0   78
2   3   4   0   103
2   3   5   0   128
2   3   6   0   146
2   3   7   0   171