bactroom <- bacteria[!duplicated(bacteria$code),c("bact","temp","run")]

bactroom$bill <- predict(bacteria.bill)[!duplicated(bacteria$code)]
bactroom$leg  <- predict(bacteria.leg) [!duplicated(bacteria$code)]
bactroom$code <- unique(bacteria$code)
bactroom$trt <-  recode(bactroom$bact,bactroom$temp)

bactroom <- bactroom[order(bactroom$code),]

bactroom.bill <- aov(bill~bact+temp,bactroom)
bactroom.leg <- aov(leg~bact+temp,bactroom)