1

変数の対数を使用してモデルを適合させています。変換された変数に関して結果を提示したいと思います。簡単な例を次に示します。

library(mgcv)
N = seq(from=1,to=10,by=.01)
a = N^2
b = 1/N
y = log(a*b*runif(length(N)))
mod = gam(y~te(a,b))
vis.gam(mod, plot.type="contour")

同じグラフを作成したいのですが、y ではなく exp(y) に基づいています。

これを行う方法はありますか?前もって感謝します。

4

1 に答える 1

0

よくあることですが、質問をした直後に解決策が思い浮かびました。

入力vis.gamして関数を取得するだけです。これをコピーしてスクリプト ファイルに貼り付けます。長い関数です。途中のどこか:

fv <- predict.gam(x, newdata = newd, se.fit = TRUE, type = type)
fv$fit = exp(fv$fit)   #INSERT THIS LINE HERE WITH WHATEVER FUNCTION YOU WANT TO MODIFY THE FITTED VALUES BY
z <- fv$fit

次に、関数名を「mod.vis.gam」などに変更し、他の関数と同様に使用します。よく働く。SE は影響を受けないことに注意してください。

于 2013-03-31T21:58:53.107 に答える