\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="") 
 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)")
 plogisatxcut <- plogis(xcut, 0, 1, log=FALSE)
 text(2, 0.015, pos = 2, paste("Pr(y = 0) = ", formatC(plogisatxcut)))
 mtext(expression(tau), 1, at = xcut)
\end{Sinput}
\end{Schunk}
