Làm thế nào để vẽ ranh giới quyết định trong R cho mô hình hồi quy logistic?


15

Tôi đã thực hiện mô hình hồi quy logistic bằng glm trong R. Tôi có hai biến độc lập. Làm thế nào tôi có thể vẽ ranh giới quyết định của mô hình của tôi trong biểu đồ phân tán của hai biến. Ví dụ: làm thế nào tôi có thể vẽ một hình như: http://onlinecferences.science.psu.edu/stat557/node/55

Cảm ơn.


2
Liên kết đến hình đã chết.
Nick Stauner

Câu trả lời:


24
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(...)
       })

văn bản thay thế

Tôi phải nhận xét rằng sự phân tách hoàn hảo xảy ra ở đây, do đó glmchức năng cung cấp cho bạn một cảnh báo. Nhưng điều đó không quan trọng ở đây vì mục đích là để minh họa cách vẽ ranh giới tuyến tính và các quan sát được tô màu theo các hiệp phương sai của chúng.


Tôi hy vọng tôi không bị lỗi thời nếu tôi sử dụng mạng :-)
suncoolsu

2
Tôi cũng hy vọng rằng nếu đây là sự cố CTNH, bạn sẽ không chỉ sao chép dán.
suncoolsu

Cảm ơn. Đây không phải là một câu hỏi CTNH và câu trả lời là hữu ích cho tôi để hiểu mô hình của tôi.
dùng2755

ồ đúng rồi :)
mpiktas

1
Ai đó có thể giải thích cho tôi logic đằng sau con dốc và đánh chặn? (liên quan đến mô hình logistic)
Fernando

22

Muốn giải quyết câu hỏi trong bình luận cho câu trả lời được chấp nhận ở trên từ Fernando: Ai đó có thể giải thích logic đằng sau con dốc và đánh chặn không?

Giả thuyết về hồi quy hậu cần có dạng:

hθ= =g(z)

g(z)z

z= =θ0+θ1x1+θ2x2

Cho rằng chúng tôi đang phân loại từ 0 đến 1, y= =1hθ0,5

θ0+θ1x1+θ2x20

ở trên là ranh giới quyết định và có thể được sắp xếp lại như sau:

x2-θ0θ2+-θ1θ2x1

y= =mx+bmb


2
Giải thích tốt kèm theo câu trả lời ở trên!
Augustin
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.