FNR in RNASeq
2
0
Entering edit mode
prat • 0
@prat-10026
Last seen 5.4 years ago

Why are false neg rates usually not seen/used for RNASeq analysis? My guess is that the TP is not known... Is it possible to still calculate/estimate the FNR in DESeq2/edgeR?

deseq2 edger • 716 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 2 days ago
United States

One way to get an idea about possible false negative rate is to conduct an RNA-seq power analysis. There are a number of Bioconductor packages for this. There are various dimensions and parameters at play, including number of samples, sequencing depth, biological variability across replicates, and of course the true effect size.

ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 4 hours ago
WEHI, Melbourne, Australia

No, it isn't possible to estimate the FNR for real data. Computing the FNR is only possible in a simulation study or with artificially constructed data for which you know all the TPs.

You could use the propTrueNull() function of the statmod package to estimate the total number of null hypotheses that are actually true, and that would give you some idea of now many false nulls have failed to reach your significance cutoff. That does rely on the p-values being uniformly distributed for truly non-DE genes, which will never be exactly true.

As Michael has said, you can get a general idea from a power analysis of what FNR might be typical for an experiment like yours, given assumptions about effect size etc, but that will not yield an estimate of FNR for any specific dataset.

ADD COMMENT

Login before adding your answer.

Traffic: 600 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6