如何解决混合效果模型,lmer
我正在尝试使用lmer创建一个混合效果模型。子图应与图嵌套。处理应与子图嵌套。因此,子图中有3到7个处理,而一个图中总是3个子图。我创建了以下模型:
model <- lmer(Depth ~ Mass + (1|Plot:SubPlot:Treatment),data=mydata)
但这给我一个错误: 错误:每个分组因子的级别数必须为
“ 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
任何想法如何进行?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。