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

Residuals:
     Min       1Q   Median       3Q      Max 
-158.748  -27.161    1.873   27.349  126.605 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) -53.30952   10.34966  -5.151 3.22e-07 ***
V043250       0.02894    0.08842   0.327   0.7435    
V041109AF    -0.38823    2.94737  -0.132   0.8952    
V043116L     -1.53779   10.46541  -0.147   0.8832    
V043116SL    25.50127   10.32322   2.470   0.0137 *  
V043116M     49.84527    9.92408   5.023 6.21e-07 ***
V043116SC    79.64971   10.23575   7.782 2.08e-14 ***
V043116C    103.70832   10.14592  10.222  < 2e-16 ***
V043116EC   111.47834   12.37917   9.005  < 2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 42.8 on 850 degrees of freedom
  (353 observations deleted due to missingness)
Multiple R-squared:  0.4094,	Adjusted R-squared:  0.4039 
F-statistic: 73.67 on 8 and 850 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
