# Trying to visualise q2-picrust2 data on R

**URL:** https://forum.qiime2.org/t/trying-to-visualise-q2-picrust2-data-on-r/25119
**Category:** PICRUSt2
**Tags:** qiime2r, picrust, r
**Created:** [January 3, 2023, 9:28am UTC](https://forum.qiime2.org/t/trying-to-visualise-q2-picrust2-data-on-r/25119 "2023-01-03T09:28:32Z")
**Posts on this page:** 3
**Page:** 1

<div class="post-metadata">

### Author: ![Johanndb](https://forum.qiime2.org/user_avatar/forum.qiime2.org/johanndb/32/17117_2.png) [@Johanndb](https://forum.qiime2.org/u/Johanndb)
#### Post date: [January 3, 2023, 9:28am UTC](https://forum.qiime2.org/t/trying-to-visualise-q2-picrust2-data-on-r/25119/1 "2023-01-03T09:28:32Z")

</div>

Hi all,

I successfully ran my data through q2-picrust2 and received the appropriate output. However, this out means nothing without a way to process the data. I am wondering if anyone has found a way to visualise their Pircust2 data through R (specifically a heat map or bar plot). Could this be accomplished using phyloseq or similar packages? How would I read the Picrust2 data into phyloseq? Lastly, (and independently of Picrust2) how would you go about creating a Venn diagram of shared sequences between all samples?

An example could be found in this paper: [Biodiversity and Activity of the Gut Microbiota across the Life History of the Insect Herbivore Spodoptera littoralis | Scientific Reports](https://www.nature.com/articles/srep29505) / [Insect-based diet, a promising nutritional source, modulates gut microbiota composition and SCFAs production in laying hens | Scientific Reports](https://www.nature.com/articles/s41598-017-16560-6)

Thank you for all the help and I hope everyone had a great festive season!

Kind regards,

Johann

---

<div class="post-metadata">

### Author: ![gregcaporaso](https://forum.qiime2.org/user_avatar/forum.qiime2.org/gregcaporaso/32/17769_2.png) [@gregcaporaso](https://forum.qiime2.org/u/gregcaporaso)
#### Post date: [January 3, 2023, 10:38pm UTC](https://forum.qiime2.org/t/trying-to-visualise-q2-picrust2-data-on-r/25119/2 "2023-01-03T22:38:39Z")

</div>

Hi @Johanndb,  
I suspect that you have a `FeatureTable[Frequency]` as output that you're interested in working with. There are a lot of tools in QIIME 2 for downstream analysis of these - e.g., `core-metrics`. Most actions that can take a QIIME 2 `FeatureTable[Frequency]` that was generated for example with q2-dada2 should work with the result of q2-picrust2 (the exception would be actions that require another input that you might not have, like a phylogenetic tree).

If you specifically want to use R for your next steps, I recommend looking into using [qiime2R](https://github.com/jbisanz/qiime2R/) - that's a third party tool that is very popular for working with QIIME 2 data in R.

As for Venn diagrams, we don't have a direct way to generate those in QIIME 2, but [this post](https://forum.qiime2.org/t/creating-venn-diagram-to-show-shared-bacterial-taxa-among-two-sample-sets/23522) provides a suggestion.

In general, I would recommend doing a test such as [`beta-group-significance`](https://docs.qiime2.org/2022.11/plugins/available/diversity/beta-group-significance/), which will tell you if your groups are significantly dissimilar from each other in their composition as that will have a clearer interpretation. You could also pair that with a Venn diagram to aid in interpretation of the diagram (Jaccard distance might be a particularly relevant metric there, as it's effectively the sum of features that fall outside the intersection of two samples in a Venn diagram of their features).

Hope this helps!

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [March 30, 2023, 10:28am UTC](https://forum.qiime2.org/t/trying-to-visualise-q2-picrust2-data-on-r/25119/4 "2023-03-30T10:28:22Z")

</div>


