#
# This is the insecticide data from Plackett and Hewlett (in my version of the book it is in Table 8.11 on EPage 159)
#
DF = data.frame(
    dose_A=c( 0.05, 0.07, 0, 0.05, 0.07 ),
    dose_B=c( 0, 0, 0.2, 0.2, 0.2 ),
    value=c( 9, 22, 5, 27, 27 )
)

# We want just zero and the highest does level for dose_A:
# 
mask = (DF$dose_A == 0) | (DF$dose_A == max(DF$dose_A))
DF = DF[mask,]

# Convert these dosages to factors:
#
DF$dose_A = as.factor( DF$dose_A )
DF$dose_B = as.factor( DF$dose_B )

fit = aov( value ~ dose_A, data=DF )
sfit = summary(fit)
print( sprintf("Classica1 ANOVA (dose_A ME): prob. type I error: %10.6f", sfit[[1]][["Pr(>F)"]][1]) )