\begin{Schunk}
\begin{Sinput}
 plot(posAffect ~ agency, data = dat, cex = 0.6, col = "gray60")
 lines(reg.mod.1.pred2[ , "fit"] ~ nd$agency, lty = 1)
 lines(reg.mod.1.pred2[ , "lwr"] ~ nd$agency, lty = 3, lwd = 2, col = "blue")
 lines(reg.mod.1.pred2[ , "upr"] ~ nd$agency, lty = 3, lwd = 2, col = "blue")
\end{Sinput}
\end{Schunk}
