R/rosenthal.R

In fsn: Rosenthal's Fail Safe Number and Related Functions

Defines functions rosenthal

Documented in rosenthal

rosenthal <- function(stat, se, alpha = 0.05) {
  k <- length(stat) ## how many studies you have
  Nr <- sum( stat/se )^2 / qnorm(1 - alpha)^2 - k
  rule <- 5 * k + 10
  res <- c(Nr, rule)
  names(res) <- c('fail-safe Nr', 'rule')
  res
}

Try the fsn package in your browser