Add a common legend to a function that generates a variable number of plots in one figure

Question

I have a function that generates two matrices from a multinomial distribution, it then produces a paired bar chart for each row, the pairs are columns within the row between