Lojistik regresyon modeli için R'de karar sınırını nasıl çizebilirim?


Yanıtlar:


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

alternatif metin

Burada mükemmel bir ayrılmanın meydana geldiğini belirtmeliyim, bu nedenle glmişlev size bir uyarı verir. Ancak burada önemli değil çünkü amaç, doğrusal sınırın nasıl çizileceğini ve ortak değişkenlerine göre renklendirilmiş gözlemleri göstermektir.


Kafes kullanırsam eski moda değilim umarım :-)
suncoolsu 13:11

2
Umarım bu bir HW problemi ise, sadece yapıştırmayı kopyalamazsınız.
suncoolsu

Teşekkürler. Bu bir HW sorusu değildir ve cevap modelimi anlamamda yardımcı olur.
user2755

oh evet sensin :)
mpiktas

1
Birisi bana eğimin ve kesişmenin arkasındaki mantığı açıklayabilir mi? (lojistik model ile ilgili)
Fernando

22

Fernando'dan yukarıdaki kabul edilen cevaba yorum olarak soruyu ele almak istedim: Birisi eğimin arkasındaki mantığı açıklayabilir ve kesişebilir mi?

Lojistik regresyon hipotezi şu şekildedir:

hθ=g(z)

burada, sigmoid fonksiyonudur ve formundadır:g(z)z

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

0 ile 1 arasında sınıflandırdığımız göz önüne alındığında, sigmoid işlevi verilen olduğunda durumlarda doğrudur:y=1hθ0.5

θ0+θ1x1+θ2x20

yukarıdaki karar sınırıdır ve aşağıdaki gibi yeniden düzenlenebilir:

x2θ0θ2+θ1θ2x1

Bu, şeklinde bir denklemdir ve o zaman ve neden kabul edilen cevapta oldukları gibi hesaplandığını görebilirsiniz.m by=mx+bmb


2
Yukarıdaki cevaba eşlik eden iyi bir açıklama!
Augustin
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.