--min_reads and --min_samples parameter in metagenomic_pipeline.py

15 views
Skip to first unread message

Vincent Darbot

unread,
Oct 25, 2021, 3:21:30 AM10/25/21
to picrust-users
Hello !

I am trying different thresholds for the --min_reads parameter in the third step of Picrust2. If I weel understood, the only difference with not putting this parameter is that the OTUs bellows the threshold will be flaged as 'RARE' in stratified output, and there is no differences in the calculation of abundances ?

I think it would be great if these OTUs could be removed for the following analysis, like the '--max-nsti' parameter. But maybe there are reasons why they are not :)

Thank's by advance for clarifications

Gavin Douglas

unread,
Oct 25, 2021, 5:57:43 PM10/25/21
to picrus...@googlegroups.com
Hey there,

If I weel understood, the only difference with not putting this parameter is that the OTUs bellows the threshold will be flaged as 'RARE' in stratified output, and there is no differences in the calculation of abundances ?

Yes - all of the OTUs/ASVs below the specificed thresholds are simply binned together as “RARE” rather than output separately.

I think it would be great if these OTUs could be removed for the following analysis, like the '--max-nsti' parameter. But maybe there are reasons why they are not :) 

Thanks for the suggestion. I think it’s better to keep them in as they aren’t necessarily erroneous sequences (which is the intention of screening based on max NSTI). However, you could remove them prior to running the pathway inference step manually by removing any rows with “RARE”.


Cheers,

Gavin

--
You received this message because you are subscribed to the Google Groups "picrust-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to picrust-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/picrust-users/7d89b2d5-816d-43c3-a454-e805fa3ce381n%40googlegroups.com.

Reply all
Reply to author
Forward
0 new messages