我想从GAM绘制单个变量/项(使用library(gam)
)。
以下是使用“自动”数据的示例:
library(ISLR)
library(gam)
# Estimate model
mygam = gam(mpg ~ s(displacement,10) + s(horsepower,10), data=Auto)
# Plot all terms (okay)
par(mfrow=c(1,2))
plot(mygam, se=T)
# Try to plot single term (problem!)
par(mfrow=c(1,1))
plot(mygam, se=T,select=1)
# -> does not work, each term is plotted separately
我认为plot.gam
已更改(更改为plot.Gam
)。也许select
选项也已更改?
我尝试了terms
(https://www.rdocumentation.org/packages/gam/versions/1.16/topics/plot.Gam),但是无法正常工作。
有趣的是,如果我使用library(mgcv)
,则select
选项有效。但我想坚持使用library(gam)
。
有人知道如何使用Hastie的library(gam)
选择单个变量/项进行绘图吗?我正在使用R版本3.5.1和gam 1.16。
谢谢!