pdf("readlpga.pdf")

lpga1 <- read.fwf("lpga1.dat", width=c(34,4,8,8,8,8,8,8,8), 
col.names=c("golfer","rounds","drivedist","frwypct","grnregpct","aveputts","avesandshot",
"sandsvpct","przperrnd"))

rounds <- matrix(lpga1$rounds,ncol=1)
drivedist <- matrix(lpga1$drivedist,ncol=1)
frwypct <- matrix(lpga1$frwypct,ncol=1)
grnregpct <- matrix(lpga1$grnregpct,ncol=1)
aveputts <- matrix(lpga1$aveputts,ncol=1)
avesandshot <- matrix(lpga1$avesandshot,ncol=1)
sandsvpct <- matrix(lpga1$sandsvpct,ncol=1)
logprzperrnd <- matrix(log(lpga1$przperrnd),ncol=1)

lpga1a <- data.frame(rounds,drivedist,frwypct,grnregpct,aveputts,avesandshot,sandsvpct,logprzperrnd)

attach(lpga1a)

mean(rounds)
mean(drivedist)
mean(frwypct)
mean(grnregpct)
mean(aveputts)
mean(avesandshot)
mean(sandsvpct)
mean(logprzperrnd)


library(leaps)

allpossreg <- regsubsets(logprzperrnd~rounds+drivedist+frwypct+grnregpct+aveputts+avesandshot+sandsvpct,nbest=5,data=lpga1)
aprout <- summary(allpossreg)


with(aprout,round(cbind(which,rsq,adjr2,cp,bic),3))



dev.off()