1

パッケージmgcvのこのスレッドからの推奨に基づいて、特定の値でGAMモデルを実行して正確な見積もりを取得しようとしています。予測関数を使用すると、結果は GAM からの推定値と一致しません。

これが私が興味を持っている私のモデルです：

mod <- gam(y~s(a, by=b), data = dat)

次を使用して結果を取得する場合:

plot(mod)

結果は、次の方法を使用して変数 a の単一のプロット ポイントを推定しようとした場合とは異なります。

pdat <- with(dat,
             data.frame(a = 0,
                        b = mean(b,na.rm = TRUE)))
pred <- predict(mod, pdat, type = "response", se.fit = TRUE)
pdat <- transform(pdat, fitted = pred$fit)
pdat <- transform(pdat, upper = fitted + (1.96 * pred$se.fit),
                  lower = fitted - (1.96 * pred$se.fit))

私が間違っているかもしれないことについて何か提案はありますか?

4

1 に答える 1

0

試行錯誤の末、私は自分の質問に対する答えを見つけました。b を 1 に置き換える必要があるため、値の平均ではなく継ぎ目があります。これにより、グラフ内の正確な値が取得されます。

おそらく、これは将来誰かを助けるでしょう。

于 2014-02-22T16:07:26.643 に答える