R的VGAM包中“恒等”函数与伽马分布的联系

Question

我正在使用VGAM包在两个伽马分布的样本之间进行回归。是否可以选择链接的标识函数？如果是，怎么做？

我正在使用这些代码行，它运行如下：

fit <- vglm(Y ~ X , family = gammaR, trace = TRUE)

但是，当我试图更改链接函数时，它不起作用，例如：

fit <- vglm(Y ˜ X, family = gammaR(link= "identitylink"), trace = TRUE)

我有以下错误消息：

unused argument (link = "identitylink")

Answer 1

我和你一样被困在了VGAM的伽玛分布中。幸运的是，我找到了一个VGAM-参考卡！

在那里读起来有点神秘，但是它列出了GLMs总是以"ff“结尾，例如"poissonff"，"binomialff"，但没有"gammaff”。然而，它是有效的！对我来说，解决办法是使用：

fit <- vglm(Y ˜ X, gammaff(link= "identitylink"))

总之："gammaff“而不是gammaR应该有效。