\begin{Schunk}
\begin{Sinput}
 mod1 <- lm(wage ~ occupation, data=dat) 
\end{Sinput}
\end{Schunk}
