Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • eEtting up DESeq 2 analysis

    I have some RNA-seq data I'm trying to do DESeq2 analysis on, and it's more complicated than what I've done before.

    Patients came in two different times and had their blood drawn and sequenced. There is some missing data, so that some patients only have data for visit 1 or 2. In between visits patients were given either a placebo or one of two treatments.

    The answer we're trying to answer is what expression levels changed in between visits in a different way between the treatment groups. Based on the vignettes and other answered questions it seems the design would be:

    ~ patient + visit + treatment + visit:treatment

    the interaction visit:treatment would be the 'difference in differences' of each treatment over the time between visits. Is this correct?

  • #2
    You can simply exclude the patients for whom you only have a single sample. They'll get ignored in the analysis anyway. Anyway, yes your model is correct and you do indeed care most about the "visit:treatment" term.

    Comment


    • #3
      Thanks for the answer!

      I did get this error however:

      Error in checkFullRank(modelMatrix) :
      the model matrix is not full rank, so the model cannot be fit as specified.
      One or more variables or interaction terms in the design formula are linear
      combinations of the others and must be removed.

      my pheno file looks like:

      sampleName visit condition patient
      1 V2 control 1
      2 V5 control 1
      3 V2 treatment 2
      4 V5 treatment 2
      5 V5 treatment 3
      6 V2 treatment 3
      7 V5 treatment 4
      8 V2 treatment 4
      9 V2 control 5
      10 V5 control 5

      Removing patients from the experimental design worked. Is there any way, or value, to preserve the patient data?

      Comment


      • #4
        Indeed, I should have foreseen that :P

        If you were to instead use "~patient+condition:visit+visit" and got rid of the "conditiontreatment:visitV2" column in the model matrix then the result would work. The original problem was that each condition is comprised of a set of patients, so you can't have patient coefficients and a "condition" coefficient (which is just the average of the patient coefficients!).

        Sorry that that's so confusing.
        Last edited by dpryan; 12-16-2015, 06:34 AM.

        Comment


        • #5
          Thank you for the help!

          I apologize, I'm not entirely clear on how to set up my model matrix based on your answer. It seems I would still need every column if I were to use "~patient+condition:visit+visit".

          Comment


          • #6
            I had a typo in my reply, I meant to remove the "conditiontreatment:visitV2" from the model matrix. That'll make it full rank,

            Comment


            • #7
              Ahh I see, I'm getting my sample table and the model matrix confused.

              So is this the correct way to use my own model matrix?

              design_string <- "~patient+condition:visit+visit"
              sample_table <- read.table(input_file, row.names = NULL, header = T, sep = ",")
              deseq_object <- DESeqDataSetFromHTSeqCount(sampleTable = sample_table,
              design = ~condition, #have to have something here
              directory = count_folder)
              mm <- model.matrix(as.formula(design_string), sample_table)
              mm <- mm[,-19] # gets rid of conditiontreatment:visitV2
              deseq_object <- DESeq(deseq_object, full=mm, betaPrior=FALSE)

              Comment


              • #8
                Something along those lines at least.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Genetic Variation in Immunogenetics and Antibody Diversity
                  by seqadmin



                  The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                  11-06-2024, 07:24 PM
                • seqadmin
                  Choosing Between NGS and qPCR
                  by seqadmin



                  Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                  10-18-2024, 07:11 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 11-08-2024, 11:09 AM
                0 responses
                136 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 11-08-2024, 06:13 AM
                0 responses
                110 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 11-01-2024, 06:09 AM
                0 responses
                67 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 10-30-2024, 05:31 AM
                0 responses
                25 views
                0 likes
                Last Post seqadmin  
                Working...
                X