\begin{Schunk}
\begin{Soutput}
Call:
lm(formula = Income ~ IQ * PSES, data = dat)

Residuals:
    Min      1Q  Median      3Q     Max 
-5.8539 -1.3270 -0.2993  1.3951  5.8935 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 3.084322   0.226037   13.64   <2e-16 ***
IQ          0.989595   0.021580   45.86   <2e-16 ***
PSES        0.990908   0.024113   41.09   <2e-16 ***
IQ:PSES     1.001278   0.002486  402.84   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 2.136 on 96 degrees of freedom
Multiple R-squared:  0.9996,	Adjusted R-squared:  0.9996 
F-statistic: 7.373e+04 on 3 and 96 DF,  p-value: < 2.2e-16
\end{Soutput}
\end{Schunk}
