我正在尝试使用lmer创建一个混合效果模型。子图应与图嵌套。处理应与子图嵌套。因此，子图中有3到7个处理，而一个图中总是3个子图。我创建了以下模型：
model <- lmer(Depth ~ Mass + (1|Plot:SubPlot:Treatment), data=mydata)
但这给我一个错误： 错误：每个分组因子的级别数必须为<观察数（问题：Plot：SubPlot：Treatment）
“ data.frame”：147磅。共有6个变量：
$情节：int 1 1 1 1 1 1 1 1 1 1 1 ...
$ SubPlot：int 1 1 1 1 1 1 1 1 2 2 2 ...
$处理：int 1 2 3 4 5 6 7 1 2 3 ...
$深度：num 0 4 4.5 5.5 6 6 6 3 3 4.5 6.5 ...
$ Mass：int 21 50 78 103 128 147 172 21 49 77 ...
以下是一些数据：
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
任何想法如何进行？