Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • DESeq2: Difference between condition+type vs. 3 conditions

    Dear all.

    I am unsure about how to use DESeq2 in the case of 3 conditions vs. 2 conditions + 2 types. Assuming I have the following design table
    Code:
              condition    type
    sample1   A            T1
    sample2   A            T1
    sample3   B            T2
    sample4   B            T2
    sample5   A            T2
    sample6   A            T2
    I am unsure about how this would be treated differently from
    Code:
              condition
    sample1   A:T1
    sample2   A:T1
    sample3   B:T2
    sample4   B:T2
    sample5   A:T2
    sample6   A:T2
    The second design table describes a 3-condition scenario.

    Now, obviously one would be interested in a detailed analysis of the counts for
    1. A:T2 vs. B:T2 (since they have the same type but a different conditions), and potentially
    2. A:T2 vs. A:T1 (since they have the same condition but different types).


    Question 1: If I reduce the problem to that of a 3-condition no-type design table, is this correctly taken into account?

    I know I would have to re-factor the columns of the 2nd matrix to reflect the correct order of fold changes that I want to calculate. So for example following re-factoring the levels as
    Code:
    levels=c("A:T2","B:T2","A:T1")
    and performing a DESeq2 analysis
    Code:
    dds<-DESeqDataSetFromMatrix(countData = countData, colData = design, design = ~ condition + type);
    dds<-DESeq(dds);
    Question 2: I could calculate the fold changes of B:T2 wrt A:T2 and A:T1 wrt A:T2, is this correct?
    I do get some issues with non-convergent dispersion fits, which I can get around if I call estimateDispersions manually with fitType="local".

    Question 3: But what happens in the case of the 1st condition+type table? I am confused as to the output of DESeq2. What role does the type play in the differential expression analysis and/or the dispersion fitting?

    Any help on this issue would be greatly appreciated.

    Regards,
    Maurits
    Last edited by mevers; 08-22-2013, 01:38 AM. Reason: Typo

  • #2
    In your first table, the type is always the same. Is this a typo? If not, I'm not sure I understand your question.

    Comment


    • #3
      Hi Simon.

      Yes, that was a silly mistake, you are absolutely right. I've changed it now in the original post. It should have read
      Code:
      type=c("T1","T1","T2","T2","T2","T2")
      Cheers,
      Maurits

      Comment


      • #4
        Question 1:

        You can technically represent it either way, although I would recommend to keep the variables separate for the following reason: if you combined the variables (as in "A:T1"), then you cannot make a clean B vs A comparison. Instead you have a B:T2 vs A:T1 comparison which mixes the effect of B vs A and T2 vs T1.

        Question 2:

        Note that fitType is also an argument for DESeq()

        Question 3:

        Both variables are used for finding fitted means (mu in the GLM formula given in the reference manual and vignette). And then the fitted means mu is used to estimate the dispersion. Dispersion is a measure of how far the counts deviate from the mu for that sample. Both variables will have fitted coefficients (betas in the GLM formula) and you can extract tests for each variable of the null hypothesis that the coefficients are equal to zero. By default the results for the last variable is provided by results(). For more, see the section in the vignette on "Multi-factor designs" and the man page for results().

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        66 views
        0 likes
        Last Post seqadmin  
        Working...
        X