在GLMM或lme中将族ID指定为随机效果

Question

我试图在R中拟合广义线性混合模型。我有很大的谱系和基因型数据。我试过这个：

m1 <- lmer(Final_sx~(1 | ID)+cohort2+cohort3+cohort4+sex,
           data=solar_new_phen,family=poisson)

它给了我留言：

  

随机效应的分组因子的级别数   等于到n，即观察次数

我在大多数教程中都看到人们使用ID作为聚类变量。在亲属关系以及lmekin函数中，它是这样的：

rand.eff=formula(paste("~1|",sub.ID))
fit <- try(lmekin(fixed=fix.eff,data=x,random = rand.eff,varlist=list(kmat)))

我应该使用个人ID或家庭ID作为随机效果吗？我很迷惑。我有血统，我认为我应该使用Family ID作为聚类变量。如果有人可以指导我对随机效应的解释很少，那就太好了。

Answer 1

因为ID可能只有一行！为此，你需要很多1个ID