Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • zeam
    Member
    • Oct 2010
    • 43

    Cuffdiff with more than two samples

    Recently,I'm using Cufflinks to process my RNA-seq data,and I have 4 tissues(more than two).In the mannual of cuffdiff,you said,if more than two are provided, Cuffdiff tests for differential expression and regulation between all pairs of samples. But when I used the cuffdiff,it only report one pairs of samples.There was only one gene_exp.diff in the output directory.Command was just as below,
    ======================================
    cuffdiff -o DEG_cuffdiff -p 20 cuffcmp.combined.gtf ./tophat_sd/accepted_hits.bam ./tophat_em/accepted_hits.bam tophat_en1/accepted_hits.bam ./tophat_en2/accepted_hits.bam
    ======================================
    What's my problem in using cuffdiff?
    Also I want to compare gene expression of duplicate gene pairs, and is it neccesery to remove the reads that mapped to multiple location in the BAM file generated by Tophat as duplicate gene pairs have high sequence similarity ?
    Thanks!
  • taozuo
    Junior Member
    • Jul 2011
    • 8

    #2
    For Question 1:

    Yes, Cuffdiff would only generate one output even though you have more than one comparisons. However, all other comparisons are also included in that file. Check the output file more carefully.

    Comment

    • zeam
      Member
      • Oct 2010
      • 43

      #3
      I have checked the outpu directory.

      In fact, I have checked the outpu directory, and there are only 11 files as follows:cds.diff , cds_exp.diff,genes.fpkm_tracking, isoforms.fpkm_tracking ,splicing.diff, tss_groups.fpkm_tracking , cds.fpkm_tracking , gene_exp.diff isoform_exp.diff , promoters.diff , tss_group_exp.diff. And I did not find any other comparison files between genes.
      I have posted this problem on the other thread ,and his met the similar issue.
      Application of sequencing to RNA analysis (RNA-Seq, whole transcriptome, SAGE, expression analysis, novel organism mining, splice variants)

      Comment

      • taozuo
        Junior Member
        • Jul 2011
        • 8

        #4
        I am not quite sure right now for your case. But I did this before and I can get all the comparisons. For example, if I open the gene_exp.diff file, I will have the q1 vs q2 first. And when I continue check the result by clicking pagedown, I will find q1 vs q3 and etc.

        Comment

        • zeam
          Member
          • Oct 2010
          • 43

          #5
          Thanks for your help

          Now I know all of results of the comparisons betweem pairs are saved in one file.
          Thanks!

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Cancer Drug Resistance: The Lingering Barrier to Rising Survival
            by SEQadmin2



            Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

            There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
            Today, 05:17 AM
          • GATTACAT
            Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
            by GATTACAT
            Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
            07-01-2026, 11:43 AM
          • SEQadmin2
            Nine Things a Sample Prep Scientist Thinks About Before Sequencing
            by SEQadmin2


            I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

            Here are nine questions we think about, in roughly the order they matter, before...
            06-18-2026, 07:11 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Yesterday, 11:05 AM
          0 responses
          7 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-02-2026, 11:08 AM
          0 responses
          28 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-30-2026, 05:37 AM
          0 responses
          28 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-26-2026, 11:10 AM
          0 responses
          27 views
          0 reactions
          Last Post SEQadmin2  
          Working...