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