I'm trying to generate boxplot from Beta diversity boxplot between two conditions from a phyloseq object.
I have a Morisita-horn distance generated by vegdist fonction :
vegdist<-vegdist(t(otu_table(physeq_genus_3a)), method="horn")#matrice distance Morisita-horn
vegdist<-as.matrix(vegdist)
Now I have a matrix and I would like to somewhat create a mean distance boxplot of the Morisita-horn index winthin two condition like the graph below.
I wanted to do a simple boxplot with ggplot but beacause it is a matrix and not a dataframe, I can't do it.
Any help or suggestions would be more than welcome.