Question: How To Get Reads Counts From Cufflins?
0
gravatar for Li, Jilong (MU-Student)
6.8 years ago by
Hi all, I got the normalized values (FPKM) from cufflinks. And I want to get relative reads counts. How can I do that? Another question: how does cufflinks handle isoform genes while calculating the reads counts? Or what papers can help me understand this? Thank you very much! Victor
rna-seq cufflinks • 3.6k views
ADD COMMENTlink modified 6.8 years ago by Jeremy Goecks2.2k • written 6.8 years ago by Li, Jilong (MU-Student)80
0
gravatar for Jeremy Goecks
6.8 years ago by
Jeremy Goecks2.2k
Jeremy Goecks2.2k wrote:
Victor, It's not clear to me what you're looking for. FPKM is a normalized read count metric where the F stands for fragment, which is a single read (or half of a paired read). Expectation maximization is used to probabilistically assign reads to isoforms. See the Cufflinks documentation for details and paper links: http://cufflinks.cbcb.umd.edu/ Best, J.
ADD COMMENTlink written 6.8 years ago by Jeremy Goecks2.2k
Dear Jeremy, Sorry, I didn't expressed my question clearly. I got the FPKM normalized values for each gene from cufflinks. And I want to get the original reads counts that were not normalized from cufflinks. Could you please tell me how to get those? Thank you very much! Victor ________________________________ To: Li, Jilong (MU-Student) Cc: galaxy-user@lists.bx.psu.edu Subject: Re: [galaxy-user] How to get reads counts from cufflins? Victor, I got the normalized values (FPKM) from cufflinks. And I want to get relative reads counts. How can I do that? It's not clear to me what you're looking for. FPKM is a normalized read count metric where the F stands for fragment, which is a single read (or half of a paired read). Another question: how does cufflinks handle isoform genes while calculating the reads counts? Or what papers can help me understand this? Expectation maximization is used to probabilistically assign reads to isoforms. See the Cufflinks documentation for details and paper links: http://cufflinks.cbcb.umd.edu/ Best, J.
ADD REPLYlink written 6.8 years ago by Li, Jilong (MU-Student)80
Reads are probabilistically assigned, so raw read counts are not available from Cufflinks. Recovering raw fragment counts could be done by reverse-engineering the FPKM value, but Cufflinks doesn't do this for you. If you choose to do this, keep in mind that Cufflinks uses an "effective transcript length". Best, J.
ADD REPLYlink written 6.8 years ago by Jeremy Goecks2.2k
Hi Victor, I also wanted to get the count reads but it did not work with Cufflinks. You can use GenomicRanges package (Bioconductor). I am also using Deseq (negative binomial distribution) to detect the differentially expressed genes from time series data. You can find it here http://bioconductor.org/packages/release/bioc/html/DESeq.html Best, Luciano Luciano Cosme PhD Candidate Texas A&M Entomology cosme@tamu.edu www.lcosme.com
ADD REPLYlink written 6.8 years ago by Luciano Cosme20
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 16.09
Traffic: 168 users visited in the last hour