R 在GLMM或lme中将族id指定为随机效果的步骤

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

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

---

m1因为您可能只有一行ID！为此，在交叉验证相关且可能有帮助的情况下，您需要对1个ID进行多个OB
。如果您试图为族估计单独的方差分量，您可能应该在公式的随机效应部分使用“族ID”。例如，如果您希望您的回答与家庭相关（可能是这样），那么这是有意义的。如果每个个体都有多个观察结果，则可以（另外）包括可能与家庭级随机效应相关的个体级随机效应。看来你没有。我会重新格式化您的问题，并在stats.stackexchange.com上询问您的问题是否更符合统计性质。
rand.eff=formula(paste("~1|",sub.ID))
fit <- try(lmekin(fixed=fix.eff,data=x,random = rand.eff,varlist=list(kmat)))