r - 在 R 中可视化 GAM 模型

Question

有谁知道如何gam很好地可视化 R 中模型的平滑分量？我真的很想可视化函数的输出之类的东西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 对象使用绘图方法，但您必须使用以下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")其他选项。