lager <- read.csv("http://www.stat.ufl.edu/~winner/data/lager_antioxidant_reg.csv", header=T) attach(lager); names(lager) # print(cbind(tpc,dsa)) lager.mod1 <- lm(dsa ~ tpc) summary(lager.mod1) anova(lager.mod1) confint(lager.mod1) plot(tpc,dsa, main="DPPH versus TPC") abline(lager.mod1)