Rの中にablineを使ってglmモデル係数をプロットする方法は？

Question

私はablineを使ってglmモデルの係数をプロットするのに苦労しています。シンプルablineを有するものをプロットする方法

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 

：GLMモデル利回り4つの係数を

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

：この単純な2Dの例を取ることができますか？

Answer 1

通常、このような二項モデルは設定されていません。通常は、1つの0 | 1応答変数があります（つまり、1つの種の中にサンプルがあるかどうかを予測します）。たぶんあなたのモデルに含まれている2つの種だけがあるので、それはまだ動作しているようです（これは3つの種をすべて含む場合はそうではありません）。私は予測のための「X」を追加しましたここ

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) 

： 秒のトリックは、個別の予測を得るためにtype="response"と、これらの値ラウンドを予測することです。色が同じであれば、予測は正しいものでした。私は予測が間違っていたところで5サンプルを数えます。