Question: export edgeR normalised OTU table
Entering edit mode
Pozdrawiam • 0
Last seen 22 days ago

Hello, I want to normalise my microbiome data (OTU table) using edgeR package and export my normalised matrix ( I want to test some normalisation methods, normalised data will serve to calculate beta diversity and visualise by PCoA). I went through bioconductor forum and biostars and I still don't know want is the proper way to export normalised data from edgeR.

So, I read my matrix and calcNorm:

d <- DGEList(counts = d, group=group)
d = calcNormFactors(d)

what happens next?:

1st option

d = estimateCommonDisp(d, verbose=TRUE)
normalised_df <- d$pseudo.alt

2nd option

normalised_df <- cpm(d,  normalized.lib.sizes = T)

3rd option

dgList <- estimateCommonDisp(dgList)
dgList <- estimateTagwiseDisp(dgList)
normalised_df <- t(t(dgList$pseudo.counts)*(dgList$samples$norm.factors))

I have an impression that every post claims something different and at this point I am very confused. Thank you!

microbiome edgeR • 99 views
Entering edit mode

Is anything wrong with the answers you were provided here ( citing the manual that 2nd option is preferred and third is discouraged? It is good practice to indicate crossposts by adding links to the other community.

Entering edit mode

I just wanted to get a 3rd opinion, as your answer was different from the first one:)


Login before adding your answer.

Similar Posts
Loading Similar Posts
Traffic: 206 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.4