使用 effects 包绘制 logit 概率

Question

我正在尝试运行一个 logit 模型并绘制一些重要预测变量的概率曲线。我正在尝试这样做 使用效果包。

df=data.frame(income=c(5,5,3,3,6,5),
              won=c(0,0,1,1,1,0),
              age=c(18,18,23,50,19,39),
              home=c(0,0,1,0,0,1))
str(df)

md1 = glm(factor(won) ~ income + age + home, 
          data=df, family=binomial(link="logit"))

summary(md1)

plot(effect("income", md1), grid=TRUE)

但我想知道如何绘制图表，以便显示基于收入（或任何其他预测变量）的获胜概率（响应）。

但是，我想做的是生成相同的图，y 轴为赢，x 轴为收入，但曲线显示年龄和家庭的概率（在单独的图中）。

在效果文档中没有看到如何执行此操作。救命！

谢谢。

Answer 1

如果我的理解正确，您需要在模型中引入交互项。例如，

md1 <- glm(factor(won) ~ income*age + home,
  data=df, family=binomial(link="logit"))

with(df, mean(age))
#[1] 27.83333
with(df, mean(age)-sd(age))
#[1] 14.33642
with(df, mean(age)+sd(age))
#[1] 41.33025

plot(effect("income*age", md1, 
  xlevels=list(income=0:10, age=c(14.3, 27.8, 41.3))), multiline=TRUE)

将向您展示收入对不同年龄值的韩元的影响。