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

Residuals:
     Min       1Q   Median       3Q      Max 
-2.10427 -0.39890  0.05395  0.44156  1.35513 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
gendermale    2.14912    0.16207  13.260  < 2e-16 ***
genderfemale  2.23329    0.16290  13.710  < 2e-16 ***
agency        0.34839    0.06226   5.596 4.24e-08 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.607 on 377 degrees of freedom
Multiple R-squared:  0.9627,	Adjusted R-squared:  0.9624 
F-statistic:  3241 on 3 and 377 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
