budding <- dotna(read.table("budding.dat",header=T))
budding <- budding[!is.na(budding$locus),]
budding$locus <- factor(budding$locus)
budding$mean <- predict(lm(db~locus,budding))

boxplot(split(budding$db,budding$locus))
tapply(budding$db,budding$locus,stem)