2

申し訳ありません。再現可能なデータセットに関する私の質問は次のとおりです。

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

2 つのプロットを 1 つの図に結合することは可能ですか?

par(mfrow=c(2,1))

これはうまくいきません。つまり、図は 2 つのグラフで別々に再現されていますが、同じ図ではありません。

4

2 に答える 2

3

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

ここに画像の説明を入力

于 2013-04-20T13:40:55.013 に答える
0

ベクトル インデックスを使用して、alleffects オブジェクトから特定の効果を参照できます。あなたの場合、私は t.effects[4:5] または同等に t.effects[c("colour:year", "colour: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])
于 2013-04-20T13:35:43.457 に答える