R data visualisation multiple boxplots showing multiple pairwise comparisons

Question

I would like to display the results of non-parametric pairwise comparisons using base R boxplots. I have multiple groups (species) and multiple variables (rates).

I w