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 分布的假设。当我在 R 的 betareg 包中为相同的数据建模时不会发生这种情况。可能导致这种差异的原因是什么？

Answer 1

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

如果你使用faimly（错别字），gam 不会抱怨并继续进行高斯运算。试试：

print (mod)

并查看它是否显示“家庭：Beta 回归”或“家庭：高斯”