线性混合效果模型未在R中完成

时间:2017-04-24 10:30:52

标签: r statistics

线性混合效果模型没有在R中完成。

我使用包lmekin来运行模型,因为我将亲属矩阵作为随机效果拟合。我使用的公式是:

  

模型< - lmekin(疾病〜分数+年龄+(1 | IID),数据=表型_df,   varlist中= kinship_mx)

phenotype_df是一个数据框,有12000个人,每个人都有疾病状态(1或0),分数和他们的年龄。 (注意:我在这里将疾病状态视为定量)。 kinship_mx是亲属矩阵,即。一个正定矩阵(12000x12000),对角线〜0.5。

当我使用~1000个人时,这个模型在大约1分钟内运行。但是,我让模型在我们的大型服务器上与所有12000个人一起运行,并且在5天后还没有完成。任何人都可以提出一个理由,为什么会发生这种情况,以及是否有任何解决方案可以使模型完成?

1 个答案:

答案 0 :(得分:0)

回答这个问题 - 他们在三周后结束了。这需要很长时间。