How to generate Beta diversity boxplot from phyloseq object?
You could convert the phyloseq object to a dataframe and just plot it using base R or ggplot2. I'm not exactly sure what the bray curtis df will be called but you can access it using the @ notation then convert with data.frame e.g.
bray_for_plot &lt;- data.frame(phylo_obj@bray)