\begin{Schunk}
\begin{Sinput}
 predictOMatic(incedmod1, predVals = list(education = "quantile"), n = 10, interval = "prediction")
\end{Sinput}
\begin{Soutput}
   education       fit        lwr      upr
1      6.380  2880.840 -4150.1994  9911.88
2      7.522  3907.285 -3084.5739 10899.14
3      8.128  4451.965 -2523.8370 11427.77
4      8.766  5025.408 -1937.0716 11987.89
5      9.708  5872.089 -1077.4777 12821.66
6     10.540  6619.902  -324.9942 13564.80
7     11.172  7187.951   242.3737 14133.53
8     12.209  8120.020  1165.4153 15074.63
9     13.620  9388.245  2405.6470 16370.84
10    14.703 10361.660  3345.4139 17377.91
11    15.970 11500.456  4431.6587 18569.25
\end{Soutput}
\end{Schunk}
