\begin{Schunk}
\begin{Sinput}
 newdf <- data.frame(income = mean(Chile$income, na.rm=TRUE), age=median(Chile$age, na.rm = TRUE), sex = "M")
 newdf$fit <- predict(m2, newdata = newdf)
 newdf
\end{Sinput}
\begin{Soutput}
    income age sex        fit
1 33.87586  36   M -0.1040288
\end{Soutput}
\end{Schunk}
