bacteria <- read.table( "bacteria.dat", header=T )

# model fits for bill and leg
bacteria.bill <- aov( bill ~ temp * bact * inoc, bacteria )
bacteria.leg <- aov( leg ~ temp * bact * inoc, bacteria )

bacteria.rcbd <- aov( bill ~ inoc + bact * temp, bacteria )