\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 
-1.0386 -0.5489 -0.0799  0.2319  4.7082 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.965354   0.655034   1.474  0.14379    
x1           0.029317   0.010552   2.778  0.00656 ** 
x2          -0.020878   0.004155  -5.025  2.3e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.949 on 97 degrees of freedom
Multiple R-squared:  0.2388,	Adjusted R-squared:  0.2231 
F-statistic: 15.21 on 2 and 97 DF,  p-value: 1.793e-06
\end{Soutput}
\end{Schunk}
