我使用VGAM
包在两个Gamma分布的样本之间进行回归。
是否可以为链接选择身份功能?如果是,怎么样?
我使用这些代码行并且可以正常工作:
fit <- vglm(Y ~ X , family = gammaR, trace = TRUE)
但是当我尝试更改链接功能时,它不起作用,例如:
fit <- vglm(Y ˜ X, family = gammaR(link= "identitylink"), trace = TRUE)
我有以下错误消息:
unused argument (link = "identitylink")
答案 0 :(得分:0)
我和VGAM中的伽玛分布一样被卡住了。幸运的是,我找到了VGAM-Reference Card!
阅读时有点神秘,但它列出的GLM总是以“ff”结尾,例如: “poissonff”,“binomialff”,但没有“gammaff”。但是,它的工作原理!我的解决方案是使用:
fit <- vglm(Y ˜ X, gammaff(link= "identitylink"))
总结:“gammaff”代替gammaR应该有效。