我有一个包含以下列的矩阵:2个固定效应,一个因变量和2个与我的方差/协方差矩阵(GRM1和GRM2)的行和列匹配的ID。在早期版本的sommer中,我使用了mmer2函数
mmer2_result <-mmer2(y〜fixed_effect1 + fixed_effect2,random =〜g(id1)+ g(id2),G = list(id1 = GRM1,id2 = GRM2),rcov =〜units,silent = TRUE,数据= data,constraint = FALSE)
如何将其映射到更新软件包中的mmer函数?
答案 0 :(得分:0)
请花一些时间阅读文档,其中有很多示例说明如何在新版本中实现更改。答案:
mmer_result <- mmer(y ~ fixed_effect1 + fixed_effect2, random=~vm(id1, Gu=GRM1)+vm(id2, Gu=GRM2),rcov=~units, silent=TRUE, data=data)
最好