\begin{Schunk}
\begin{Sinput}
 termplot(mod2, terms=c("occupationf"), partial=T, se=T)
\end{Sinput}
\end{Schunk}
