¿Cómo trazar el límite de decisión en R para el modelo de regresión logística?

Hice un modelo de regresión logística usando glm en R. Tengo dos variables independientes. ¿Cómo puedo trazar el límite de decisión de mi modelo en el diagrama de dispersión de las dos variables? Por ejemplo, ¿cómo puedo trazar una figura como: http://onlinecourses.science.psu.edu/stat557/node/55

Gracias.

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

Debo señalar que aquí se produce una separación perfecta, por lo tanto, la glmfunción le da una advertencia. Pero eso no es importante aquí ya que el propósito es ilustrar cómo dibujar el límite lineal y las observaciones coloreadas de acuerdo con sus covariables.

Quería abordar la pregunta en comentario a la respuesta aceptada de Fernando: ¿Alguien puede explicar la lógica detrás de la pendiente e interceptar?

La hipótesis de la regresión logística toma la forma de:

$$h_{\theta} = g(z)$$

$g(z)$$z$

$$z = \theta_{0} + \theta_{1}x_{1} + \theta_{2}x_{2}$$

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

$$\theta_{0} + \theta_{1}x_{1} + \theta_{2}x_{2} \geq 0$$

lo anterior es el límite de decisión y se puede reorganizar como:

$$x_{2} \geq \frac{-\theta_{0}}{\theta_{2}} + \frac{-\theta_{1}}{\theta_{2}}x_{1}$$

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