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)
}

Try the IHW package in your browser

Any scripts or data that you put into this service are public.

IHW documentation built on Nov. 8, 2020, 7:44 p.m.