\begin{Schunk}
\begin{Sinput}
 newdf <- data.frame(income=c(mean(Chile$income, na.rm = T)),  age = c(55), sex = "F")
 predict(m1, newdata = newdf)
\end{Sinput}
\begin{Soutput}
        1 
0.1957864 
\end{Soutput}
\end{Schunk}
