Nothing
meanz <-
function(p, log.p = FALSE) {
keep <- (p > 0) & (p < 1)
invalid <- sum(1L * keep) < 2
if(invalid) {
warning("Must have at least two valid p values")
res <- list(z = NA_real_, p = NA_real_,
validp = p[keep])
} else {
if(sum(1L * keep) != length(p)) {
warning("Some studies omitted")
}
zvals <- (qnorm(p[keep], lower.tail = FALSE))
zp <- mean(zvals) / (sd(zvals) / sqrt(length(p[keep])))
res <- list(z = zp, p = pnorm(zp, lower.tail = FALSE,
log.p = log.p), validp = p[keep])
}
class(res) <- c("meanz", "metap")
res
}
print.meanz <- function(x, ...) {
cat("meanz = ", x$z, "p = ", x$p, "\n")
invisible(x)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.