R 使用effects软件包绘制logit概率图

我试图运行一个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轴为收入，但曲线显示年龄和家庭的概率（在单独的图中）

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

---

谢谢。

如果我理解正确，您需要在模型中引入一个交互术语。比如说,

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)

---

md1对，这不是问题。我试图在x轴上绘制收入曲线，在y轴上绘制收入曲线，以及代表不同预测值的曲线。因此，预测的家庭问题将被绘制在一个图表上，其中x轴在x轴上有收入，y轴上有赢家。我编辑了我的回答，因为我没有得到你想要的。希望这个答案更好。。。！在stats.stackexchange.com上可能更好，因为它看起来更像是一个统计问题，而不是一个
R
编程问题。问题是交互项可能与响应变量无关，因此与统计模型无关。我试图找出如何在R中做到这一点，但并没有真正关注统计数据，但也许我应该把这个问题转移到那里。