\begin{Schunk}
\begin{Sinput}
 lm1 <- lm(y ~ x1 + x2, data=dat)
 summary(lm1)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = y ~ x1 + x2, data = dat)

Residuals:
       Min         1Q     Median         3Q        Max 
-133325385  -81548444  -29144773   26306378 1114657874 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept) -366588914   85770267  -4.274 4.53e-05 ***
x1            19986500    3681228   5.429 4.25e-07 ***
x2              621443    4685999   0.133    0.895    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 160200000 on 96 degrees of freedom
  (1 observation deleted due to missingness)
Multiple R-squared:  0.2384,	Adjusted R-squared:  0.2225 
F-statistic: 15.03 on 2 and 96 DF,  p-value: 2.102e-06
\end{Soutput}
\end{Schunk}
