\begin{Schunk}
\begin{Sinput}
 mod1 <- lm(PosAFF1 ~ Agency1 + Intrin1, data = dat)
 summary(mod1)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = PosAFF1 ~ Agency1 + Intrin1, data = dat)

Residuals:
     Min       1Q   Median       3Q      Max 
-2.17693 -0.35887  0.05093  0.55123  1.12382 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  2.05775    0.18125  11.353  < 2e-16 ***
Agency1      0.18594    0.06581   2.825  0.00497 ** 
Intrin1      0.20778    0.04414   4.707 3.54e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.6363 on 377 degrees of freedom
Multiple R-squared:  0.09673,	Adjusted R-squared:  0.09194 
F-statistic: 20.19 on 2 and 377 DF,  p-value: 4.693e-09
\end{Soutput}
\end{Schunk}
