h1Proportion <- function (pv, method = "storey", lambda = 0.5)
{
method = match.arg(method)
pi1 = 1 - mean(pv > lambda, na.rm = TRUE)/(1 - lambda)
if (pi1 < 0) {
warning(paste("estimated pi1 =", round(pi1, digits = 4),
"set to 0.0"))
pi1 = 0
}
if (pi1 > 1) {
warning(paste("estimated pi1 =", round(pi1, digits = 4),
"set to 1.0"))
pi1 = 1
}
return(pi1)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.