如何将三个独立的交互图​​组合成一个图？

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（“color：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])