addmod <- read.table("tukadd.dat",header=T)
addmod$tau <- factor(addmod$tau)
addmod$beta <- factor(addmod$beta)
addmod.fit <- aov(resp~beta+tau,addmod)