First off, thanks for providing the code! It really helps with understanding QMP! This is more of a question than an issue:
Given that flow cytometry cell count data can have rather large standard deviations (eg., a median sd of 6.2e9 for frozen feces in Sup Table 5 in Vandeputte et al., 2017, it seems prudent to include this error in the QMP calculations. Otherwise, this possibly large source of analytical error is ignored. Any suggestions on how to alter your R code to include error propagation of the cell count error in the QMP calculations? Moreover, any recommendations on how one should incorporate that error into the downstream statistical analyses?
First off, thanks for providing the code! It really helps with understanding QMP! This is more of a question than an issue:
Given that flow cytometry cell count data can have rather large standard deviations (eg., a median sd of 6.2e9 for frozen feces in Sup Table 5 in Vandeputte et al., 2017, it seems prudent to include this error in the QMP calculations. Otherwise, this possibly large source of analytical error is ignored. Any suggestions on how to alter your R code to include error propagation of the cell count error in the QMP calculations? Moreover, any recommendations on how one should incorporate that error into the downstream statistical analyses?