> pdf("E:\\sta4211\\KNNL_TA27_07.pdf")
> 
> shoe <- read.table("http://www.stat.ufl.edu/~winner/sta4211/datasets/CH27TA07.txt",header=F,
+     col.names=c("sales","mkt","advcmp","tprd"))
> 
> attach(shoe)
> 
> a <- max(advcmp)
> b <- max(tprd)
> s <- max(mkt)
> mkt1 <- (advcmp-1)*s + mkt    # Gives each market a distinct level 
> 
> mkt1 <- factor(mkt1)
> advcmp <- factor(advcmp)
> tprd <- factor(tprd)
> 
> par(mfrow=c(1,2))
> interaction.plot(tprd[advcmp==1],mkt[advcmp==1],sales[advcmp==1],lty=1,pch=1)
> interaction.plot(tprd[advcmp==2],mkt[advcmp==2],sales[advcmp==2],lty=1,pch=1)
> 
> shoe.aov1 <- aov(sales ~ advcmp + advcmp/mkt1 + tprd + advcmp:tprd)
> summary(shoe.aov1)
            Df  Sum Sq Mean Sq  F value    Pr(>F)    
advcmp       1  168151  168151 469.7380 2.762e-13 ***
tprd         2   67073   33537  93.6862 1.468e-09 ***
advcmp:mkt1  8 1833681  229210 640.3113 < 2.2e-16 ***
advcmp:tprd  2     391     196   0.5468    0.5892    
Residuals   16    5727     358                       
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
> 
> shoe.aov2 <- aov(sales ~ advcmp*tprd + Error(mkt1))
> summary(shoe.aov2)

Error: mkt1
          Df  Sum Sq Mean Sq F value Pr(>F)
advcmp     1  168151  168151  0.7336 0.4166
Residuals  8 1833681  229210               

Error: Within
            Df Sum Sq Mean Sq F value    Pr(>F)    
tprd         2  67073   33537 93.6862 1.468e-09 ***
advcmp:tprd  2    391     196  0.5468    0.5892    
Residuals   16   5727     358                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
>