\begin{Schunk}
\begin{Sinput}
 mod2 <- lm(th.bk ~ V043250 + V041109A, data = anes1)
 summary(mod2)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = th.bk ~ V043250 + V041109A, data = anes1)

Residuals:
     Min       1Q   Median       3Q      Max 
-113.174  -42.222   -2.782   42.478  107.164 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)  
(Intercept) -3.08501    4.83135  -0.639   0.5232  
V043250      0.19128    0.09166   2.087   0.0371 *
V041109AF   -7.71339    3.12318  -2.470   0.0137 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 53.77 on 1188 degrees of freedom
  (21 observations deleted due to missingness)
Multiple R-squared:  0.008467,	Adjusted R-squared:  0.006798 
F-statistic: 5.072 on 2 and 1188 DF,  p-value: 0.006404
\end{Soutput}
\end{Schunk}
