*2.4. Data Management and Analysis*

The data were exported as a Microsoft Excel 2013 worksheet into R Statistical Software (v4.1.2; R Core Team 2021) for statistical analysis. District vaccination coverage was weighted using completeness of the data entry according to regional data completeness in 2020. Data on regional and district data entry completeness reported in the DHIS-2 database (in percentages) was downloaded and cleaned, and the formulas below were used to calculate the adjusted vaccine coverages.

$$Adjusted\,DTP-1\,coverage = \frac{\text{Regional data completeness}}{\text{District data completes}} \times \,DTP-1\,coverage \tag{1}$$

$$Adjusted\,DTP-3\,coverage = \frac{\text{Regional data completeness}}{\text{District data completes}} \times \,DTP-3\,coverage \tag{2}$$

Using the Chi-square test, we compared the adjusted DTP-1 and DTP-3 coverages in 2019 (pre-pandemic) and 2020 (pandemic). The 2020 data were considered the observed outcome, and the 2019 data reflected the expected outcome in the analysis. A *p*-value < 0.05 was considered statistically significant.
