R/pi0_estimators.R

Defines functions weighted_storey_pi0

weighted_storey_pi0 <- function(pvalues, weights, tau=0.5, m = length(pvalues)){
    w_inf <- max(weights)
    num <- w_inf + sum( weights * (pvalues > tau))
    num/m/(1-tau)
}
nignatiadis/IHW documentation built on Aug. 22, 2023, 2:11 p.m.