r - 在 GLMM 或 lme 中将 Family id 指定为随机效应

Question

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

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

它给了我信息：

随机效应的分组因子的水平数 等于n，即观察数

我在大多数教程中看到人们使用 ID 作为聚类变量。在 kinship 包以及 lmekin 函数中是这样的：

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

我应该使用个人身份证还是家庭身份证作为随机效应？我很困惑。我有家谱，我认为我应该使用 Family ID 作为聚类变量。如果有人可以指导我对随机效应的解释很少，那就太好了。

Answer 1

因为您可能只有一行 ID ！为此，您需要为 1 个 ID 提供许多 obs