R/bFDR.R

In MRS: Multi-Resolution Scanning for Cross-Sample and Cross-Group Differences

bFDR <- function(ans_mrs, threshold) {
  PMAPs <- ans_mrs$RepresentativeTree$AltProbs
  selected <- PMAPs > threshold
  if (!any(selected)) {
    return(0)
  }
  mean(1 - PMAPs[selected])
}