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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.