sink(file='Estimation-Output.txt') 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) ############################################ #Fit the simple linear regression using lm fit.RP<-lm(y~x) #Summary of inference summary(fit.RP) sink()