@@ -270,24 +270,25 @@ current and past COVID-19 case rates for all states. That is, we can make
270
270
forecasts on the archive, ` x ` , and compare those to forecasts on the latest
271
271
data, ` x_latest ` using the same general set-up as above. For version-aware
272
272
forecasting, note that ` x ` is fed into ` epix_slide() ` , while for version-unaware
273
- forecasting, ` x_latest ` is fed into ` epi_slide() ` . #%% update to include
274
- percent_cli after that issue is fixed?
273
+ forecasting, ` x_latest ` is fed into ` epi_slide() ` .
275
274
276
275
``` {r make-ar-kweek-asof}
277
276
k_week_ahead_as_of <- function(ahead = 7, version_aware = TRUE) {
278
277
if (version_aware) {
279
278
epix_slide(
280
279
x,
281
- ~ arx_forecaster(.x, "case_rate", "case_rate",
280
+ ~ arx_forecaster(.x, "case_rate", c( "case_rate", "percent_cli") ,
282
281
args_list = arx_args_list(ahead = ahead)) %>%
283
282
extract2("predictions") %>%
284
- select(-c(geo_value, time_value) ),
283
+ select(-time_value),
285
284
before = 120 - 1,
286
285
ref_time_values = fc_time_values,
287
286
new_col_name = "fc") %>%
288
- mutate(engine_type = "lm", version_aware = version_aware)
287
+ mutate(engine_type = "lm", version_aware = version_aware) %>%
288
+ rename(geo_value = fc_geo_value)
289
289
} else {
290
- k_week_ahead(x_latest, "case_rate", "case_rate", ahead, linear_reg()) %>%
290
+ k_week_ahead(x_latest, "case_rate", c("case_rate", "percent_cli"),
291
+ ahead, linear_reg()) %>%
291
292
mutate(version_aware = version_aware)
292
293
}
293
294
}
0 commit comments