tukmod <- data(tukmod)
tukmod$tau <- factor(tukmod$tau)
tukmod$beta <- factor(tukmod$beta)
tukmod.fit <- aov(resp~beta+tau,tukmod)