\begin{Schunk}
\begin{Sinput}
 m6 <- lm(statusquo ~ income * region + sex, data = Chile)
 plotSlopes(m6, plotx = "income", modx = "sex", plotLegend = FALSE, ylim = magRange(Chile$statusquo, c(1,1.3)), llwd = c(2,2), cex = 1.0, lwd = 0.5)
 legend("topleft", c("Female", "Male"), col = c(1, 2), lty = 1:2, bg = "white") 
\end{Sinput}
\end{Schunk}
