線性混合效應模型未在R中完成

Question

我在R中未完成線性混合效果模型時遇到問題。

我使用包lmekin來運行模型，因為我將親屬關系矩陣擬合為隨機效果。 我使用的公式是：

模型<-lmekin（疾病〜得分+年齡+（1 | IID），數據=表型_df，varlist =血統_mx）

phenotype_df是一個數據框，包含12000個個體，每個個體的疾病狀態（1或0），得分和年齡。 （注意：我在這里將疾病狀況視為定量的）。 kinship_mx是一個親屬矩陣，即。 對角線〜0.5的正定矩陣（12000x12000）。

當我使用約1000個人時，此模型將在1分鍾左右運行。 但是，我保留了該模型以在大型服務器上與所有12000個人一起運行，並且在5天后還沒有完成。 誰能提出一個原因，說明為什么會發生這種情況以及是否有任何解決方案可以使模型完成？

Answer 1

只是為了回答這個問題-他們在三周后就完成了。 只需要很長時間。