tests/testthat/test-powerSignificance.R

In florafauna/RStest: Design and Analysis of Replication Studies

## library(testthat)
## sapply(list.files("../../R", pattern='\\.R$', full.names = TRUE), source)

context("powerSignificance")

test_that("numeric test for powerSignificance(): 1", {
    expect_equal(object = powerSignificance(zo = qnorm(p = 1 - 0.05/2),
                                            c = 1, level = 0.05,
                                            alternative = "two.sided"),
                 expected = 0.5,
                 tol = 0.0001)
    expect_equal(object = powerSignificance(zo = qnorm(p = 1 - 0.0056/2),
                                            c = 1, level = 0.05,
                                            alternative = "two.sided"),
                 expected = 0.791,
                 tol = 0.001)
    
})


test_that("numeric test for powerSignificance(): 2", {
    zo <- seq(-1, 1, 1)
    apply_grid <- expand.grid(priors = c("conditional", "predictive", "EB"),
                              c = c(0.5, 2),
                              h = c(0, 1),
                              alt = c("two.sided", "one.sided"),
                              shrinkage = c(0, 0.5),
                              stringsAsFactors = FALSE)
    out <- lapply(X=seq_len(nrow(apply_grid)), FUN=function(i){
        powerSignificance(zo = zo,
                          c = apply_grid$c[i],
                          level = 0.05,
                          designPrior = apply_grid$priors[i],
                          alternative = apply_grid$alt[i],
                          h = apply_grid$h[i],
                          shrinkage = apply_grid$shrinkage[i])
    })
    expect_equal(out,
                 list(c(0.105128839609805, 0.025, 0.105128839609805),
                      c(0.153164887002923, 0.0547655919261876,
                        0.153164887002923),
                      c(0.025, 0.025, 0.025),
                      c(0.29261875345421, 0.025, 0.29261875345421),
                      c(0.376346981956041, 0.128904320268035, 0.376346981956041),
                      c(0.025, 0.025, 0.025),
                      c(0.105128839609805, 0.025, 0.105128839609805),
                      c(0.214070597849871, 0.107563528680747, 0.214070597849871),
                      c(0.0547655919261876, 0.0547655919261876, 0.0547655919261876),
                      c(0.29261875345421, 0.025, 0.29261875345421),
                      c(0.418288337676037, 0.229408341649887, 0.418288337676037),
                      c(0.128904320268035, 0.128904320268035, 0.128904320268035),
                      c(0.174187261617932, 0.05, 0.174187261617932),
                      c(0.221937205879032, 0.08963317592404, 0.221937205879032),
                      c(0.05, 0.05, 0.05),
                      c(0.40879721979387, 0.05, 0.40879721979387),
                      c(0.447033381345704, 0.171143362995444, 0.447033381345704),
                      c(0.05, 0.05, 0.05),
                      c(0.174187261617932, 0.05, 0.174187261617932),
                      c(0.276562743130718, 0.149101020851907, 0.276562743130718),
                      c(0.08963317592404, 0.08963317592404, 0.08963317592404),
                      c(0.40879721979387, 0.05, 0.40879721979387),
                      c(0.465266702002727, 0.26707081412311, 0.465266702002727),
                      c(0.171143362995444, 0.171143362995444, 0.171143362995444),
                      c(0.0540918623618536, 0.025, 0.0540918623618536),
                      c(0.0948227125608762, 0.0547655919261876, 0.0948227125608762),
                      c(0.025, 0.025, 0.025),
                      c(0.105128839609805, 0.025, 0.105128839609805),
                      c(0.234736296915959, 0.128904320268035, 0.234736296915959),
                      c(0.025, 0.025, 0.025),
                      c(0.0540918623618536, 0.025, 0.0540918623618536),
                      c(0.154818677461375, 0.107563528680747, 0.154818677461375),
                      c(0.0547655919261876, 0.0547655919261876,
                        0.0547655919261876),
                      c(0.105128839609805, 0.025, 0.105128839609805),
                      c(0.317915585520223, 0.229408341649887, 0.317915585520223),
                      c(0.128904320268035, 0.128904320268035, 0.128904320268035),
                      c(0.0982997943579869, 0.05, 0.0982997943579869),
                      c(0.145863130961703, 0.08963317592404, 0.145863130961703),
                      c(0.05, 0.05, 0.05),
                      c(0.174187261617932, 0.05, 0.174187261617932),
                      c(0.294113061454605, 0.171143362995444, 0.294113061454605),
                      c(0.05, 0.05, 0.05),
                      c(0.0982997943579869, 0.05, 0.0982997943579869),
                      c(0.207052810036395, 0.149101020851907, 0.207052810036395),
                      c(0.08963317592404, 0.08963317592404, 0.08963317592404),
                      c(0.174187261617932, 0.05, 0.174187261617932),
                      c(0.361506457143659, 0.26707081412311, 0.361506457143659),
                      c(0.171143362995444, 0.171143362995444, 0.171143362995444)))
})