具有高级绘图功能的多个绘图，尤其是 plot.rqs()

Question

我正在尝试将两个回归摘要与一个居中的标题并排绘制。每个回归摘要由plot.rqs() 生成，总共有 9 个图。

我已经尝试过使用par(mfrow=c(1,2))，但正如我从 Paul Murrel (2006) 的书中了解到的那样，像 plot.rqs() 或 pairs() 这样的高级函数会在绘制之前保存图形状态，然后再恢复一次图形状态已完成，因此对par() 或layout() 的抢先呼叫无法帮助我。 plot.rqs() 也没有“面板”功能。

看来实现结果的唯一办法就是修改plot.rqs()函数，得到一个新函数，比如说modified.plot.rqs()，然后运行

par(mfrow=c(1,2))
modified.plot.rqs(summary(fit1))
modified.plot.rqs(summary(fit2))
par(mfrow=c(1,1))

从那里我也许可以弄清楚如何使用layout() 为图像添加整体标题。有谁知道如何创建可以以这种方式使用的modified.plot.rqs() 函数？

谢谢

Answer 1

您可以按如下方式修补函数： 使用dput 和capture.output 检索 函数的代码，作为字符串； 根据需要更改它（在这里，我只是替换每次出现的par 一个什么都不做的函数）； 最后评估结果以生成新函数。

library(quantreg)
a <- capture.output(dput(plot.summary.rqs))
b <- gsub("^\\s*par\\(", "nop(", a)
nop <- function(...) {}
my.plot.summary.rqs <- eval(parse(text=b))

Answer 2

首先我们生成一个示例对象 fm 。然后我们复制plot.rqs 并在副本上使用trace 以在顶部插入par <- list，从而有效地取消函数内par 的任何使用。然后我们对plot.summary.rqs 做同样的事情。最后我们用我们自己的par来测试它：

library(quantreg)
example(plot.rqs) # fm to use in example

# plot.rqs
plot.rqs <- quantreg::plot.rqs
trace("plot.rqs", quote(par <- list), print = FALSE)

# plot.summary.rqs
plot.summary.rqs <- quantreg::plot.summary.rqs
trace("plot.summary.rqs", quote(par <- list), print = FALSE)

# test it out
op <- par(mfrow = c(2, 2))

plot(summary(fm))
plot(fm)
title("My Plots", outer = TRUE, line = -1)

par(op)

编辑：添加plot.summary.rqs。