Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Deseq2 multifactor analysis

    Hello,

    I have a RNAseq experiment where I have 2 conditions and 3 replicates each and replicates are taken from different individuals. So my design is as follows:

    Sample Condition Individual
    S1 treated A
    S2 untreated A
    S3 treated B
    S4 untreated B
    S5 treated C
    S6 untreated C

    Here is my code to analyse the data using deseq2:
    Code:
    library('DESeq2')
    directory<-"counts"
    sampleFilesN1_N2 <- grep("*N[1-2]",list.files(directory),value=TRUE)
    sampleConditionN1_N2<-c("treated","untreated","treated","untreated","treated","untreated")
    sampleIndividual<-c("A","A","B","B","C","C")
    sampleTable1<-data.frame(sampleName=sampleFilesN1_N2, fileName=sampleFilesN1_N2, condition=sampleConditionN1_N2, individual=sampleIndividual)
    ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable1, directory = directory,design= ~ condition + individual)
    ddsHTSeq$condition <- relevel(ddsHTSeq$condition, "untreated")
    dds<-DESeq(ddsHTSeq)
    res<-results(dds)
    res<-res[order(res$padj),]
    write.csv(as.data.frame(res),file="DeSEQ2_treated_vs_untreated.csv")
    However, when I run these and look at results, following message:
    Code:
    "log2 fold change (MAP): individual C vs A
    Wald test p-value: individual C vs A"
    So, its only giving me results after comparing C vs A. Am I missing anything here?

    Thanks for any help.

  • #2
    results() will just give you one of the possible comparisons (the last one in the model matrix, by default). See help(results) for how to get the results you want.

    Comment


    • #3
      Originally posted by dpryan View Post
      results() will just give you one of the possible comparisons (the last one in the model matrix, by default). See help(results) for how to get the results you want.
      Thanks for the hint, that's really helpful. So, I understand if I want to get comparison of conditions considering the pairing of samples into account, I should give

      Code:
      res <- results(dds, contrast=c("condition","treated","untreated"))
      Is this correct?

      Thanks again.
      Last edited by genomica; 09-24-2015, 07:41 AM.

      Comment


      • #4
        Yup and I see you figured out the error in your earlier attempt

        Comment


        • #5
          Thanks very much Devon. Yes, I did figure out the first error

          I am a bit paranoid as I am doing this analysis for the first time and hence wondered if you could confirm if my experimental design is ok.

          Code:
          ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable1, directory = directory,[B]design= ~ condition + individual[/B])
          or do I need:
          Code:
          ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable1, directory = directory,[B]design= ~ condition + individual +condition:individual[/B])
          Many thanks again.

          Comment


          • #6
            Go for the first one, the second one ends up directly comparing individual samples and therefore isn't reliable.

            Comment


            • #7
              Thanks a lot for your help Devon.

              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
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X