hotel.energy <- read.csv("http://www.stat.ufl.edu/~winner/data/hotel_energy.csv") attach(hotel.energy); names(hotel.energy) enrgcons <- enrgcons/1000000 he.mod1 <- lm(enrgcons ~ area + age + effrooms) summary(he.mod1) anova(he.mod1) he.mod2 <- lm(enrgcons ~ area) summary(he.mod2) anova(he.mod2) anova(he.mod2, he.mod1)