compute_holm <- function(statistic = NULL,pvalue){
n <- length(pvalue)
i <- seq_len(n)
o <- order(pvalue)
ro <- order(o)
pvalue <- pmin(1, cummax((n - i + 1L) * pvalue[o]))[ro]
out=list(main = cbind(statistic = statistic,pvalue = pvalue))
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.