Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • jorgebm
    Member
    • Feb 2010
    • 18

    Outliers management using DESeq

    Hello,


    Due my lack in statistical skills I would like you to ask for advice:

    If I understood it properly, DESeq manages outliers at "gene" level. I mean, large variance makes a gene "harder" to be differential expressed (specially using "sharingmode=maximun") ; so, it prevents false positives.

    But what happens If having some samples; one of this samples is quite diferent compared with the remaining ones of its condition. This sample is an outlier and induces to DESeq to avoid several otherwise (excluding this sample) truly diferential expresed genes.

    As I understand It's a must to perform a preocedure to exclude sample outliers from conditions If we (are lucky) to have several samples for each condition. doesn't it?

    In my case, for an RNA-Seq experiment, I've modestly large number of samples (in my case 40 replicates for condition A and 36 replicates for condition B). What procedure do you recomend to prepare sample sets (filter outliers) to feed DESeq?.

    Thanks in advance
  • Simon Anders
    Senior Member
    • Feb 2010
    • 995

    #2
    With such an unusually large number of samples, you can indeed afford to remove outlier samples that may compromise your detection power. Cluster your samples using the variance stabilizing transformation, as described in the DESeq vignette to identify outlier samples. Furthermore, once you have more than, say, 10 or 15 degrees of freedom for variance estimation, you might be better off not using the fitted data at all (sharingMode="gene-est-only").

    Comment

    • jorgebm
      Member
      • Feb 2010
      • 18

      #3
      Thank you, Simon.

      I'll try what you say.

      Comment

      • Wolfgang Huber
        Senior Member
        • Aug 2009
        • 109

        #4
        Hi Jorge

        the arrayQualityMetrics package might provide some help in rationalising the hunt for outliers (it does not fully automate it, in the end there is still a human judgement call). See e.g. http://permalink.gmane.org/gmane.sci...onductor/40146 for some hints.

        Best wishes
        Wolfgang
        Wolfgang Huber
        EMBL

        Comment

        • jorgebm
          Member
          • Feb 2010
          • 18

          #5
          Thank you, Wolfgang.

          Just now I'm busy with another (priorized) task but soon I'll back yo this issue...probably with a lot of dubts

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Pathogen Surveillance with Advanced Genomic Tools
            by seqadmin




            The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
            03-24-2025, 11:48 AM
          • seqadmin
            New Genomics Tools and Methods Shared at AGBT 2025
            by seqadmin


            This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

            The Headliner
            The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
            03-03-2025, 01:39 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-20-2025, 05:03 AM
          0 responses
          49 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-19-2025, 07:27 AM
          0 responses
          57 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-18-2025, 12:50 PM
          0 responses
          49 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-03-2025, 01:15 PM
          0 responses
          200 views
          0 reactions
          Last Post seqadmin  
          Working...