\begin{Schunk}
\begin{Sinput}
 ## Data must be subdivided by groups
 ps71 <- plotSlopes(reg.mod.7, plotx = "agency", modx = "ethnicity", interval = "confidence")
 ps71$newdata$posAffect <- ps71$newdata$fit
 p1 <- ggplot(dat, aes(x = agency, y = posAffect)) + geom_point(shape = 1, alpha = 0.5) + 
  facet_wrap( ~ ethnicity, ncol = 2) + geom_line(data = ps71$newdata, color = "blue") +
  geom_ribbon(data = ps71$newdata, aes(ymin = lwr, ymax = upr), fill = "pink", alpha = 0.5) +
  theme_bw() + 
  theme(strip.background = element_rect(color="darkgoldenrod4", fill="lightgoldenrod"))
 p1
\end{Sinput}
\end{Schunk}
