THSD <- TukeyHSD(anova) DT::datatable(THSD$Type) %>% formatRound(columns=colnames(THSD$Type),digits=3)