2

我有以下问题。运行有序的 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))

在后一种情况下，两个图的位置很好，但是一旦我添加第二个图（或者更好，它被覆盖），第一个图仍然消失。任何如何防止这种行为的建议将不胜感激。

4

2 回答 2

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()不起作用的原因是因为这个包使用的是基于网格系统的点阵图形，与基础图形不兼容。既不会也不会对网格图形产生任何影响。par()layout

于 2013-03-05T15:28:06.133 回答
1

这似乎有效：

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"))

编辑：太晚了...

于 2013-03-05T15:34:15.177 回答