\begin{Schunk}
\begin{Sinput}
 p1 <- p1 + geom_ribbon(data = ps31$newdata, 
              aes(ymin = lwr, ymax = upr), 
  			fill = "pink", alpha = 0.5)
 p1
\end{Sinput}
\end{Schunk}
