Featurecounts normalize by gene length
1
1
Entering edit mode
2.9 years ago
David ▴ 210

Hello,

I have run featureCounts on aligned bam files (DNA from bacterial genomes).

I have a table with counts per sample as expected. However i would like to normalize by gene length since genes of different length can have more reads mapped then genes that are shorter.

I was wondering what would be the best way to do that ?? Any tool or script ? Once i have normalized by gene length i would like to feed it to DESEq2

Thanks

featureCount tpm deseq2 • 3.0k views
ADD COMMENT
1
Entering edit mode
2.9 years ago
ATpoint 55k

DESeq2 takes care of all necessary steps internally. There is no need to manipulate the data prior to creating your DESeq2 object. Raw counts are expected. Please read the manual.

ADD COMMENT
0
Entering edit mode

Hi , I understand deseq will normalize by library size but prior to that you need to normalize by gene length. Am i correct?

ADD REPLY
1
Entering edit mode

Not really. There is no need to account for gene length because gene length is identical for all conditions. See also as a starting point for further online searches: https://support.bioconductor.org/p/67131/ and DESeq2 normalisation: is the size of the gene taken into account?

ADD REPLY
0
Entering edit mode

Ok thanks, So ranking genes by expression in this context should not be possible.

ADD REPLY
0
Entering edit mode

How do you mean it's not possible? DESeq will do a far more sophisticated normalization than you can think off, and you can get the normalized counts. Then you can rank whatever you want. And now, before we help you any further, you are going to read all the links ATpoint gave to you.

ADD REPLY

Login before adding your answer.

Traffic: 2824 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6