r - margins.plot: 'which' 引数を使用して、プロットに含めるマージンを選択します

Question

ロジスティック回帰に基づいて r の限界効果をプロットしようとしています。例えば：

data <- mtcars
mod <- glm(am ~ cyl + hp + wt + mpg, family = binomial, data = data)

library(margins)
marg <- margins(mod, atmeans = TRUE)
summary(marg)

margins plotコマンドを実行できます。

plot(marg)

すべての IV の限界効果と信頼区間をプロットします。cylプロットに含めたいのは、hp関心のある説明変数のみです。r のドキュメントによると、これは文字ベクトルを取る 'which' 引数を使用して実現できます。ただし、ドキュメントには、この引数の使用方法は記載されていません。「どの」引数を使用して margins.plot に選択した限界効果のみをプロットするように依頼する方法を知っている人はいますか? 残念ながら、上にリンクされているマージン プロットのヘルプ ページには例がありません。

プロット イメージ

Answer 1

variablesプロットする前に、関数内のオプションで目的の変数を指定できますmargins()。

mod <- glm(am ~ cyl + hp + wt + mpg, family=binomial, data=mtcars)

library(margins)
marg <- margins(mod, variables=c("cyl", "hp"))

plot(marg)

与えます: