Связь функции "идентичности" с гамма-распределением в пакете VGAM R
Я использую VGAM
пакет, чтобы сделать регрессию между двумя выборками, которые распределены по гамме. Можно ли выбрать функцию идентификации для ссылки? Если да, то как?
Я использую эти строки кода, и это работает:
fit <- vglm(Y ~ X , family = gammaR, trace = TRUE)
Но когда я пытаюсь изменить функцию ссылки, она не работает, например:
fit <- vglm(Y ˜ X, family = gammaR(link= "identitylink"), trace = TRUE)
У меня есть следующее сообщение об ошибке:
unused argument (link = "identitylink")
1 ответ
Как и вы, я застрял с гамма-дистрибутивом в VGAM. К счастью, я нашел справочную карту VGAM!
Там немного загадочно читать, но в нем перечислены GLM, всегда заканчивающиеся на "ff", например, "poissonff", "binomialff", но без "gammaff". Тем не менее, это работает! Решение для меня было использовать:
fit <- vglm(Y ˜ X, gammaff(link= "identitylink"))
Подводя итог: "гаммафф" вместо гаммаР должен работать.