tests/testthat/_snaps/enw_nowcast_summary.md

enw_nowcast_summary() can extract the summarised nowcast as expected

Code
  round_numerics(nowcast[, c("rhat", "ess_bulk", "ess_tail") := NULL])
Output
      reference_date report_date .group max_confirm location age_group confirm
   1:     2021-08-03  2021-08-22      1         149       DE       00+     149
   2:     2021-08-04  2021-08-22      1         166       DE       00+     166
   3:     2021-08-05  2021-08-22      1         133       DE       00+     133
   4:     2021-08-06  2021-08-22      1         137       DE       00+     137
   5:     2021-08-07  2021-08-22      1         139       DE       00+     139
   6:     2021-08-08  2021-08-22      1          97       DE       00+      97
   7:     2021-08-09  2021-08-22      1          58       DE       00+      58
   8:     2021-08-10  2021-08-22      1         175       DE       00+     175
   9:     2021-08-11  2021-08-22      1         233       DE       00+     233
  10:     2021-08-12  2021-08-22      1         237       DE       00+     237
  11:     2021-08-13  2021-08-22      1         204       DE       00+     204
  12:     2021-08-14  2021-08-22      1         189       DE       00+     189
  13:     2021-08-15  2021-08-22      1         125       DE       00+     125
  14:     2021-08-16  2021-08-22      1          98       DE       00+      98
  15:     2021-08-17  2021-08-22      1         242       DE       00+     242
  16:     2021-08-18  2021-08-22      1         223       DE       00+     223
  17:     2021-08-19  2021-08-22      1         202       DE       00+     202
  18:     2021-08-20  2021-08-22      1         171       DE       00+     171
  19:     2021-08-21  2021-08-22      1         112       DE       00+     112
  20:     2021-08-22  2021-08-22      1          45       DE       00+      45
      cum_prop_reported delay prop_reported mean median  sd mad  q5 q20 q35 q50
   1:                 1    19             0  149    149   0   0 149 149 149 149
   2:                 1    18             0  167    167   1   1 166 166 167 167
   3:                 1    17             0  136    136   2   1 133 134 135 136
   4:                 1    16             0  141    141   2   2 138 139 140 141
   5:                 1    15             0  146    146   3   3 141 143 144 146
   6:                 1    14             0  104    103   3   3  99 101 102 103
   7:                 1    13             0   63     62   2   1  59  61  62  62
   8:                 1    12             0  185    185   4   4 180 182 183 185
   9:                 1    11             0  256    255   6   6 246 250 253 255
  10:                 1    10             0  267    267   8   7 255 261 264 267
  11:                 1     9             0  236    236   8   9 223 229 232 236
  12:                 1     8             0  231    230  10  10 216 222 226 230
  13:                 1     7             0  165    165  10  10 150 157 160 165
  14:                 1     6             0  129    128   9   9 116 122 126 128
  15:                 1     5             0  293    292  12  12 275 283 288 292
  16:                 1     4             0  293    292  16  15 270 280 286 292
  17:                 1     3             0  291    290  19  19 263 274 284 290
  18:                 1     2             0  296    292  29  28 255 271 282 292
  19:                 1     1             0  309    304  51  47 237 267 285 304
  20:                 1     0             1  381    371 104  96 241 294 332 371
      q65 q80 q95
   1: 149 149 149
   2: 168 168 170
   3: 136 137 139
   4: 142 143 145
   5: 147 148 151
   6: 105 106 109
   7:  63  65  67
   8: 186 188 192
   9: 258 261 267
  10: 270 274 281
  11: 239 243 252
  12: 234 239 249
  13: 168 173 183
  14: 132 136 144
  15: 297 303 313
  16: 298 306 320
  17: 298 306 325
  18: 304 319 348
  19: 323 349 406
  20: 405 459 569

enw_nowcast_summary() can extract summarised nowcast with custom quantiles

Code
  round_numerics(nowcast[, c("rhat", "ess_bulk", "ess_tail") := NULL])
Output
      reference_date report_date .group max_confirm location age_group confirm
   1:     2021-08-03  2021-08-22      1         149       DE       00+     149
   2:     2021-08-04  2021-08-22      1         166       DE       00+     166
   3:     2021-08-05  2021-08-22      1         133       DE       00+     133
   4:     2021-08-06  2021-08-22      1         137       DE       00+     137
   5:     2021-08-07  2021-08-22      1         139       DE       00+     139
   6:     2021-08-08  2021-08-22      1          97       DE       00+      97
   7:     2021-08-09  2021-08-22      1          58       DE       00+      58
   8:     2021-08-10  2021-08-22      1         175       DE       00+     175
   9:     2021-08-11  2021-08-22      1         233       DE       00+     233
  10:     2021-08-12  2021-08-22      1         237       DE       00+     237
  11:     2021-08-13  2021-08-22      1         204       DE       00+     204
  12:     2021-08-14  2021-08-22      1         189       DE       00+     189
  13:     2021-08-15  2021-08-22      1         125       DE       00+     125
  14:     2021-08-16  2021-08-22      1          98       DE       00+      98
  15:     2021-08-17  2021-08-22      1         242       DE       00+     242
  16:     2021-08-18  2021-08-22      1         223       DE       00+     223
  17:     2021-08-19  2021-08-22      1         202       DE       00+     202
  18:     2021-08-20  2021-08-22      1         171       DE       00+     171
  19:     2021-08-21  2021-08-22      1         112       DE       00+     112
  20:     2021-08-22  2021-08-22      1          45       DE       00+      45
      cum_prop_reported delay prop_reported mean median  sd mad  q5 q50 q95
   1:                 1    19             0  149    149   0   0 149 149 149
   2:                 1    18             0  167    167   1   1 166 167 170
   3:                 1    17             0  136    136   2   1 133 136 139
   4:                 1    16             0  141    141   2   2 138 141 145
   5:                 1    15             0  146    146   3   3 141 146 151
   6:                 1    14             0  104    103   3   3  99 103 109
   7:                 1    13             0   63     62   2   1  59  62  67
   8:                 1    12             0  185    185   4   4 180 185 192
   9:                 1    11             0  256    255   6   6 246 255 267
  10:                 1    10             0  267    267   8   7 255 267 281
  11:                 1     9             0  236    236   8   9 223 236 252
  12:                 1     8             0  231    230  10  10 216 230 249
  13:                 1     7             0  165    165  10  10 150 165 183
  14:                 1     6             0  129    128   9   9 116 128 144
  15:                 1     5             0  293    292  12  12 275 292 313
  16:                 1     4             0  293    292  16  15 270 292 320
  17:                 1     3             0  291    290  19  19 263 290 325
  18:                 1     2             0  296    292  29  28 255 292 348
  19:                 1     1             0  309    304  51  47 237 304 406
  20:                 1     0             1  381    371 104  96 241 371 569

enw_nowcast_summary() can extract the summarised nowcast as expected when a delay smaller than specified was modelled

Code
  round_numerics(nowcast[, c("rhat", "ess_bulk", "ess_tail") := NULL])
Output
      reference_date report_date .group max_confirm location age_group confirm
   1:     2021-08-01  2021-08-20      1          77       DE       00+      77
   2:     2021-08-02  2021-08-21      1          59       DE       00+      59
   3:     2021-08-03  2021-08-22      1         149       DE       00+     149
   4:     2021-08-04  2021-08-22      1         166       DE       00+     166
   5:     2021-08-05  2021-08-22      1         133       DE       00+     133
   6:     2021-08-06  2021-08-22      1         137       DE       00+     137
   7:     2021-08-07  2021-08-22      1         139       DE       00+     139
   8:     2021-08-08  2021-08-22      1          97       DE       00+      97
   9:     2021-08-09  2021-08-22      1          58       DE       00+      58
  10:     2021-08-10  2021-08-22      1         175       DE       00+     175
  11:     2021-08-11  2021-08-22      1         233       DE       00+     233
  12:     2021-08-12  2021-08-22      1         237       DE       00+     237
  13:     2021-08-13  2021-08-22      1         204       DE       00+     204
  14:     2021-08-14  2021-08-22      1         189       DE       00+     189
  15:     2021-08-15  2021-08-22      1         125       DE       00+     125
  16:     2021-08-16  2021-08-22      1          98       DE       00+      98
  17:     2021-08-17  2021-08-22      1         242       DE       00+     242
  18:     2021-08-18  2021-08-22      1         223       DE       00+     223
  19:     2021-08-19  2021-08-22      1         202       DE       00+     202
  20:     2021-08-20  2021-08-22      1         171       DE       00+     171
  21:     2021-08-21  2021-08-22      1         112       DE       00+     112
  22:     2021-08-22  2021-08-22      1          45       DE       00+      45
      reference_date report_date .group max_confirm location age_group confirm
      cum_prop_reported delay prop_reported mean median  sd mad  q5 q20 q35 q50
   1:                 1    19             0   77     77   0   0  77  77  77  77
   2:                 1    19             0   59     59   0   0  59  59  59  59
   3:                 1    19             0  149    149   0   0 149 149 149 149
   4:                 1    18             0  167    167   1   1 166 166 167 167
   5:                 1    17             0  136    136   2   1 133 134 135 136
   6:                 1    16             0  141    141   2   2 138 139 140 141
   7:                 1    15             0  146    146   3   3 141 143 144 146
   8:                 1    14             0  104    103   3   3  99 101 102 103
   9:                 1    13             0   63     62   2   1  59  61  62  62
  10:                 1    12             0  185    185   4   4 180 182 183 185
  11:                 1    11             0  256    255   6   6 246 250 253 255
  12:                 1    10             0  267    267   8   7 255 261 264 267
  13:                 1     9             0  236    236   8   9 223 229 232 236
  14:                 1     8             0  231    230  10  10 216 222 226 230
  15:                 1     7             0  165    165  10  10 150 157 160 165
  16:                 1     6             0  129    128   9   9 116 122 126 128
  17:                 1     5             0  293    292  12  12 275 283 288 292
  18:                 1     4             0  293    292  16  15 270 280 286 292
  19:                 1     3             0  291    290  19  19 263 274 284 290
  20:                 1     2             0  296    292  29  28 255 271 282 292
  21:                 1     1             0  309    304  51  47 237 267 285 304
  22:                 1     0             1  381    371 104  96 241 294 332 371
      cum_prop_reported delay prop_reported mean median  sd mad  q5 q20 q35 q50
      q65 q80 q95
   1:  77  77  77
   2:  59  59  59
   3: 149 149 149
   4: 168 168 170
   5: 136 137 139
   6: 142 143 145
   7: 147 148 151
   8: 105 106 109
   9:  63  65  67
  10: 186 188 192
  11: 258 261 267
  12: 270 274 281
  13: 239 243 252
  14: 234 239 249
  15: 168 173 183
  16: 132 136 144
  17: 297 303 313
  18: 298 306 320
  19: 298 306 325
  20: 304 319 348
  21: 323 349 406
  22: 405 459 569
      q65 q80 q95


epiforecasts/epinowcast documentation built on Feb. 3, 2025, 4:17 p.m.