Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • narges
    Member
    • Aug 2012
    • 29

    DESeq raw variance

    Variance in DESeq is the sum of shot noise and raw variance. The raw variance is the smooth function of estimated read counts and experimental condition. In the corresponding paper, it is mentioned that this function aim to pool the genes with the same strength to estimate the variance(due to the more probable few number of replicates).
    I am not sure if I have understand it correctly. Does this means that DESeq groups n genes (i1, i2, ...in) from a common condition with similar expression levels and then calculate the variance for them? So, instead of calculating the variance for a specific gene, it calculates the variance for a group of similar gene and associate that value for all of them?
    And how does it decide which genes have similar strength (? similar expression levels)?
    Then can we conclude that DESeq does not calculate the gene-specific variance but group of genes-specific variance?
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    Well, variance is the sum of technical (shot noise) and biological sources, though perhaps you mean the latter by "raw".

    The process of fitting the smooth curve to the data is essentially the "pooling information" step, since variances are then shrunken toward it (generally you use distance from this line as a penalty and perform maximum-likelihood estimation of the variance with that, though I recall DESeq2 also treats genes with variance >3sigma (or something like that, it's in their paper) from the expected differently). So no, DESeq2 doesn't group genes from a common condition with similar expression levels and then calculates variance from that. It calculates per-gene variance and then shrinks that toward the expected value.

    Comment

    • narges
      Member
      • Aug 2012
      • 29

      #3
      Originally posted by dpryan View Post
      Well, variance is the sum of technical (shot noise) and biological sources, though perhaps you mean the latter by "raw".
      Yes, I meant the biological oriented variance.
      Thanks!
      Is it true that point "Px" in smooth curve is the mean of count values for gene gene x in the condition p? So, it forms the curve based on the means for each gene of a condition and next shrunken the variances to their corresponding point on the mean-specific smooth curve?

      Comment

      • dpryan
        Devon Ryan
        • Jul 2011
        • 3478

        #4
        No, that's not strictly true. It will typically be the case that there's only 1 gene with a mean expression of a certain value, so that would mean that the curve must go through that mean, which will typically not be the case. The curve is fit to the means, as you mentioned (at least that's my recollection), so the remainder of what you wrote looks correct.

        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
        50 views
        0 reactions
        Last Post seqadmin  
        Started by seqadmin, 03-03-2025, 01:15 PM
        0 responses
        201 views
        0 reactions
        Last Post seqadmin  
        Working...