tests/testthat/test-prepare_mdi.R

In healthiar: Quantifying and Monetizing Health Impacts Attributable to Exposure

testthat::test_that("results correct", { # with verbose = FALSE

  exdat_prepare_mdi <- base::readRDS(testthat::test_path("data", "exdat_get_mdi.rds"))

  testthat::expect_equal(
    object =
      healthiar::prepare_mdi(
        geo_id_micro = exdat_prepare_mdi$id,
        edu = exdat_prepare_mdi$edu,
        unemployed = exdat_prepare_mdi$unemployed,
        single_parent = exdat_prepare_mdi$single_parent,
        pop_change = exdat_prepare_mdi$pop_change,
        no_heating = exdat_prepare_mdi$no_heating,
        n_quantile = 10,
        verbose = FALSE)$mdi_main$MDI[1:50] |> base::round(digits = 7), # Extract the MDI of the first 50 geo units
    expect = c(0.2117721, 0.4319924, 0.1847750, 0.3787937, 0.3121354, 0.2565185, 0.2245822,
               0.2140148, 0.2656597, 0.3566141, 0.1746261, 0.2108900, 0.2221277, 0.2226641,
               0.2478800, 0.2513628, 0.3595167, 0.2321278, 0.2879502, 0.2962355, 0.1866685,
               0.3168898, 0.2928375, 0.2658067, 0.2646959, 0.2306425, 0.2273271, 0.1892787,
               0.3941645, 0.2268015, 0.2973753, 0.1853536, 0.2850227, 0.2642978, 0.2867885,
               0.2775215, 0.3538885, 0.2835787, 0.2554690, 0.2423936, 0.3548941, 0.2564889,
               0.3041343, 0.1991622, 0.2355921, 0.2671123, 0.2708377, 0.2348108, 0.2902528,
               0.3112479) # Results on 2025-04-15; Sciensano results
  )
})

## COMMENTED OUT SO NOTHING IS PRINTED TO CONSOLE DURING OVERALL TESTING
# testthat::test_that("results correct", { # with verbose = TRUE
#
#   # exdat_prepare_mdi <- base::readRDS(testthat::test_path("data", "exdat_get_mdi.rds"))
#
#   testthat::expect_equal(
#     object =
#       healthiar::prepare_mdi(
#         geo_id_micro = exdat_prepare_mdi$id,
#         edu = exdat_prepare_mdi$edu,
#         unemployed = exdat_prepare_mdi$unemployed,
#         single_parent = exdat_prepare_mdi$single_parent,
#         pop_change = exdat_prepare_mdi$pop_change,
#         no_heating = exdat_prepare_mdi$no_heating,
#         n_quantile = 10,
#         verbose = TRUE)$mdi_main$MDI[1:50] |> base::round(digits = 7), # Extract the MDI of the first 50 geo units
#     expect = c(0.2117721, 0.4319924, 0.1847750, 0.3787937, 0.3121354, 0.2565185, 0.2245822,
#                0.2140148, 0.2656597, 0.3566141, 0.1746261, 0.2108900, 0.2221277, 0.2226641,
#                0.2478800, 0.2513628, 0.3595167, 0.2321278, 0.2879502, 0.2962355, 0.1866685,
#                0.3168898, 0.2928375, 0.2658067, 0.2646959, 0.2306425, 0.2273271, 0.1892787,
#                0.3941645, 0.2268015, 0.2973753, 0.1853536, 0.2850227, 0.2642978, 0.2867885,
#                0.2775215, 0.3538885, 0.2835787, 0.2554690, 0.2423936, 0.3548941, 0.2564889,
#                0.3041343, 0.1991622, 0.2355921, 0.2671123, 0.2708377, 0.2348108, 0.2902528,
#                0.3112479) # Results on 2025-04-15; Sciensano results
#   )
# })