I am writing on this forum because the phyloseq github issues page is not being answered and I need help with a very specific formating question. Ultimately, I want to have a grid of graphs like this where the column headers (where asia, europe, and north america are) is body site, sex, and pregnancy status, while my 3 row headers are the alpha diversity metrics (shannon, observed, and Inverse Simpson). I want to use the plot_richness function with phyloseq and create diversity metrics for each of my categories.