text: user7427: I am working with profiled metagenomic taxonomy abundance data. I want to generate beta diversity (bray-curtis) boxplot from a phyloseq object where two groups (control and test) will be shown. Something like this figure:
  
  How can I make it? Can anyone share any code or tutorial?
  Thanks,
  dpc
  
  user220: 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 <- data.frame(phylo_obj@bray)
  
  
original_text: 7427: <p>I am working with profiled metagenomic taxonomy abundance data. I want to generate beta diversity (bray-curtis) boxplot from a phyloseq object where two groups (control and test) will be shown. Something like this figure:</p>
  <p><a href="https://i.stack.imgur.com/pgauk.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/pgauk.png" alt="enter image description here" /></a></p>
  <p>How can I make it? Can anyone share any code or tutorial?</p>
  <p>Thanks,
  dpc</p>
   220: <p>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.</p>
  <pre><code>bray_for_plot &lt;- data.frame(phylo_obj@bray)
  </code></pre>
  