I am trying to figure out what is the estimated variance (i.e. the estimated "error") of residuals around a fitted line.
> summary(model)
Call:
lm(formula = fecundity ~ Organic)
Residuals:
Min 1Q Median 3Q Max
-2.2909 -1.6439 -0.4606 1.5121 3.7273
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 47.6667 1.4907 31.98 9.97e-10 ***
Organic -8.6788 0.4805 -18.06 9.06e-08 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2.182 on 8 degrees of freedom
Multiple R-squared: 0.9761, Adjusted R-squared: 0.9731
F-statistic: 326.2 on 1 and 8 DF, p-value: 9.063e-08
> anova(model)
Analysis of Variance Table
Response: fecundity
Df Sum Sq Mean Sq F value Pr(>F)
Organic 1 1553.5 1553.50 326.22 9.063e-08 ***
Residuals 8 38.1 4.76
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1