\begin{Schunk}
\begin{Sinput}
 Y <- 30 + 2*X + rnorm(100, 0, 50)
 plot(Y ~ X)
 summary(lm(Y ~ X))
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = Y ~ X)

Residuals:
     Min       1Q   Median       3Q      Max 
-104.755  -40.437   -9.674   32.133  152.864 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  28.9183     5.4696   5.287 7.55e-07 ***
X             1.6869     0.5418   3.114  0.00242 ** 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 54.64 on 98 degrees of freedom
Multiple R-squared:  0.09003,	Adjusted R-squared:  0.08075 
F-statistic: 9.696 on 1 and 98 DF,  p-value: 0.002421
\end{Soutput}
\end{Schunk}
