Announcement

Collapse
No announcement yet.

What happened to per-condition and gene-est-only in DESeq2?

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • What happened to per-condition and gene-est-only in DESeq2?

    In DESeq, there was an option to estimate the gene dispersions within each condition instead of across all conditions, and per gene instead of sharing dispersion estimates across genes.

    However, in DESeq2, I can't find these options anywhere -- how to specify these tweaks to the dispersion estimates now?

    Thank you!

  • #2
    In DESeq2 we only estimate a single dispersion value per gene. See the methods description in the paper:

    http://genomebiology.com/2014/15/12/550

    Comment


    • #3
      Thanks, Michael.

      So it is clear to me now that the dispersion values are now per-gene, so gene-est-only no longer applies.

      Quoting from the paper:

      Many methods for differential expression analysis of RNA-seq data perform such information sharing across genes for variance (or, equivalently, dispersion) estimation. edgeR[2],[3] moderates the dispersion estimate for each gene toward a common estimate across all genes, or toward a local estimate from genes with similar expression strength, using a weighted conditional likelihood. Our DESeq method [4] detects and corrects dispersion estimates that are too low through modeling of the dependence of the dispersion on the average expression strength over all samples.

      What about estimating the dispersions within each condition instead of across all samples? (per-condition)?

      Thanks.

      Comment

      Working...
      X