使用visreg控制刻面列数

Question

虽然我喜欢使用package visreg来可视化我的回归，但有一件事我还不能控制:切面时的列数。有关示例，请参阅下面的系数曲线广义加性回归：

library(dplyr)
library(mgcv)
library(visreg)
data(airquality)
test <- gam(
    Ozone ~ s(Temp, by = Month),
    data = airquality %>% mutate(Month = as.factor(Month))
)

如果我这样做了

visreg(test, xvar = "Temp", by = "Month", gg = TRUE)

我得到了一个1行5列的系数曲线。

​

​

有趣的是，如果我去掉gg = TRUE，它变成了2行。但无论是哪种情况，我都希望能够在刻面时控制列数和行数。到目前为止，通过操作visreg的省略号参数或直接尝试操作生成的ggplot对象，我都没有成功。

例如，如果我想用3行2列的gg = TRUE做visreg，我最好的机会是什么-或者有没有推荐的其他包？

Answer 1

您只需修改ggplot对象并按常规方式添加facet_wrap即可：

p <- visreg(test, xvar = "Temp", by = "Month", gg = TRUE)
p + facet_wrap(vars(Month), nrow = 3)

​

实际上，您不需要先创建p，这会产生相同的结果：

visreg(test, xvar = "Temp", by = "Month", gg = TRUE) +
  facet_wrap(vars(Month), nrow = 3)