kenton <- read.table("http://www.stat.ufl.edu/~winner/sta4211/datasets/CH16TA04.txt",header=F, col.names=c("package","store","cases","X1","X2","X3")) attach(kenton) kenton.reg1 <- lm(cases ~ X1+X2+X3) summary(kenton.reg1) anova(kenton.reg1) X1w <- X1 X2w <- X2 X3w <- c(rep(0,10),rep(1,4),rep(-0.8,5)) kenton.reg2 <- lm(cases ~ X1w+X2w+X3w) summary(kenton.reg2) anova(kenton.reg2) X1cm <- c(rep(1,5),rep(0,5),rep(0,4),rep(0,5)) X2cm <- c(rep(0,5),rep(1,5),rep(0,4),rep(0,5)) X3cm <- c(rep(0,5),rep(0,5),rep(1,4),rep(0,5)) X4cm <- c(rep(0,5),rep(0,5),rep(0,4),rep(1,5)) kenton.reg3 <- lm(cases ~ X1cm+X2cm+X3cm+X4cm-1) summary(kenton.reg3) anova(kenton.reg3) X1cmeq <- rep(1,19) kenton.reg4 <- lm(cases ~ X1cmeq-1) summary(kenton.reg4) anova(kenton.reg4) anova(kenton.reg4,kenton.reg3)