data(BrendonSmall)
model = aov(Calories ~ Instructor, data=BrendonSmall)
TUK = TukeyHSD(model, "Instructor", ordered = TRUE)
### Convert the TukeyHSD output to a standard data frame
TUK = as.data.frame(TUK$Instructor)
names(TUK) = gsub(" ", ".", names(TUK))
HSD = data.frame(Comparison=row.names(TUK),
diff=TUK$diff, lwr=TUK$lwr, lwr=TUK$lwr, p.adj=TUK$p.adj)
HSD
cldList(p.adj ~ Comparison, data = HSD,
threshold = 0.05,
remove.space=FALSE)
Run the code above in your browser using DataLab