tests/testthat/test_ci_1pop_bern_success_number_success.R

In statBasics: Basic Functions to Statistical Methods Course

testthat::test_that("Number of sucess and confidence level 90% - type = 'two.sided'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90)$lower_ci, qnorm(0.05) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90)$upper_ci, qnorm(0.95) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90)$lower_ci, qnorm(0.05) / (2 * sqrt(100)) + 50 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90)$upper_ci, qnorm(0.95) / (2 * sqrt(100)) + 50 / 100)
})

testthat::test_that("Number of sucess and confidence level 90% - type = 'left'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90, type = 'left')$lower_ci, qnorm(0.1) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90, type = 'left')$upper_ci, 1)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90, type = 'left')$lower_ci, qnorm(0.1) / (2 * sqrt(100)) + 50 / 100 )
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90, type = 'left')$upper_ci, 1)
})

testthat::test_that("Number of sucess and confidence level 90% - type = 'right'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.90, type = 'right')$upper_ci, qnorm(0.9) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.90, type = 'right')$upper_ci, qnorm(0.9) / (2 * sqrt(100)) + 50 / 100)
})

#################################################################################################
# two.sided, and number of sucess
# conf_level = 95%

testthat::test_that("Number of sucess and confidence level 95% - type = 'two.sided'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95)$lower_ci, qnorm(0.025) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95)$upper_ci, qnorm(0.975) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95)$lower_ci, qnorm(0.025) / (2 * sqrt(100)) + 50 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95)$upper_ci, qnorm(0.975) / (2 * sqrt(100)) + 50 / 100)
})

testthat::test_that("Number of sucess and confidence level 95% - type = 'left'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95, type = 'left')$lower_ci, qnorm(0.05) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95, type = 'left')$upper_ci, 1)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95, type = 'left')$lower_ci, qnorm(0.05) / (2 * sqrt(100)) + 50 / 100 )
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95, type = 'left')$upper_ci, 1)
})

testthat::test_that("Number of sucess and confidence level 95% - type = 'right'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.95, type = 'right')$upper_ci, qnorm(0.95) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.95, type = 'right')$upper_ci, qnorm(0.95) / (2 * sqrt(100)) + 50 / 100)
})

#################################################################################################
# two.sided, and number of sucess
# conf_level = 99%

testthat::test_that("Number of sucess and confidence level 99% - type = 'two.sided'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99)$lower_ci, qnorm(0.005) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99)$upper_ci, qnorm(0.995) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99)$lower_ci, qnorm(0.005) / (2 * sqrt(100)) + 50 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99)$upper_ci, qnorm(0.995) / (2 * sqrt(100)) + 50 / 100)
})

testthat::test_that("Number of sucess and confidence level 99% - type = 'left'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99, type = 'left')$lower_ci, qnorm(0.01) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99, type = 'left')$upper_ci, 1)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99, type = 'left')$lower_ci, qnorm(0.01) / (2 * sqrt(100)) + 50 / 100 )
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99, type = 'left')$upper_ci, 1)
})

testthat::test_that("Number of sucess and confidence level 99% - type = 'right'", {
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(20, 100, conf_level = 0.99, type = 'right')$upper_ci, qnorm(0.99) / (2 * sqrt(100)) + 20 / 100)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99, type = 'right')$lower_ci, 0)
    testthat::expect_equal(ci_1pop_bern(50, 100, conf_level = 0.99, type = 'right')$upper_ci, qnorm(0.99) / (2 * sqrt(100)) + 50 / 100)
})