Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Comparative Time Series Analysis of RNA-Seq data?

    I need help with comparative time-series analysis of RNA-Seq Data. I have an untreated sample, then 3 treatments at day 3, the same 3 treatments at day 6 and then at day 9. I have tag-counts(generated from HT-Seq counts) for all of these samples across time. Could you please suggest some tools to get Differentially expressed genes. I do not wish to perform pair-wise comparisons.

    Thank you so much

    Ashu

  • #2
    DESeq(2), edgeR, limma/voom, etc. Anything that can handle GLMs will work.

    Comment


    • #3
      Great, than you so much.

      Could you or anybody explain this to me with an example? I am new to this and its tough to understand what functions/methods to use.

      Originally posted by dpryan View Post
      DESeq(2), edgeR, limma/voom, etc. Anything that can handle GLMs will work.

      Comment


      • #4
        I expect there's an example of this in one of the vignettes/user guides for a couple of the packages I listed. In general, you just want to fit your data with a GLM of the form "counts ~ time*treatment", assuming that there can be a meaningful time:treatment interaction in your experiment (otherwise, just exchange a plus sign for the asterisk).

        Comment


        • #5
          The pdf(http://bioconductor.org/packages/rel.../doc/DESeq.pdf) from bioconductor does this ->
          > fit1 = fitNbinomGLMs( cdsFull, count ~ libType + condition )
          > fit0 = fitNbinomGLMs( cdsFull, count ~ libType )
          > pvalsGLM = nbinomGLMTest( fit1, fit0 )
          > padjGLM = p.adjust( pvalsGLM, method="BH" )

          I dont understand what are they trying to do here. Below is what I tried with my data

          >glm_design

          samples libType time treatment

          Samp_1 single-end d0 Untreated
          Samp_2 single-end d3 treatment1
          Samp_3 single-end d6 treatment1
          Samp_4 single-end d9 treatment1
          Samp_5 single-end d3 treatment2
          Samp_6 single-end d6 treatment2
          Samp_7 single-end d9 treatment2
          Samp_8 single-end d3 treatment3
          Samp_9 single-end d6 treatment3
          Samp_10 single-end d9 treatment3

          glm_cds = newCountDataSet(combined_file , glm_design)

          glm_cds = estimateSizeFactors(glm_cds)

          # I dont have replicates
          glm_cds = estimateDispersions(glm_cds, method="blind", sharingMode="fit-only" )

          fit1 = fitNbinomGLMs( glm_cds, count ~ time*treatment)

          fit0 = fitNbinomGLMs( glm_cds, count ~ time)

          pvalsGLM = nbinomGLMTest( fit1, fit0 )

          padjGLM = p.adjust( pvalsGLM, method="BH" )


          IS THIS WHAT I AM SUPPOSED TO DO?

          unique(padjGLM)
          [1] 1 NA

          padj values do not seem right??
          *************************************

          Originally posted by dpryan View Post
          I expect there's an example of this in one of the vignettes/user guides for a couple of the packages I listed. In general, you just want to fit your data with a GLM of the form "counts ~ time*treatment", assuming that there can be a meaningful time:treatment interaction in your experiment (otherwise, just exchange a plus sign for the asterisk).
          Last edited by ashuchawla; 07-23-2013, 10:13 AM.

          Comment


          • #6
            # I dont have replicates
            glm_cds = estimateDispersions(glm_cds, method="blind", sharingMode="fit-only" )

            fit1 = fitNbinomGLMs( glm_cds, count ~ time*treatment)
            fit0 = fitNbinomGLMs( glm_cds, count ~ time)
            Well, you sort of have replicates, just not full replicates. You might try the default dispersion method and see how well that fits. For fit0, try "fit- <- fitNbinomGLMs(glm_cds, count~time+time:treatment)". Otherwise, you getting output for things that can vary by both treatment and a time:treatment interaction, which you indicated not wanting.

            BTW, an adjusted p-value of 1 is normal and to be expected. That all of them are 1 is too bad, though doing as I suggested above my change that. You can also try doing some independent filtering (see the genefilter package and the accompanying paper in PNAS), which might improve things further.

            Comment


            • #7
              > glm_cds = estimateDispersions(glm_cds)
              Error in .local(object, ...) :
              None of your conditions is replicated. Use method='blind' to estimate across conditions, or 'pooled-CR', if you have crossed factors.



              Originally posted by dpryan View Post
              Well, you sort of have replicates, just not full replicates. You might try the default dispersion method and see how well that fits. For fit0, try "fit- <- fitNbinomGLMs(glm_cds, count~time+time:treatment)". Otherwise, you getting output for things that can vary by both treatment and a time:treatment interaction, which you indicated not wanting.

              BTW, an adjusted p-value of 1 is normal and to be expected. That all of them are 1 is too bad, though doing as I suggested above my change that. You can also try doing some independent filtering (see the genefilter package and the accompanying paper in PNAS), which might improve things further.

              Comment


              • #8
                Ah, right, you might try pooled-CR.

                Comment


                • #9
                  Thank you so much for your help.

                  This is what it tells me when I try pooled-CR

                  glm_cds = estimateDispersions(glm_cds, method="pooled-CR", sharingMode="fit-only" , modelFormula=count~time+time:treatment)
                  Error in FUN(newX[, i], ...) :
                  No residual degrees of freedom. Most likely the design is lacking sufficient replication.

                  Originally posted by dpryan View Post
                  Ah, right, you might try pooled-CR.

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Recent Developments in Metagenomics
                    by seqadmin





                    Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                    09-23-2024, 06:35 AM
                  • seqadmin
                    Understanding Genetic Influence on Infectious Disease
                    by seqadmin




                    During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                    Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                    09-09-2024, 10:59 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 10-02-2024, 04:51 AM
                  0 responses
                  11 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 10-01-2024, 07:10 AM
                  0 responses
                  19 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 09-30-2024, 08:33 AM
                  0 responses
                  24 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 09-26-2024, 12:57 PM
                  0 responses
                  17 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X