高斯GAM（MGCV程序包）中因变量的方差估计？

Question

考虑以下代码：

library(mgcv)

set.seed(123)
X = runif(300, 0, 1)
set.seed(123)
Y = X^3 + 2*X^2 + 1 + rnorm(300)

model = gam(Y~s(X), family= gaussian)

因此model是高斯广义加性模型（GAM）。如何在Y中找到因变量（model）的估计方差？

更新：在广义加性模型中，当族为高斯时，尺度参数等于Y的方差。因此，我认为我可以使用summary(model)$scale，它实际上提供了比例参数估计，但也可以认为与Y的方差估计相等。

Answer 1

您可以通过拟合模型的sig2组件直接从模型对象获取此信息：

> summary(model)$scale
[1] 0.9006256
> model$sig2
[1] 0.9006256

scale.estimated组件还会告诉您这是由模型估算还是由提供：

> model$scale.estimated
[1] TRUE