r - R での GAM モデルの視覚化

Question

gamRでモデルの滑らかなコンポーネントを非常によく視覚化する方法を知っている人はいますか? 関数の出力のようなものを視覚化したいと思いますvisreg。以下のこのコードは私の問題を示しています

library(gam)

f=function(v){exp(v)}
n=100
x=runif(n)
t=runif(n)
y=x+f(t)+rnorm(n, sd=0.1)

fit=gam(y~x+s(t))

plot(t,y)
lines(t,as.numeric(fit$smooth))

#want something more like
library(visreg)
visreg(fit)

Answer 1

gam オブジェクトには plotting メソッドを使用できますが、次のdataパラメーターを使用する必要がありますgam。

library(gam)

f <- function(v){exp(v)}
n <- 100
x <- runif(n)
t <- runif(n)
y <- x+f(t)+rnorm(n, sd=0.1)
DF <- data.frame(y, x, t)

fit <- gam(y~x+s(t), data = DF)
layout(t(1:2))
plot(fit, se=TRUE)

help("plot.gam")その他のオプションについては、 を参照してください。