r - 如何将三个单独的交互图组合成一个图形？

Question

对此表示歉意。这是我对可重现数据集的问题：

library(effects)
data(Arrests)

Arrests$year <- as.factor(Arrests$year)

arrests.mod <- glm(released ~ employed + citizen + checks + colour*year + 
                   colour*age, family=binomial, data=Arrests)

t.effects  <- allEffects(arrests.mod)

plot(t.effects, "colour:year")
plot(t.effects, "colour:age")

是否可以将这两个图组合成一个图形？

par(mfrow=c(2,1)) 

这行不通。即这些数字分别在两个图表中再现，但不在同一个图中。

Answer 1

par(mfrow=c(2,1))不适用于网格图。它仅适用于基本图形。您可以使用gridExtra来排列格子图。

library(gridExtra)
p1 <- plot(t.effects, "colour:year")
p2 <- plot(t.effects, "colour:age")
grid.arrange(latticeGrob(p1),
            latticeGrob(p2))

Answer 2

您可以使用矢量索引从您的 alleffects 对象中引用特定效果。在你的情况下，我相信 t.effects[4:5] 或者，等效地， t.effects[c("colour:year", "color:age")]

library(effects)
data(Arrests)

Arrests$year <- as.factor(Arrests$year)
arrests.mod <- glm(released ~ employed + citizen + checks + colour*year + colour*age, family=binomial, data=Arrests)
t.effects <- allEffects(arrests.mod)

plot(t.effects[4:5])