gam模型擬合值中的Beta族大於1且小於0。這是怎么回事？ （mgcv）

Question

我使用R中的mgcv軟件包的gam函數將gam擬合為間隔（0,1）上的數據。我的模型代碼如下所示：

mod <- gam(y ~ x1 + x2 + s(latitude, longitude), faimly=betar(link='logit'), data = data)

模型擬合得很好，但是當我繪制擬合值與觀察值時，看起來像這樣：

plot(data$y ~ fitted(mod), ylab='observed',xlab='fitted')

顯然，該模型的擬合值大於1且小於0。這不應該發生。 它違反了beta分布的假設。 當我在betareg包中為R建模相同的數據時，不會發生這種情況。什么可能導致此差異？

Answer 1

mod <- gam(y ~ x1 + x2 + s(latitude, longitude), faimly=betar(link='logit'), data = data)

如果您faimly使用（typo），則gam不會抱怨，而是繼續進行高斯faimly 。 嘗試：

print (mod)

並查看是否顯示“家庭：Beta回歸”或“家庭：高斯”