I had a first go at trying to reproduce the #2019nCoV #R0 results of Imai et al (Fig 1, Report 3 https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news--wuhan-coronavirus/ …), using code from @JulRiou @C_Althaus and didn't have much luck. Can anyone see where I'm going wrong?https://gist.github.com/sdwfrost/b3d0c4cbff7ff460562162affceffb17 …
-
No need to be embarrassed, but yes, there's lots of reallocations! Pre-allocating a fixed vector of max_cases times would be much faster (as I'm sure you know), but introduces more bookkeeping (and indexing issues cause sooo many bugs). I'll update the ref in my Google doc.
Yes, the bookkeeping and indexing was exactly where I didn’t want to go to.