我正在尝试使用伽马分布拟合混合效果模型。最基本的模型有一个固定预测器和一个随机效应。无论我指定哪个链接(我已经尝试过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,找不到另一个遇到此错误的人。我的语法不正确吗?
感谢您的帮助。