R 中的廣義線性模型 (GLM)

Question

我有一個響應變量 (A)，我從數據 (X) 轉換 (logA) 和預測變量 (B)，它們都是連續的。 如何使用 R 中的廣義加性模型 (GAM) 檢查兩個變量之間的線性。我使用以下代碼

model <- gamlss(logA ~ pb(B) , data = X, trace = F)

但我不確定，當 logA 在 GLM 中連續時，我可以在代碼中添加"family=Poisson"嗎？ 對此有何想法？

提前致謝

Answer 1

如果您的因變量是計數變量，則可以使用family=PO()而不進行對數轉換。 使用family=PO()已應用日志鏈接來轉換變量。 請參閱gamlss 系列的幫助頁面以及計數回歸第 2.1 節的小插圖。

所以它會像：

library(gamlss)
fit = gamlss(gear ~ pb(mpg),data=mtcars,family=PO())

您可以看到預測是對數轉換的，您需要采用指數：

with(mtcars,plot(mpg,gear))
points(mtcars$mpg,exp(predict(fit,what="mu")),col="blue",pch=20)