\begin{Schunk}
\begin{Sinput}
 modlm1 <- lm(posAffect ~ -1 + gender/agency, data = dat)
 summary(modlm1, signif.stars = FALSE)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = posAffect ~ -1 + gender/agency, data = dat)

Residuals:
     Min       1Q   Median       3Q      Max 
-2.10608 -0.40401  0.02606  0.44460  1.32508 

Coefficients:
                    Estimate Std. Error t value Pr(>|t|)    
gendermale           1.88976    0.22056   8.568 2.75e-16 ***
genderfemale         2.51356    0.22955  10.950  < 2e-16 ***
gendermale:agency    0.45182    0.08624   5.239 2.69e-07 ***
genderfemale:agency  0.23701    0.08950   2.648  0.00843 ** 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.6054 on 376 degrees of freedom
Multiple R-squared:  0.963,	Adjusted R-squared:  0.9626 
F-statistic:  2444 on 4 and 376 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
