We developed two R packages for this analysis: ‘rightTruncation’ [19] and ‘covidseir’ [34]. ‘rightTruncation’ performs maximum likelihood estimates of delay distributions accounting for right truncation. We used this approach to estimate the time between symptom onset and reporting in BC and New Zealand. The ‘covidseir’ package facilitates the SEIR model fitting of case-count data in Stan.