วิธีการพล็อตขอบเขตการตัดสินใจใน R สำหรับรูปแบบการถดถอยโลจิสติก?

ฉันสร้างโมเดลการถดถอยโลจิสติกโดยใช้ glm ใน R. ฉันมีตัวแปรอิสระสองตัว ฉันจะพล็อตขอบเขตการตัดสินใจของโมเดลของฉันในพล็อตกระจายของสองตัวแปรได้อย่างไร ตัวอย่างเช่นฉันจะพล็อตร่างเช่น: http://onlinecourses.science.psu.edu/stat557/node/55

ขอบคุณ

set.seed(1234)

x1 <- rnorm(20, 1, 2)
x2 <- rnorm(20)

y <- sign(-1 - 2 * x1 + 4 * x2 )

y[ y == -1] <- 0

df <- cbind.data.frame( y, x1, x2)

mdl <- glm( y ~ . , data = df , family=binomial)

slope <- coef(mdl)[2]/(-coef(mdl)[3])
intercept <- coef(mdl)[1]/(-coef(mdl)[3]) 

library(lattice)
xyplot( x2 ~ x1 , data = df, groups = y,
   panel=function(...){
       panel.xyplot(...)
       panel.abline(intercept , slope)
       panel.grid(...)
       })

ฉันต้องสังเกตว่าการแยกที่สมบูรณ์แบบเกิดขึ้นที่นี่ดังนั้นglmฟังก์ชั่นจึงเตือนคุณ แต่นั่นไม่สำคัญที่นี่เพราะจุดประสงค์คือเพื่อแสดงให้เห็นถึงวิธีการวาดขอบเขตเชิงเส้นและการสังเกตสีตาม covariates ของพวกเขา

ต้องการที่จะตอบคำถามในความคิดเห็นต่อคำตอบที่ได้รับการยอมรับจากเฟอร์นันโด: ใครสามารถอธิบายตรรกะที่อยู่เบื้องหลังความชันและการสกัดกั้นได้หรือไม่?

สมมติฐานสำหรับการถดถอยโลจิสติกส์ในรูปแบบของ:

$g(z)$$z$

$y = 1$$h_{\theta} \geq 0.5$

ข้างต้นเป็นขอบเขตการตัดสินใจและสามารถจัดใหม่เป็น:

$y = mx + b$$m$$b$