@@ -1605,7 +1605,7 @@ knitr::kable(nowcast_check)
1605
1605
#| eval: true
1606
1606
1607
1607
all_nowcast_dates = nchs_archive$DT |>
1608
- filter(time_value >= as.Date("2021 -01-10 ")) |>
1608
+ filter(time_value >= as.Date("2022 -01-01 ")) |>
1609
1609
distinct(time_value) |>
1610
1610
pull(time_value)
1611
1611
```
@@ -2036,7 +2036,6 @@ The first value here is a version of $Y_{\text{2020-11-30}}$ as it was reported
2036
2036
2037
2037
``` {r regression-trial-nowcast-predictors-data}
2038
2038
#| echo: true
2039
- #| eval: false
2040
2039
predictors <- predictor_descriptions |>
2041
2040
pmap(function(varname, lag_days, predictor_name) {
2042
2041
get_predictor_training_data(nchs_ca_past_archive, varname, lag_days, predictor_name)
@@ -2069,7 +2068,7 @@ reported at our trial nowcast date $t$.
2069
2068
2070
2069
``` {r regression-trial-nowcast-fit}
2071
2070
#| echo: true
2072
- #| eval: false
2071
+
2073
2072
training_test <- full_join(predictors, target, by = c("geo_value", "time_value"))
2074
2073
2075
2074
training <- training_test |> drop_na()
@@ -2092,7 +2091,7 @@ pred
2092
2091
2093
2092
``` {r regression-trial-nowcast-pred}
2094
2093
#| echo: true
2095
- #| eval: false
2094
+
2096
2095
pred
2097
2096
```
2098
2097
@@ -2296,7 +2295,7 @@ reg4_settings <- c(reg2_settings, median = TRUE)
2296
2295
```
2297
2296
2298
2297
``` {r regression-run-nowcasts-backtesting}
2299
- #| eval: false
2298
+
2300
2299
hhs_hosp_archive <- pub_covidcast(
2301
2300
source = "hhs",
2302
2301
signals = "confirmed_admissions_covid_1d_7dav",
@@ -2346,7 +2345,6 @@ reg4_nowcasts <- hosp_mort_archive |>
2346
2345
## Data wrangling
2347
2346
2348
2347
``` {r regression-nowcast-wrangling}
2349
- #| eval: false
2350
2348
ratio_nowcasts_archive <- nowcasts |>
2351
2349
filter(geo_value == "ca") |>
2352
2350
rename(time_value = target_date,
@@ -2383,7 +2381,6 @@ nowcast_comparison <-
2383
2381
2384
2382
``` {r regression-nowcast-plot-linreg}
2385
2383
#| fig-width: 7
2386
- #| eval: false
2387
2384
nowcast_comparison |>
2388
2385
filter(target_date >= min(all_nowcast_dates) - 35,
2389
2386
!(Nowcaster %in% c("QuantReg model", "QuantReg + hosp"))) |>
@@ -2399,7 +2396,6 @@ nowcast_comparison |>
2399
2396
2400
2397
``` {r regression-nowcast-plot-quantreg}
2401
2398
#| fig-width: 7
2402
- #| eval: false
2403
2399
2404
2400
nowcast_comparison |>
2405
2401
filter(target_date >= min(all_nowcast_dates) - 35,
@@ -2416,7 +2412,6 @@ nowcast_comparison |>
2416
2412
2417
2413
``` {r regression-nowcast-eval-comparison}
2418
2414
#| echo: true
2419
- #| eval: false
2420
2415
2421
2416
n_models <- length(unique(nowcast_comparison$Nowcaster))
2422
2417
nowcast_comparison |>
0 commit comments