correct damping term

epiforecasts · Oct 11, 2023 · 544997e · 544997e
1 parent ef3e2ad
commit 544997e
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/README.Rmd b/README.Rmd
@@ -144,7 +144,7 @@ estimates <- epinow(
   delays = delay_opts(incubation_period + reporting_delay),
   rt = rt_opts(
     prior = list(mean = 2, sd = 0.2),
-    incidence_feedback = list(mean = 0, sd = 0.0001)
+    incidence_feedback = list(mean = 0.01, sd = 0.01)
   ),
   stan = stan_opts(cores = 4, control = list(adapt_delta = 0.99)),
   verbose = interactive()

diff --git a/inst/stan/functions/infections.stan b/inst/stan/functions/infections.stan
@@ -48,7 +48,7 @@ vector generate_infections(vector oR, int uot, vector gt_rev_pmf,
   for (s in 1:ot) {
     infectiousness[s] = update_infectiousness(infections, gt_rev_pmf, uot, s);
     if (if_present) {
-      infections[s + uot] = exp(
+      R[s] = exp(
         log(R[s]) - incidence_feedback[1] .* infectiousness[s]
       );
     }