shapiro.test(A$income)
shapiro.test(A$urban)

A1<-lm(urban~income,data=Anscombe)
summary(A1)
scatterplot(urban~income,data=Anscombe,smooth=F)

Call:
lm(formula = urban ~ income, data = Anscombe)

Residuals:
    Min      1Q  Median      3Q     Max
-351.06  -51.07   -4.37   81.75  220.13

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept) 67.04877   91.99292   0.729     0.47
income       0.18524    0.02811   6.590 2.87e-08 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 111.3 on 49 degrees of freedom
Multiple R-squared:  0.4699,    Adjusted R-squared:  0.459
F-statistic: 43.43 on 1 and 49 DF,  p-value: 2.866e-08