Normalization factors in metagenomeSeq
2
0
Entering edit mode
jovel_juan ▴ 10
@jovel_juan-7129
Last seen 12 months ago

I have two questions:

QUESTION 1. In metagenomeSeq:

# Calculate proper percentile to normalize data

percentile = cumNormStatFast(genusObj, pFlag = FALSE, rel = 0.1)

percentile

genusObj = cumNorm(genusObj, p = percentile)

genusObj

If I substitute p = percentile by p = 1 , will metagenomeSeq perform total sum normalization?

QUESTION 2:

Is it possible to conduct DESeq and TMM normalization using metagenomeSeq?

microbiome • 1.3k views
1
Entering edit mode
@joseph-nathaniel-paulson-6442
Last seen 4.6 years ago
United States

Hi Jovel, Thanks for the interest. First, an answer to your first question: Yes, setting the p = 1 will result in total sum normalization.

To answer the second question: Yes, it is possible to use any other scaling method, please see: A: normFactors when aggregating in metagenomeSeq

pData(obj@expSummary$expSummary)$normFactors = values

0
Entering edit mode
jovel_juan ▴ 10
@jovel_juan-7129
Last seen 12 months ago

Thanks for the prompt response. It was really helpful!