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

          • mylaser
            Reply to Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
            by mylaser
            Kheloyar – Everything You Need to Know About Kheloyaar Login and Kheoyar Id
            If you are looking for an online gaming platform that offers a user-friendly experience, Kheloyar has become a name that many users search for. Whether you're interested in creating a new account, accessing your dashboard through Kheloyaar Login, or learning how to obtain a Kheoyar Id, understanding the platform's features and account process is essential.
            This guide explains everything you need to know about...
            Today, 01:13 AM
          • SEQadmin2
            Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
            by SEQadmin2



            Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
            ...
            07-09-2026, 11:10 AM
          • 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...
            07-08-2026, 05:17 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 07-09-2026, 10:04 AM
          0 responses
          16 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-08-2026, 10:08 AM
          0 responses
          10 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-07-2026, 11:05 AM
          0 responses
          22 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 07-02-2026, 11:08 AM
          0 responses
          31 views
          0 reactions
          Last Post SEQadmin2  
          Working...