\begin{Schunk}
\begin{Sinput}
 dat$religion2 <- combineLevels(dat$religion, c("cath", "roman"), "cath")
\end{Sinput}
\begin{Soutput}
The original levels cath jewish muslem other prot roman 
have been replaced by jewish muslem other prot cath 
\end{Soutput}
\end{Schunk}
