gam模型中的Beta族拟合值大于1且小于0.什么事情发生了？ （mgcv）

Question

我使用R中gam包的mgcv函数将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（拼写错误），gam似乎没有抱怨并继续前进并执行高斯。尝试：

print (mod)

看看它是否说“家庭：Beta回归”或“家庭：高斯”