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

Residuals:
    Min      1Q  Median      3Q     Max 
-63.072 -10.687  -0.023  11.442  47.739 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 4.707e+01  9.577e-01   49.15   <2e-16 ***
instcst     1.491e-03  9.003e-05   16.56   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 16.85 on 1176 degrees of freedom
  (123 observations deleted due to missingness)
Multiple R-squared:  0.1891,	Adjusted R-squared:  0.1884 
F-statistic: 274.3 on 1 and 1176 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
