將「身份」函數鏈接到VGAM包中的伽瑪分佈R

我正在使用VGAM包對Gamma分佈的兩個樣本進行迴歸。是否可以選擇鏈接的身份識別功能？如果是，如何？將「身份」函數鏈接到VGAM包中的伽瑪分佈R

我使用這些代碼行，它的工作原理：

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

但是當我嘗試改變鏈接功能，這是行不通的，如：

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

我有以下錯誤信息：

unused argument (link = "identitylink")

2015-01-08 Didi86

我很難理解這將如何工作。伽馬函數只在正值上定義（因爲它的典型鏈接函數是'log'，所以這是有意義的）。如果鏈接是「身份」，算法應該如何加重負偏差？ –

它不會對觀察到的數據做出負面預測。問題解決了。 – Dason

感謝您的回答。我沒有提到我在Y的價值觀是非常積極的。對不起，我沒有得到你的觀點，因爲我不是VGLM的家庭成員。我認爲它會起作用，因爲當我使用帶有Gamma分佈和標識鏈接的相同變量的R的glm包時，我估計Y | X分佈的均值沒有問題，但是使用VGLM包時，不可能具有比例和比率參數。 – Didi86

我和你一樣卡在VGAM中的gamma分佈。幸運的是，我發現了一個VGAM-Reference Card！

這裏有點含糊不清，但它列出了總是以「ff」結尾的GLM，例如，「poissonff」，「binomialff」，但不包含「gammaff」。但是，它的工作原理！對我來說，解決方案是使用：

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

綜上所述：不是gammaR 「gammaff」應該工作。

2016-06-15 14:54:21 Jens

回答