In each province, the incidence rate for each AESI will be calculated as the number of cases per 100,000 person-years at risk during the study entire study period, excluding the COVID-19 pandemic period, in order to avoid its influence on the baseline rates of these outcomes, and 95% confidence intervals will be estimated using the Poisson distribution. Using data from the same period as the overall rates, incidence rates will be stratified by the following subgroups: calendar year, calendar quarter, time period relative the COVID-19 pandemic (pre-: 2015-2019; during: 2020-2022; and post-: 2023 to latest available data), sex, age group (pediatric cohort: 0-12 and 13-17; adult cohort: 18-29, 30-39, and so on to ≥ 80 years), and sex and age group. All incidence rates will be further age and sex standardized to the Canadian and European populations by the direct method.
For each incident AESI cohort, demographic and clinical characteristics of cases will be summarized. To contextualize these characteristics, each case will be matched 1:1 to individuals without the AESIs on age, sex, and calendar time of the event. Characteristics between the AESI and matched cohorts will be compared using standardized mean differences.
Three sensitivity analyses will be conducted. First, as in the DARWIN EU study, different case definitions for selected AESIs (encephalitis and ischemic stroke) will be applied to assess the impact of varying sensitivity and specificity. Second, rates will be estimated in Manitoba and Ontario using data only from the post-pandemic period, as this may better represent the current baseline rates of these AESIs. Third, a 90-day washout window for AESIs will be used in Manitoba and Ontario for closer alignment with the DARWIN EU study. If necessary, these sensitivity analyses will be conducted in the remaining provinces.
Results will be reported for each province separately and pooled using random-effects meta-analyses to provide national baseline rates.