如何用R中的abline绘制glm模型系数？

Question

我正在努力使用abline绘制glm模型的cofficients。让我们来看看这个简单的2D例子：

d <- iris[51:150, c(3:4,5)]
d[,3] <- factor(d[,3])
plot(d[,1:2], col=d[,3])

glm模型产生4个系数：

m <- glm(formula = Species~Petal.Length*Petal.Width, data = d, family = "binomial")
m$coefficients
# (Intercept)             Petal.Length              Petal.Width Petal.Length:Petal.Width 
#  -131.23813                 22.93553                 63.63527                -10.63606 

如何用简单的abline绘制那些？

Answer 1

二项式模型通常不会像这样设置。您通常会有一个0 | 1响应变量（即预测单个物种中的样本）。也许是因为你的模型中只包含2个物种，它似乎仍然有效（当包括所有3个物种时，情况并非如此）。 第二个技巧是预测type="response"并舍入这些值以获得离散预测：

d$pred <- factor(levels(d[,3])[round(predict(m, type="response"))+1])
plot(d[,1:2], col=d[,3])
points(d[,1:2], col=d$pred, pch=4)

在这里我添加了一个＆＃34; X＆＃34;对于预测。如果颜色相同，那么预测是正确的。我计算了预测不正确的5个样本。