Skip to content

Commit

Permalink
differences for PR #138
Browse files Browse the repository at this point in the history
  • Loading branch information
actions-user committed Mar 5, 2024
1 parent 53e987e commit 140d168
Show file tree
Hide file tree
Showing 17 changed files with 109 additions and 2,483 deletions.
87 changes: 0 additions & 87 deletions config.yaml

This file was deleted.

53 changes: 30 additions & 23 deletions create-forecast.md
Original file line number Diff line number Diff line change
Expand Up @@ -61,10 +61,10 @@ estimates <- epinow(
```

```{.output}
WARN [2024-02-27 02:04:29] epinow: There were 3 divergent transitions after warmup. See
WARN [2024-03-05 01:48:51] epinow: There were 3 divergent transitions after warmup. See
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
to find out why this is a problem and how to eliminate them. -
WARN [2024-02-27 02:04:29] epinow: Examine the pairs() plot to diagnose sampling problems
WARN [2024-03-05 01:48:51] epinow: Examine the pairs() plot to diagnose sampling problems
-
```

Expand Down Expand Up @@ -131,17 +131,28 @@ estimates <- epinow(
rt = rt_opts(prior = list(mean = rt_log_mean, sd = rt_log_sd)),
obs = obs_opts(scale = obs_scale)
)
```

```{.output}
WARN [2024-03-05 01:55:24] epinow: There were 5 divergent transitions after warmup. See
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
to find out why this is a problem and how to eliminate them. -
WARN [2024-03-05 01:55:24] epinow: Examine the pairs() plot to diagnose sampling problems
-
```

```r
summary(estimates)
```

```{.output}
measure estimate
<char> <char>
1: New confirmed cases by infection date 18077 (9954 -- 32284)
2: Expected change in daily cases Likely decreasing
3: Effective reproduction no. 0.89 (0.56 -- 1.3)
4: Rate of growth -0.015 (-0.065 -- 0.041)
5: Doubling/halving time (days) -46 (17 -- -11)
measure estimate
<char> <char>
1: New confirmed cases by infection date 17939 (9853 -- 32383)
2: Expected change in daily cases Likely decreasing
3: Effective reproduction no. 0.89 (0.56 -- 1.3)
4: Rate of growth -0.015 (-0.066 -- 0.04)
5: Doubling/halving time (days) -46 (17 -- -11)
```


Expand Down Expand Up @@ -235,10 +246,6 @@ deaths_forecast <- forecast_secondary(
plot(deaths_forecast)
```

```{.warning}
Warning: Removed 30 rows containing missing values (`position_stack()`).
```

<img src="fig/create-forecast-rendered-unnamed-chunk-11-1.png" style="display: block; margin: auto;" />

The plot shows the forecast secondary observations (deaths) over the dates which we have recorded cases for.
Expand Down Expand Up @@ -328,10 +335,10 @@ ebola_estimates <- epinow(
```

```{.output}
WARN [2024-02-27 02:13:59] epinow: There were 18 divergent transitions after warmup. See
WARN [2024-03-05 01:57:16] epinow: There were 29 divergent transitions after warmup. See
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
to find out why this is a problem and how to eliminate them. -
WARN [2024-02-27 02:13:59] epinow: Examine the pairs() plot to diagnose sampling problems
WARN [2024-03-05 01:57:16] epinow: Examine the pairs() plot to diagnose sampling problems
-
```

Expand All @@ -340,16 +347,16 @@ summary(ebola_estimates)
```

```{.output}
measure estimate
<char> <char>
1: New confirmed cases by infection date 101 (49 -- 270)
2: Expected change in daily cases Increasing
3: Effective reproduction no. 1.7 (1.1 -- 3)
4: Rate of growth 0.043 (0.0058 -- 0.094)
5: Doubling/halving time (days) 16 (7.3 -- 120)
measure estimate
<char> <char>
1: New confirmed cases by infection date 101 (48 -- 261)
2: Expected change in daily cases Increasing
3: Effective reproduction no. 1.7 (1.1 -- 2.9)
4: Rate of growth 0.042 (0.0048 -- 0.09)
5: Doubling/halving time (days) 17 (7.7 -- 140)
```

The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.7 (1.1 -- 3). The exponential growth rate of case numbers is 0.043 (0.0058 -- 0.094).
The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.7 (1.1 -- 2.9). The exponential growth rate of case numbers is 0.042 (0.0048 -- 0.09).

Visualize the estimates:

Expand Down
Binary file modified fig/create-forecast-rendered-unnamed-chunk-11-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/create-forecast-rendered-unnamed-chunk-3-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/create-forecast-rendered-unnamed-chunk-9-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/introduction-rendered-unnamed-chunk-5-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/model-choices-rendered-unnamed-chunk-3-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/model-choices-rendered-unnamed-chunk-4-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/quantify-transmissibility-rendered-unnamed-chunk-15-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/quantify-transmissibility-rendered-unnamed-chunk-16-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/quantify-transmissibility-rendered-unnamed-chunk-19-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified fig/simulating-transmission-rendered-plot-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
5 changes: 1 addition & 4 deletions introduction.md
Original file line number Diff line number Diff line change
Expand Up @@ -135,12 +135,9 @@ epinow_estimates <- epinow(
```

```{.output}
WARN [2024-02-27 01:05:38] epinow: Bulk Effective Samples Size (ESS) is too low, indicating posterior means and medians may be unreliable.
WARN [2024-03-05 00:49:40] epinow: Bulk Effective Samples Size (ESS) is too low, indicating posterior means and medians may be unreliable.
Running the chains for more iterations may help. See
https://mc-stan.org/misc/warnings.html#bulk-ess -
WARN [2024-02-27 01:05:39] epinow: Tail Effective Samples Size (ESS) is too low, indicating posterior variances and tail quantiles may be unreliable.
Running the chains for more iterations may help. See
https://mc-stan.org/misc/warnings.html#tail-ess -
```

As an output, we get the time-varying (or [effective](../learners/reference.md#effectiverepro)) reproduction number, as well as the cases by date of report and date of infection:
Expand Down
36 changes: 18 additions & 18 deletions md5sum.txt
Original file line number Diff line number Diff line change
@@ -1,19 +1,19 @@
"file" "checksum" "built" "date"
"CODE_OF_CONDUCT.md" "549f00b0992a7743c2bc16ea6ce3db57" "site/built/CODE_OF_CONDUCT.md" "2024-02-27"
"LICENSE.md" "14377518ee654005a18cf28549eb30e3" "site/built/LICENSE.md" "2024-02-27"
"config.yaml" "d11647b34aa74192704ce8944e83b0d8" "site/built/config.yaml" "2024-02-27"
"index.md" "adfca1a79e0106ee8b3f7731d0678b59" "site/built/index.md" "2024-02-27"
"links.md" "8184cf4149eafbf03ce8da8ff0778c14" "site/built/links.md" "2024-02-27"
"episodes/introduction.Rmd" "1861ec911e99b7feab3b65e0dbd91d41" "site/built/introduction.md" "2024-02-27"
"episodes/read-delays.Rmd" "2c7422e552f8e16c28449877182c0704" "site/built/read-delays.md" "2024-02-27"
"episodes/quantify-transmissibility.Rmd" "f507b32d4af315d95acc3f3f5a60c202" "site/built/quantify-transmissibility.md" "2024-02-27"
"episodes/create-forecast.Rmd" "3fdb0f4223f220d061431f18846bf1f4" "site/built/create-forecast.md" "2024-02-27"
"episodes/simulating-transmission.Rmd" "3f006aa142099aaeaa4688c064bd5d4a" "site/built/simulating-transmission.md" "2024-02-27"
"episodes/model-choices.Rmd" "775bfba024103820482bc56907d86f14" "site/built/model-choices.md" "2024-02-27"
"episodes/modelling-interventions.Rmd" "4955dd207b237ad495b601315bb3e8df" "site/built/modelling-interventions.md" "2024-02-27"
"episodes/compare-interventions.Rmd" "1f04f6fa8980d60273755c7faf20e70c" "site/built/compare-interventions.md" "2024-02-27"
"instructors/instructor-notes.md" "ca3834a1b0f9e70c4702aa7a367a6bb5" "site/built/instructor-notes.md" "2024-02-27"
"learners/reference.md" "3c3b0a87a98f91231d47f5a2e24ad27d" "site/built/reference.md" "2024-02-27"
"learners/setup.md" "9030a5d0d8e636cb89ee8d48e864469e" "site/built/setup.md" "2024-02-27"
"profiles/learner-profiles.md" "31b503c4b5bd1f0960ada730eca4a25e" "site/built/learner-profiles.md" "2024-02-27"
"renv/profiles/lesson-requirements/renv.lock" "3c56dface88f3db2ac611fca64c6e7c6" "site/built/renv.lock" "2024-02-27"
"CODE_OF_CONDUCT.md" "549f00b0992a7743c2bc16ea6ce3db57" "site/built/CODE_OF_CONDUCT.md" "2024-03-05"
"LICENSE.md" "14377518ee654005a18cf28549eb30e3" "site/built/LICENSE.md" "2024-03-05"
"config.yaml" "d11647b34aa74192704ce8944e83b0d8" "site/built/config.yaml" "2024-03-05"
"index.md" "adfca1a79e0106ee8b3f7731d0678b59" "site/built/index.md" "2024-03-05"
"links.md" "8184cf4149eafbf03ce8da8ff0778c14" "site/built/links.md" "2024-03-05"
"episodes/introduction.Rmd" "1861ec911e99b7feab3b65e0dbd91d41" "site/built/introduction.md" "2024-03-05"
"episodes/read-delays.Rmd" "2c7422e552f8e16c28449877182c0704" "site/built/read-delays.md" "2024-03-05"
"episodes/quantify-transmissibility.Rmd" "f507b32d4af315d95acc3f3f5a60c202" "site/built/quantify-transmissibility.md" "2024-03-05"
"episodes/create-forecast.Rmd" "3fdb0f4223f220d061431f18846bf1f4" "site/built/create-forecast.md" "2024-03-05"
"episodes/simulating-transmission.Rmd" "3f006aa142099aaeaa4688c064bd5d4a" "site/built/simulating-transmission.md" "2024-03-05"
"episodes/model-choices.Rmd" "775bfba024103820482bc56907d86f14" "site/built/model-choices.md" "2024-03-05"
"episodes/modelling-interventions.Rmd" "4955dd207b237ad495b601315bb3e8df" "site/built/modelling-interventions.md" "2024-03-05"
"episodes/compare-interventions.Rmd" "1f04f6fa8980d60273755c7faf20e70c" "site/built/compare-interventions.md" "2024-03-05"
"instructors/instructor-notes.md" "ca3834a1b0f9e70c4702aa7a367a6bb5" "site/built/instructor-notes.md" "2024-03-05"
"learners/reference.md" "3c3b0a87a98f91231d47f5a2e24ad27d" "site/built/reference.md" "2024-03-05"
"learners/setup.md" "9030a5d0d8e636cb89ee8d48e864469e" "site/built/setup.md" "2024-03-05"
"profiles/learner-profiles.md" "31b503c4b5bd1f0960ada730eca4a25e" "site/built/learner-profiles.md" "2024-03-05"
"renv/profiles/lesson-requirements/renv.lock" "91eb6875f1fb9811cc55bce90ae063db" "site/built/renv.lock" "2024-03-05"
Loading

0 comments on commit 140d168

Please sign in to comment.