报错:
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)

方差组分结果:

Error: grouping factors must have > 1 sampled level_数据

这样就搞定了!