\begin{Schunk}
\begin{Sinput}
 m6 <- lm(th.bush.kerry ~ V043250 * V041109A, data=mydta1)
 plotSlopes(m6, plotx = "V043250", modx = "V041109A", plotPoints = FALSE, plotLegend = FALSE,  xlab = "Age", ylab = "Bush-Kerry Thermometer", ylim = c(-100,150), col = c("black","gray60"), llwd=c(2,3))
 legend("topleft", legend = levels(mydta1$V041109A), col=c("black","gray60"), lty=c(1,2), lwd=c(2,3))
\end{Sinput}
\end{Schunk}
