pairwise_stats performs a provided function pairwise between
all combinations of groups. The first two arguments of the function passed
to pairwise_stats must accept vectors of values as inputs. These vectors
should correspond to the values for group A and group B, respectively.
The function takes as input a data.frame or tibble, the column names of
grouping variables, the column name for a variable of interest, and a
function.