#r
#r
Вопрос:
У меня есть матрица со следующими столбцами: 2 фиксированных эффекта, зависимая переменная и 2 идентификатора, которые соответствуют строкам и столбцам моей матрицы дисперсии / ковариации (GRM1 и GRM2). В предыдущей версии sommer я использовал функцию mmer2
mmer2_result <- mmer2 (y ~ fixed_effect1 fixed_effect2, random=~g (id1) g (id2), G=список (id1= GRM1, id2 = GRM2), rcov = ~ единиц, молчание = ИСТИНА, данные =data, ограничение = FALSE)
Как бы мне сопоставить это с функцией mmer в обновленном пакете?
Ответ №1:
Пожалуйста, уделите некоторое время чтению документации, там есть множество примеров, показывающих, как изменения реализованы в новой версии. Ответ:
mmer_result <- mmer(y ~ fixed_effect1 fixed_effect2, random=~vm(id1, Gu=GRM1) vm(id2, Gu=GRM2),rcov=~units, silent=TRUE, data=data)
Лучшие,