报错:
Error: grouping factors must have > 1 sampled level
报错代码:
分析一年多点数据,需要考虑地点和地点内区组,都用混合线性模型的随机因子。
> ## 建模
> blup <- lmer(First~(1|Line)+(1|Loc)+(1|Rep%in%Loc)+(1|Line:Loc), data=phe)
Error: grouping factors must have > 1 sampled level
这里%in%
换为/
。
正确的代码如下:
> ## 建模
> blup <- lmer(First~(1|Line)+(1|Loc)+(1|Loc/Rep)+(1|Line:Loc), data=phe)
> summary(blup)
方差组分结果:
这样就搞定了!