Hsu's multiple comparison for the best method

R

dht <- glht(aov.model, linfct = mcp( site = "Dunnett" ), + alternative = "two.sided")
summary(dht)