Output of featureCounts() is not integer
0
1
Entering edit mode
4.6 years ago
Francis ▴ 20

Hi, I using featureCount() to get the reads count for next analysis from sorted SAM file, and I find the read count is not integer. It is not feasible because I want to use DESeq2 to analyse the different genes.

So how to get integer by featureCount() or normalize the output to integer?

thanks.

RNA-Seq • 3.5k views
ADD COMMENT
3
Entering edit mode

As a side comment for you and for future readers on this topic, the fractional counts usually mean that they are expected counts and not raw counts. Simply rounding them up won't be the best choice, because raw and expected counts are different numbers.

The best approach, if you want to use DESeq2 afterwards, is to output raw counts (which will be therefore integers).

ADD REPLY
0
Entering edit mode

From the link above.

The help page for featureCounts() says

"fraction:

logical indicating if fractional counts will be produced for multi-mapping reads. If TRUE, a fractional count, 1/n, will be generated for each reported alignment of a multi-mapping read, where n is the total number of alignments reported for that read. countMultiMappingReads must be set to TRUE when fraction is TRUE."

ADD REPLY
0
Entering edit mode

So those fractions are indeed expected based on distributing multimapping reads, they're not exact nor raw.

ADD REPLY
0
Entering edit mode

Hi,

I do not run with fraction=TURE & countMultiMappingReads=TURE.

I use

featureCounts -T 8 -g gene_id -p -a Yeast.ercc.gff -o count.txt Yeast.sam

Now the output file is not integer. So if I use fraction and countMultiMappingReads, this should be fraction?

The following comment also said "Simply rounding them up won't be the best choice" and I think so. But I don`t have a good idea to output the integer reads count except rounding tnem up.

thank u for your repply.

ADD REPLY
1
Entering edit mode

Maybe a workaround: htseq-count outputs integer counts.

ADD REPLY
0
Entering edit mode

so does featureCounts with default settings

ADD REPLY
1
Entering edit mode

Did you run featureCounts with (fraction=TRUE)? Then you could round the fractions as shown in this thread.

ADD REPLY

Login before adding your answer.

Traffic: 2953 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