Open metabolin opened 2 years ago
Hello, may I ask a question? In the script B. Calculating Rt using Maximum likelihood estimation. You give:
mutate(imported_delay = imported/delay_precalculation(t), domestic_delay = domestic/delay_precalculation(t)) %>% mutate(total=imported_delay+domestic_delay)
Here, if I understand correctly, you are using the equation (2) in [Tsuzuki et al 2017] (doi:10.2807/1560-7917.ES.2017.22.46.17-00710). But why you divide by time-delay distribution instead of multiplying it like what equation(2) does?
Hello, may I ask a question? In the script B. Calculating Rt using Maximum likelihood estimation. You give:
Here, if I understand correctly, you are using the equation (2) in [Tsuzuki et al 2017] (doi:10.2807/1560-7917.ES.2017.22.46.17-00710). But why you divide by time-delay distribution instead of multiplying it like what equation(2) does?