Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • RockChalkJayhawk
    Senior Member
    • Mar 2009
    • 192

    Calculating Sample Size Estimates for RNA Sequencing Data

    I know this was picked up by the Newsbot, but I would like to mention that we created an R package on bioconductorthat can calculate sample size, power, effect sizes, etc. The paper is located here and the link to the Bioconductor package is here

    Feel free to let me know if you have any issues.
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    Out of curiosity, how does RNAseqPower compare to Scotty when the same pilot data is used for both? Regardless, that's a nice tool to have!

    Comment

    • RockChalkJayhawk
      Senior Member
      • Mar 2009
      • 192

      #3
      I hadn't actually seen it before, but it is a very good looking web app! The power estimates should be similar, but won't be identical for 2 reasons:

      1. Scotty assumes that variance follows a lognormal distribution. I think this is a valid and logical assumption. RNASeqPower uses the measured variance from the data instead.

      2. Scotty pseudo-randomly selects 200 values from a 2 dimensional matrix (variance and depth) then uses that mean as the power of the dataset. The way we use RNASeqPower is to say that we want to know the minimum requirements to reach the desired power. We need to know two things, how many genes do I want to detect and how much variance to I want to allow. If you only want to detect the top 90% of genes, then you take the 10th percentile of gene counts for your pilot (all other genes will have more than this count). You do the same for variance except you select the 90th percentile of variance. Using the minimum read count and the maximum allowable variance, that is the power you have in your experiment. Granted most genes will have higher power since they will have lower variance and higher counts. So in the end, it comes down to how you define what power actually means.

      RNASeqPower can also give you the power for each gene in a dataset.

      Comment

      • dpryan
        Devon Ryan
        • Jul 2011
        • 3478

        #4
        Cool, thanks!

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
          by SEQadmin2


          Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

          The systematic characterization of the human proteome has
          ...
          07-20-2026, 11:48 AM
        • SEQadmin2
          Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
          by SEQadmin2



          Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
          ...
          07-09-2026, 11:10 AM
        • SEQadmin2
          Cancer Drug Resistance: The Lingering Barrier to Rising Survival
          by SEQadmin2



          Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

          There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
          07-08-2026, 05:17 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Today, 12:17 PM
        0 responses
        9 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, Yesterday, 11:41 AM
        0 responses
        10 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 07-20-2026, 11:10 AM
        0 responses
        23 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 07-13-2026, 10:26 AM
        0 responses
        37 views
        0 reactions
        Last Post SEQadmin2  
        Working...