线性混合效果模型没有在R中完成。
我使用包lmekin来运行模型,因为我将亲属矩阵作为随机效果拟合。我使用的公式是:
模型< - lmekin(疾病〜分数+年龄+(1 | IID),数据=表型_df, varlist中= kinship_mx)
phenotype_df
是一个数据框,有12000个人,每个人都有疾病状态(1或0),分数和他们的年龄。 (注意:我在这里将疾病状态视为定量)。
kinship_mx
是亲属矩阵,即。一个正定矩阵(12000x12000),对角线〜0.5。
当我使用~1000个人时,这个模型在大约1分钟内运行。但是,我让模型在我们的大型服务器上与所有12000个人一起运行,并且在5天后还没有完成。任何人都可以提出一个理由,为什么会发生这种情况,以及是否有任何解决方案可以使模型完成?
答案 0 :(得分:0)
回答这个问题 - 他们在三周后结束了。这需要很长时间。