nfl2007 <- read.csv("http://www.stat.ufl.edu/~winner/data/nfl2007teamscore.csv", header=T) attach(nfl2007); names(nfl2007) nfl.mod1 <- lm(obsScore ~ vegasScore + homeTeam + vS_hT) summary(nfl.mod1) anova(nfl.mod1) drop1(nfl.mod1,test="F") Y.test <- obsScore - vegasScore nfl.mod2 <- lm(obsScore ~ -1, offset=vegasScore) summary(nfl.mod2) anova(nfl.mod2,nfl.mod1)