This function performs Tukey HSD to a data frame of dimension scores where more
than one corpus have been anaylsed. If using the result of dtag_tbl
,
the data frame needs to be changed to a long format, e.g. by using the
pivot_longer(tbl , contains("Dimension"), names_to = "dimension",
names_transform = list(dimension = as.factor))
, with an extra column to identify
the source corpus.