# read data and fit model
bactroom <- read.table( "bactroom.dat", header=T )

bactroom.full <- aov( bill ~ temp * bact, bactroom )

bactroom.bill <- aov( bill ~ temp * bact * inoc, bactroom )
bactroom.badd <- aov( bill ~ temp + bact + inoc, bactroom )