Dear qiime2 forum,
I am planning to use the q2-fmt plugin on my data, but I would like to combine it with q2-boots so I can perform Rarefaction if possible.
My thought was that I could take the inputs from q2-boots and sub them into q2-fmt, and that I could create an "average feature table" from the multiple feature tables created through rarefaction in q2-boots by using q2-feature-table merge by using the average
as the overlap method.
However, as I recognize the tool was not designed for this I am concerned that there are plugin features within it that would not be compatible, or even worse that using the q2-boots outputs instead could jeopardize the outputs of q2-fmt. As such, I would greatly appreciate your input on:
a) the feasibility of this (would this be possible, is this not recommended?) and
b) if there is a better way of combining the two than I have outlined below:
Specifically, I would like to use the following plugin features:
-
Distance to Baseline / Donor - input seems to be a diversity metric
- Would it be possible to simply input a diversity metric generated from q2 boots instead?
-
Proportional Engraftment of Donor Features (PEDF) - input is rarefied feature table
- Would it be possible to input the merged "average rarefaction feature table" instead?
-
Permutation Test of PEDF - input is rarefied table but also requires information on sampling depth
- Would it be possible to input the merged "average rarefaction feature table" instead?
- Out of curiosity, why is sampling depth information used for in this function?
-
Feature Engraftment input is taxonomy collapsed, relative frequency feature table
- What is the advantage of using relative frequency feature table rather than a Rarefaction corrected feature table? For consistency, would it be possible to use the "average rarefaction feature table" instead?
Thank you so much in advance for your help!!
kindest regards,
Zoë
p.s. upon thinking about this more, I wonder if the best way for PEDF would be to create it like how q2-boots does -- by calculating a value for each of the eg. 100 subsampled tables, and then making an average from that. Would there be any way that this would be possible? Thank you so much!!