The dates displayed for an article provide information on when various publication milestones were reached at the journal that has published the article activities on preceding journals at which the article was previously under consideration are not shown (for instance submission All content on this site: Copyright © 2025 Elsevier B.V. Metrics details vaccination has been consolidated as one of the most cost effective and successful public health interventions and a right of every human being This study aimed to assess the spatial dynamics of the vaccine coverage (VC) rate of children aged < 1 year per municipality in the Brazilian Northeast at 2016 and 2017 This is a mixed-type ecological study that use a Public domain data Health Information Vaccine doses were obtained from the Information System of the Brazilian National Immunization Program and live births from the Brazilian Information System of Live Births of the Brazilian Unified Health System Descriptive analysis of the coverage of all the vaccines for each year of the study was conducted and Mann–Whitney U test was used to compare VC between the study years Chi-squared test was used to evaluate the association between the years and VC Spatial distribution was analyzed according to both each study year and vaccine and presented as thematic maps Spatial autocorrelation was analyzed using Moran’s Global and Local statistics In the spatial analysis of the studied vaccines the Global Moran’s Index did not show any spatial autocorrelation (p > 0.05) but the Local Moran’s Index showed some municipalities and a positive influence on neighboring municipalities (p < 0.05) most municipalities with low VC were concentrated in the Mata Paraibano region negatively influencing their neighbors (p < 0.05) Uneven geographic regions and clusters of low VC for children aged < 1 year in the State of Paraíba were spatially visualized Health policy makers and planners need to urgently devise and coordinate an action plan directed at each state’s regions to fulfill the vaccination calendar thereby reversing the vulnerability of this age group which is at a higher risk of diseases preventable by vaccination the objective of this study is to use a geographic information system (GIS) and spatial analysis techniques to analyze the spatial dynamics of the VC rates and each vaccine’s spatial autocorrelation in children under the age of 1 year per municipality in the northeast Brazilian State in 2016 and 2017 This is a mixed-type ecological study with the municipality of residence as the analysis unit [9] This research is part of a broader project titled Spatial Analysis of Children’s Vaccine Coverage and its Relationship with Socioeconomic and Health Characteristics in Brazil funded by the Bill and Melinda Gates Foundation and the Brazilian National Council for Scientific and Technological Development (CNPq) Data were collected from the Department of Informatics of the Brazilian Unified Health System and publishing health information nationwide data from the Information System of the National Immunization Program (SI–PNI) and Information System on Live Births (SINASC) were used The applied doses of each vaccine for calculating VC were obtained from SI-PNI and data on live births were obtained from SINASC both collected by year and municipality of residence VC was calculated by antigen administered to children aged under 1 year as a fraction using either the number of doses applied (for single-dose vaccines) or the number of last doses applied (for multiple-dose vaccines) as the numerator and the live-birth data per year and municipality as the denominator The vaccines considered in the calculation of VC are part of the Brazilian National Vaccination Calendar for children aged under 1 year: Bacillus Calmette–Guérin (BCG) at birth; vaccine against Hepatitis B (HepB) at birth and at 2 and 6 months; vaccine against meningococcus type C (MnCc) at 3 and 7 months; vaccine against Diphtheria–Tetanus–Pertussis (DTP) vaccine against Haemophilus influenzae type B (HiB) The following formula was used to calculate the coverage of each vaccine: The PNI agreed upon the target of 90% of the population below the age of 1 year to be immunized with the BCG and rota vaccines Attaining these coverage rates is considered to be adequate and signifies that the target has been met A descriptive analysis of all the study variables was performed Quantitative variables were analyzed in terms of their central tendency and dispersion values Qualitative variables were analyzed in terms of their absolute and relative values and Chi-squared test was used to evaluate the association between the years and the categorized VCs USA) was used for the statistical analysis of the study data thematic maps were built from the calculation of VC for each vaccine by municipality and year using thematic cartography of Geographic Information System (GIS) The maps were built in QGIS 3.10 (QGIS Development Group) a free and open-source GIS that allows for the analysis of georeferenced data According to the definitions of thematic representation meaning above the target; and two light colors which translate into below target and far below target This study used the neighborhood matrix “w” of first order which signified that municipalities sharing a common physical border are considered as neighbors The connected regions are believed to interact more than unconnected ones and these connections are represented by matrix wherein a value of 1 signifies a common border and 0 signifies no border The spatial autocorrelation analysis was initially performed using the Global Moran’s Index (I) which identified the State of Paraíba as a unique study area and allowed a general measurement of the spatial association as it does not allow for the detailed analysis of spatial patterns the identification of the spatial correlations between the municipalities the Local Index of Spatial Association (LISA) was used to identify areas with values of similar attribute (clusters) and they were visualized through Moran’s Map which indicated a municipality with high VC and a positive influence on its neighbors (i.e. neighbors with high VC); “low-low” was identified in blue and indicated a municipality with low VC and a negative influence on its neighbors (i.e. municipalities with low VC); “high-low” was identified in green and indicated the municipalities with high VC surrounded by low VC municipalities; “low–high,” indicated in light blue included the municipalities with low VC surrounded with cities with high VC; and nonsignificant was indicated in white and included municipalities with no statistically significant spatial autocorrelation (p ≤ 0.05) and Moran’s Map analyses were performed using the free and open-source software R Studio (R Development Core Team 2019) with the spatial autocorrelation and LISA tools In this study, the VC data, one of the PNI performance indicators, were analyzed for 2016 and 2017 in the 223 municipalities of the State of Paraíba (Table 1) Paraiba’s total CV was 50.10% and in 2017 it was 70.08% Spatial analysis using the spatial autocorrelation indicator Global Moran’s Index found no statistically significant spatial autocorrelation in the two years (Table 2) When calculating the Spatial Association Index it is possible to perform the local spatial autocorrelation visualize the municipalities with similarities and form clusters of high and low VC which consequently positively or negatively influence their neighbors Vaccine coverage for Bacillus Calmette–Guérin (BCG), Diphtheria–Tetanus–Pertussis (DTP), hepatitis B, and Haemophilus influenzae type B (HiB) in 2016 and 2017, State of Paraíba, Brazil Moran’s Map of the vaccines given to children under 1 year of age some agglomerations of Sertão Paraibano municipalities were seen including one cluster toward the north and another in the center of the region with high VC and high similarity and a third cluster toward the east with low VC and a negative influence on their neighbors except for DTP in one municipality of the Mata Paraibano which showed low VC and a negative influence on its neighbors Figure 4 shows that the polio, rota, and pneumo vaccines had completely different spatial patterns in the two years analyzed. These vaccines must be simultaneously administered with the application of DTP, HiB, and HepB, according to the child vaccination calendar recommended by the Brazilian PNI. Moran’s Map of the vaccines given to children aged under 1 year João Pessoa state capital city showed a lower VC and a negative influence on its neighbors a few municipalities showed spatial autocorrelation (p < 0.05) with high VC and a positive influence on their neighbors A cluster of municipalities can also be seen in the northwest and the south of the Mata Paraibana region and in the southeast of the Agreste Paraibano region with low rota VC and a negative influence on their neighbors the municipalities with spatial autocorrelation (p < 0.05) included a high VC cluster toward the north of Sertão Paraibano but not influence its neighbors A few municipalities in the Mata Paraibana and Agreste regions showed low rota VC and a negative influence on their neighbors showed high pneumo VC and a positive influence on their neighbors A cluster of municipalities in the Mata Paraibana region was observed with low pneumo VC and a negative influence on their neighbors 13 municipalities showed statistical significance (p < 0.05) but only two municipalities showed high VC and a positive influence on their neighbors both in the north of the Sertão Paraibano region one in the Agreste Paraibano region and another in the Sertão Paraibano region showed low pneumo VC and a negative influence on their neighbors a few municipalities with positive spatial autocorrelation (p < 0.05) were seen and only two municipalities with high MnCc VC and a positive influence on their neighbors both in the Sertão Paraibano region The Mata Paraibana showed spatial autocorrelation (p < 0.05) with low MnCc VC and a negative influence on their neighbors but none had high MnCc VC with a positive influence on its neighbors a considerable number of municipalities with low or very low VC according to the PNI targets were observed for all the vaccines analyzed in this study 49.3% and 46.7% of the municipalities had VC < 50% in 2016 and 2017 and only 15.2% of the municipalities had adequate VC in both the years When analyzing the spatial distribution of all the evaluated vaccines a large number of municipalities with no statistical significance were found further corroborating that VC in the Paraíba region is far below the targets set by PNI From all the results of the spatial analysis the Sertão Paraibano region had the most municipalities with adequate VC that region showed the largest number of spatial clusters with high VC and high similarity presented clusters of low VC and a negative influence on neighboring municipalities which analyzed areas with decreased VC in Brazil from 2006 to 2016 it was shown that BCG and polio had the lowest VC rates in 2016 The study also demonstrated the heterogeneous spatial distribution of the drop in VC among the country’s various regions the study’s 2016 results agreed with those of the present study which reported VC to be below target in both the years the 2016 polio VC data in the above study are not similar with that of the present study it was observed that the Sertão Paraibano region had the largest number of municipalities showed high VC and high similarity with their neighbors we have been observed a low VC clusters negatively influencing the neighboring municipalities identifying areas with low VC is important to avoid epidemics of diseases preventable by vaccination and the implementation of immunization strategies no studies could be found that analyzed VC in children aged less than 1 year using Moran’s Index to measure spatial autocorrelation It is important to highlight the possibilities presented by this study including the formation of a database that allowed for a more accurate VC calculation and visualization of the spatial distribution of “vaccination clusters” through spatial analysis assessing VC through spatial analysis is an innovative approach in health care as few Brazilian and international studies have addressed this topic that is relevant to the worldwide public health To resolve possible inconsistencies in the data regarding the doses administered in vaccination facilities VC for each vaccine was calculated by extracting the SI–PNI data by the place of residence We declare that the data used are from the public domain health (link: https://datasus.saude.gov.br) and were obtained according to the criteria of good research practice and ethical precepts Assessment report of the Global Vaccine Action Plan - Strategic Advisory Group of Experts on Immunization Programa Nacional de Imunizações - vacinação Classificação de risco de transmissão de doenças imunopreveníveis a partir de indicadores de coberturas vacinais nos municípios brasileiros Departamento de Análise em Saúde e Vigilância de Doenças não Transmissíveis Avaliação dos indicadores de desempenho da vacinação do Programa Nacional de Imunizações e os desafios para elevar as coberturas vacinais no Brasil Saúde Brasil 2019: uma análise da situação de saúde com enfoque nas doenças imunopreveníveis e na imunização Abordagens espaciais na saúde pública: Série Capacitação e Atualização em Geoprocessamento em Saúde Ecologic studies in epidemiology: concepts Resolução no 510 de 7 de abril de 2016 que trata das especificidades éticas das pesquisas nas ciências humanas e sociais e de outras que utilizam metodologias próprias dessas áreas CONEP. Conselho Nacional de Saúde.RESOLUÇÃO No 580, DE 22 DE MARÇO DE 2018. Disponível em: https://conselho.saude.gov.br/resolucoes/2018/Reso580.pdf CONEP. Conselho Nacional de Saúde.RESOLUÇÃO No 466, DE 12 DE DEZEMBRO DE 2012. Disponível em: https://bvsms.saude.gov.br/bvs/saudelegis/cns/2013/res0466_12_12_2012.html Instituto Brasileiro de Geografia e Estatística (IBGE) Mapas: bases e referenciais: bases cartográficas: malhas digitais: municipal Sistemas e ciência da informação geográfica [recurso  eletrônico]/ Paul A et al.]; revisão técnica: Heinrich Hasenack Editado também como livro impresso em 2013 Sistemas e Ciência da Informação Geográfica (p Barata RB, de Almeida Ribeiro MCS, de Moraes JC, Flannery B. Socioeconomic inequalities and vaccination coverage: results of an immunisation coverage survey in 27 Brazilian capitals, 2007–2008. J Epidemiol Community Heal. 2012;66:934–41. https://doi.org/10.1136/JECH-2011-200341 Ferreira VLR, Waldman EA, Rodrigues LC, Martineli E, Costa ÂA, Inenami M, et al. Avaliação de coberturas vacinais de crianças em uma cidade de médio porte (Brasil) utilizando registro informatizado de imunização. Cad Saude Publica. 2018;34:e00184317. https://doi.org/10.1590/0102-311X00184317 Nóvoa T d’Avila, Cordovil VR, Pantoja GM, Ribeiro MES, Cunha AC dos S, Benjamin AIM, et al. Cobertura vacinal do programa nacional de imunizações (PNI)/Vacinal coverage of the national immunization program (PNI). Brazilian J Heal Rev. 2020;3:7863–73. https://doi.org/10.34119/BJHRV3N4-053 Oliveira G da S, Bitencourt EL, Amaral PFF, Vaz GP, Júnior PMR, Costa SB da. Cobertura vacinal: uma análise comparativa entre os estados da Região Norte do Brasil. Rev Patol do Tocantins. 2020;7:14–7. https://doi.org/10.20873/UFT.2446-6492.2020V7N1P14 Queiroz LLC, Monteiro SG, Mochel EG, Veras MA de SM, de Sousa FGM, Bezerra ML de M, et al. Cobertura vacinal do esquema básico para o primeiro ano de vida nas capitais do Nordeste brasileiro. Cad Saude Publica. 2013;29:294–302. https://doi.org/10.1590/S0102-311X2013000200016 Khan J, Shil A, Prakash R. Exploring the spatial heterogeneity in different doses of vaccination coverage in India. PLoS ONE. 2018;13:e0207209. https://doi.org/10.1371/JOURNAL.PONE.0207209 Vyas P, Kim D, Adams A. Understanding Spatial and Contextual Factors Influencing Intraregional Differences in Child Vaccination Coverage in Bangladesh. Asia-Pacific J Public Heal. 2019;31:51–60. https://doi.org/10.1177/1010539518813604 Joy TM, George S, Paul N, Renjini BA, Rakesh PS, Sreedevi A. Assessment of vaccine coverage and associated factors among children in urban agglomerations of Kochi, Kerala. India J Fam Med Prim Care. 2019;8:91. https://doi.org/10.4103/JFMPC.JFMPC_276_18 Yourkavitch J, Burgert-Brucker C, Assaf S, Delgado S. Using geographical analysis to identify child health inequality in sub-Saharan Africa. PLoS ONE. 2018;13:e0201870. https://doi.org/10.1371/JOURNAL.PONE.0201870 Applying an equity lens in the Decade of Vaccines Spatial analysis of vaccination coverage in Paraíba in 2016 Download references This study was funded by the Bill & Melinda Gates Foundation National Council for Scientific and Technological Development (CNPq) and Brazilian Ministry of Health in the so-called Grand Challenges Explorations – Brazil: Data Science To Improve Maternal and Child Health in Brazil Number of CNPq No 443790/2018–3; Fundação Bill & Melinda Gates OPP 1202115 Catholic University of Santos (Universidade Católica de Santos – Programa de Pós- Graduação strictu senso em Saúde Coletiva) Ysabely de Aguiar Pontes Pamplona & Lourdes Conceição Martins Hospital das Clínicas - Faculdade de Medicina da Universidade de São Paulo Nairmara Soares Pimentel Cunha & Sylvia Costa Lima Fahrat design of the work; the acquisition of database interpretation of data; and in writing the manuscript contributed with the analysis and interpretation of the data substantial review of the manuscript and interpretation of the data contributed with the analysis and interpretation of the data and substantial revision of the manuscript interpretation of data; and revised the manuscript substantial contributions to the conception interpretation of data and revised the manuscript All authors read and approved the final manuscript We are informed that all the methods were performed in accordance with the relevant guidelines and regulation The authors declare that they have no competing interests Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations unless otherwise stated in a credit line to the data Download citation DOI: https://doi.org/10.1186/s12889-022-13589-9 Anyone you share the following link with will be able to read this content: a shareable link is not currently available for this article