Оценка дисперсии зависимой переменной в гауссовой 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