\begin{Schunk}
\begin{Sinput}
 dat$occupationf2[dat$occupationf2 %in% occlev[5]] <- occlev[1] 
 levels(dat$occupationf2)[1] <- "manag/prof"
 dat$occupationf2 <- dat$occupationf2[, drop=T]
\end{Sinput}
\end{Schunk}
