I know how to run the basic skbio permanova code in python:
import skbio as skb
skperm=skb.stats.distance.permanova(dm, host_md_test, column='Host1', permutations=999)
But I only get stats on the overall test. Is there a way to show the pairwise comparisons like in qiime2?
I am trying to keep this in python because I want to run a for loop on many permanovas.
My other option will be to go into R and use adonis's pairwise option.
Any info would be appreciated!!