R/F_calcWeights.R

In reconsi: Resampling Collapsed Null Distributions for Simultaneous Inference

#' Obtain weights as posterior probabilities to calculate the consensus null
#'
#' @param logDensPerm A matrix with B rows  of logged density estimates of the B
#'    permutation distributions, and p columns for the p observed test
#'    statistics
#' @param fdr A vector of local false discovery rates for the observed tests
#'  statistics of length p
#'
#' @return A vector of weights of length B
calcWeights = function(logDensPerm, fdr){
    weights = exp(stabExp(colSums(logDensPerm*fdr)))
    weights/sum(weights)
}