- Open Access
Local government responses for COVID-19 management in the Philippines
BMC Public Health volume 21, Article number: 1711 (2021)
Responses of subnational government units are crucial in the containment of the spread of pathogens in a country. To mitigate the impact of the COVID-19 pandemic, the Philippine national government through its Inter-Agency Task Force on Emerging Infectious Diseases outlined different quarantine measures wherein each level has a corresponding degree of rigidity from keeping only the essential businesses open to allowing all establishments to operate at a certain capacity. Other measures also involve prohibiting individuals at a certain age bracket from going outside of their homes. The local government units (LGUs)–municipalities and provinces–can adopt any of these measures depending on the extent of the pandemic in their locality. The purpose is to keep the number of infections and mortality at bay while minimizing the economic impact of the pandemic. Some LGUs have demonstrated a remarkable response to the COVID-19 pandemic. The purpose of this study is to identify notable non-pharmaceutical interventions of these outlying LGUs in the country using quantitative methods.
Data were taken from public databases such as Philippine Department of Health, Philippine Statistics Authority Census, and Google Community Mobility Reports. These are normalized using Z-transform. For each locality, infection and mortality data (dataset Y) were compared to the economic, health, and demographic data (dataset X) using Euclidean metric d=(x−y)2, where x∈X and y∈Y. If a data pair (x,y) exceeds, by two standard deviations, the mean of the Euclidean metric values between the sets X and Y, the pair is assumed to be a ‘good’ outlier.
Our results showed that cluster of cities and provinces in Central Luzon (Region III), CALABARZON (Region IV-A), the National Capital Region (NCR), and Central Visayas (Region VII) are the ‘good’ outliers with respect to factors such as working population, population density, ICU beds, doctors on quarantine, number of frontliners and gross regional domestic product. Among metropolitan cities, Davao was a ‘good’ outlier with respect to demographic factors.
Strict border control, early implementation of lockdowns, establishment of quarantine facilities, effective communication to the public, and monitoring efforts were the defining factors that helped these LGUs curtail the harm that was brought by the pandemic. If these policies are to be standardized, it would help any country’s preparedness for future health emergencies.
Since the emergence of the COVID-19 pandemic, the number of cases have already reached 82 million worldwide at the end of 2020. In the Philippines, the number of cases exceeded 473,000. As countries around the world face the continuing threat of the COVID-19 pandemic, national governments and health ministries formulate, implement and revise health policies and standards based on recommendations by world health organization (WHO), experiences of other countries, and on-the-ground experiences. Early health measures were primarily aimed at preventing and reducing transmission in populations at risk. These measures differ in scale and speed among countries, as some countries have more resources and are more prepared in terms of healthcare capacity and availability of stringent policies [1, 2].
During the first months of the pandemic, several countries struggled to find tolerable, if not the most effective, measures to ‘flatten’ the COVID-19 epidemic curve so that health facilities will not be overwhelmed [3, 4]. In responding to the threat of the pandemic, public health policies included epidemiological and socio-economic factors. The success or failure of these policies exposed the strengths or weaknesses of governments as well as the range of inequalities in the society [5, 6].
As national governments implemented large-scale ‘blanket’ policies to control the pandemic, local government units (LGUs) have to consider granular policies as well as real-time interventions to address differences in the local COVID-19 transmission dynamics due to heterogeneity and diversity in communities. Some policies in place, such as voluntary physical distancing, wearing of face masks and face shields, mass testing, and school closures, could be effective in one locality but not in another [7–9]. Subnational governments like LGUs are confronted with a health crisis that have economic, social and fiscal impact. While urban areas have been hot spots of the COVID-19 pandemic, there are health facilities that are already well in placed as compared to less developed and deprived rural communities . The importance of local narratives in addressing subnational concerns are apparent from published experiences in the United States , China [12, 13], and India .
In the Philippines, the Inter-Agency Task Force on Emerging Infectious Diseases (IATF) was convened by the national government in January 2020 to monitor a viral outbreak in Wuhan, China. The first case of local transmission of COVID-19 was confirmed on March 7, 2020. Following this, on March 8, the entire country was placed under a State of Public Health Emergency. By March 25, the IATF released a National Action Plan to control the spread of COVID-19. A community quarantine was initially put in place for the national capital region (NCR) starting March 13, 2020 and it was expanded to the whole island of Luzon by March 17. The initial quarantine was extended up to April 30 [5, 15]. Several quarantine protocols were then implemented based on evaluation of IATF:
Community Quarantine (CQ) refers to restrictions in mobility between quarantined areas.
In Enhanced Community Quarantine (ECQ), strict home quarantine is implemented and movement of residents is limited to access essential goods and services. Public transportation is suspended. Only economic activities related to essential and utility services are allowed. There is heightened presence of uniformed personnel to enforce community quarantine protocols.
Modified Enhanced Community Quarantine (MECQ) is implemented as a transition phase between ECQ and GCQ. Strict home quarantine and suspension of public transportation are still in place. Mobility restrictions are relaxed for work-related activities. Government offices operates under a skeleton workforce. Manufacturing facilities are allowed to operate with up to 50% of the workforce. Transportation services are only allowed for essential goods and services.
In General Community Quarantine (GCQ), individuals from less susceptible age groups and without health risks are allowed to move within quarantined zones. Public transportation can operate at reduced vehicle capacity observing physical distancing. Government offices may be at full work capacity or under alternative work arrangements. Up to 50% of the workforce in industries (except for leisure and amusement) are allowed to work.
Modified General Community Quarantine (MGCQ) refers to the transition phase between GCQ and the New Normal. All persons are allowed outside their residences. Socio-economic activities are allowed with minimum public health standard.
LGUs are tasked to adopt, coordinate, and implement guidelines concerning COVID-19 in accordance with provincial and local quarantine protocols released by the national government .
In this study, we identified economic and demographic factors that are correlated with epidemiological metrics related to COVID-19, specifically to the number of infected cases and number of deaths [17, 18]. At the regional, provincial, and city levels, we investigated the localities that differ with the other localities, and determined the possible reasons why they are outliers compared to the average practices of the others.
We categorized the data into economic, health, and demographic components (See Table 1). In the economic setting, we considered the number of people employed and the number of work hours. The number of health facilities provides an insight into the health system of a locality. Population and population density, as well as age distribution and mobility, were used as the demographic indicators. The data (as of November 10, 2020) from these seven factors were analyzed and compared to the number of deaths and cumulative cases in cities, provinces or regions in the Philippines to determine the outlier.
The Philippine government’s administrative structure and the availability of the data affected its range for each factor. Regional data were obtained for the economic component. For the health and demographic components, data from cities and provinces were retrieved from the sources. Due to the NCR exhibiting the highest figures in all key components, an investigation was conducted to identify an outlier among its cities. The z-transform
where x is the actual data, μ is the mean and σ is the standard deviation were applied to normalize the dataset. Two sets of normalized data X and Y were compared by assigning to each pair (x,y), where x∈X and y∈Y, its Euclidean metric d given by d=(x−y)2. Here, the Y’s are the number of COVID-19 cases and deaths, and X’s are the other demographic indicators. Since 95% of the data fall within two standard deviations from the mean, this will be the threshold in determining an outlier. This means that if a data pair (x,y) exceeds, by two standard deviations, the mean of the Euclidean metric values between the sets X and Y, the pair is assumed to be an outlier.
To identify a good outlier, a bias computation was performed. In this procedure, Y represents the normalized data set for the number of deaths or the number of cases while X represents the normalized data set for every factor that were considered in this study. The bias is computed using the metric
for all x in X and y in Y. To categorize a city, province, or region as a good outlier, the bias corresponding to this locality must exceed two standard deviations from the mean of all the bias computations between the sets X and Y.
Results and discussion
The data used were the reported COVID-19 cases and deaths in the Philippines as of November 10, 2020 which is 240 days since community lockdowns were implemented in the country. Figure 1 shows the different lockdowns implemented per province since March 15. It can be seen that ECQ was implemented in Luzon and major cities in the country in the first few weeks since March 15, and slowly eased into either GCQ or MGCQ as time progressed. By August, the most stringent lockdown was MECQ in the National Capital Region (NCR) and some nearby provinces. Places under MECQ on September were Iloilo City, Bacolod City, and Lanao del Sur, with the last province as the lone community to be placed under MECQ the month after. By November 1, 2020, communities were either placed under GCQ or MGCQ.
Comparison of economic, health, and demographic components and COVID-19 parameters
The economic, health and demographic components were compared to COVID-19 cases and deaths. These comparisons were done for different community levels (regional, provincial, city/metropolitan) (See Tables 2, 3, and 4). Figure 2 summarizes the correlation of components to COVID-19 cases and deaths at the regional level. In all components, correlations with other parameters to both COVID-19 cases and deaths are close. Every component except Residential Mobility and GRDP have slightly higher correlation coefficient for COVID-19 cases as compared to COVID-19 deaths.
Among the components, the number of ICU beds component has the highest correlation with COVID-19 parameters. This makes sense as this is one of the first-degree measures of COVID-19 transmission. Population density comes in second, followed by mean hours worked and working population, which are all related to how developed the region is economy-wise. Regions having larger population density also have a huge working population and longer working hours . Thus, having a huge population density implies high chance of having contact with each other [25, 26]. Another component with high correlation to the cases and deaths is the number of doctors on quarantine, which can be looked at two ways; (i) huge infection rate in the region which is the reason the doctors got exposed or are on quarantine, and (ii) lots of doctors on quarantine which resulted to less frontliners taking care of the infected individuals. All definitions of mobility and the GDP are not strongly correlated to any of the COVID-19 measures.
In each data set, outliers were identified depending on their distance from the mean. For simplicity, we denote components that are compared with COVID-19 cases by (C) and with COVID-19 deaths by (D). The summary of outliers among regions in the Philippines is shown in Figs. 3 and 4. Data is classified according to groups of component. In each outlier region, non-pharmaceutical interventions (NPI) implemented and their timing are identified.
Region III is an outlier in terms of working population (C) and the number of ICU beds (C) (see Fig. 5 and Table 5). This means that considering the working population of the region, the number of COVID-19 infections are better than that of other regions. Same goes with the number of ICU beds in relation to COVID-19 deaths. Region III is comprised of Aurora, Bataan, Nueva Ecija, Pampanga, Tarlac, Zambales, and Bulacan. This good performance might be attributed to their performance especially on their programs against COVID-19. As early as March 2020, the region had been under a community lockdown together with other regions in Luzon. Being the closest to NCR, Bulacan has been the most likely to have high number of COVID-19 cases in the region. But the province responded by opening infection control centers which offer free healthcare, meals, and rooms for moderate-severe COVID-19 patients . They have also implemented strict monitoring of entry-exit borders, organization of provincial task force and incident command center, establishment of provincial quarantine facilities for returning overseas Filipino workers, mandated municipal quarantine facilities for asymptomatic cases, and mass testing, among others . Most of which have been proven effective in reducing the number of COVID-19 cases and deaths .
Region IV-A is an outlier in terms of population and working population (D) and doctors on quarantine (D) (see Fig. 5 and Table 5). Considering their population and working population, the COVID-19 death statistics show better results compared to other regions. Same goes with the number of doctors in the region which are in quarantine in relation to the reported COVID-19 deaths. This shows that the region is doing well in terms of decreasing the COVID-19 fatalities compared to other regions in terms of populations and doctors on quarantine. Region IV-A is comprised of Batangas, Cavite, Laguna, Quezon, and Rizal. Same with Region III, they have been under the community lockdown since March of last year. Provinces of the region such as Rizal have been proactive in responding to the epidemic as they have already suspended classes and distributed face masks even before the nationwide lockdown . Despite being hit by natural calamities, the region still continue ramping up the response to the pandemic through cash assistance, first aid kits, and spreading awareness .
An interesting result is that NCR, the center of the country and the most densely populated, is a good outlier in terms of GRDP (C) and GRDP (D). Cities in the region launched various programs in order to combat the disease. They have launched mass testings with Quezon City, Taguig City, and Caloocan City starting as early as April 2020. Pasig City started an on-the-go market called Jeepalengke. Navotas, Malabon, and Caloocan recorded the lowest attack rate of the virus. Caloocan city had good strategies for zoning, isolation and even in finding ways to be more effective and efficient. Other programs also include color-coded quarantine pass, and quarantine bands. It is also possible that NCR may just have a very high GRDP compared to other regions. A breakdown of the outliers within NCR can be seen in Fig. 8.
Region VII is also an outlier in terms of population density (D) and frontliners (D) (see Fig. 6 and Table 5). This means that given the population density and the number of frontliners in the region, their COVID-related deaths in the region is better than the rest of the country. This region consists of four provinces (Cebu, Bohol, Negros Oriental, and Siquijor) and three highly urbanized cities (Cebu City, Lapu-Lapu City, and Mandaue City), referred to as metropolitan Cebu. This significant decline may be explained by how the local government responded after they were placed in stricter community quarantine measures despite the rest of the country easing in to more lenient measures. Due to the longer and stricter quarantine in Cebu, the lockdown had a greater impact here than in other areas where restrictions were eased earlier . Dumaguete was one of the destinations of the first COVID case in the Philippines , their local government was able to keep infections at bay early on. Siquijor was also COVID-19-free for 6 months . The compounded efforts of the different provinces in the region can account for the region being identified as an outlier.
Among the metropolitan cities, Davao came out as a good outlier in terms of population (C) and working population (C) (see Figs. 7, 9, and Table 5). This result may be attributed to their early campaign on consistent communication of COVID-19-related concerns to the public . They were also able to set up transportation for essential workers early on .
This study identified outliers in each data group and determined the NPIs implemented in the locality. Economic, health and demographic components were used to identify these outliers. For the regional data, three regions in Luzon and one in Visayas were identified as outliers. Apart from the minimum IATF recommended NPIs, various NPIs were implemented by different regions in containing the spread of COVID-19 in their areas. Some of these NPIs were also implemented in other localities yet these other localities did not come out as outliers. This means that one practice cannot be the sole explanation in determining an outlier. The compounding effects of practices and their timing of implementation are seen to have influenced the results. A deeper analysis of daily data for different trends in the epidemic curve is considered for future research.
Correlation tables, outliers and community quarantine timeline
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, Ren R, Leung KSM, Lau EHY, Wong JY, Xing X, Xiang N, Wu Y, Li C, Chen Q, Li D, Liu T, Zhao J, Liu M, Tu W, Chen C, Jin L, Yang R, Wang Q, Zhou S, Wang R, Liu H, Luo Y, Liu Y, Shao G, Li H, Tao Z, Yang Y, Deng Z, Liu B, Ma Z, Zhang Y, Shi G, Lam TTY, Wu JT, Gao GF, Phil D, Cowling BJ, Yang B, Leung GM, Feng Z. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N Engl J Med. 2020; 382(13):1199–207.
Hsiang S, Allen D, Annan-Phan S, Bell K, Bolliger I, Chong T, Druckenmiller H, Huang LY, Hultgren A, Krasovich E, Lau P, Lee J, Rolf E, Tseng J, Wu T. The effect of large-scale anti-contagion policies on the covid-19 pandemic. Nature. 2020; 584:262–67.
Anderson R, Heesterbeek JAP, Klinkenberg D, Hollingsworth T. Comment how will country-based mitigation measures influence the course of the covid-19 epidemic?Lancet. 2020; 395. https://doi.org/10.1016/S0140-6736(20)30567-5.
Buhat CA, Torres M, Olave Y, Gavina MK, Felix E, Gamilla G, Verano KV, Babierra A, Rabajante J. A mathematical model of covid-19 transmission between frontliners and the general public. Netw Model Anal Health Inform Bioinforma. 2021; 10. https://doi.org/10.1007/s13721-021-00295-6.
Ocampo L, Yamagishic K. Modeling the lockdown relaxation protocols of the philippine government in response to the covid-19 pandemic: an intuitionistic fuzzy dematel analysis. Socioecon Plann Sci. 2020; 72. https://doi.org/10.1016/j.seps.2020.100911.
Weible C, Nohrstedt D, Cairney P, Carter D, Crow D, Durnová A, Heikkila T, Ingold K, McConnell A, Stone D. Covid-19 and the policy sciences: initial reactions and perspectives. Policy Sci. 2020; 53:225–41. https://doi.org/10.1007/s11077-020-09381-4.
Wibbens PD, Koo WW-Y, McGahan AM. Which covid policies are most effective? a bayesian analysis of covid-19 by jurisdiction. PLoS ONE. 2020. https://doi.org/10.1371/journal.pone.0244177.
Mintrom M, O’Connor R. The importance of policy narrative: effective government responses to covid-19. Policy Des Pract. 2020; 3(3):205–27. https://doi.org/10.1080/25741292.2020.1813358.
Chin T, Kahn R, Li R, Chen JT, Krieger N, Buckee CO, Balsari S, Kiang MV. Us-county level variation in intersecting individual, household and community characteristics relevant to covid-19 and planning an equitable response: a cross-sectional analysis. BMJ Open. 2020; 10(9). https://doi.org/10.1136/bmjopen-2020-039886.
OECD. The territorial impact of COVID-19: managing the crisis across levels of government. 2020. https://www.oecd.org/coronavirus/policy-responses/the-territorial-impact-of-covid-19-managing-the-crisis-across-levels-of-government-d3e314e1/#biblio-d1e5202. Accessed 20 Feb 2007.
White ER, Hébert-Dufresne L. State-level variation of initial covid-19 dynamics in the united states. PLoS ONE. 2020; 15. https://doi.org/10.1371/journal.pone.0240648.
Lin S, Huang J, He Z, Zhan D. Which measures are effective in containing covid-19? — empirical research based on prevention and control cases in China. medRxiv. 2020. https://doi.org/10.1101/2020.03.28.20046110. https://www.medrxiv.org/content/early/2020/03/30/2020.03.28.20046110.full.pdf.
Mei C. Policy style, consistency and the effectiveness of the policy mix in China’s fight against covid-19. Policy Soc. 2020; 39(3):309–25. https://doi.org/10.1080/14494035.2020.1787627. http://arxiv.org/abs/https: //doi.org/10.1080/14494035.2020.1787627.
Dutta A, Fischer HW. The local governance of covid-19: disease prevention and social security in rural india. World Dev. 2021; 138:105234. https://doi.org/10.1016/j.worlddev.2020.105234.
Vallejo BM, Ong RAC. Policy responses and government science advice for the covid 19 pandemic in the philippines: january to april 2020. Prog Disaster Sci. 2020; 7:100115. https://doi.org/10.1016/j.pdisas.2020.100115.
Inter-Agency Task Force for the Management of Emerging Infectious Diseases. Omnibus guidelines on the implementation of community quarantine in the Philippines. 2020. https://doh.gov.ph/node/27640. Accessed 20 Feb 2020.
Roy S, Ghosh P. Factors affecting covid-19 infected and death rates inform lockdown-related policymaking. PloS ONE. 2020; 15(10):0241165. https://doi.org/10.1371/journal.pone.0241165.
Pullano G, Valdano E, Scarpa N, Rubrichi S, Colizza V. Evaluating the effect of demographic factors, socioeconomic factors, and risk aversion on mobility during the covid-19 epidemic in france under lockdown: a population-based study. Lancet Digit Health. 2020; 2(12):638–49.
Department of Health. COVID-19 tracker. 2020. https://doh.gov.ph/covid19tracker. Accessed 25 Nov 2020.
Authority PS. Philippine population density (based on the 2015 census of population). 2020. https://psa.gov.ph/content/philippine-population-density-based-2015-census-population. Accessed 11 Apr 2020.
Google. COVID-19 community mobility report. 2020; https://www.google.com/covid19/mobility?hl=en. Accessed 25 Nov 2020.
Authority PS. Labor force survey. 2020. https://psa.gov.ph/statistics/survey/labor-and-employment/labor-force-survey?fbclid=IwAR0a5GS7XtRgRmBwAcGl9wGwNhptqnSBm-SNVr69cm8sCVd9wVmcoKHRCdU. Accessed 11 Apr 2020.
Authority PS. https://psa.gov.ph/grdp/tables?fbclid=IwAR3dKvo3B5eauY7KcWQG4VXbuiCrzFHO4b-f1k5Od76ccAlYxUimUIaqs94. Accessed 11 Apr 2020. 2020.
Peterson E. The role of population in economic growth. SAGE Open. 2017; 7:215824401773609. https://doi.org/10.1177/2158244017736094.
Buhat CA, Duero JC, Felix E, Rabajante J, Mamplata J. Optimal allocation of covid-19 test kits among accredited testing centers in the philippines. J Healthc Inform Res. 2021; 5. https://doi.org/10.1007/s41666-020-00081-5.
Hamidi S, Sabouri S, Ewing R. Does density aggravate the covid-19 pandemic?: early findings and lessons for planners. J Am Plan Assoc. 2020; 86:1–15. https://doi.org/10.1080/01944363.2020.1777891.
Philippine News Agency. Bulacan shares anti-COVID-19 best practices. 2020. https://mb.com.ph/2020/08/16/bulacan-shares-anti-covid-19-best-practices/. Accessed Mar 2020.
Buhat CA, Villanueva SK. Determining the effectiveness of practicing non-pharmaceutical interventions in improving virus control in a pandemic using agent-based modelling. Math Appl Sci Eng. 2020; 1:423–38. https://doi.org/10.5206/mase/10876.
Hallare K. Cainta, Rizal suspends classes, distributes face masks over coronavirus threat. 2020. https://newsinfo.inquirer.net/1238217/cainta-rizal-suspends-classes-distributes-face-masks-over-coronavirus-threat. Accessed Mar 2020.
Relief International. Responding to COVID-19 in the Aftermath of Volcanic Eruption. 2020. https://www.ri.org/projects/responding-to-covid-19-in-the-aftermath-of-volcanic-eruption/. Accessed Mar 2020.
Macasero R. Averting disaster: how Cebu City flattened its curve. 2020. https://www.rappler.com/newsbreak/explainers/how-cebu-city-flattened-covid-19-curve/. Accessed Mar 2020.
Edrada EM, Lopez EB, Villarama JB, Salva-Villarama EP, Dagoc BF, Smith C, Sayo AR, Verona JA, Trifalgar-Arches J, Lazaro J, Balinas EGM, Telan EFO, Roy L, Galon M, Florida CHN, Ukawa T, Villaneuva AMG, Saito N, Nepomuceno JR, Ariyoshi K, Carlos C, Nicolasor AD, Solante RM. First covid-19 infections in the philippines: a case report. Trop Med Health. 2020; 48(30). https://doi.org/10.1186/s41182-020-00218-7.
Macasero R. Coronavirus-free for 6 months, Siquijor reports first 2 cases. 2020. https://www.rappler.com/nation/siquijor-coronavirus-cases-august-2-2020. Accessed Mar 2020.
Davao City. Mayor Sara, disaster radio journeying with dabawenyos. 2020. https://www.davaocity.gov.ph/disaster-risk-reduction-mitigation/mayor-sara-disaster-radio-journeying-with-dabawenyos. Accessed Mar 2020.
Davao City. Davao city free rides to serve GCQ-allowed workers. 2020. https://www.davaocity.gov.ph/transportation-planning-traffic-management/davao-city-free-rides-to-serve-gcq-allowed-workers/. Accessed Mar 2020.
JFR is supported by the Abdus Salam International Centre for Theoretical Physics Associateship Scheme.
This research is funded by the UP System through the UP Resilience Institute.
Ethics approval and consent to participate
Not applicable. We used secondary data. These are from the public database of the Philippine Department of Health (https://www.doh.gov.ph/covid19tracker) and Philippine Statistics Authority Census (https://psa.gov.ph)
Consent for publication
The authors declare that they have no competing interests.
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.
About this article
Cite this article
S. Talabis, D.A., Babierra, A.L., H. Buhat, C.A. et al. Local government responses for COVID-19 management in the Philippines. BMC Public Health 21, 1711 (2021). https://doi.org/10.1186/s12889-021-11746-0
- Local government
- Quantitative methods