我有以下问题。运行有序 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))
在后一种情况下,两个图的位置都很好,但是一旦我添加第二个图(或者更好,它被覆盖),第一个图仍然消失。任何有关如何防止这种行为的建议将不胜感激。
仔细阅读
?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
都不会对网格图形产生任何影响。
这似乎有效:
plot(eff,col=1,row=2,ncol=1,nrow=2,style="stacked",rug=F,
key.args=list(space="right"),more=T)
plot(eff2,col=1,row=1,ncol=1,nrow=2,style="stacked",rug=F,
key.args=list(space="right"))
编辑:太迟了...
gridExtra 包为此提供了简洁的语法
library(gridExtra)
plot1 <- plot(eff, grid=FALSE)
plot2 <- plot(eff2, grid=FALSE)
grid.arrange(plot1, plot2, ncol=2)