Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • ekimmike
    Member
    • Apr 2012
    • 14

    with or w/o pseudocounts in DESeq?

    Hi,

    I'd like to ask about your opinion on using pseudocounts in DE analysis, as some people say they're absolutely necessary some disagree ...
    I've read all the vignettes - still no unclear, hence decided to test it on my data

    I did comparison between DEseq, edgeR, BaySeq and cuffDiff on RNAseq data (these are tophat genes, HTseq tables, FDR = 10% for all)

    please see the file


    In case of condition 1 it seems clear DESeq performs best, but in condition 2 it appears most affected by pseudocounts detecting about 1600 in case 1, 2000 in case2 with only 1000 overlap

    any suggestions?
  • Gordon Smyth
    Member
    • Apr 2011
    • 91

    #2
    I suggest that you use the packages that you mention according to the instructions that come with them. None of the package permit you to do ad hoc transformations of the counts prior to analysis.

    BTW, the term 'pseudocount' has a technical meaning for the edgeR package, but it part of the internal computations, not something that is manipulated by a user. It is not a meaningful concept for any of the other packages.

    Comment

    • Simon Anders
      Senior Member
      • Feb 2010
      • 995

      #3
      I really wonder where this idea that "pseudocounts are absolutely necessary" comes from.

      The DESeq vignette states "The count values must be raw counts of sequencing reads. This is important for DESeq’s statistical model to hold, as only the actual counts allow assessing the measurement precision correctly. Hence, please do do not supply other quantities, such as (rounded) normalized counts, or counts of covered base pairs – this will only lead to nonsensical results."

      I guess in the next release I will put the first sentence in red boldface. Every other week somebody comes up with a new idea what to try as input. So, again: If the numbers in your count tables are not counts of reads, they are not suitable. Pseudocounts are not counts of reads, but counts of nothing.

      So, kudos for your work of testing all these methods with pseudocounts, but why did you bother?

      Comment

      • ekimmike
        Member
        • Apr 2012
        • 14

        #4
        Originally posted by Simon Anders View Post
        I guess in the next release I will put the first sentence in red boldface.
        I know, I know ... I'll move this discussion to bioconductor forum, and try to explain better

        Originally posted by Simon Anders View Post
        So, kudos for your work of testing all these methods with pseudocounts, but why did you bother?
        mostly because I switched from nonparametric tools
        Last edited by ekimmike; 10-04-2012, 01:22 PM.

        Comment

        Latest Articles

        Collapse

        • 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
        • seqadmin
          Investigating the Gut Microbiome Through Diet and Spatial Biology
          by seqadmin




          The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
          02-24-2025, 06:31 AM

        ad_right_rmr

        Collapse

        News

        Collapse

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