test_that("Optimal_bias works for additive adjustment method", {
skip_on_cran()
expect_equal(optimal_bias(w = 0.3, hr1 = 0.69, hr2 = 0.88,
id1 = 210, id2 = 420,
d2min = 20, d2max = 100, stepd2 = 5,
hrgomin = 0.7, hrgomax = 0.9, stephrgo = 0.05,
adj = "additive",
lambdamin = 0.6, lambdamax = 1, steplambda = 0.05,
alphaCImin = 0.25, alphaCImax = 0.5, stepalphaCI = 0.025,
alpha = 0.025, beta = 0.1, xi2 = 0.7, xi3 = 0.7,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE, num_cl = 2)$u, 859.71)
})
test_that("Optimal_bias works for multiplicative method", {
skip_on_cran()
expect_equal(optimal_bias(w = 0.3, hr1 = 0.69, hr2 = 0.88,
id1 = 210, id2 = 420,
d2min = 20, d2max = 100, stepd2 = 5,
hrgomin = 0.7, hrgomax = 0.9, stephrgo = 0.05,
adj = "multiplicative",
lambdamin = 0.6, lambdamax = 1, steplambda = 0.05,
alphaCImin = 0.025, alphaCImax = 0.5, stepalphaCI = 0.025,
alpha = 0.025, beta = 0.1, xi2 = 0.7, xi3 = 0.7,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = FALSE, num_cl = 2)$u, 98.28)
})
test_that("Optimal_bias works for both methods", {
skip_on_cran()
expect_equal(optimal_bias(w = 0.3, hr1 = 0.69, hr2 = 0.88,
id1 = 210, id2 = 420,
d2min = 50, d2max = 150, stepd2 = 5,
hrgomin = 0.7, hrgomax = 0.9, stephrgo = 0.05,
adj = "both",
lambdamin = 0.6, lambdamax = 1, steplambda = 0.05,
alphaCImin = 0.025, alphaCImax = 0.5, stepalphaCI = 0.025,
alpha = 0.025, beta = 0.1, xi2 = 0.7, xi3 = 0.7,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE, num_cl = 2)$u, c(972.87,922.63))
})
test_that("Optimal_bias works for method `all`", {
skip_on_cran()
expect_equal(optimal_bias(w = 0.3, hr1 = 0.69, hr2 = 0.88,
id1 = 210, id2 = 420,
d2min = 50, d2max = 150, stepd2 = 5,
hrgomin = 0.7, hrgomax = 0.9, stephrgo = 0.05,
adj = "all",
lambdamin = 0.6, lambdamax = 1, steplambda = 0.05,
alphaCImin = 0.25, alphaCImax = 0.5, stepalphaCI = 0.025,
alpha = 0.025, beta = 0.1, xi2 = 0.7, xi3 = 0.7,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE, num_cl = 2)$u, c(972.87,922.63,966.59,921.11))
})
test_that("Optimal_bias works for method `all` with prior distributions", {
skip_on_cran()
expect_equal(optimal_bias(w = 0.3, hr1 = 0.69, hr2 = 0.88,
id1 = 210, id2 = 420,
d2min = 50, d2max = 150, stepd2 = 5,
hrgomin = 0.7, hrgomax = 0.9, stephrgo = 0.05,
adj = "all",
lambdamin = 0.6, lambdamax = 1, steplambda = 0.05,
alphaCImin = 0.25, alphaCImax = 0.5, stepalphaCI = 0.025,
alpha = 0.025, beta = 0.1, xi2 = 0.7, xi3 = 0.7,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = FALSE, num_cl = 2)$u, c(98.28,75.95,96.69,77.40))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.