\begin{Schunk}
\begin{Sinput}
 reg.mod.12 <- lm(posAffect ~ agency * intMotiv * extMotiv, data = dat)
 summary(reg.mod.12)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = posAffect ~ agency * intMotiv * extMotiv, data = dat)

Residuals:
     Min       1Q   Median       3Q      Max 
-1.88373 -0.36143  0.03298  0.40975  1.14755 

Coefficients:
                         Estimate Std. Error t value
(Intercept)                0.2944     2.1812   0.135
agency                     1.0015     0.8970   1.117
intMotiv                   0.9717     0.6850   1.419
extMotiv                   1.1748     1.4165   0.829
agency:intMotiv           -0.3254     0.2685  -1.212
agency:extMotiv           -0.5861     0.5547  -1.057
intMotiv:extMotiv         -0.5460     0.4473  -1.221
agency:intMotiv:extMotiv   0.2386     0.1680   1.420
                         Pr(>|t|)
(Intercept)                 0.893
agency                      0.265
intMotiv                    0.157
extMotiv                    0.407
agency:intMotiv             0.226
agency:extMotiv             0.291
intMotiv:extMotiv           0.223
agency:intMotiv:extMotiv    0.156

Residual standard error: 0.5861 on 372 degrees of freedom
Multiple R-squared:  0.1548,	Adjusted R-squared:  0.1389 
F-statistic: 9.735 on 7 and 372 DF,  p-value: 3.806e-11
\end{Soutput}
\end{Schunk}
