R/z2p.R

In ReplicationSuccess: Design and Analysis of Replication Studies

z2p <- function(z, 
                alternative = "two.sided"){
    
    # vectorize function in all arguments
    pV <- mapply(FUN = function(z, alternative) {
        if (!is.numeric(z))
            stop("z must be numeric")
        if (!(alternative %in% c("less", "greater", "two.sided", "one.sided")))
            stop('alternative must be either "less", "greater", "two.sided", or "one.sided"')
        
        if (alternative == "two.sided")
            p <- 2*pnorm(abs(z), lower.tail = FALSE)
        
        if (alternative == "less")
            p <- pnorm(q = z, lower.tail = TRUE)
        
        if (alternative == "greater" | alternative == "one.sided")
            p <- pnorm(q = z, lower.tail = FALSE)
        
        return(p)    
    }, z, alternative)
    return(pV)
}