Skip to content

Commit

Permalink
remove naomi-spectrum comparison check from navigator code
Browse files Browse the repository at this point in the history
  • Loading branch information
rtesra committed Dec 11, 2024
1 parent 25e6565 commit d07a38f
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 70 deletions.
65 changes: 2 additions & 63 deletions R/model.R
Original file line number Diff line number Diff line change
Expand Up @@ -570,7 +570,7 @@ naomi_model_frame <- function(area_merged,
spec_unaware_untreated_prop_t4 = unaware_untreated_prop,
asfr_t4 = asfr,
frr_plhiv_t4 = frr_plhiv,
frr_already_art_t4 = frr_already_art
frr_already_art_t4 = frr_already_art
),
by = c("spectrum_region_code", "sex", "age_group")
) %>%
Expand Down Expand Up @@ -881,65 +881,6 @@ select_naomi_data <- function(

stopifnot(methods::is(naomi_mf, "naomi_mf"))

## Check anc_testing and art_number against Spectrum inputs.
## Return NA if spec_program_data not provided
anc_testing_spectrum_aligned <- NA
art_number_spectrum_aligned <- NA

if (!is.null(spec_program_data)) {
stopifnot(methods::is(spec_program_data, "spec_program_data"))

if (!is.null(anc_testing)) {

anc_merged <- anc_testing %>%
dplyr::left_join(
dplyr::select(naomi_mf$mf_areas, area_id, spectrum_region_code),
by = "area_id"
) %>%
tidyr::pivot_longer(dplyr::starts_with("anc"),
names_to = "indicator",
values_to = "value_naomi") %>%
dplyr::count(spectrum_region_code, year, indicator,
wt = value_naomi, name = "value_naomi") %>%
dplyr::inner_join(
spec_program_data$anc_testing %>%
dplyr::rename("value_spectrum" = "value"),
by = c("spectrum_region_code", "indicator", "year")
)

anc_testing_spectrum_aligned <- all(anc_merged$value_naomi == anc_merged$value_spectrum)

} else {
## If no ANC testing data, return TRUE
anc_testing_spectrum_aligned <- TRUE
}

if (!is.null(art_number)) {

art_merged <- art_number %>%
dplyr::left_join(
dplyr::select(naomi_mf$mf_areas, area_id, spectrum_region_code),
by = "area_id"
) %>%
dplyr::count(spectrum_region_code, sex, age_group, calendar_quarter,
wt = art_current, name = "art_current_naomi") %>%
dplyr::inner_join(
spec_program_data$art_dec31 %>%
dplyr::mutate(
calendar_quarter = paste0("CY", year, "Q4"),
year = NULL
),
by = c("spectrum_region_code", "sex", "age_group", "calendar_quarter")
)

art_number_spectrum_aligned <- all(art_merged$art_current_naomi == art_merged$art_dec31)

} else {
## If no ANC testing data, return TRUE
art_number_spectrum_aligned <- TRUE
}
}

common_surveys <- intersect(artcov_survey_ids, vls_survey_ids)
if (length(common_surveys)) {
stop(t_("ART_COV_AND_VLS_SAME_SURVEY",
Expand Down Expand Up @@ -1101,9 +1042,7 @@ select_naomi_data <- function(
artnum_calendar_quarter_t1 = artnum_calendar_quarter_t1,
artnum_calendar_quarter_t2 = artnum_calendar_quarter_t2,
anc_prev_year_t1 = anc_artcov_year_t1,
anc_prev_year_t2 = anc_artcov_year_t2,
art_number_spectrum_aligned = art_number_spectrum_aligned,
anc_testing_spectrum_aligned = anc_testing_spectrum_aligned)
anc_prev_year_t2 = anc_artcov_year_t2)

naomi_mf$data_options <- data_options

Expand Down
4 changes: 0 additions & 4 deletions R/unaids-navigator.R
Original file line number Diff line number Diff line change
Expand Up @@ -91,10 +91,6 @@ write_navigator_checklist <- function(naomi_output,

if (!is.null(data_options)) {

## Compare aggregated naomi inputs to national Spectrum totals
v$TrueFalse[v$NaomiCheckPermPrimKey == "ART_is_Spectrum"] <- data_options$art_number_spectrum_aligned
v$TrueFalse[v$NaomiCheckPermPrimKey == "ANC_is_Spectrum"] <- data_options$anc_testing_spectrum_aligned

## Check that all surveys used are from most recent quarter available
most_recent_survey_available <- max(data_options$prev_survey_available_quarters)
v$TrueFalse[v$NaomiCheckPermPrimKey == "Opt_recent_survey_only"] <-
Expand Down
4 changes: 1 addition & 3 deletions tests/testthat/test-outputs.R
Original file line number Diff line number Diff line change
Expand Up @@ -165,7 +165,7 @@ test_that("subset_output_package() saves expected output package", {


## Test that can **drop** selected indicators (rather than keep)

sub_drop_file <- tempfile(fileext = ".zip")

expect_warning(
Expand Down Expand Up @@ -328,8 +328,6 @@ test_that("navigator checklist returns expected results", {

adj_output$fit$data_options$prev_survey_ids <- "DEMO2020PHIA"
adj_output$fit$data_options$prev_survey_quarters <- "CY2020Q3"
adj_output$fit$data_options$art_number_spectrum_aligned <- TRUE
adj_output$fit$data_options$anc_testing_spectrum_aligned <- TRUE

adj_output$fit$calibration_options$spectrum_population_calibration <- "subnational"
adj_output$fit$calibration_options$spectrum_artnum_calibration_level <- "subnational"
Expand Down

0 comments on commit d07a38f

Please sign in to comment.