在R中绘制GLM模型

我正在做一个线性回归模型，所以现在我有

mylogit <- glm(Gender ~ Drug*Environment + Ethnicity + Age, data = mydata, family = "binomial")    

---

但这只是将它分成两个图，每个图中有两条线，并且没有在上面绘制单独的数据点

与此未经测试的代码大致相同。（如果需要测试代码，请提供数据。）

---

您可以根据需要添加其他组，并注意predict.glm需要种族变量的有效值。

这里需要澄清，

Gender

是您的目标变量吗？我不太使用glm，所以我可能会遗漏一些东西。是的，性别是因变量。嗯，例如，如果这些是数据（忘记种族）。啊，对不起，我不知道该怎么做！性别=（M，M，F，F，F，F，M，F，M，F，M）药物=（Y，Y，Y，N，N，N，N，N，N，N）环境=（H，L，H，H，H，H，L，L，L，L，L，L，L，L，L，L，M）年龄=（20，33，55，45，46，47，23，31，43，38）你应该把dput（mydata）的结果作为编辑而不是评论发布。上面的向量不是R码，也没有种族。

plot(effect("Gender*Drug*Environment*Age", mylogit, xlevels=list(Age=20:60)), + multiline=TRUE, rug=FALSE)

plot(x=seq(20, 65, by=5), predict( mlogit, 
                    newdata=data.frame(Drug="Y", Environment="H", Ethnicity="White",
                   Age=seq(20, 65, by=5) ), type="response" ) )
lines(x=seq(20, 65, by=5), predict( mlogit, 
                   newdata=data.frame(Drug="N", Environment="L", Ethnicity="White",
                   Age=seq(20, 65, by=5) ), type="response" ) )