Skip to main content

Frequency, inequalities and spatial distribution of oral health services utilization in Peruvian children under twelve years of age: a population-based comparative analysis of the years 2017 and 2021

Abstract

Background

Oral diseases pose a significant public health challenge among Peruvian children. However, oral health services utilization among them is marked by inequalities, which may have been exacerbated by the COVID-19 pandemic. The objective was to compare the frequency, inequalities, determinants of inequality, and spatial distribution in oral health services utilization in Peruvian children under 12 years of age in 2017 and 2021.

Methods

Comparative secondary data analysis from the Demographic and Family Health Survey (ENDES) for the years 2017 (38,787 minors) and 2021 (36,729 minors). Age-adjusted prevalence ratios (PR) with 95% confidence intervals (CI) were calculated to assess the change in oral health services utilization in the last 12 months between the years 2017 and 2021, stratifying by demographic and socioeconomic covariates. Inequality was assessed by decomposing the Erreygers Concentration Index (ECI) and calculating the percentage contribution to inequality of the independent variables. Spatial analysis was performed using spatial autocorrelation methods, Hot Spot Analysis, Cluster and Outlier analysis, Inverse Distance Weighting, ordinary Kriging, and Kulldorff analysis.

Results

The probability of Peruvian children under 12 years of age using oral health services in the last 12 months was reduced by 45% in 2021 compared to 2017. The ECI identified a significant reduction in oral health services utilization inequalities at the national level (Diff: -0.0963; p < 0.001). The main contributor to inequality was higher tutor education level (55.2% in 2017 and 82.7% in 2021). In the comparison of spatial distribution, there was a greater dispersion of the conglomerates in which the use of oral health services is concentrated in 2021.

Conclusions

The frequency of oral health services utilization in the Peruvian children under 12 years of age was halved between 2017 and 2021. This problem is transversal to the entire population at the demographic and socioeconomic level. The key factor contributing to inequalities in the utilization of oral health services was the higher educational attainment of caregivers or guardians. Despite the improvement observed in inequalities and spatial distribution of the concentration of oral health services utilization, it is necessary to keep monitoring these patterns to guide decision-making.

Peer Review reports

Background

Oral diseases are one of the main non-communicable diseases worldwide, the most frequent being dental caries, periodontal disease, edentulism and oral cancer [1]. According to the Global Oral Health Status Report, almost half of the population (approximately 3.5 billion people) suffer from oral diseases throughout their lives. Likewise, 3 out of 4 people affected by these pathologies reside in low- and middle-income countries (LMIC) [2]. Additionally, since 1990 total disability-adjusted life years due to oral conditions have increased 64% to 16.9 million in 2015 [3]. Furthermore, LMIC have the lowest coverage and the highest socioeconomic inequality in terms of oral health care, which limits the opportunities for timely preventive and recuperative care [4, 5].

Children, especially those living in LMIC, have the highest prevalence of oral diseases in the world [1, 6, 7]. In Peru, in 2013, in the population of children affiliated with the Health Social Security (contributory insurance applicable to formal workers and their dependents), the prevalence of dental caries was 79.8% in children aged 3 to 5 years, and 90.4% in those of 12 years of age [8]. Likewise, according to the Global Oral Health Status Report, in 2019, the prevalence of untreated caries was estimated between 41.4 and 45.8% in deciduous teeth in children aged 1 to 9 years, and between 35.6 and 40.6% in permanent teeth of individuals 5 or more years of age in Peru [2]. In this sense, oral diseases in the population of Peruvian children are an important public health problem that needs urgent action on the part of the government, mainly to guarantee adequate access to free preventive and recuperative oral health treatments.

Since 2009, universal health insurance has been implemented in Peru, and given that this includes oral health care, increases in the use of oral health services and decreases in related inequality were evidenced in the period 2004–2017 [9]. Also, according to the Demographic and Family Health Survey (ENDES, by its acronym in Spanish), there is evidence of an increase over the years in the frequency of oral health services utilization in children under 12 years of age, taking into account the last six months (2014: 27.6% vs. 2019: 31.0%). However, during the COVID-19 pandemic, a decrease was reported in this percentage (2020: 19.6%) [10, 11]. Likewise, it has previously been reported that Peruvian children under 12 years of age belonging to the wealthiest quintiles, had a higher probability of using oral health services in 2019, a scenario that may have been exacerbated by the COVID-19 pandemic, increasing inequalities in oral health [12]. Finally, the pandemic has hit Peru particularly hard as it is a developing country with a limited capacity to respond to this type of event, insufficient access to basic services, high rates of poverty and unemployment [13], and spending on health that represented only 5.22% of the Gross Domestic Product in 2019, being below other countries in the Latin American and the Caribbean region and the regional average (7.88%) [14].

To date, no study assessing inequalities in oral health services utilization before and after the pandemic or the factors contributing to inequalities in this area has been conducted in Peruvian infants. Therefore, the objective of this study was to evaluate the frequency, inequalities, and determinants of inequality in the use of oral health services in children under 12 years of age, based on a comparative analysis of information from the ENDES for the years 2017 and 2021. Additionally, we carried out a comparative spatial analysis to project the distribution and clusters of oral health services utilization in the same years. Our results present an updated overview of oral health services utilization among Peruvian infants, offering a valuable foundation for the development of policies and interventions aimed at enhancing their oral health. Moreover, our study underscores the importance of obtaining more comprehensive, nationally representative data to facilitate further research on this topic.

Methods

Study context

In 2021, Peru had a population of approximately 33 million inhabitants. The territorial surface of Peru is 1,285,215.60 km2. The country’s landscape is diverse and geographically divided into 24 departments and one constitutional province (Callao). These administrative divisions are further categorized into three natural regions: the coast, the highlands, and the jungle [15, 16]. The coastal region is comprised by the departments of Lima, La Libertad, Piura, Ica, Lambayeque, Callao, Tumbes, Tacna, and Moquegua. The highland departments include Cusco, Pasco, Ayacucho, Huánuco, Junín, Cajamarca, Arequipa, Ancash, Apurímac, Huancavelica, and Puno. Finally, the jungle departments are San Martín, Ucayali, Loreto, Madre de Dios, and Amazonas [17].

The coastal region of Peru stretches along the Pacific Ocean and is home to nearly 60% of the country’s total population [18]. This region is known for its arid and desert-like landscapes, as well as for housing the majority of urban areas, including cities such as Lima, the capital of Peru. The coastal region is the country’s economic hub, foresting a diverse range of industries, international trade, and dynamic urban development. The highlands, dominated by the majestic Andes mountain range, traverse the heart of the country, defining the central part of Peru [18]. This region is distinguished by its rugged terrain, encompassing lofty plateaus at high altitudes. Moreover, it is home to nearly 30% of the Peruvian population, including numerous indigenous communities. The jungle, also known as the Amazon rainforest, covers a significant portion of Peru’s territory, primarily in the eastern part of the country. This region is renowned for its great biodiversity, dense vegetation, and ecosystems [18].

In terms of socioeconomic status, during 2020, more than 75% of individuals in rural areas experienced multidimensional poverty, whereas in urban areas, only 25% of individuals were affected [19]. Regarding the natural regions in the same year, monetary poverty in the coastal region was 25.9%, which was below the national average of 30.1% [19]. Additionally, the highlands were found to be a poorer region compared to the jungle (37.4% vs. 31.0%). However, when considering multidimensional poverty, the situation is reversed (highlands 17.3% vs. jungle 38.5%) [19]. This indicates that individuals in the jungle region may have a higher level of affordability compared to their counterparts in the highlands, but they also experience a greater degree of non-monetary deprivations [19].

Study design and data source

This was a cross-sectional study that used information from the ENDES for the years 2017 and 2021 for comparative purposes. The ENDES is a nationally representative survey carried out annually in Peru by the National Institute of Statistics and Informatics (INEI, by its acronym in Spanish), which collects information on demographic dynamics and the health status of mothers and children under 5 years of age, associated factors with communicable and non-communicable diseases, as well as access to diagnostic and treatment services for all children under 12 years of age [20, 21].

Sampling and data collection process

The ENDES has a complex, two-stage, probabilistic, balanced, independent sampling, stratified at the departmental level by urban and rural areas. The sampling frame is obtained from the statistical and cartographic information provided by the National Population and Housing Censuses and updated cartographic material from ENDES. The primary sampling unit at the urban level is the conglomerates and private dwellings, and, at the rural level, the rural census areas and private dwellings. The research unit for this study comprises individuals residing in both in urban and rural areas, who have spent the night prior to the application of the survey in the selected dwelling. Further specifications on the survey methodology are available in the ENDES technical reports [20, 21].

After removing the records with incomplete information on the variables of interest for the study, the final sample consisted of 38,787 and 36,729 Peruvian children under 12 years of age in the years 2017 and 2021, respectively (See Fig. 1).

Fig. 1
figure 1

Flow chart of participants included in the study

Measures

Outcome

The dependent variable was oral health services utilization in the last 12 months (This variable is categorical dichotomous (yes/no), constructed from the variables: QS803 (Have you ever been attended in a dental service by a dentist?); QS804U (Unit of time—How long ago was the last time you were seen in a dental service?), and QS804C (How long ago was the last service).

Independent variable for the analysis of inequalities

Characterization of the socioeconomic level was carried out based on the home well-being index (HV271). This variable is measured continuously and is calculated in the ENDES database and follows the methodology established by the Demographic and Health Surveys Program, which is based on the characteristics of the home and household assets.

Covariates

The following covariates were included based on previous studies, the PROGRESS framework [22], and the availability of survey variables: age (until 2 years/3 to 5 years/6 to 11 years); sex (women/men); area of residence (urban/rural); wealth index (richest/richer/middle/poorer/poorest); natural region (coast/highlands/jungle); health insurance (yes/does not have, does not know), yes is applicable when children are covered by at least one of the following types of health insurance: Integral Health Insurance (a publicly funded system primarily for the impoverished), Social Security Health Insurance (for formal employed workers and their dependents), Armed Forces and Police Health System (for employees of these institutions and their dependents), or Private insurance; mother's language (native/no native); tutor educational level (higher/secondary/until primary); and department of residence (the 24 departments and 1 constitutional province that make up the Peruvian administrative political division).

Variables for the spatial analysis

The variables used for the spatial analysis were the latitude and longitude of the conglomerate in which the minor’s home was located. For each of these points, the number of minors with and without oral health services utilization in the last 12 months was counted. In ENDES 2017 and 2021, geolocation measurement is given by taking points in the global positioning system measured with a tablet one meter from the main door of the informant’s home. More information about the measurement of the geolocation process in the ENDES 2017 and 2021 can be found in the survey manual.

Statistical analysis

Stata 17.0 software was used for the statistical analysis (StataCorp, College Station, TX, USA) and Fig. 2 was created in Python programming language. All the analyses included the configuration of the ENDES complex sampling design, using the svy command.

A descriptive analysis of the characteristics of the sample was carried out, as well as the distribution of the outcome according to such characteristics, based on the report of absolute and relative frequencies for both years of the survey. Additionally, age-adjusted prevalence ratios (PR) were calculated (except for age when it was evaluated as an independent variable). The PRs and their 95% confidence intervals (CI) were calculated to assess the change in oral health services utilization in the last 12 months between the years 2017 and 2021, stratifying by demographic and socioeconomic covariates. In all cases, p values were reported, where the null hypothesis was the absence of significant differences in the frequency of the outcome between the years 2017 and 2021 for each category of the variables of interest.

The analysis of inequalities in our study employs the Erreygers Concentration Index (ECI), which was calculated using the ‘conindex’ command in Stata [23]. The ECI is a widely used measure that assesses the concentration of a certain variable of interest, such as oral health services utilization, relative to a socioeconomic variable within a population. Specifically, we utilize the ECI to evaluate the extent of inequality in oral health services utilization among Peruvian children under twelve years of age. To interpret the ECI, it is important to note that it varies in a range from -1 to 1. A value equal to 0 indicates perfect equality, indicating that oral health services utilization is distributed proportionally among all socioeconomic strata of the population. A negative value denotes a concentration of oral health services utilization among individuals with lower socioeconomic status. Conversely, a positive value indicates a concentration among individuals with higher socioeconomic status, implying inequality in the opposite direction. Additionally, a decomposition analysis of the ECI was conducted using a generalized linear model in order to establish the percentage of contribution of the independent variables to the inequality of the occurrence of the event studied. Likewise, the elasticity, ECI, and absolute and percentage contribution to inequality are reported for each independent variable. More details about the methodological aspects and interpretation of these estimates can be found in O'Donnell et al. [24].

Spatial analysis was conducted in accordance with the methodology outlined in a previous study [25]. The spatial autocorrelation analysis (Global Moran Index), Hot Spot Analysis (Getis-Ord Gi*), and Cluster and Outlier Analysis (Anselin Local Moran’s I) were performed using the ArcGIS Desktop version 10.5 geographic information system software (ESRI Inc., Redlands, CA, USA) [26, 27]. A total of 3163 and 3023 clusters were utilized for the years 2017 and 2021 respectively. In our study, the Moran’s Global Index assessed the overall pattern and trend of the cases that used oral health services to determine if they were clustered, scattered, or random. The Global Moran Index ranges between -1 and + 1, with positive value indicating spatial clustering, 0 indicating a randomly distributed pattern, and negative values indicating a sparse pattern. To ensure a continuous spatial distribution with values at any point in the Peruvian territory, a spatial analysis based on interpolation was carried out. The deterministic Inverse Distance Weighting (IDW) method and ordinary Kriging geostatistical method with a spherical variogram, fitted to the data with parameters obtained in R-3.6.3, were used [28, 29]. Interpolation involves predicting values at unsampled spatial locations based on the closest sampled data, resulting in estimates of oral health services utilization on a 1 × 1 km spatial resolution grid [30]. Lastly, using the SaTScan V10.1 software (Martin Kulldorff, Boston, MA, USA), a Kulldorff spatial exploration analysis (Bernoulli model) was used to identify clusters with a high concentration of oral health services utilization [31]. This analysis uses the creation of a circular window that scans the area of influence of the study, where the radius of the circle depends on the maximum size of the population considered for the clusters. A maximum of 25% was set, that is, a group can contain a maximum of 25% of the total population with cases of oral health services utilization. The results are reported as a log-likelihood ratio (LLR) and p-value, where it is interpreted that the risk within the window is different from that outside the window. The p-value is calculated through the simulation model of Monte Carlo with 999 replicates; a low p-value suggests that the identified cluster is significant and probably not the result of chance. If the LLR value is positive and high it indicates a possible spatial concentration of cases [32]. The SaTScan output for statistically significant clusters includes the location of the center of the scan window, the radius of the scan window, and the p-value of the cluster [32]. Statistically significant clusters were defined as those with a p < 0.05.

Ethical considerations

The conduct of the study did not require the approval of an ethics committee as it was an analysis of aggregated secondary data that is in the public domain and does not allow the identification of the evaluated participants.

Results

A total of 38,787 and 36,729 children under 12 years of age were included in the years 2017 and 2021, respectively. In both years (2017 and 2021), a higher frequency of participants was identified in the age group of 6 to 11 years (50.2%; 51.3%), male sex (50.7%; 50.9%), the urban area of residence (74.1%; 77.1%), the poorer wealth quintile (24.1%; 23.5%), the coastal region (56.3%; 57.8%), those with health insurance (83.8%; 88.7%), children with mothers with non-native languages (81.6%; 83.8%) and those with mothers who have secondary education level (45.9%; 47.1%) (See Table 1).

Table 1 Characteristics of the Peruvian children under 12 years included in the ENDES 2017 and 2021

The probability of oral health services utilization in the last 12 months in the population of Peruvian children under 12 years of age was reduced by 45% in 2021 compared to 2017. A reduction in the probability of oral health services utilization was observed between the years compared for all the categories of the characteristics evaluated (See Table 2). The greatest reduction in the probability of oral health services utilization (55%) was identified in children of the middle wealth quintile, those living in the coastal region (50%), and children with mothers who have a secondary (50%) or primary education level (50%) (See Table 2).

Table 2 Frequencies and prevalence ratios of access to oral health services during the last year in Peruvian children under 12 years of age between 2017 and 2021

The ECI identified significant reductions of oral health services utilization inequality at the population level in Peruvian children under 12 years of age (Diff: -0.0963; p < 0.001). Significant reductions in the ECI were also identified in all age groups, in males and females, in the area of rural residence, the highlands and jungle region, the insured and uninsured, native and non-native mother tongue, and the secondary and primary educational levels (See Table 3).

Table 3 Changes in inequalities measured by means of the concentration index of access to oral health services during the last year in Peruvian children under 12 years of age between 2017 and 2021

The main contributor to inequality in both years was the higher education level. In 2017, the higher education contribution percentage was 55.2%, with rural residency as the second largest contributor at 24.7%. In 2021, the majority of the inequality was explained by higher education, comprising 82.7%, while the age group of 3 to 5 years age group contributed 4.8% (See Table 4).

Table 4 Decomposition analysis of inequalities of access to oral health services during the last year in Peruvian children under 12 years of age in 2017 and 2021

At the departmental level, it was detected that the greatest percentage reduction in the use of oral health services between 2021 and 2017 occurred in Lambayeque, with 61%, followed by Moquegua, with 55%, and both, Ica and La Libertad with 54%. On the other hand, the departments that experienced the lowest percentage reductions were Loreto, with 13%, Amazonas, with 21%, and San Martín and Puno, with 28% (See Fig. 2).

Fig. 2
figure 2

Departmental 2017–2021 differences in oral health services utilization in Peruvian children under 12 years of age

Spatial autocorrelation of the clusters in the two years was positive and statistically significant in 2017 (Moran’s Index 0.0649, p < 0.001 and z-score 9.93) and 2021 (Moran’s Index 0.0878, p < 0.001 and z-score 16.21). In 2017, clusters with a high concentration of use of oral health services were identified in Loreto, as well as in departments of the central and southern highlands, such as Huánuco, Pasco, Junín, Huancavelica, Ayacucho, Apurímac, Cusco and Arequipa and the departments of Lima and Ica. In 2021, the departments of Huánuco, Huancavelica, Apurimac, Ayacucho and Cusco maintained conglomerates with a high concentration of oral health services. Likewise, unlike in 2017, clusters with a high concentration were identified in the departments of Puno, Tacna, Madre de Dios, Amazonas, San Martín, and Ancash (See Fig. 3). The analyses of the spatial distribution of oral health services utilization (IDW interpolation, spherical Kriging interpolation and SaTScan analysis) showed consistent patterns of cluster grouping with a high concentration of use of oral health services, in comparison with the Moran local analysis and Hot Spot analysis (See Figs. 4, 5, and 6). When spatially comparing oral health services utilization in the study years, a greater dispersion of the conglomerates in which the use of oral health services was concentrated was found in 2021. According to SaTScan analysis, the statistically significant clusters increased from 7 to 17 between 2017 and 2021, respectively. According to the Hot spot analysis, the department of Amazonas and Madre de Dios had a particular change in the increase in oral health services utilization in 2021 compared to 2017.

Fig. 3
figure 3

Spatial distribution of oral health services utilization using Local Moran and Hot Spot analysis. A Cluster and outlier analysis, 2017. B Cluster and outlier analysis, 2021. C Hot spot analysis Getis Ord Gi, 2017. D Hot spot analysis Getis Ord Gi, 2021

Fig. 4
figure 4

Spatial distribution of oral health services utilization using spatial interpolation IDW. A Interpolation IDW, 2017. B Interpolation IDW, 2021

Fig. 5
figure 5

Spatial distribution of oral health services utilization using spatial interpolation ordinary spherical Kriging. A Interpolation ordinary spherical Kriging, 2017. B Interpolation ordinary spherical Kriging, 2021

Fig. 6
figure 6

Spatial distribution of oral health services utilization using spatial SaTScan analysis. A Spatial SatScan analysis, 2017. B Spatial SatScan analysis, 2021

Discussion

In 2021, the probability of oral health services utilization by Peruvian children under 12 years of age decreased by 45% compared to 2017. A similar trend of reduced oral health services utilization has been identified in some ecological studies, encompassing the majority of municipalities in Brazil as well [33, 34]. This problem could be due to the fact that during the COVID-19 pandemic the World Health Organization issued a series of recommendations to reduce the application of aerosol-generating procedures in the context of oral health care as a measure to prevent the spread of the virus [35]. Consequently, at the national level, the Peruvian Ministry of Health restricted face-to-face care in cases of oral urgencies or emergencies in the public and private spheres, limiting the opportunities for primary prevention in the dental office [36]. Another contributing factor was the economic and labor uncertainty, which affected many families throughout Latin America during the pandemic [37]. This may have forced the prioritization of expenses (e.g., food, housing, and transportation) and the use of health insurance, leaving oral health aside, even after vaccination and the lifting of care restrictions. As a result of unemployment, the reduction in family income and the increase in food prices also brought with it a scenario of food insecurity, which affected Peru to a greater extent compared to other Latin American countries [38, 39]. Likewise, due to the fear of contagion, especially due to the high rate of excess mortality in Peru [40], it is possible that both patients and dentists postponed treatments or preventive visits [41,42,43,44].

Our results show that at the national level, there was a reduction in inequality in oral health services utilization between 2017 and 2021. This could be due to the increase in public health coverage since the implementation of universal health insurance in 2009, which is publicly funded with no cost to the poor [45]. Likewise, greater coverage of conditional transfer programs over the years should be highlighted, such as JUNTOS, which reaches districts with a high percentage of poverty [46], as well as the efforts to integrate the promotion of oral health in programs such as the control of growth and development (CRED, by its acronym in Spanish) of girls and boys under five years of age, implemented in 2017, and of mandatory application in public and private health establishments nationwide [47]. All these measures may have generated, from different fronts, a positive impact on the inequalities in oral health services utilization among different population groups. Nevertheless, the relatively diminished inequality in oral health services utilization between the study years appears contradictory, considering the economic repercussions of the COVID-19 pandemic, particularly on the socioeconomically disadvantaged segment of the population. This paradoxical phenomenon might be attributed to changes in patterns of oral health services utilization, resulting from pandemic-related restrictions. For instance, a study utilizing data from the Public Health System of Brazil reported a nationwide decline in the total number of oral health procedures between March and August of 2020, with more pronounced reductions observed in both individual and collective preventive procedures [33]. Unfortunately, the ENDES lacks information regarding the purpose of dental visits, impeding a specific evaluation. Nonetheless, we hypothesize that individuals with higher socioeconomic status may have sought oral health services less frequently for preventive or cosmetic purposes, while those with lower socioeconomic standing, who may have poorer overall health, may have sought dental care due to genuine emergencies or urgent needs more frequently. Indeed, another study conducted in Brazil, reported greater reductions in urgent dental care during the most stringent period of the pandemic in municipalities with higher human development indices compared to 2019 [34]. In the Peruvian context, the pandemic scenario could have created a temporary and deceptive reduction in inequalities in oral health services utilization. Further research is warranted to investigate and substantiate our hypothesis.

The higher educational level of the legal guardian was the main contributing factor to the inequality in oral health services utilization in children under 12 years of age. This could be explained in that, with more education, guardians have greater knowledge and understanding of the importance of oral health care for the general health of their children [48, 49]. Additionally, guardians with a higher level of education often have higher incomes, health insurance, and better financial security, allowing them to pay for quality dental services for their children [50]. These differences in knowledge, economic capacity, and access to information can result in a gap in oral health services utilization among children whose parents have a higher level of education. Therefore, it is important to address these inequalities through strategies, such as parent education programs through schools or the health system, which could help improve their knowledge and understanding of the importance of oral health care in their minors and how to access oral health services in the national territory [51].

On the other hand, the largest percentage reductions in oral health services utilization occurred in the coastal departments (Lambayeque, Moquegua, Ica and La Libertad) and the smallest reduction were observed in the jungle (Loreto, Amazonas and San Martín). This may be because the coastal regions had higher mortality figures during the COVID-19 pandemic compared to the jungle regions [52], and this may have generated fear of exposure to SARS-CoV-2 in residents of the coast, and consequently, a reduction in health care, including oral health care. In this regard, although a decrease in oral health services utilization has been observed across all departments, it is imperative to prioritize actions and allocate resources to implement and strengthen strategies that promote the utilization of oral health services, particularly in the departments most significantly affected. On the other hand, in the spatial analysis, it was evidenced that, compared to 2017, a greater number of regions were found with conglomerates with a high concentration of oral health services utilization in the year 2021, denoting a greater dispersion of the conglomerates in which oral health services utilization in the last year under evaluation was concentrated. This could be considered favorable since the concentration of use would no longer only occur in some areas of the national territory. On the other hand, this may be related to the percentage reductions in oral health services utilization identified in our study, which in 2017 may have been concentrated in some regions, mainly Lima, Ica and the central highlands, and, given that, during the pandemic, care was reduced to urgencies or emergencies throughout the territory, a dispersion on the use of services can be seen, probably because routine and/or preventive face-to-face care was not carried out during confinement. However, these results should be considered exploratory, and further spatial investigations and monitoring of concentration patterns are warranted in the future to inform decision-making processes.

The main limitation of our study is its cross-sectional design, which prevents the establishment of cause-and-effect relationships. In addition, there is the possible introduction of recall bias since the outcome of interest was assessed through self-reported data in the past. In addition, there is the possibility of social desirability bias, as respondents could have modified their answers to fit the standards of their peers. Also, this study performed an analysis of secondary data, limiting the scope of the evaluation to the variables available in the survey. As a result, it hindered the assessment of more insightful aspects of oral health services utilization patterns, such as the intensity of utilization or the specific reasons for dental attendance. As for the spatial analysis, the results were obtained according to the ENDES sampling conglomerates in which the respondents’ dwelling was located, since the ENDES reports the geographic coordinates in this way to preserve the identity of the participants, which could limit the accuracy of the spatial distribution. Despite these limitations, this study provides information on the frequency, inequalities, and spatial distribution on oral health services utilization in Peruvian children under 12 years of age, using the most up-to-date information available, and calculating estimates with national representativeness. Understanding of Peru’s socio-economic and geographic context is crucial for comprehending the challenges and opportunities related to access to oral health services. By acknowledging the country’s diverse regions, including the coast, highlands, and jungle, researchers and policymakers can develop context-specific approaches to improve oral health outcomes and ensure equitable access to oral healthcare throughout Peru.

Conclusions

Our findings show a worrying decrease in the frequency of oral health services utilization in the population of Peruvian children under 12 years of age between 2017 and 2021. The reduction was generalized in all the categories of characteristics evaluated, indicating a problem that affects all population groups. It is important that the Peruvian government reflects on these results and take measures to improve oral health services utilization in the child population, especially in the most vulnerable groups, such as economically disadvantaged children, who continue to be affected by inequalities in this area to date. Therefore, it is essential to monitor and strengthen the existing upstream strategies designed to increase oral health service utilization in Peruvian children. Additionally, we identified a reduction in inequalities regarding oral health services utilization and an improvement in the spatial distribution of this outcome at the national level. However, these patterns are likely influenced by pandemic-related restrictions, and further research is necessary to gain a better understanding of these paradoxical trends in dental attendance after the pandemic, which can inform decision-making processes. It is highly recommended that future studies thoroughly investigate the impact of the pandemic and inequalities considering the reasons for oral health services utilization as the main outcome. In light of this, we urge the responsible entities in Peru tasked with providing nationally representative data to prioritize the collection of more comprehensive information for oral health research.

Availability of data and materials

The data supporting the findings of this study are publicly available from the National Institute for Statistics and Informatics (INEI) website (https://proyectos.inei.gob.pe/microdatos/). The databases can be obtained by entering the “Query by surveys” tab, selecting in “Survey” the option “Demographic and Family Health Survey – ENDES”, for the years 2017 and 2021.

Abbreviations

LMIC:

Low and middle-income countries

ENDES:

Demographic and Family Health Survey

INEI:

National Institute of Statistics and Informatics

PR:

Prevalence ratios

CI:

Confidence interval

ECI:

Erreygers concentration index

IDW:

Inverse distance weighting

LLR:

Log-likelihood ratio

References

  1. Peres MA, Macpherson LMD, Weyant RJ, Daly B, Venturelli R, Mathur MR, et al. Oral diseases: a global public health challenge. Lancet. 2019;394:249–60.

    Article  PubMed  Google Scholar 

  2. World Health Organization. Global oral health status report: towards universal health coverage for oral health by 2030. Geneva: World Health Organization; 2022.

    Google Scholar 

  3. Kassebaum NJ, Smith AGC, Bernabé E, Fleming TD, Reynolds AE, Vos T, et al. Global, Regional, and National Prevalence, Incidence, and Disability-Adjusted Life Years for Oral Conditions for 195 Countries, 1990–2015: A Systematic Analysis for the Global Burden of Diseases, Injuries, and Risk Factors. J Dent Res. 2017;96:380–7.

    Article  PubMed  Google Scholar 

  4. Singh A, Antunes JLF, Peres MA. Socio-Economic Inequalities in Oral Health. In: Peres MA, Antunes JLF, Watt RG, editors. Oral Epidemiology. Cham: Springer International Publishing; 2021. p. 279–94.

    Chapter  Google Scholar 

  5. World Health Organization. Oral health. 2022. https://www.who.int/news-room/fact-sheets/detail/oral-health. Accessed 24 Feb 2023.

  6. Golikeri SS, Grenfell J, Kim D, Pae C. Pediatric Oral Diseases. Dent Clin North Am. 2020;64:229–40.

    Article  PubMed  Google Scholar 

  7. Vasireddy D, Sathiyakumar T, Mondal S, Sur S. Socioeconomic Factors Associated With the Risk and Prevalence of Dental Caries and Dental Treatment Trends in Children: A Cross-Sectional Analysis of National Survey of Children’s Health (NSCH) Data, 2016–2019. Cureus. 2021. https://doi.org/10.7759/cureus.19184.

    Article  PubMed  PubMed Central  Google Scholar 

  8. León FAO. Perfil epidemiológico de salud bucal en niños atendidos en el Seguro Social del Perú. Revista Odontología Pediátrica. 2014;13:94–103.

    Google Scholar 

  9. Hernández-Vásquez A, Bendezu-Quispe G, Azañedo D, Santero M. Use of oral health care services in Peru: trends of socio-economic inequalities before and after the implementation of Universal Health Assurance. BMC Oral Health. 2019;19:39.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Instituto Nacional de Estadística e Informática. Perú: Enfermedades no transmisibles y transmisibles, 2020. INEI; 2021. https://proyectos.inei.gob.pe/endes/2020/SALUD/ENFERMEDADES_ENDES_2020.pdf. Accessed 24 Feb 2023.

  11. Instituto Nacional de Estadística e Informática. Perú: Enfermedades no transmisibles y transmisibles, 2019. INEI; 2020. https://proyectos.inei.gob.pe/endes/2019/SALUD/ENFERMEDADES_ENDES_2019.pdf. Accessed 24 Feb 2023.

  12. Torres-Mantilla JD, Newball-Noriega EE. Factors associated with the use of oral health services in Peruvian children under the age of 12 years. Clin Exp Dent Res. 2023;9:230–9.

    Article  PubMed  Google Scholar 

  13. Schwalb A, Seas C. The COVID-19 Pandemic in Peru: What Went Wrong? Am J Trop Med Hyg. 2021;104:1176–8.

    Article  PubMed  PubMed Central  Google Scholar 

  14. The World Bank. Current health expenditure (% of GDP) - Peru. 2022. https://data.worldbank.org/indicator/SH.XPD.CHEX.GD.ZS?locations=PE&name_desc=false. Accessed 1 Mar 2023.

  15. Carrillo-Larco RM, Guzman-Vilca WC, Leon-Velarde F, Bernabe-Ortiz A, Jimenez MM, Penny ME, et al. Peru – Progress in health and sciences in 200 years of independence. Lancet Reg Health Am. 2022;7:100148.

    PubMed  Google Scholar 

  16. Instituto Nacional de Estadística e Informática. Perú: Proyecciones de Población por Departamento, Provincia y Distrito, 2018 - 2020. 2020.

    Google Scholar 

  17. Instituto Nacional de Estadística e Informática. Mapa de pobreza provincial y distrital 2013. 2015.

    Google Scholar 

  18. Instituto Nacional de Estadística e Informática. Perú: Perfil Sociodemográfico. Informe Nacional. 2018.

    Google Scholar 

  19. Banco Central de Reserva del Perú. Perú: Historia de dos pobrezas. 2022.

    Google Scholar 

  20. Instituto Nacional de Estadística e Informática. Ficha Técnica ENDES 2017. 2018. https://proyectos.inei.gob.pe/endes/recursos/endes_2017_ficha_tecnica.pdf.

  21. Instituto Nacional de Estadística e Informática. Ficha Técnica ENDES 2021. 2022. https://proyectos.inei.gob.pe/endes/2021/DOCUMENTOS/FICHA_TECNICA_ENDES_2021.pdf.

  22. O’Neill J, Tabish H, Welch V, Petticrew M, Pottie K, Clarke M, et al. Applying an equity lens to interventions: using PROGRESS ensures consideration of socially stratifying factors to illuminate inequities in health. J Clin Epidemiol. 2014;67:56–64.

    Article  PubMed  Google Scholar 

  23. O’Donnell O, O’Neill S, Van Ourti T, Walsh B. Conindex: Estimation of Concentration Indices. Stand Genomic Sci. 2016;16:112–38.

    Google Scholar 

  24. O’Donnell O, World Bank, editors. Analyzing health equity using household survey data: a guide to techniques and their implementation. Washington: World Ban; 2008.

  25. Hernández-Vásquez A, Vargas-Fernández R, Turpo Cayo EY. Determinants, inequalities and spatial patterns of diarrhea in the Peruvian under-five population: Findings from nationally representative survey data. Front Public Health. 2023;11:1170670.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wong DWS, Lee J, Lee J. Statistical analysis of geographic information with ArcView GIS and ArcGIS. Fully rev. & updated. Hoboken: John Wiley & Sons, Inc; 2005.

    Google Scholar 

  27. Anselin L. Local Indicators of Spatial Association-LISA. Geogr Anal. 2010;27:93–115.

    Article  Google Scholar 

  28. Lu GY, Wong DW. An adaptive inverse-distance weighting spatial interpolation technique. Comput Geosci. 2008;34:1044–55.

    Article  Google Scholar 

  29. Webster R, Oliver MA. Geostatistics for environmental scientists. 2nd ed. Chichester: John Wiley & sons; 2007.

    Book  Google Scholar 

  30. Cressie NAC. Statistics for spatial data. Revised. Hoboken, NJ: John Wiley & Sons, Inc; 2015.

    Google Scholar 

  31. Kulldorff M. A spatial scan statistic. Commun Stat Theory Methods. 1997;26:1481–96.

    Article  Google Scholar 

  32. Wangdi K, Penjor K, Tobgyal, Lawpoolsri S, Price RN, Gething PW, et al. Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination. IJERPH 2021;18:5553.

  33. Chisini LA, Costa FDS, Sartori LRM, Corrêa MB, D’Avila OP, Demarco FF. COVID-19 Pandemic impact on Brazil’s Public Dental System. Braz Oral Res. 2021;35:e082.

    Article  PubMed  Google Scholar 

  34. Silva HGE, Móra PMPK, Zajkowski LA, Celeste RK, Scarparo RK. Urgent dental care in the Brazilian public health system: learning lessons from the COVID-19 pandemic for future situations. Cad Saúde Pública. 2022;38:e00013122.

    Article  PubMed  Google Scholar 

  35. World Health Organization. Maintaining essential health services: operational guidance for the COVID-19 context: interim guidance; 1 June 2020. World Health Organization; 2020. https://iris.who.int/bitstream/handle/10665/332240/WHO-2019-nCoV-essential_health_services-2020.2-eng.pdf?sequence=1. Accessed 1 Mar 2023.

  36. Ministerio de Salud. Directiva Sanitaria N° 100/MINSA/2020/DGIESP. 2020.

    Google Scholar 

  37. Economic Commission for Latin America and the Caribbean (ECLAC). The sociodemographic impacts of the COVID-19 pandemic in Latin America and the Caribbean. ECLAC; 2022. https://repositorio.cepal.org/bitstream/handle/11362/47923/1/S2200158_en.pdf. Accessed 1 Feb 2023.

  38. Hernández-Vásquez A, Visconti-Lopez FJ, Chacón-Torrico H, Azañedo D. COVID-19 and Food Insecurity in Latin America and the Caribbean. J Hunger Environ Nutr. 2022;18:372–79.

  39. Hernández-Vásquez A, Visconti-Lopez FJ, Vargas-Fernández R. Factors Associated with Food Insecurity in Latin America and the Caribbean Countries: A Cross-Sectional Analysis of 13 Countries. Nutrients. 2022;14:3190.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Financial Times. Coronavirus tracker: the latest figures as countries fight the Covid-19 resurgence. Financial Times; 2021. https://www.ft.com/content/a2901ce8-5eb7-4633-b89c-cbdf5b386938. Accessed 28 Feb 2023.

  41. Campus G, Diaz Betancourt M, Cagetti M, Giacaman R, Manton D, Douglas G, et al. The COVID-19 pandemic and its global effects on dental practice. An International survey. J Dent. 2021;114:103749.

    Article  Google Scholar 

  42. Moraes RR, Cuevas-Suárez CE, Escalante-Otárola WG, Fernández MR, Dávila-Sánchez A, Grau-Grullon P, et al. A multi-country survey on the impact of COVID-19 on dental practice and dentists’ feelings in Latin America. BMC Health Serv Res. 2022;22:393.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Taravati S, Farokhnia M. The impact of parents’ fears of the SARS-CoV-2 pandemic on children’s visits to pediatric dental offices in Iran. J Family Med Prim Care. 2022;11:6472.

    Article  PubMed  PubMed Central  Google Scholar 

  44. González-Olmo MJ, Delgado-Ramos B, Ortega-Martínez AR, Romero-Maroto M, Carrillo-Díaz M. Fear of COVID-19 in Madrid. Will patients avoid dental care? Int Dent J. 2022;72:76–82.

    Article  PubMed  Google Scholar 

  45. Mezones-Holguin E, Amaya E, Bellido-Boza L, Mougenot B, Murillo JP, Villegas-Ortega J, et al. Cobertura de aseguramiento en salud: el caso peruano desde la Ley de Aseguramiento Universal. Rev Peru Med Exp Salud Publica. 2019;36:196.

    Article  PubMed  Google Scholar 

  46. Ministerio de Desarrollo e Inclusión Social. InfoJUNTOS. InfoJUNTOS - Nacional. 2023. http://www2.juntos.gob.pe/infojuntos/indexe.html. Accessed 2 Mar 2023.

  47. Ministerio de Salud. Norma Técnica Sanitaria N° 137-MINSA/2017/DGIESP. 2017.

    Google Scholar 

  48. Marquez-Arrico C, Almerich-Silla J, Montiel-Company J. Oral health knowledge in relation to educational level in an adult population in Spain. J Clin Exp Dent. 2019;11:e1143–150.

  49. Saldūnaitė K, Bendoraitienė EA, Slabšinskienė E, Vasiliauskienė I, Andruškevičienė V, Zūbienė J. The role of parental education and socioeconomic status in dental caries prevention among Lithuanian children. Medicina. 2014;50:156–61.

    Article  PubMed  Google Scholar 

  50. Chen L, Hong J, Xiong D, Zhang L, Li Y, Huang S, et al. Are parents’ education levels associated with either their oral health knowledge or their children’s oral health behaviors? A survey of 8446 families in Wuhan. BMC Oral Health. 2020;20:203.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Nakre P, Harikiran A. Effectiveness of oral health education programs: A systematic review. J Int Soc Prevent Communit Dent. 2013;3:103.

    Article  Google Scholar 

  52. Flores López MG, Soto Tarazona A, De La Cruz-Vargas JA. Regional distribution of COVID-19 mortality in Peru. RFMH. 2021;21:326–34.

    Article  Google Scholar 

Download references

Acknowledgements

The authors are grateful to Donna Pringle for reviewing the language and style of the manuscript.

Funding

Nil.

Author information

Authors and Affiliations

Authors

Contributions

DA and AHV came up with the research idea, designed the study, collected, processed and analyzed the data. AHV and EYTC performed the spatial distribution analysis. DA, AHV, FJVL and EYTC participated in the interpretation of the data, wrote the manuscript and approved the final version.

Corresponding author

Correspondence to Akram Hernández-Vásquez.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Azañedo, D., Hernández-Vásquez, A., Visconti-Lopez, F.J. et al. Frequency, inequalities and spatial distribution of oral health services utilization in Peruvian children under twelve years of age: a population-based comparative analysis of the years 2017 and 2021. BMC Oral Health 23, 778 (2023). https://doi.org/10.1186/s12903-023-03528-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12903-023-03528-4

Keywords