Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 4 conditions in DESeq2

    I have four conditions A, B, C, and D that each have three HTSeq count files.

    I want to use DESeq2 to compare A to B, A to C, A to D, B to C, B to D etc.

    With DESeq, I would bring all the 12 files into a countDataSet and then I could specify which of the conditions to compare.

    With DESeq2, I can bring all the 12 files into a DESeqDataSet, but I can not figure out how to do specify the particular comparisons.

  • #2
    It's easy. Use the contrast variable of the results function. Comparing all pairs of contrasts is not a good idea, though. Why don't you have a better biological hypothesis ?

    Comment


    • #3
      We have wild type condition A and wild type condition B as well as treated condition A and treated condition B. So, we would like to see changes in gene expression between condition A and condition B, as well as changes in gene expression between wild type and treated.

      Will contrast provide us with this ? To me it seems like the analysis is done by the DESeq function and the results function just organizes the output of the DESeq function.

      Comment


      • #4
        Yes, a contrast allows you to perform pretty much any comparison you want. Note that the comparisons to A don't require contrasts, since those are just fit coefficients.

        Comment


        • #5
          hi Matt,

          For lots more information about what you want to do, check out the man page for ?results and section 3.2 Contrasts in: vignette("DESeq2")

          Comment


          • #6
            Thank you very much for your answers Dario, Devon and Michael !

            Comment


            • #7
              > "Note that the comparisons to A don't require contrasts, since those are just fit coefficients."

              Safest always to use contrast=c("condition","B","A"), because when we have expanded model matrices for symmetric shrinkage, then the coefficients are not B_vs_A, but simply A, B, etc. This will be clear from resultsNames(dds)
              Last edited by Michael Love; 01-22-2015, 12:24 PM.

              Comment


              • #8
                OK. Thank you for clarifying that.
                (And also for making DESeq2 !)

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Recent Innovations in Spatial Biology
                  by seqadmin


                  Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.

                  3D Genomics
                  While spatial biology often involves studying proteins and RNAs in their...
                  Yesterday, 07:30 PM
                • seqadmin
                  Advancing Precision Medicine for Rare Diseases in Children
                  by seqadmin




                  Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
                  12-16-2024, 07:57 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 12-30-2024, 01:35 PM
                0 responses
                21 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 12-17-2024, 10:28 AM
                0 responses
                41 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 12-13-2024, 08:24 AM
                0 responses
                55 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 12-12-2024, 07:41 AM
                0 responses
                40 views
                0 likes
                Last Post seqadmin  
                Working...
                X