私は私のモデルのコードは次のようになりますR.にmgcv
パッケージのgam
機能を使用して間隔上のデータへのGAMの(0,1)をフィッティングしています:gamモデルのベータファミリは、1より大きく0より小さい値に適合します。 (mgcv)
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未満である。これは、HAに想定されていない:私は観測値対フィットをプロットし、それはこのようになりますppen。これは、ベータ版配布の前提条件に違反します。 Rのbetareg
パッケージの同じデータをモデリングすると、この不一致が発生する可能性があります。
どのように「フィット」を決定していますか? 'predict'を使っているのなら、引数' type = "response" 'を使いたいでしょう。そうでなければ、範囲外になります。 'fitted'関数を使用している場合は、わかりません。 – Wayne
@Wayne私は 'fitted(mod)'を使って近似値を得ています。プロットは 'plot(y〜fitted(mod))'で行われます。 – colin