# BearsReg Example -- bearsreg.r source code # Read data from input file. raw = read.fwf(file="bears85.txt", widths=c(3, 20, 3, 1, 1, 3, 4, 3, 3, 14), header=FALSE, skip=2, col.names=c("JerseyNumber", "Name", "Pos", "HtFt", "dash", "HtIn", "WtLbs", "Age", "YrsInPros", "College")) # Convert height and weight to metric. height = (raw$HtFt * 12 + raw$HtIn) * 0.0254 weight = raw$WtLbs * 0.45359237 # Create data frame. bears = data.frame(Height=height, Weight=weight) # Print bears data frame. print("Bears Dataset") print(bears) # Obtain simple regression model information. srm = lm(Weight ~ Height, data=bears) print("Regression Summary of Bears Dataset:") print(summary(srm)) # Get residuals and predicted values. r = residuals(srm) p = fitted(srm) pdf("bearsreg.pdf") # Create residual plot. plot(p, r, main="Residual Plot for Bears Dataset", xlab="Predicted Values", ylab="Residuals") # Create normal plot. qqnorm(r, main="Normal Plot for Bears Dataset") dev.off( )