Ошибка в спецификации гамма-ссылки для lmer
Я пытаюсь согласовать модель со смешанными эффектами с гамма-распределением. Самая базовая модель имеет один фиксированный предиктор и 1 случайный эффект. Независимо от того, какую ссылку я указываю (я пробовал log, identity и inverse), я получаю следующую ошибку. Мои реальные данные имеют нули в Y, но даже когда я применяю смоделированные данные только с положительным Y, как показано ниже, они выдают ту же ошибку.
mockdf = data.frame(y = rnorm(100,77,6.5), x1 = sample(letters,100,replace = T), x2 = seq(1900,1999,1))
mod = lmer(y ~ (1|x1) + x2, family = gamma(link = 'identity'), na.action = na.exclude, data = mockdf)
Error in gamma(link = "identity") :
supplied argument name 'link' does not match 'x'
Я искал через SO и не смог найти другого человека, который столкнулся с этой ошибкой. Мой синтаксис неверен?
Спасибо за вашу помощь.