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

Residuals:
     Min       1Q   Median       3Q      Max 
-2.11727 -0.38900  0.04804  0.44363  1.38301 

Coefficients:
                         Estimate Std. Error t value
(Intercept)               2.47790    0.47151   5.255
agency                    0.15378    0.18164   0.847
ethnicityBlack           -0.59417    0.82482  -0.720
ethnicityHispanic        -0.20871    0.59406  -0.351
ethnicityWhite           -0.43142    0.51081  -0.845
genderfemale              0.09997    0.06343   1.576
agency:ethnicityBlack     0.29965    0.33543   0.893
agency:ethnicityHispanic  0.13503    0.22878   0.590
agency:ethnicityWhite     0.24453    0.19813   1.234
                         Pr(>|t|)    
(Intercept)               2.5e-07 ***
agency                      0.398    
ethnicityBlack              0.472    
ethnicityHispanic           0.726    
ethnicityWhite              0.399    
genderfemale                0.116    
agency:ethnicityBlack       0.372    
agency:ethnicityHispanic    0.555    
agency:ethnicityWhite       0.218    
---
Signif. codes:  
0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.6078 on 371 degrees of freedom
Multiple R-squared:  0.0933,	Adjusted R-squared:  0.07375 
F-statistic: 4.772 on 8 and 371 DF,  p-value: 1.347e-05
\end{Soutput}
\end{Schunk}
