【发布时间】:2019-03-06 20:57:28
【问题描述】:
我有一个包含以下列的矩阵:2 个固定效应、一个因变量和 2 个与我的方差/协方差矩阵(GRM1 和 GRM2)的行和列匹配的 ID。在之前版本的 sommer 中,我使用了 mmer2 函数
mmer2_result
如何将其映射到更新包中的 mmer 函数?
【问题讨论】:
标签: r
我有一个包含以下列的矩阵:2 个固定效应、一个因变量和 2 个与我的方差/协方差矩阵(GRM1 和 GRM2)的行和列匹配的 ID。在之前版本的 sommer 中,我使用了 mmer2 函数
mmer2_result
如何将其映射到更新包中的 mmer 函数?
【问题讨论】:
标签: r
请花点时间阅读文档,有大量示例展示了如何在新版本中实现更改。答案:
mmer_result <- mmer(y ~ fixed_effect1 + fixed_effect2, random=~vm(id1, Gu=GRM1)+vm(id2, Gu=GRM2),rcov=~units, silent=TRUE, data=data)
最好的,
【讨论】: