\begin{Schunk}
\begin{Sinput}
 bp1 <- predict(mod1, interval = "confidence")
 head(bp1)
\end{Sinput}
\begin{Soutput}
       fit      lwr      upr
1 66.10013 64.88272 67.31754
3 50.18203 48.56150 51.80256
5 50.18203 48.56150 51.80256
6 66.10013 64.88272 67.31754
7 50.18203 48.56150 51.80256
8 50.18203 48.56150 51.80256
\end{Soutput}
\end{Schunk}
