提取用于在 R 中绘制 GAM 参数组件图的值

Question

我执行了一个包含连续平滑项和分类变量的 GAM。我已经使用 plot(mod,residuals=T,all.terms=T,pages=1) 绘制了模型 (mod)。这会生成两个平滑参数以及参数化参数的图。我想提取用于制作这些图的值，以便我可以重新制作它们并使它们看起来更好。如果我将绘图保存在一个对象中，这将为我提供平滑项所需的一切，但不包含有关参数组件的任何信息：plot.mod=plot(mod,residuals=T,all.terms=T,select=0)。但是我看不到参数组件默认绘图的数字来自哪里。还有办法提取这些吗？

这是我到目前为止所做的一个可重现的例子

library(mgcv)

# create some data
data=data.frame(response=c(10,12,8,9,3,4,5,5,4,5,4,5,4,1),pred1=c(9,8,8,9,6,7,6,4,3,4,2,3,3,1),pred2=as.factor(c("A","C","B","B","A","A","C","B","C","A","C","B","A","B")),pred3=c(1,6,3,4,8,6,4,5,7,10,11,3,12,1))

# run the GAM
mod <- gam(response ~ s(pred1,k=8) + pred2 + s(pred3,k=5), data=data, family=gaussian(), method="REML")

# the default plot
plot(mod,residuals=T,all.terms=T,pages=1)

# save values in an object. But this only saves the smooth terms.
plot.mod=plot(mod,residuals=T,all.terms=T,select=0)

# How can I extract the values used to plot the parametric term?

我试图提取数据制作的情节：

Answer 1

根据 plot.gam 文档，termplot 用于参数项，因此

plot.para <- termplot(mod, se = TRUE, plot = FALSE)

将该绘图保存到列表中。
格式和其他的不一样，但是数据是有的

提取用于在 R 中绘制 GAM 参数组件图的值

Extract values used to make plot for parametric component of GAM in R

r

plot

gam

mgcv