\begin{Schunk}
\begin{Sinput}
 x <- seq (-7, 7, by=0.05)
 logis10 <- dlogis(x, location = 0, scale = 1, log = FALSE)
 plot(x, logis10, type="l", ylim=c(0,0.35), xlab=expression(epsilon[i]), ylab="Probability Density", main="") 
 legend("topleft", c(expression(paste("Logistic(", location==0, " , ",
                                       scale == 1, ")"))), lty=1:2, col=1:2, cex = 0.8) 
 title("Logistic variable")
 xcut <- 2
 xlt <- x < xcut
 xmax <- max(x[xlt])
 xseq <- c(x[xlt], rev(x[xlt])) 
 yseq <- c(rep(0, sum(xlt)), rev(logis10[xlt])) 
 polygon(xseq, yseq, col = "gray80")
 text(2.05, 0.015, pos = 4, "Pr(y = 1)")
 text(-1.5, 0.015, pos = 2, "Pr(y = 0)")
\end{Sinput}
\end{Schunk}
