This function is merely a thin wrapper around stats::TukeyHSD with tidying done
by broom::tidy and optional formatting via pixiedust::sprinkle.
Its input is not a aov model like in the original function, but instead the
aov model is fit internally based on the arguments given. This is meant to
enable a consistent usage between the tadaa_pairwise-functions.