Using alternative LFC shrinkage methods for computing p-values in DESeq2
Entering edit mode
zefrieira • 0
Last seen 3.5 years ago

In this post, Mike Love explained the reasoning for moving LFC shrinkage into its own step using the lfcShrink function, decoupling testing from shrinkage.

Is there any way to use apeglm-shrinked LFC to compute p-values (similar to using betaPrior=TRUE, but with the new estimator)?

I know that you can get significance levels when you use lfcThreshold with apeglm. Is there a reason that this is only possible when you incorporate lfcThreshold into the null hypothesis?


deseq2 apeglm statistics • 1.5k views
Entering edit mode
Last seen 16 hours ago
United States

For apeglm and ashr we compute s-values if you set svalue = TRUE, or for apeglm if you set a threshold. See the apeglm publication for details on the s-value. We found this to make more sense for summarizing over a set of posterior distributions than p-values.

Entering edit mode

Thank you very much!

Is there a reason you recommend using a stricter threshold for s-values? Is this something you observed empirically?

Also, are there any plans to unify the interface for getting significance levels based on shrunken LFC? Right now, there are two different interfaces:

  • If the user wants to get padj values based on the classical shrinkage, they method must use betaPrior=TRUE in the DESeq function.
  • If the user wants s-values based on ashr or apeglm-shrunken LFCs, they must use betaPrior=FALSE in the DESeq function and then set svalue=TRUE in lfcShrink.

Nevertheless, thanks for your help. From what I've seen so for, it seems like apeglm is giving me more biologically relevant results.

Entering edit mode

Yes, for example, see here:

A very large p-value is 1. A very large s-value is 0.5 (probability of wrong sign is equal to random guessing). In general, they represent fairly distinct conceptual amounts though, so I'd just think about what bound you want on wrong sign LFCs.

Also, along the same lines, no we want be harmonizing these interfaces. betaPrior=TRUE is for all intents deprecated. I didn't remove it altogether for backward compatibility. But it is not recommended (e.g. see the performance of the Normal prior vs apeglm and ashr in the apeglm paper).

Also, the s-value are the more meaningful statistic than p-values for sets of posterior distributions, of which the shrunken LFC represent the mode.


Login before adding your answer.

Traffic: 1168 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6