Question: Inconsistency Between Cufflinks And Cuffdiff. This Is Killing Me.
5.6 years ago by
Jia Meng • 10
Jia Meng • 10 wrote:
<<problem>> The cufflinks and cuffdiff results are not consistent with each other. This is killing me. Does it make sense? << Obseration >> We have 3 control samples and 3 treated sample. For many genes, their FPKM in cufflinks and cuffdiff are far from consistent. In cufflinks result, for a gene’s FPKM are control group (sample: 1,2,3)：0， 0， 4.8 treated group(sample: 1,2,3)：0， 0， 6.0 In cuffdiff, the estimated FPKM are control group： 12.6 treated group：2.0 <<method>> Use ucsc gene annotation gtf file, mm9, downloaded from UCSC table database Use cufflinks on each individual sample. Cufflinks: galaxy mirror at cistrome, minimal count:10, no quantile normalization, use gtf as reference, no background correction Use cufflinks on treated groups (3 biological replicates) and control groups (3 biological replicates) Cuffdiff: galaxy mirror at cistrome, minimal count:10, no normalization, use gtf as reference, no background correction <<additional comments="">> Cufflinks returns 55350 transcripts, while cuffdiff return 55418 transcripts, even though they use the same gene annotation gtf file. For the 6 cufflinks results (corresponding to 6 samples), the transcript ids are all the same, but the order are not, <<question>> Does it make sense? Or did I do anything wrong?
ADD COMMENT • link •