Skip to content

Commit daf69d3

Browse files
committed
differences for PR #187
1 parent 551bfe1 commit daf69d3

25 files changed

+114
-6352
lines changed

config.yaml

Lines changed: 0 additions & 86 deletions
This file was deleted.

create-forecast.md

Lines changed: 41 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -140,6 +140,14 @@ estimates <- EpiNow2::epinow(
140140
)
141141
```
142142

143+
``` output
144+
WARN [2025-09-16 12:37:04] epinow: There were 2 divergent transitions after warmup. See
145+
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
146+
to find out why this is a problem and how to eliminate them. -
147+
WARN [2025-09-16 12:37:04] epinow: Examine the pairs() plot to diagnose sampling problems
148+
-
149+
```
150+
143151
::::::::::::::::::::::::::::::::: callout
144152

145153
### Do not wait for this to complete!
@@ -196,18 +204,36 @@ estimates <- EpiNow2::epinow(
196204
# Add observation model
197205
obs = EpiNow2::obs_opts(scale = obs_scale)
198206
)
207+
```
208+
209+
``` output
210+
WARN [2025-09-16 12:42:08] epinow: There were 1 divergent transitions after warmup. See
211+
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
212+
to find out why this is a problem and how to eliminate them. -
213+
WARN [2025-09-16 12:42:08] epinow: There were 500 transitions after warmup that exceeded the maximum treedepth. Increase max_treedepth above 12. See
214+
https://mc-stan.org/misc/warnings.html#maximum-treedepth-exceeded -
215+
WARN [2025-09-16 12:42:08] epinow: Examine the pairs() plot to diagnose sampling problems
216+
-
217+
WARN [2025-09-16 12:42:09] epinow: Bulk Effective Samples Size (ESS) is too low, indicating posterior means and medians may be unreliable.
218+
Running the chains for more iterations may help. See
219+
https://mc-stan.org/misc/warnings.html#bulk-ess -
220+
WARN [2025-09-16 12:42:10] epinow: Tail Effective Samples Size (ESS) is too low, indicating posterior variances and tail quantiles may be unreliable.
221+
Running the chains for more iterations may help. See
222+
https://mc-stan.org/misc/warnings.html#tail-ess -
223+
```
199224

225+
``` r
200226
base::summary(estimates)
201227
```
202228

203229
``` output
204-
measure estimate
205-
<char> <char>
206-
1: New infections per day 20091 (13574 -- 30897)
207-
2: Expected change in reports Likely decreasing
208-
3: Effective reproduction no. 0.97 (0.77 -- 1.2)
209-
4: Rate of growth -0.012 (-0.09 -- 0.071)
210-
5: Doubling/halving time (days) -58 (9.7 -- -7.7)
230+
measure estimate
231+
<char> <char>
232+
1: New infections per day 20148 (13173 -- 30177)
233+
2: Expected change in reports Stable
234+
3: Effective reproduction no. 0.97 (0.77 -- 1.2)
235+
4: Rate of growth -0.0099 (-0.09 -- 0.067)
236+
5: Doubling/halving time (days) -70 (10 -- -7.7)
211237
```
212238

213239

@@ -512,34 +538,21 @@ ebola_estimates <- EpiNow2::epinow(
512538
# horizon needs to be 14 days to create two week forecast (default is 7 days)
513539
forecast = EpiNow2::forecast_opts(horizon = 14)
514540
)
515-
```
516541

517-
``` output
518-
WARN [2025-08-26 01:46:34] epinow: There were 6 divergent transitions after warmup. See
519-
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
520-
to find out why this is a problem and how to eliminate them. -
521-
WARN [2025-08-26 01:46:34] epinow: Examine the pairs() plot to diagnose sampling problems
522-
-
523-
WARN [2025-08-26 01:46:38] epinow: Tail Effective Samples Size (ESS) is too low, indicating posterior variances and tail quantiles may be unreliable.
524-
Running the chains for more iterations may help. See
525-
https://mc-stan.org/misc/warnings.html#tail-ess -
526-
```
527-
528-
``` r
529542
summary(ebola_estimates)
530543
```
531544

532545
``` output
533-
measure estimate
534-
<char> <char>
535-
1: New infections per day 92 (49 -- 185)
536-
2: Expected change in reports Increasing
537-
3: Effective reproduction no. 1.6 (1.2 -- 2.4)
538-
4: Rate of growth 0.04 (0.0038 -- 0.083)
539-
5: Doubling/halving time (days) 17 (8.3 -- 180)
546+
measure estimate
547+
<char> <char>
548+
1: New infections per day 91 (48 -- 189)
549+
2: Expected change in reports Increasing
550+
3: Effective reproduction no. 1.6 (1.2 -- 2.4)
551+
4: Rate of growth 0.041 (0.0036 -- 0.087)
552+
5: Doubling/halving time (days) 17 (8 -- 190)
540553
```
541554

542-
The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.6 (1.2 -- 2.4). The exponential growth rate of case numbers is 0.04 (0.0038 -- 0.083).
555+
The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.6 (1.2 -- 2.4). The exponential growth rate of case numbers is 0.041 (0.0036 -- 0.087).
543556

544557
Visualize the estimates:
545558

delays-functions.md

Lines changed: 4 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -161,8 +161,8 @@ generate(covid_serialint, times = 10)
161161
```
162162

163163
``` output
164-
[1] 3.133504 6.369584 2.651628 10.017361 1.305227 2.536725 2.996175
165-
[8] 7.880128 6.562764 5.961404
164+
[1] 3.492157 2.822130 5.519427 3.704806 2.092871 1.705613 4.284163 3.869586
165+
[9] 3.035544 2.764683
166166
```
167167

168168
::::::::: instructor
@@ -588,12 +588,10 @@ epinow_estimates_cgi <- epinow(
588588
```
589589

590590
``` output
591-
WARN [2025-08-26 01:35:59] epinow: There were 1 divergent transitions after warmup. See
591+
WARN [2025-09-16 12:31:55] epinow: There were 1 divergent transitions after warmup. See
592592
https://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup
593593
to find out why this is a problem and how to eliminate them. -
594-
WARN [2025-08-26 01:35:59] epinow: There were 375 transitions after warmup that exceeded the maximum treedepth. Increase max_treedepth above 12. See
595-
https://mc-stan.org/misc/warnings.html#maximum-treedepth-exceeded -
596-
WARN [2025-08-26 01:35:59] epinow: Examine the pairs() plot to diagnose sampling problems
594+
WARN [2025-09-16 12:31:55] epinow: Examine the pairs() plot to diagnose sampling problems
597595
-
598596
```
599597

-5 Bytes
Loading
-47 Bytes
Loading
537 Bytes
Loading
-23 Bytes
Loading
476 Bytes
Loading
879 Bytes
Loading
-496 Bytes
Loading

0 commit comments

Comments
 (0)