说,我有一个看起来像这样的 GAM:
# Load library
library(mgcv)
# Load data
data(mtcars)
# Model for mpg
mpg.gam <- gam(mpg ~ s(hp) + s(wt), data = mtcars)
Run Code Online (Sandbox Code Playgroud)
现在,我想使用ggplot2. 所以,我plot.gam用来生成我需要的所有信息,如下所示:
foo <- plot(mpg.gam)
Run Code Online (Sandbox Code Playgroud)
这也会产生一个不需要的数字。(是的,我意识到我在抱怨绘图函数绘制了一些东西......)visreg以相同的方式使用时,我只是指定plot = FALSE抑制数字,但plot.gam似乎没有这个选项。我的第一个想法是可能invisible会完成这项工作(例如,invisible(foo <- plot(mpg.gam))),但这似乎不起作用。是否有一种简单的方法可以在不将不需要的图形输出到文件的情况下执行此操作?
好的,所以我终于在发布这个 5 分钟后想通了。尽管默认行为是绘制所有项,但可以选择绘制哪个项(例如,select = 1是第一项,select = 2是第二项)。但是,如果我使用select = 0它不会绘制任何内容并且不会给出错误,但会返回完全相同的信息。一探究竟:
# Load library
library(mgcv)
# Load data
data(mtcars)
# Model for mpg
mpg.gam <- gam(mpg ~ s(hp) + s(wt), data = mtcars)
# Produces figures for all terms
foo1 <- plot(mpg.gam)
# Doesn't produce figures
foo2 <- plot(mpg.gam, select = 0)
# Compare objects
identical(foo1, foo2)
[1] TRUE
Run Code Online (Sandbox Code Playgroud)
好礼!
| 归档时间: |
|
| 查看次数: |
325 次 |
| 最近记录: |