\begin{Schunk}
\begin{Sinput}
 x <- runif(500, -10,100)
 y <- -0.8 + 0.079 * x
 prob <- 1/(1+ exp(-y))
 yobs <- rbinom(500, p=prob, size=1)
 yobsf <- factor(yobs, labels=c("No","Yes"))
 plot(x, as.numeric(yobsf)-1, xlab="Socio Economic Status", ylab="Will Vote for Republican", axes=F, ylim=c(-0.3,1.3))
 axis(1)
 axis(2, at=c(0,1), labels=c("No","Yes"))
\end{Sinput}
\end{Schunk}
