Question: What's the difference of FPKM generated by htseq and stringtie?
0
gravatar for hxlei613
16 months ago by
hxlei61390
Bangladesh
hxlei61390 wrote:

Hi I get FPKM for gene (not transcript) from stringtie by these parameter: -e -G -A. I also use htseq (default parameter )to generate counts matrix for the same gene sets and calculate fpkm using R. The gtf I used in this two method is the same. However the result fpkm is not the same. It seems that the overall fpkm generated by htseq is lower than stingtie. I try to use stringtie to get counts matrix for each gene but stringtie seems not generate count for gene level. I have searched the difference between them but I didn't find any answers. Thank you very much for helping me!

stringtie fpkm htseq • 778 views
ADD COMMENTlink modified 16 months ago by Devon Ryan94k • written 16 months ago by hxlei61390
2
gravatar for Devon Ryan
16 months ago by
Devon Ryan94k
Freiburg, Germany
Devon Ryan94k wrote:

StringTie is able to use multimappers, which htseq will ignore them. That is the most probably cause for the difference you're seeing.

ADD COMMENTlink written 16 months ago by Devon Ryan94k

That is helpful. However I also find some gene whose counts are not zero in all samples and the fpkm I calculate is not zero (for example the average fpkm of such genes can be 300 or 3000) but fpkm generated by stringtie is zero. What's possible reason for that?

ADD REPLYlink modified 15 months ago • written 15 months ago by hxlei61390

Without further information it's impossible for anyone to know.

ADD REPLYlink written 15 months ago by Devon Ryan94k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1678 users visited in the last hour