我正在尝试以一个居中标题并排绘制两个回归摘要。每个回归摘要由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()
函数?谢谢
最佳答案
您可以按如下所示修补功能:
使用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))