rock1 <- read.csv("http://www.stat.ufl.edu/~winner/data/tunnel_rock_igneous.csv", header=T) attach(rock1); names(rock1) plot(rock1) rock.mod1 <- lm(CAI ~ QC + UCS + BTS + Brittle1 + Brittle2 + Brittle3) summary(rock.mod1) anova(rock.mod1)