differences for PR #187

Author: actions-user

config.yaml

@@ -196,18 +196,27 @@ estimates <- EpiNow2::epinow(
   # Add observation model
   obs = EpiNow2::obs_opts(scale = obs_scale)
 )
+```
 
+``` output
+WARN [2025-10-13 17:36:52] epinow: There were 1 transitions after warmup that exceeded the maximum treedepth. Increase max_treedepth above 12. See
+https://mc-stan.org/misc/warnings.html#maximum-treedepth-exceeded - 
+WARN [2025-10-13 17:36:52] epinow: Examine the pairs() plot to diagnose sampling problems
+ - 
+```
+
+``` r
 base::summary(estimates)
 ```
 
 ``` output
                         measure                 estimate
                          <char>                   <char>
-1:       New infections per day   20133 (13440 -- 30167)
+1:       New infections per day   19961 (13183 -- 30101)
 2:   Expected change in reports        Likely decreasing
 3:   Effective reproduction no.       0.96 (0.76 -- 1.2)
-4:               Rate of growth -0.013 (-0.095 -- 0.067)
-5: Doubling/halving time (days)         -55 (10 -- -7.3)
+4:               Rate of growth -0.013 (-0.096 -- 0.069)
+5: Doubling/halving time (days)         -55 (10 -- -7.2)
 ```
 
 
@@ -291,12 +300,6 @@ estimate_cases_to_deaths <- EpiNow2::estimate_secondary(
 )
 ```
 
-``` output
-WARN [2025-09-25 13:36:12] estimate_secondary (chain: 1): 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 - 
-```
-
 
 ::::::::::::::::::::::::::::::::::::: callout
 
@@ -527,21 +530,34 @@ ebola_estimates <- EpiNow2::epinow(
   # horizon needs to be 14 days to create two week forecast (default is 7 days)
   forecast = EpiNow2::forecast_opts(horizon = 14)
 )
+```
+
+``` output
+WARN [2025-10-13 17:37:56] epinow: There were 1 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 [2025-10-13 17:37:56] epinow: Examine the pairs() plot to diagnose sampling problems
+ - 
+WARN [2025-10-13 17:37:59] 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 - 
+```
 
+``` r
 summary(ebola_estimates)
 ```
 
 ``` output
                         measure                estimate
                          <char>                  <char>
-1:       New infections per day          91 (47 -- 190)
+1:       New infections per day          92 (50 -- 184)
 2:   Expected change in reports              Increasing
-3:   Effective reproduction no.        1.6 (1.1 -- 2.4)
-4:               Rate of growth 0.041 (0.0027 -- 0.088)
-5: Doubling/halving time (days)         17 (7.9 -- 260)
+3:   Effective reproduction no.        1.6 (1.2 -- 2.3)
+4:               Rate of growth 0.042 (0.0063 -- 0.083)
+5: Doubling/halving time (days)         17 (8.4 -- 110)
 ```
 
-The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.6 (1.1 -- 2.4). The exponential growth rate of case numbers is 0.041 (0.0027 -- 0.088).
+The effective reproduction number $R_t$ estimate (on the last date of the data) is 1.6 (1.2 -- 2.3). The exponential growth rate of case numbers is 0.042 (0.0063 -- 0.083).
 
 Visualize the estimates:
 

create-forecast.md

@@ -167,8 +167,8 @@ epiparameter::generate(covid_serialint, times = 10)
 ```
 
 ``` output
- [1]  3.270725  1.934659  3.429970  4.610215  1.973068  1.984880  6.728371
- [8]  2.991406 11.102558  9.452121
+ [1] 5.703076 4.331661 4.620432 2.113233 1.803581 4.448651 7.457367 3.663820
+ [9] 4.105711 2.722164
 ```
 
 ::::::::: instructor
@@ -762,27 +762,7 @@ epinow_estimates_cgi <- EpiNow2::epinow(
   generation_time = EpiNow2::generation_time_opts(covid_serial_interval),
   delays = EpiNow2::delay_opts(covid_incubation_time)
 )
-```
-
-``` output
-WARN [2025-10-10 18:10:57] epinow: There were 336 transitions after warmup that exceeded the maximum treedepth. Increase max_treedepth above 12. See
-https://mc-stan.org/misc/warnings.html#maximum-treedepth-exceeded - 
-WARN [2025-10-10 18:10:57] epinow: There were 1 chains where the estimated Bayesian Fraction of Missing Information was low. See
-https://mc-stan.org/misc/warnings.html#bfmi-low - 
-WARN [2025-10-10 18:10:57] epinow: Examine the pairs() plot to diagnose sampling problems
- - 
-WARN [2025-10-10 18:10:58] epinow: The largest R-hat is NA, indicating chains have not mixed.
-Running the chains for more iterations may help. See
-https://mc-stan.org/misc/warnings.html#r-hat - 
-WARN [2025-10-10 18:10:59] 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 [2025-10-10 18:11:00] 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 - 
-```
 
-``` r
 base::plot(epinow_estimates_cgi)
 ```
 
@@ -934,10 +914,10 @@ epinow_estimates_egi <- EpiNow2::epinow(
 ```
 
 ``` output
-WARN [2025-10-10 18:11:52] epinow: There were 3 divergent transitions after warmup. See
+WARN [2025-10-13 17:28:16] epinow: There were 4 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 [2025-10-10 18:11:52] epinow: Examine the pairs() plot to diagnose sampling problems
+WARN [2025-10-13 17:28:16] epinow: Examine the pairs() plot to diagnose sampling problems
  - 
 ```