\begin{Schunk}
\begin{Sinput}
 reg.mod.1.pm <- predictOMatic(reg.mod.1, predVals = c("agency"), n = 5)
 reg.mod.1.pm
\end{Sinput}
\begin{Soutput}
    agency      fit
1 1.000000 2.537395
2 2.175525 2.947789
3 2.499983 3.061062
4 2.832975 3.177315
5 4.000000 3.584741
\end{Soutput}
\end{Schunk}
