R/levelSceptical.R

In ReplicationSuccess: Design and Analysis of Replication Studies

levelSceptical <- function(level, 
                           alternative = "one.sided", 
                           type = "golden"){
    if (!(type %in% c("nominal", "liberal", "controlled", "golden")))
        stop('type must be either "nominal", "liberal", "controlled", or "golden"')
    if (!(alternative %in% c("one.sided", "two.sided")))
        stop('type must be either "one.sided" or "two.sided"')
    
    if(type == "nominal")
        res <- level
    
    if(type == "liberal")
        res <- pIntrinsic(p = level, alternative = alternative, type = "Held")
    
    if(type == "controlled"){
        if (alternative == "two.sided") {
            t1 <- level^2 ## level is two-sided significance level
            # t1 <- (2*level)^2 ## level is a one-sided significance level
            res <- 2*(1 - pnorm(q = qnorm(p = 1 - t1/2)/2))
        } 
        if (alternative == "one.sided") {
            # t1 <- level*(level/2) ## level is a two-sided significance level
            t1 <- 2*level^2 ## level is a one-sided significance level
            res <- 1 - pnorm(q = qnorm(p = 1 - t1)/2)
        }
        if (alternative == "less") {
            ## TODO implement this
        }
        if (alternative == "greater") {
            ## TODO implement this
        }
        # if(alternative == "two.sided") level <- level/2
        # pI2 <- pIntrinsic(p = 2*level^2, alternative = "one.sided", type = "Held")
        # res <- pIntrinsic(p = pI2, alternative = "one.sided", type = "Held")
        # if(alternative == "two.sided") res <- 2*res
    }
    
    if(type == "golden"){
        res <- pIntrinsic(p = level, alternative = alternative, type = "Matthews")
    }
    
    return(res)
}