tests/testthat/test_overview.R

In rcalinjageman/esci: Estimation Statistics with Confidence Intervals

test_that("overview to esci examples", {

  # Setup
  mydata <- esci::data_halagappa

  # Check - same result as in esci with 95% CI
  res <- overview(
    means = mydata$Mean,
    sds = mydata$SD,
    ns = mydata$n,
    assume_equal_variance = TRUE,
    outcome_variable_name = "% time near target",
    grouping_variable_levels = as.character(mydata$Groups),
    grouping_variable_name = "Diet"
  )

  # Should return mean_LL of 32.32519301	26.72519301	36.02519301	28.22519301	24.72519301	33.12519301
  # Should return mean_UL of 42.67480699	37.07480699	46.37480699	38.57480699	35.07480699	43.47480699


  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_LL"], 32.32519301)
  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_UL"], 42.67480699)

  testthat::expect_equal(res[res$grouping_variable_level == "AFree10", "mean_LL"], 26.72519301)
  testthat::expect_equal(res[res$grouping_variable_level == "AFree10", "mean_UL"], 37.07480699)

  testthat::expect_equal(res[res$grouping_variable_level == "ADiet10", "mean_LL"], 36.02519301)
  testthat::expect_equal(res[res$grouping_variable_level == "ADiet10", "mean_UL"], 46.37480699)

  testthat::expect_equal(res[res$grouping_variable_level == "NFree17", "mean_LL"], 28.22519301)
  testthat::expect_equal(res[res$grouping_variable_level == "NFree17", "mean_UL"], 38.57480699)

  testthat::expect_equal(res[res$grouping_variable_level == "AFree17", "mean_LL"], 24.72519301)
  testthat::expect_equal(res[res$grouping_variable_level == "AFree17", "mean_UL"], 35.07480699)

  testthat::expect_equal(res[res$grouping_variable_level == "ADiet17", "mean_LL"], 33.12519301)
  testthat::expect_equal(res[res$grouping_variable_level == "ADiet17", "mean_UL"], 43.47480699)


  # Check - same result as in esci with 99% CI
  res <- overview(
    means = mydata$Mean,
    sds = mydata$SD,
    ns = mydata$n,
    assume_equal_variance = TRUE,
    outcome_variable_name = "% time near target",
    grouping_variable_levels = as.character(mydata$Groups),
    grouping_variable_name = "Diet",
    conf_level = .99
  )

  # Should return mean_LL of 30.65450169	25.05450169	34.35450169	26.55450169	23.05450169	31.45450169
  # Should return mean_UL of 44.34549831	38.74549831	48.04549831	40.24549831	36.74549831	45.14549831

  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_LL"], 30.65450169)
  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_UL"], 44.34549831)

  testthat::expect_equal(res[res$grouping_variable_level == "AFree10", "mean_LL"], 25.05450169)
  testthat::expect_equal(res[res$grouping_variable_level == "AFree10", "mean_UL"], 38.74549831)

  testthat::expect_equal(res[res$grouping_variable_level == "ADiet10", "mean_LL"], 34.35450169)
  testthat::expect_equal(res[res$grouping_variable_level == "ADiet10", "mean_UL"], 48.04549831)

  testthat::expect_equal(res[res$grouping_variable_level == "NFree17", "mean_LL"], 26.55450169)
  testthat::expect_equal(res[res$grouping_variable_level == "NFree17", "mean_UL"], 40.24549831)

  testthat::expect_equal(res[res$grouping_variable_level == "AFree17", "mean_LL"], 23.05450169)
  testthat::expect_equal(res[res$grouping_variable_level == "AFree17", "mean_UL"], 36.74549831)

  testthat::expect_equal(res[res$grouping_variable_level == "ADiet17", "mean_LL"], 31.45450169)
  testthat::expect_equal(res[res$grouping_variable_level == "ADiet17", "mean_UL"], 45.14549831)


  #Ok for just one group
  res <- overview(
    means = mydata$Mean[1],
    sds = mydata$SD[1],
    ns = mydata$n[1],
    assume_equal_variance = TRUE,
    outcome_variable_name = "% time near target",
    grouping_variable_levels = as.character(mydata$Groups)[1],
    grouping_variable_name = "Diet"
  )

  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_LL"], 32.68015428)
  testthat::expect_equal(res[res$grouping_variable_level == "NFree10", "mean_UL"], 42.31984572)


  # data frame
  mydata <- esci::data_penlaptop1

  res <- overview(
    data = mydata,
    outcome_variable = transcription,
    grouping_variable = condition,
    assume_equal_variance = FALSE
  )

  testthat::expect_equal(res[res$grouping_variable_level == "Pen", "mean"], 8.81176, tolerance = 1e-5)
  testthat::expect_equal(res[res$grouping_variable_level == "Pen", "mean_LL"], 7.15464, tolerance = 1e-5)
  testthat::expect_equal(res[res$grouping_variable_level == "Pen", "mean_UL"], 10.4689, tolerance = 1e-5)

  testthat::expect_equal(res[res$grouping_variable_level == "Laptop", "mean"], 14.5194, tolerance = 1e-5)
  testthat::expect_equal(res[res$grouping_variable_level == "Laptop", "mean_LL"], 11.847, tolerance = 1e-5)
  testthat::expect_equal(res[res$grouping_variable_level == "Laptop", "mean_UL"], 17.1917, tolerance = 1e-5)


})