\begin{Schunk}
\begin{Sinput}
  dat$grpfac4 <- with(dat,  factor(grpfac2, levels=c("Eastern","Western","Midwest")))
  with(dat, table(grpfac4, grpfac2))
\end{Sinput}
\begin{Soutput}
         grpfac2
grpfac4   Western Midwest Eastern
  Eastern       0       0       6
  Western       7       0       0
  Midwest       0       6       0
\end{Soutput}
\end{Schunk}
