\begin{Schunk}
\begin{Sinput}
 mod4pmo <- predictOMatic(mod2, predVals = list(z = levels(dat$z), xnum1 = range(dat$xnum1)))
 mod4pmo
\end{Sinput}
\begin{Soutput}
      xnum1 z        fit
1 -3.128437 A -0.9176073
2 -3.128437 B -1.8681696
3 -3.128437 C -0.9025288
4 -3.128437 D -2.0413029
5  2.168475 A  0.6174936
6  2.168475 B  1.7448742
7  2.168475 C  1.2437233
8  2.168475 D  0.9373036
\end{Soutput}
\end{Schunk}
