\begin{Schunk}
\begin{Sinput}
 mod2 <- lm(y ~ xnum1 * z, data = dat)
 mod2sum <- summary(mod2)
 coef(mod2sum)
\end{Sinput}
\begin{Soutput}
               Estimate Std. Error    t value   Pr(>|t|)
(Intercept) -0.01095327  0.1581058 -0.0692781 0.94491862
xnum1        0.28981053  0.1681968  1.7230442 0.08824016
zB           0.27670253  0.2417953  1.1443668 0.25544018
zC           0.37603376  0.2569725  1.4633227 0.14678733
zD          -0.27113912  0.2088686 -1.2981327 0.19748670
xnum1:zB     0.39229326  0.2801136  1.4004794 0.16473350
xnum1:zC     0.11537877  0.3911294  0.2949887 0.76866715
xnum1:zD     0.27251832  0.2047313  1.3311026 0.18644453
\end{Soutput}
\end{Schunk}
