Patients were followed up for a maximum of 14 days from their earliest COVID-19-positive specimen until the hospital admission or emergency care attendance date. Patients were censored at the date of death if this occurred without a previous hospital attendance event within the 14-day period.
In the primary analysis, stratified Cox regression was used to estimate hazard ratios (HRs) of the hospitalisation outcomes (hospital admission or emergency care attendance) for patients with the delta variant compared with patients with the alpha variant. Strata were created by intersecting the likely confounders. Additional potential confounders were included using main effects. Linear main effects terms for age and calendar date were used to adjust for residual confounding after stratification.
In the secondary analysis, the HRs of the hospitalisation outcomes by variant were estimated by vaccination status. The base models were refitted with an interaction term between variant and vaccination. Due to low numbers of patients with COVID-19 who had been vaccinated, and consequently low numbers within some vaccination categories, vaccination status was grouped into two categories: unvaccinated or less than 21 days since the first vaccination dose; and 21 days or more since the first vaccination dose, with or without the second dose.
In additional analyses, the proportional hazards assumption of the Cox regression model was graphically assessed using Schoenfeld residual plots and formally tested using the Schoenfeld test. Post-evaluations of the relative magnitudes of the confounders' contribution to the adjusted HRs were done by sequentially adding the adjustment variables in the order of the percentage change in the adjusted HRs for patients with the delta variant versus the alpha variant. To assess the impact of stratification versus regression modelling on the HRs and 95% CIs, the primary model was refitted with each stratification variable instead included as a regression variable.
HRs were assessed for sensitivity to stratification by alternative region or calendar period covariates, confounding due to recent international travel or symptomatic status subgroups, or to the precise outcome definitions. Details are shown in the appendix (p 8).
Data were prepared using Stata version 15.1. Statistical analyses were done in R version 4.1.0.
Free full text: Click here