data.source<-"http://www.math.mcgill.ca/dstephens/Regression/Data/2-1-RocketProp.csv"

RocketProp<-read.csv(file=data.source)

names(RocketProp)<-c('i','Strength','Age')

x<-RocketProp$Age
y<-RocketProp$Strength
n<-length(x)
summary(fit.RP)

anova(fit.RP)

anova(fit.RP)[,'Df']
anova(fit.RP)[,'Sum Sq']

sum(anova(fit.RP)[,'Sum Sq'])-(n-1)*var(y)