我正在尝试以一个居中标题并排绘制两个回归摘要。每个回归摘要由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()函数?

谢谢

最佳答案

您可以按如下所示修补功能:
使用dputcapture.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))

08-25 07:06