# 1. OBS Observation index # 2. ATHLETE =1 if athlete # 3. BLACK =1 if black # 4. COLGPA GPA after fall semester # 5. FEMALE =1 if female # 6. HSIZE size graduating class, 100s # 7. HSIZESQ hsize^2 # 8. HSPERC 100*(hsrank/hssize) # 9. HSRANK rank in graduating class # 10.SAT combined SAT score # 11.TOTHRS total hours through fall semest # 12.VERBMATH verbal/math SAT score # 13.WHITE =1 if white data = read.table("http://hedibert.org/wp-content/uploads/2014/04/gpa2.txt",header=TRUE) n = nrow(data) attach(data) par(mfrow=c(2,4)) plot(ATHLETE,COLGPA) abline(lm(COLGPA~ATHLETE),col=2) title(paste("R2=",round(summary(lm(COLGPA~ATHLETE))$r.sq,4),sep="")) plot(WHITE,COLGPA) abline(lm(COLGPA~WHITE),col=2) title(paste("R2=",round(summary(lm(COLGPA~WHITE))$r.sq,4),sep="")) plot(FEMALE,COLGPA) abline(lm(COLGPA~FEMALE),col=2) title(paste("R2=",round(summary(lm(COLGPA~FEMALE))$r.sq,4),sep="")) plot(HSIZE,COLGPA) abline(lm(COLGPA~HSIZE),col=2) title(paste("R2=",round(summary(lm(COLGPA~HSIZE))$r.sq,4),sep="")) plot(HSRANK,COLGPA) abline(lm(COLGPA~HSRANK),col=2) title(paste("R2=",round(summary(lm(COLGPA~HSRANK))$r.sq,4),sep="")) plot(HSPERC,COLGPA) abline(lm(COLGPA~HSPERC),col=2) title(paste("R2=",round(summary(lm(COLGPA~HSPERC))$r.sq,4),sep="")) plot(TOTHRS,COLGPA) abline(lm(COLGPA~TOTHRS),col=2) title(paste("R2=",round(summary(lm(COLGPA~TOTHRS))$r.sq,4),sep="")) plot(SAT,COLGPA) abline(lm(COLGPA~SAT),col=2) title(paste("R2=",round(summary(lm(COLGPA~SAT))$r.sq,4),sep=""))