r - 结合 R 中使用效果包创建的两个图

Question

我有以下问题。运行有序的 logit 模型后，我想要 R 的效果包来可视化结果。这很好用，我对两个自变量这样做了，然后我试着把这两个图结合起来。但是，这似乎不起作用。我在这里提供了一个可复制的小例子，这样您就可以自己查看我的问题:

library(car)
data(Chile)
mod <- polr(vote ~ age + log(income), data=Chile)

eff <- effect("log(income)", mod)
plot1 <- plot(eff, style="stacked",rug=F, key.args=list(space="right"))

eff2 <- effect("age", mod)
plot2 <- plot(eff2, style="stacked",rug=F, key.args=list(space="right"))

我现在可以独立打印这两个图，但是当我尝试将它们一起绘制时，第一个图被覆盖了。我尝试设置 par(mfrow=c(2,1))，但没有用。接下来我尝试了以下操作:

print(plot1, position=c(0, .5, 1, 1), more=T)
print(plot2, position=c(0,0, 1, .5))

在后一种情况下，两个地 block 的位置都很好，但是在我添加第二个地 block 后，第一个地 block 仍然消失了(或者更好的是，它被覆盖了)。任何有关如何防止此行为的建议都将不胜感激。

Answer 1

阅读 ?print.eff 的一长串参数，我们发现有一些参数可以做到这一点:

plot(eff, style="stacked",rug=F, key.args=list(space="right"),
     row = 1,col = 1,nrow = 1,ncol = 2,more = TRUE)
plot(eff2, style="stacked",rug=F, key.args=list(space="right"),
     row = 1,col = 2,nrow = 1,ncol = 2)

par() 没有工作的原因是因为这个包使用了lattice 图形，它基于grid 系统，这与基本图形不兼容。 par() 和 layout 都不会对网格图形产生任何影响。