Skip to content

Commit c7c4c60

Browse files
committed
substituting clearly bad data
1 parent 5f7892b commit c7c4c60

File tree

3 files changed

+18
-2
lines changed

3 files changed

+18
-2
lines changed

R/utils.R

Lines changed: 14 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -153,7 +153,20 @@ get_exclusions <- function(
153153
}
154154
return("")
155155
}
156-
156+
data_substitutions <- function(dataset, disease, forecast_generation_date) {
157+
disease <- "flu"
158+
forecast_generation_date <- as.Date("2025-01-08")
159+
substitutions <- readr::read_csv(
160+
glue::glue("{disease}_data_substitutions.csv"),
161+
comment = "#",
162+
show_col_types = FALSE) %>%
163+
filter(forecast_date == forecast_generation_date) %>%
164+
select(-forecast_date) %>%
165+
rename(new_value = value)
166+
dataset %>% left_join(substitutions) %>%
167+
mutate(value = ifelse(!is.na(new_value), new_value, value)) %>%
168+
select(-new_value)
169+
}
157170
parse_prod_weights <- function(filename = here::here("covid_geo_exclusions.csv"),
158171
gen_forecast_date) {
159172
all_states <- c(

flu_data_substitutions.csv

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,2 @@
1+
geo_value, forecast_date, time_value, value
2+
sd, 2025-01-08, 2025-01-01, 43

scripts/flu_hosp_prod.R

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -96,7 +96,8 @@ rlang::list2(
9696
time_value = time_value - 3
9797
) %>%
9898
filter(version == max(version)) %>%
99-
select(-version)
99+
select(-version) %>%
100+
data_substitutions(disease = "flu")
100101
},
101102
cue = tar_cue(mode = "always")
102103
),

0 commit comments

Comments
 (0)