adjusted scaling in xds_cohort: eir is the daily eir

dd-harp · Sep 26, 2024 · 33b8a3e · 33b8a3e
1 parent d267a06
commit 33b8a3e
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/R/xde_cohort.R b/R/xde_cohort.R
@@ -36,7 +36,7 @@ xds_solve_cohort = function(pars, bday=0, A = 10, da = 10){
 
   stats::integrate(F_eir, 0, 365, bday=0)$val -> scale
 
-  pars$EIRpar$scale = scale
+  pars$EIRpar$scale = scale*365
 
   F_eir <- with(pars$EIRpar, function(age, bday){
     eir/scale*F_season(age+bday)*F_trend(age+bday)*F_age(age)

diff --git a/R/xds_setup.R b/R/xds_setup.R
@@ -441,7 +441,7 @@ xds_setup_cohort = function(eir=1,
 
   stats::integrate(F_eir, 0, 365, bday=0)$val -> scale
 
-  pars$EIRpar$scale = scale
+  pars$EIRpar$scale = scale*365
 
   F_eir <- with(pars$EIRpar, function(age, bday){
     eir/scale*F_season(age+bday)*F_trend(age+bday)*F_age(age)