\begin{Schunk}
\begin{Sinput}
 mod1 <- lm(y2 ~ x1*x3, data = dat)
 ps1 <- plotSlopes(mod1, plotx = "x1", modx = "x3")
\end{Sinput}
\end{Schunk}
