\begin{Schunk}
\begin{Sinput}
 table(dat$religion2, dat$religion, dnn = c("religion2", "religion"))
\end{Sinput}
\begin{Soutput}
         religion
religion2 cath jewish muslem other prot roman
   jewish    0     87      0     0    0     0
   muslem    0      0     94     0    0     0
   other     0      0      0   294    0     0
   prot      0      0      0     0  169     0
   cath    177      0      0     0    0   105
\end{Soutput}
\end{Schunk}
