\begin{Schunk}
\begin{Sinput}
 sexlev <- levels(Chile$sex)
 newdf <- data.frame(income=c(mean(Chile$income, na.rm = T)),  age = c(55), sex = sexlev[1])
 newdf$fit <- predict(m1, newdata = newdf)
\end{Sinput}
\end{Schunk}
