Hi All, I have a featurecount output row counts and need to convert them to FPKM. I knew about the countToFPKM package, however, the description of the package says its for paired end RNA seq data. My data are single end. Could I still use this package? and if not could you please suggest an easy way to do it?
The reason I need to convert from row counts to FPKM is to perform ssGSEA analyses on gene pattern. Thanks.
Mohamed
1 answer
Once you have a matrix of counts per feature, there is no way to tell if the counts for the feature came from paired-end data or single end data. The single or double endedness of your DNA Fragments is accounted for in the generation of an integer count representing a feature. Thus, your count matrix can be converted to FPKM whether it was generated from single or paired-end data. FPKM is used to represent abundance regardless of end type, and is identical to RPKM for this same reason. RPKM simply suffers from a historical prematurity in that "Reads" actually represent "Fragments", whether sequenced by one end or two. One caveat: the package, and the method it implements (from Trapnell et al., 2010), uses the average DNA fragment length for a given library in the calculation. This can be determined easily in paired end data, but for single end data you usually have to estimate the library insert size (usually an educated guess).
Log in to answer this question.
-
How to convert HTSeq raw read counts to FPKMs?
written by s.singh •
-
CPM and FPKM for machine learning
written by bioinfo
-
Question regarding TCGA and TCGAbiolinks --> RNAseq data as HTseq-Counts, FPKM, FPKM-UQ
written by mario.red8976
-
How to export STAR alignment output files
written by salehm •
-
ssGSEA analysis workflow and interpretation
written by salehm •
-
Download data from the TCGA for gene expression analyses in R
written by Lucy •
-
ssGSEA 2.0 R Package?
written by ivanarg2
-
blastn for paired end reads
written by badribio
-
RPKM to FPKM conversion
written by kanwarjag
-
Using Bayespeak for Paired-end reads
written by Nathan •
Does the FPKM metric apply to SE data? What would be the difference between RPKM and FPKM in the context of SE data?
It sounds that yes the FPKM needs paired end. So I may need to convert my single end data to RPKM, any easy idea how I can do that? Thanks.
What have you tried? Have you searched online for "get RPKM from raw counts"?
I was thinking of using countToFPKM package but it says it works for paired end reads, I am not sure if I could also use it for single end reads.
So you searched for "get RPKM from counts" and the only thing you see is the countToFPKM package?
Is featurecounts going to do FPKM right? Don't you need something like Kallisto or RSEM that understands how the different transcripts are proportioned?