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
कार्य आपको चेतावनी देता है। लेकिन यहाँ यह महत्वपूर्ण नहीं है क्योंकि उद्देश्य यह बताना है कि रेखीय सीमा और उनके कोवरिएट्स के अनुसार रंगीन टिप्पणियों को कैसे आकर्षित किया जाए।