tests/testthat/test_ci_1pop_chisq_sd.R

In statBasics: Basic Functions to Statistical Methods Course

size <- 1000
mean_pop <- 3
sd_pop <- 2
x <- rnorm(size, mean = mean_pop, sd = sd_pop)

###############################################################################
# type = 'two.sided'

lower <- function(conf_level, type = 'two.sided') {
  if (type == 'two.sided') {
    suppressWarnings({
      return(unname(EnvStats::varTest(x, conf.level = conf_level)$conf.int[1]))
    })
  } else if (type == 'left') {
    suppressWarnings({
      return(unname(EnvStats::varTest(x, conf.level = conf_level, alternative = "greater")$conf.int[1]))
    })
  } else if (type == 'right') {
    return(0)
  }
}

upper <- function(conf_level, type = 'two.sided') {
  if (type == 'two.sided') {
    suppressWarnings({
      return(unname(EnvStats::varTest(x, conf.level = conf_level)$conf.int[2]))
    })
  } else if (type == 'left') {
    return(Inf)
  } else if (type == 'right') {
    suppressWarnings({
      return(unname(EnvStats::varTest(x, conf.level = conf_level, alternative = 'less')$conf.int[2]))
    })
  }
}

testthat::test_that("Comparing, type = 'two.sided', conf_level = '90%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90)$lower_ci, lower(0.90))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90)$upper_ci, upper(0.90))
})

testthat::test_that("Comparing, type = 'two.sided', conf_level = '95%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95)$lower_ci, lower(0.95))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95)$upper_ci, upper(0.95))
})

testthat::test_that("Comparing, type = 'two.sided', conf_level = '99%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99)$lower_ci, lower(0.99))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99)$upper_ci, upper(0.99))
})

###############################################################################
# type = 'left'

testthat::test_that("Comparing, type = 'left', conf_level = '90%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90, type = 'left')$lower_ci, lower(0.90, type = 'left'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90, type = 'left')$upper_ci, upper(0.90, type = 'left'))
})

testthat::test_that("Comparing, type = 'two.sides', conf_level = '95%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95, type = 'left')$lower_ci, lower(0.95, type = 'left'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95, type = 'left')$upper_ci, upper(0.95, type = 'left'))
})

testthat::test_that("Comparing, type = 'two.sides', conf_level = '99%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99, type = 'left')$lower_ci, lower(0.99, type = 'left'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99, type = 'left')$upper_ci, upper(0.99, type = 'left'))
})

###############################################################################
# type = 'right'

testthat::test_that("Comparing, type = 'right', conf_level = '90%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90, type = 'right')$lower_ci, lower(0.90, type = 'right'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.90, type = 'right')$upper_ci, upper(0.90, type = 'right'))
})

testthat::test_that("Comparing, type = 'right', conf_level = '95%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95, type = 'right')$lower_ci, lower(0.95, type = 'right'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.95, type = 'right')$upper_ci, upper(0.95, type = 'right'))
})

testthat::test_that("Comparing, type = 'right', conf_level = '99%'", {
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99, type = 'right')$lower_ci, lower(0.99, type = 'right'))
  testthat::expect_equal(ci_1pop_norm(x, parameter = 'variance', conf_level = 0.99, type = 'right')$upper_ci, upper(0.99, type = 'right'))
})