# NOT RUN {
library(TukeyC)
data(sorghum)
av <- aov(y ~ r/bl + x,
data=sorghum$dfm)
tk <- TukeyC(av,
which='x',
sig.level=0.05)
summary(tk)
plot(tk)
# }
Run the code above in your browser using DataLab