Соммер, отображающий функцию mmer2 на mmer

Question

У меня есть матрица со следующими столбцами: 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 = ~ единиц, молчание = ИСТИНА, данные = данные, ограничение = ЛОЖЬ) </p>

Как мне сопоставить это с функцией mmer в обновленном пакете?

Answer 1

Пожалуйста, найдите время, чтобы прочитать документацию, есть множество примеров, показывающих, как изменения реализованы в новой версии.Ответ:

mmer_result <- mmer(y ~ fixed_effect1 + fixed_effect2, random=~vm(id1, Gu=GRM1)+vm(id2, Gu=GRM2),rcov=~units, silent=TRUE, data=data)

Лучший,