\begin{Schunk}
\begin{Sinput}
 mod1 <- lm (gradrate ~ schtype, data = dat)
 summary(mod1)
\end{Sinput}
\begin{Soutput}
Call:
lm(formula = gradrate ~ schtype, data = dat)

Residuals:
   Min     1Q Median     3Q    Max 
-51.10 -11.18  -0.10  11.90  49.82 

Coefficients:
               Estimate Std. Error t value Pr(>|t|)    
(Intercept)      50.182      0.826   60.75   <2e-16 ***
schtypeprivate   15.918      1.033   15.41   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 17.21 on 1201 degrees of freedom
  (98 observations deleted due to missingness)
Multiple R-squared:  0.1651,	Adjusted R-squared:  0.1644 
F-statistic: 237.4 on 1 and 1201 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
