tests/testthat/test-datasets-communes.R

In vdes2020/eneRgyVD: Municipal statistics for climate-related policies, Canton de Vaud (Switzerland)

# These tests aim at checking that supplied communes are valid,
#  and that the sum of all communes equals the `Canton de Vaud` entry.


test_that("Check valid `commune` variable in every dataset (choices_canton_communes in utils_helpers.R)",

          {
            # Define vector of clean choices (requires unlisting)
            values_canton_communes <- unlist(choices_canton_communes, use.names = FALSE)


            purrr::map(energy_datasets, \(df) expect_in(unique(df$commune), values_canton_communes))
            purrr::map(mobility_datasets, \(df) expect_in(unique(df$commune), values_canton_communes))
            purrr::map(adaptation_datasets, \(df) expect_in(unique(df$commune), values_canton_communes))


          })

test_that("Check missing 'commune' in every dataset so we explicitely plot/table communes without values.
          For good reasons (complexity/confusion) we allow some datasets to not pass this test.",
          {
            # Define vector of clean choices (requires unlisting)
            values_canton_communes <- unlist(choices_canton_communes, use.names = FALSE)


            # for energy : we allow subsidies_* to remain with missing communes, too many combinations etc. to handle otherwise before the app
            purrr::map(energy_datasets[!stringr::str_detect(names(energy_datasets), pattern = "subsidies")], \(df) expect_in(values_canton_communes, unique(df$commune)))
            purrr::map(mobility_datasets, \(df) expect_in(values_canton_communes, unique(df$commune)))
            purrr::map(adaptation_datasets, \(df) expect_in(values_canton_communes, unique(df$commune)))

          })


test_that("Check for each dataset that the value for 'Canton de Vaud' equals the sum of all other `commune` entries",
          {

            # !!!! ELEC_PROD not valid because of (Confidentiel) inserts !!!! #


            # expect_equal(label = "The sum of 'Canton de Vaud' (GWh)",
            #              expected.label = "the sum of all communes. This is expected as long as there's filter for N<3 installations in the raw-data -> .rda script",
            #              sum(elec_prod$production[elec_prod$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
            #              sum(elec_prod$production[elec_prod$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            # )

            expect_equal(
              sum(energy_datasets$elec_cons$consommation[energy_datasets$elec_cons$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$elec_cons$consommation[energy_datasets$elec_cons$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$ng_cons$consommation[energy_datasets$ng_cons$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$ng_cons$consommation[energy_datasets$ng_cons$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$regener_cons_ae_aff$consommation[regener_cons_ae_aff$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$regener_cons_ae_aff$consommation[regener_cons_ae_aff$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$regener_cons_ae_aff$consommation[regener_cons_ae_aff$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$regener_cons_ae_aff$consommation[regener_cons_ae_aff$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$regener_cons_ae_use$consommation[regener_cons_ae_use$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$regener_cons_ae_use$consommation[regener_cons_ae_use$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$regener_needs$besoins[regener_needs$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$regener_needs$besoins[regener_needs$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            # subsidies datasets ----

            expect_equal(
              sum(energy_datasets$subsidies_by_building$N_EGID[subsidies_by_building$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$subsidies_by_building$N_EGID[subsidies_by_building$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$subsidies_by_building$SRE[subsidies_by_building$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$subsidies_by_building$SRE[subsidies_by_building$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )

            expect_equal(
              sum(energy_datasets$subsidies_by_measure$nombre[subsidies_by_measure$commune == "Canton de Vaud"], na.rm = TRUE)/1e6,
              sum(energy_datasets$subsidies_by_measure$nombre[subsidies_by_measure$commune != "Canton de Vaud"], na.rm = TRUE)/1e6
            )


          })