Nothing
test_that("Output of function 'powerSignificance' stays the same.", {
zo <- seq(-1, 1, 1)
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(grid)),
FUN = function(i) {
powerSignificance(
zo = zo,
c = grid$c[i],
level = 0.05,
designPrior = grid$priors[i],
alternative = grid$alt[i],
h = grid$h[i],
shrinkage = 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)
)
)
})
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)
})
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.