Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Primary care physicians and infant mortality: Evidence from Brazil

  • Letícia Xander Russo,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Economics, State University of Maringá, Maringá Paraná, Brazil

  • Anthony Scott ,

    Roles Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

    a.scott@unimelb.edu.au

    Affiliation Melbourne Institute: Applied Economic and Social Research, The University of Melbourne, Melbourne, Victoria, Australia

  • Peter Sivey,

    Roles Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation School of Economics Finance and Marketing, RMIT University, Melbourne, Victoria, Australia

  • Joilson Dias

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

    Affiliation Department of Economics, State University of Maringá, Maringá Paraná, Brazil

Abstract

Primary health care has been recognized as a critical strategy for improving population health in developing countries. This paper investigates the effect of primary care physicians on the infant mortality rate in Brazil using a dynamic panel data approach. This method accounts for the endogeneity problem and the persistence of infant mortality over time. The empirical analysis uses an eight-year panel of municipalities between 2005 and 2012. The results indicate that primary care physician supply contributed to the decline of infant mortality in Brazil. An increase of one primary care physician per 10,000 population was associated with 7.08 fewer infant deaths per 10,000 live births. This suggests that, in addition to other determinants, primary care physicians can play an important role in accounting for the reduction of infant mortality rates.

Introduction

Health care systems with a strong primary care orientation are more likely to provide better population health and more equitable health outcomes [1]. Primary care is generally defined as dealing with a wide range of health problems in the community, health promotion, disease prevention and rehabilitation as set out in Alma-Ata [2]. Although studies show the importance of primary care to health outcomes, evidence of the effect of changes in the level of human resources for health is limited. The majority of previous studies investigating the association between primary care physicians and health outcomes are cross-sectional and derived from sources that do not allow testing of the effect of changes in primary care on changes in outcomes over time. To our knowledge, there is only two longitudinal study that analyzed the relationship between primary care physicians (PCPs) and mortality rates at disaggregated level over time. Basu et al. [3] found that an increase in primary care physician density was associated with lower mortality in the United States. Greater PCP density showed a higher effect on life expectancy than specialist physician density, 51.5 day and 19.2 day increase in life expectancy, respectively. Aakvik and Holmas [4] did not find a statistically significant relationship between the level of PCP supply and mortality in Norway. However, Norway has less inequality than the United States and Brazil, and so there may be little room for improvement.

Using longitudinal data, but at state level, Shi et al. [5] found that increased PCP supply led to reductions in infant mortality in US. They investigated the relationship between primary care and infant mortality in 50 US states for eleven years (1985–95) and control for Gini coefficient, racial composition, rural areas, unemployment and education.

In this paper, we investigate the effect of changes in primary care physician supply on improvements in infant mortality rates (IMR) in Brazilian municipalities between 2005 and 2012. Infant mortality is one of the most used indicators of population health. Moreover, infant mortality is sensitive to a country’s health care system as well as socioeconomic conditions [6].

There are some particular channels through which primary care may influence infant mortality during both the neonatal and post-neonatal period. In 2012, neonatal mortality (deaths between 0–27 days) represented 69.3% of overall infant mortality (under-1 year) in Brazil. Neonatal deaths are associated with the quality of prenatal care [7,8], which allows early identification and intervention to minimize damage to child health. Quality of prenatal care is influenced by access to health care facilities, numbers of prenatal consultations and qualified health workers including physicians [9].

Post-neonatal mortality (deaths between 28–364 days) for ambulatory care-sensitive conditions (ACSC) corresponded to 22% of all post-neonatal mortality in 2012, with a higher incidence of lobar pneumonia, intestinal infectious, acute bronchiolitis and malnutrition (DATASUS/Brazilian Ministry of Health). Shi et al. [5] pointed out two mechanisms through which primary care is associated with lower post-neonatal mortality. First, primary care can improve mother’s and child’s health and better manage conditions affecting post-neonatal age (for example, by providing information about infections and safety). Second, primary care can deal with a series of maternal risk factors linked to infant deaths, such as smoking, alcoholism, poor weight gain during pregnancy, sexually transmitted diseases and poor nutrition.

In Brazil, as in other developing countries, primary health care has been recognized as a central strategy for improving population health, as well as a key instrument to reduce child mortality and achieve the fourth of the Millennium Development Goals (MDG) set by United Nations. The MDG-4 aimed to reduce the under-5 mortality rate by two-thirds for the period from 1990 to 2015. Though the global goal for child mortality rate was not reached, significant achievement has been made in the world (global under-5 mortality has declined more than half) [10].

Brazil showed substantial progress in reducing child mortality over the period, falling by 73 per cent from 1990 to 2015. The number of under-5 deaths was 60.8 deaths per 1,000 live births in 1990 and 16.4 deaths per 1,000 live births in 2015 [11]. However, despite the reduction in the infant mortality rate in Brazil, the country still has relatively high rates of infant deaths, more than 2.4 times when compared with the average for high income countries.

We contribute to the literature through our innovative use of data on primary care physicians from the Brazilian Ministry of Health, which provides annual data on human resources for health for all municipalities since 2005, and a robust methodological approach, dynamic panel data analysis. Dynamic panel data models have also been used in previous studies for analyzing the determinants of infant mortality [12,13]. This method accounts for endogeneity and unobserved municipality-specific effects. Endogeneity may lead to an inconsistent estimator and generate misleading results [14]. This issue has been addressed by only few studies investigating the association between primary care and health outcomes [4,15,16]. In addition, the dynamic panel model improves the accuracy of estimates when the dependent variable is persistent over time. In our case, such persistence could be due to health habits, mother’s behaviours, and unobserved characteristics of the environment which are associated with child health. Furthermore, previous findings of inter-family heterogeneity that drives a positive correlation in infant mortality across siblings [17] suggest that a dynamic model is appropriate.

Our findings from the dynamic panel model demonstrate the effect of primary care physician supply on reducing infant mortality, even after including a range of important determinants of child health. The estimated results strengthen the case for the role of primary care in improving health in developing countries.

Primary care in Brazil

The Brazilian health system (Unified Health System—SUS) was founded in 1988 on the principles of universality and equality, as well as decentralized management and popular participation [18]. The SUS is funded primarily through taxes and offers free health services, from medical consultations to organ transplants.

In the 1990s, after the creation of SUS, efforts focused on universal access to primary care. The government adopted measures aimed at strengthening primary care incorporating principles such as first contact, continuity, comprehensiveness and coordination—as outlined by Starfield [19]. The implementation of the Family Health Program (now called Family Health Strategy—FHS) by the Ministry of Health in 1994, was a strategy to reorganize health care and establish a new dynamic in the organization of health services and preventive services [20].

Although a broad expansion of the FHS has occurred, the traditional model of primary care still remains. Physicians and staff in both settings are salaried and funded by government, but the models differ in the organization of care, in the composition of the teams and territorial definition of coverage areas.

Traditional healthcare centers include physicians (general practitioners, pediatricians and obstetrician-gynecologists), nurses, nurse assistants and dentists who respond reactively to patient demand. FHS units are composed of a team as specified by the Ministry of Health comprising a primary care physician, nurse, nurse assistant, community health workers, and sometimes dental health professionals. FHS teams are organized geographically and they are guided by the principles of responsibility for the health of the local population through proactive needs assessment. They are designed to bring primary care closer to the communities, including periodic visits to households and consultations at home as well as in healthcare centers [21]. The FHS is a model focused on preventive health care with a multi-professional team, while the traditional model prioritizes curative and individual measures. Currently, primary care is provided by both the FHS and traditional models.

Data and methods

Measures

The definition of primary care physicians can be broad and in the US can include family and general practitioners, general internists, and general pediatricians [22]. In this paper, primary care physicians are measured using the sum of the hours worked per week of FHS physicians, community physicians, general practitioners (GPs), general pediatricians and gynecologists/obstetricians who work in healthcare centers (healthcare centers include Centro de Saúde/Unidade Básica de Saúde, Posto de Saúde and Centro de Apoio a Saúde da Família—CASF). Healthcare centers are designed to provide primary health care to the population and they concentrate primary health teams. Physicians working in hospitals are not included. Our measure of primary care physician density is defined as the total number of PCPs per 10,000 population in each municipality. The total number of PCPs is calculated by dividing the total hours worked per week by 40 to arrive at a “full-time equivalent” measure of PCP supply. This measure allows more accuracy to compare the availability of physicians among municipalities, since physicians have diverse workloads in Brazil.

Based on previous literature [23,24], a set of variables known to be associated with infant mortality was included in the analyses as covariates: other FHS team members per 10,000 population, hospital beds per 10,000 population (public and private), percentage of population with supplemental health insurance (private), percentage of female illiteracy, Gross Domestic Product (GDP) per capita and four measures of household facilities (percentage of household with trash collection, available electricity, access to piped water and percentage of household with sewage).

The data are obtained from several databases. Infant mortality and hospital beds are from the DATASUS database at the Brazilian Ministry of Health. Primary care physician and other FHS team are collected from TabWin—program developed by DATASUS/Ministry of Health. Data on supplemental health insurance is taken from National Regulatory Agency for Private Health Insurance and Plans (ANS). Data on GDP and population are obtained from the Brazilian Institute of Geography and Statistics (IBGE). Data on household facilities (trash collection, sewage, available electricity and access to piper water) are taken from National Census Bureau. The data on female illiteracy is based on the Annual List of Social Information (RAIS) from the Ministry for Labor and Employment (MTE); this data is only available for female workers in the formal sector. All data were extracted from publicly sources from links listed in S1 Appendix. The data are also available from authors on request.

Other FHS team members include nurses and nurse assistants. Nurses of FHS are included as covariates based on their relevant role for primary care and infant health [25,26]. The total number of other FHS team members is calculated by dividing the total hours worked per week by 40.

Linear interpolation is used to estimate annual values of the four household facilities variables for the period 2005 to 2009. The estimates are based on the National Census (2000 and 2010). For the posterior period, 2011 and 2012, linear extrapolation is used based on the same function. This technique has been utilized to deal with the lack of annual data from Census [15,27]. Final results are stable when these variables are excluded. Annual values for all other variables are available from the database.

Descriptive statistics and data

Descriptive statistics of each variable over the whole period are shown in Table 1.

thumbnail
Table 1. Descriptive statistics for Brazilian municipalities (2005–2012).

https://doi.org/10.1371/journal.pone.0217614.t001

Table 2 reports the means of the variables for each time period. The infant mortality rate decreased 20% over period, while all of the household facilities recorded an improvement over the entire period. The number of primary care physicians per 10,000 population increased from 3.0 in 2005 to 3.7 in 2012. After the implementation of FHS in 1994 the number of FHS physicians and other FHS team members has increased, while the number of physicians of traditional models has decreased over the period.

thumbnail
Table 2. Variable means for Brazilian municipalities, by year (2005–2012).

https://doi.org/10.1371/journal.pone.0217614.t002

We use data from the Brazilian Ministry of Health (DATASUS—Department of the SUS), and follow a balanced panel of 5,563 municipalities in Brazil over a period of 8 years from 2005–2012 for the main analysis. Data on PCPs was available at municipal level only after 2005.

In addition, we performed sensitivity analyses for alternative samples. The reason for these alternative estimations is that the quality of mortality data may vary across regions of Brazil because of inadequate death registration systems, especially in the North and Northeast regions. These regional differences in death reporting may lead to the underestimation of the infant mortality rate especially in poorer areas and those with limited access to health care. These issues may not be as concerning for longitudinal data analyses. Nevertheless, the model estimates may be affected by improvements in the reporting system over time though this would be captured by the year fixed effects [28]. To ensure that the results from the main model are consistent and valid, we employed a method to reduce measurement error in the municipality-level mortality data. This method consists of analyzing the quality of information on births and deaths for all Brazilian municipalities and including only those municipalities who met a previously criteria defined by Frias et al [28].

Frias et al. [28] analyzed the adequacy of the reporting system using municipal indicators constructed with data from Ministry of Health and developed a method for selecting data from only the most reliable municipalities which has since been used by previous studies [27,29]. The Frias et al. [28] method validates the information about births and deaths at municipality level using five indicators. i) Age standardized mortality rate: values above 5 per 1,000 inhabitants indicate adequate death reporting. ii) Average annual deviation from the three-year period mean of the age standardized mortality rate: deviations of less than 20% indicate reliable information. iii) Ratio of reported live births to estimated live births: rates above 0.85 represent reliable information. iv) Average annual deviation from the three-year period mean of the birth rate: deviations below 20% indicate it was adequately reported. v) Proportion of ill-defined deaths (the use of imprecise terms and dubious expressions that prevent the identification and coding of cause of death): values below 20.7% for municipalities with population of less than 50,000 and values below 16.2% for municipalities with population of more than 50,000 indicate adequate death reporting [30].

The sensitivity analyses were performed for two alternative estimation samples: (A) Frias et al. [28] criteria in a balanced panel (5,110 observations from 730 municipalities) and (B) Frias et al. [28] criteria in an unbalanced panel (14,095 observations from 3,850 municipalities).

Method

The empirical analysis of this study uses the System Generalized Method of Moments (GMM) estimator developed by Arellano and Bover [31] and Blundell and Bond [32] for dynamic models. Consider the dynamic equation represented below (1) where imit is the infant mortality rate of municipality i at year t; imit-1 denotes the lagged infant mortality rate and it is included to capture persistence in infant mortality; PCPit is the density of primary care physicians; Xit contains a set of explanatory variables; λ is the time trend in the infant mortality rate of all municipalities; ηi is the unobserved municipality-specific effect; and εit denotes the error term.

The estimation of Eq (1) using pooled OLS will potentially lead to biased and inconsistent results as the orthogonality condition is violated [33]. This is a well-known problem of pooled OLS, as the correlation between ηi and imit−1 or Xit is very likely to be non-zero. The fixed effects model accounts for the municipality-specific effects but requires an assumption of strict exogeneity, which never holds in models with a lagged dependent variable [14].

We follow Arellano and Bond [34], who suggest an estimator using the Generalized Method of Moments that includes the lagged levels of the explanatory variables as instruments. To eliminate the fixed municipality specific effect, differenced-GMM estimates the Eq (1) in first differences (2) where Δimit = imitimit−1. Since Δimit−1 and Δεit are correlated by definition, the OLS estimate is inconsistent, hence the need for an instrumental variables approach. Blundell and Bond [32] pointed out that the GMM estimator obtained after first differencing is likely to perform poorly when the number of time periods is moderately small and the autoregressive parameter is moderately large. In such a case lagged levels of the series provide weak instruments for first differences. To cope with this problem, Arellano and Bover (1995) and Blundell and Bond [32] proposed a system of regressions in differences and levels, where the instruments used for the levels are lagged first-differences of the explanatory variables. System GMM can handle endogenous regressors, fixed individual effects and unobservable heteroskedasticity [35].

One concern in using System GMM is instrument proliferation. The number of instruments tends to easily grow with the number of time periods, making some asymptotic results and weakened specification tests. We follow the approach of Roodman [36], who suggests that the instrumental variables can be reduced by collapsing instruments and limiting lag depth amounts.

Results

Eq (1) is estimated with four different estimation methods (pooled OLS, random effects, fixed effects and System GMM). The regressions include a full set of covariates informed by the previous literature and control for time trends. All regressions use robust standard errors.

S2 Appendix reports the results of pooled OLS, random effects and fixed effects (within). Our estimates consistently show the lagged dependent variable (the infant mortality rate) is statistically significant, indicating a low degree of persistence in the infant mortality. Column (1) and (2) show the results of pooled OLS, column (3) and (4) report the results of random effects and column (5) and (6) the results of fixed effects. The coefficient on the number of PCPs is statistically significant using all estimation methods, except using fixed effects estimation after including time trends (column (6)).

However, as discussed previously and addressed by Nickell (1981), in dynamic panel data models, fixed effects estimation, as well as pooled OLS and random effects, may be biased and inconsistent. To deal with endogenous regressors, among other concerns, the Eq (1) is estimated using System GMM.

Table 3 reports the System GMM estimations. The table presents different specifications to test the sensitivity of the results to including additional controls. The Hansen/Sargan test of over-identifying restrictions show the null hypothesis cannot be rejected in the System GMM, consistent with the hypothesis that the instruments are valid for the estimation. In addition, we also test for second-order residual serial correlation and there is no evidence for significant AR(2) serial correlation of the residuals. To correct for downwards biased standard errors in finite samples, we employ a System GMM two-step estimator implementing the correction developed by Windmeijer [37]. In addition, we follow an empirical strategy to deal with instrument proliferation in the System GMM. Following the approach of Roodman [35] we collapsed the instruments and restricted the number of lags used in the model.

thumbnail
Table 3. Estimation results for infant mortality using System GMM.

https://doi.org/10.1371/journal.pone.0217614.t003

All specifications of the model find a negative and statistically significant coefficient on PCPs of a similar order of magnitude (ranging from -5.479 to -7.082). The first column shows the regression results including only PCP and time trends. Column (2) includes other FHS team, column (3) adds hospital beds and supplemental health insurance, column (4) contains female illiteracy and GDP per capita.

The columns 5–8 include household facilities separately. The trash collection, available electricity and piped water coefficients are negative and statistically significant. This indicates the negative association between household facilities and the infant mortality rate. The statistical insignificance of the electricity coefficient and the positive effect of sewage coefficient in column (9), which included all covariates, may be related to the presence of multicollinearity between the household variables.

In the most complete specification (column (9)), an increase of one primary care physician per 10,000 is associated with a 7.08 fewer infant deaths per 10,000 live births or 4.5% of the mean infant mortality rate.

In order to test the individual effect of the traditional and FHS models of primary care, the main specification, with all controls, was estimated considering PCPs in both models separately (S3 Appendix). The results show that only the coefficient on FHS physicians (-10.05) was statistically significant.

In addition, we test for a non-linear relationship between PCPs and infant mortality, where the effect of PCPs may at first increase and then decline with each additional PCP. The main model was re-estimated including a squared term and both variables were found to statistically significant with a U-shaped relationship (Table 4).

thumbnail
Table 4. Estimation results for a nonlinear specification.

https://doi.org/10.1371/journal.pone.0217614.t004

Based on these results (Table 4), we can estimate the approximate level of primary care physicians that minimizes infant mortality rates. According to our data the turning point is approximately at 6.7 PCPs per 10,000 population. The mean value across municipalities in 2012 was 3.4 PCPs per 10,000 population (S.D 1.93), suggesting that in the range of the estimation sample further increases in the number of PCPs could continue to reduce infant mortality but with diminishing returns to the supply of PCPs. The predicted relationship between the density of primary care physicians and infant mortality is shown in Fig 1.

thumbnail
Fig 1. Non-linear relationship between infant mortality rate and density of primary care physicians.

https://doi.org/10.1371/journal.pone.0217614.g001

A placebo test was also conducted using the mortality rate for all ages from external causes (ICD-10 codess V00-Y99) as an alternative dependent variable. External causes include mortality by transport accident, falls, assault, drowning, among others and hence, the number of primary care physicians should have no impact as there is no plausible mechanism through which PCPs could influence this measure of mortality. The coefficient on PCPs is found to be near zero and statistically insignificant (S4 Appendix). This result supports the validity of our empirical strategy.

Finally, since the results are based on all Brazilian municipalities, we performed a sensitivity analysis. Table 5 reports the results for the main model (column (9) in Table 3) using alternative estimation samples. Besides the full sample, Eq (1) was estimated using two different estimation samples: Frias et al. [28] criteria in a balanced panel (column A) and Frias et al. [28] criteria in an unbalanced panel (column B). The coefficient of primary care physicians was found to be statistically significant, but to have a larger order of magnitude in both samples. The other FHS team members coefficient also had a larger magnitude in these sub-samples when compared with the full sample. A larger magnitude of the coefficients could be explained by inadequate records of deaths and births over the period in the full sample; this follows from the method of selecting the restricted sample which involves excluding municipalities with a high standard deviation in mortality rates and municipalities with very low reported mortality rates. Although with a variation in magnitude, these results confirm that the negative effect of primary care physicians on infant mortality is found in all samples.

thumbnail
Table 5. Estimation results for sensitivity analyses on alternative estimation samples.

https://doi.org/10.1371/journal.pone.0217614.t005

Discussion

This study assesses the relationship between primary care physicians and infant mortality at the municipal level in Brazil using panel data from 2005 to 2012. Our findings support the results of previous studies which have emphasized the role of primary care in improving health in developing countries while addressing concerns that much of the previous literature suffer, such as their weak methodological approach, including inadequate covariates, and the lack of longitudinal sources which allows to evaluate the effect of primary care on outcomes over time [38]. In this study, we control for a wide range of covariates and use a robust dynamic panel model to account for the endogeneity problem.

Our results show a negative effect of the density of primary care physicians on infant mortality in Brazil. In particular, our results are consistent with two similar previous studies from Brazil but using data from an earlier time period, though we find larger effects. Sousa et al. [9] used data between 1991 and 2000, of 4,267 Minimum Comparable Areas (MCA), and found that an increase of one physician (including non-PCPs) per 1,000 population reduces neonatal mortality by 2.3%. Macinko et al. [39], using microregional data from 1999 to 2004, found that a 10% increase in population coverage of the FHS reduced infant mortality by 0.45%. Our study finds an increase of one PCP per 10,000 reduces infant mortality by 7.08 per 10,000 live births, or 4.5% of the mean infant mortality rate. Extrapolating our results, an increase of 10% in primary care physicians (0.37 PCPs per 10,000 population) predicts a reduction in infant deaths by 765 per year. Sensitivity analysis, which includes only those municipalities with high quality of information [28], shows that this effect can be even stronger.

In addition, two other studies, using different methods, also support our results. Rocha & Soares [40], using a difference-in-difference estimator and data from 1994 to 2004, found that municipalities which implemented the FHS for three years were associated with a 6.5% reduction in IMR and 20% reduction in municipalities with 8 year of FHS implantation. Aquino et al. [27], using panel data with a negative binomial response, reported a robust association between FHS coverage and infant mortality reduction from 1996 to 2004.

Hospital beds and other FHS team members are not found to be statistically significant in explaining infant mortality in the preferred specification. A high correlation was noted between PCPs and other FHS team members. However, the PCP coefficient was -6.40 after omitting the other FHS team members variable in column (9), showing no influence on its statistical significance and small influence on its magnitude. Previous studies for Brazil that use hospital beds and nurses density as independent variables to account for infant mortality have not reached a consensus. Macinko et al. [23] found both variables negatively associated with infant mortality, but only hospital beds was statistically significant. Sousa et al. [9] showed nurse density had an adverse and significant effect on infant mortality. By contrast, Macinko et al. [39] recorded that the hospital beds coefficient is positive and significant, but the study did not includes nurses.

The results show that municipalities with higher GDP per capita were associated with lower infant mortality rates. The association between income and health outcomes is well established in the literature. This relationship is found using different measurements of income, including GPD per capita [41]. Preston [42] argues that national income is probably the best indicator of living standards in a region, as well as representing a large range of factors that can influence mortality.

Kassouf [43] found that, in the same way as income, supplemental health insurance coverage is associated with better health status. Although The Brazilian health system—Unified Health System (SUS)—has universal access, around 24.4% of population was covered by supplemental health insurance in 2012 [44]. There are great inequalities in the use of health services in Brazil. Cambota & Rocha [45] showed that income, education and supplemental health insurance coverage contribute to increase the inequalities in these services. There are evidence that individuals who have supplemental health insurance visit a physician on average around 25% more than individual with no supplemental health insurance coverage; for hospital admissions the differential ranges from 8% to 15% [46]. We found a negative effect of supplemental health insurance coverage on infant mortality.

Our results are robust to including time trends and a range of important determinants of infant mortality rate such household facilities. Trash collection, piped water and available electricity are found relevant in accounting for changes in infant mortality rate in developing countries [4751]. Electricity access facilitates more hygienic preparation of food and sterilization [52]. Piped water can cause a significant reduction in the infant mortality rate mainly by reducing the incidence of diarrhea [53]. Diarrhea is more prevalent in the developing countries due the lack of safe drinking water, sanitation and hygiene and it was responsible for 1.5 deaths among children under-5 worldwide in 2004 [54].

In additional tests, we included the Bolsa Familia Program (BFP) as a control variable. The literature has indicated the relevance of this conditional cash transfer program for infant/child mortality [55,56]. Our results remained consistent and showed small variations in the magnitude of the effect of PCPs after controlling for the BFP. However, the size of the BFP coefficient was sensitive to the control of time trends.

The effect of FHS PCPs was tested separately from the ‘traditional’ PCPs (S3 Appendix). The results showed that only the FHS PCPs was associated with lower infant mortality. The result confirms the relevance of the FHS for infant health, as reported in previous studies [23,27,39,40]. The coefficient of traditional PCPs (general practitioners, pediatricians and obstetrician-gynecologists) was not statistically significant, maybe because traditional physicians have increasingly been replaced by the FHS in recent years across Brazilians municipalities. The reduction of the traditional model of primary care, and expansion of FHS, has been based on a national strategy for reorganization of primary care [20].

In addition to the density and coverage of primary health care team, the literature has also indicated the relevance of consolidating primary care. There is evidence that municipalities with a consolidated FHS recorded a further reduction in infant mortality rates than municipalities without FHS coverage, incipient coverage and intermediate coverage [27]. This finding was also reported by Rasella, et al. [57] for neonatal, postneonatal and children mortality.

Finally, our extensive robustness checks including alternative model specifications and a placebo regression strengthen the validity of our findings. The study’s limitations include the quality of mortality data, which vary across region of Brazil. However, overall these additional estimations improves the results’ validity since we find similar results.

The above discussion highlighted that studies using data from earlier time periods produced results of a different magnitude but a similar direction of effect. Since 2012, when our data ended, there have been significant changes in both the economy and in PCP policy that could potentially lead to different effects than we have observed here. First, there was a significant economic downturn with a significant recession in 2015 and 2016. If this differentially affected municipalities, then the effect of PCP supply on infant mortality could be different post-2012. Second, alongside the economic crisis, fiscal austerity policies have recently been implemented in Brazil. These measures will reduce public expenditure as a percentage of GDP over 20 years, substantially affecting the healthcare system [58]. Studies have pointed out that the major impacts of this policy may be felt by the poorest and the most vulnerable [59]. In the long-term, child morbidity and mortality may also be considerably affected [60].

The Brazilian government introduced the More Doctors Program in 2013 and although evaluations of this program have reported a positive effect on some health outcomes [61], they found little evidence of an effect on infant mortality [62]. This could either be due to diminishing returns and/or the effects of the economic downturn on other factors that influence infant mortality. This highlights that the effectiveness of PCPs may change over time and may vary with economic cycles.

Conclusions

Brazilian health policy has placed substantial importance on primary care in past decades, potentially providing lessons for other developing countries. The objective of this study has been to investigate infant mortality and the level of primary care physicians, including a range of important determinants of infant health. The results indicate a negative association between primary care physicians and infant mortality rates in Brazil. Moreover, a reduction of illiteracy among woman and a rise in the private health insurance coverage, GDP per capita, percentage of household with trash collection and piped water has also contributed to the decline of infant mortality.

Primary health care is regarded as a key component for building a strong healthcare system that ensures positive health outcomes, quality in healthcare and universal access. Brazil’s experience demonstrates that increases of primary care physician supply represents one strategy to reduce infant deaths and remains a potential instrument to achieve international goals related to child health.

Supporting information

S2 Appendix. Estimation results for infant mortality using different estimation techniques.

https://doi.org/10.1371/journal.pone.0217614.s002

(DOCX)

S3 Appendix. Estimation results for traditional and FHS physicians.

https://doi.org/10.1371/journal.pone.0217614.s003

(DOCX)

S4 Appendix. Estimation results for external causes mortality.

https://doi.org/10.1371/journal.pone.0217614.s004

(DOCX)

Acknowledgments

This work was conducted at the University of Melbourne during a scholarship supported by the Coordination of Improvement of Higher Level Personnel (CAPES), Ministry of Education of Brazil.

References

  1. 1. Macinko J, Starfield B, Shi L. The contribution of primary care systems to health outcomes within Organization for Economic Cooperation and Development (OECD) countries, 1970–1998. Health Serv Res. 2003 Jun;38(3):831–865. pmid:12822915
  2. 2. WHO Declaration of Alma-Ata: International Conference on Primary Health Care, Alma-Ata. WHO—World Health Organization. 1978. 6 p.
  3. 3. Basu S, Berkowitz SA, Phillips RL, Bitton A, Landon BE, Phillips RS. Association of Primary Care Physician Supply With Population Mortality in the United States, 2005–2015. JAMA Intern Med. 2019 Feb 18.
  4. 4. Aakvik A, Holmas TH. Access to primary health care and health outcomes: the relationships between GP characteristics and mortality rates. J Health Econ. 2006 Nov;25(6):1139–1153. pmid:16675052
  5. 5. Shi L, Macinko J, Starfield B, Xu J, Regan J, et al. Primary care, infant mortality, and low birth weight in the states of the USA. J Epidemiol Community Health. 2004 May;58(5):374–380. pmid:15082734
  6. 6. MacDorman MF, Mathews TJ. The challenge of infant mortality: have we reached a plateau? Public Health Rep. 2009 Sep-Oct;124(5):670–681. pmid:19753945
  7. 7. Nascimento RM, Leite AJ, Almeida NM, Almeida PC, Silva CF. [Determinants of neonatal mortality: a case-control study in Fortaleza, Ceara State, Brazil]. Cad Saude Publica. 2012 Mar;28(3):559–572. pmid:22415188
  8. 8. Ribeiro AM, Guimaraes MJ, Lima Mde C, Sarinho SW, Coutinho SB. [Risk factors for neonatal mortality among children with low birth weight]. Rev Saude Publica. 2009 Apr;43(2):246–255. pmid:19225697
  9. 9. Sousa A, Dal Poz MR, Boschi-Pinto C. Reducing inequities in neonatal mortality through adequate supply of health workers: evidence from newborn health in Brazil. PLoS One. 2013;8(9):e74772. pmid:24073222
  10. 10. WHO The Millennium Development Goals Report 2015. New York: WHO—World Health Organization.2015. pp.
  11. 11. WHO WHO Data. WHO—World Health Organization.2017. pp.
  12. 12. Farahani M, Subramanian SV, Canning D. The effect of changes in health sector resources on infant mortality in the short-run and the long-run: a longitudinal econometric analysis. Soc Sci Med. 2009 Jun;68(11):1918–1925. pmid:19362762
  13. 13. Mishra P, Newhouse D. Does health aid matter? J Health Econ. 2009 Jul;28(4):855–872. pmid:19616864
  14. 14. Wooldridge JrM Econometric Analysis of Cross Section and Panel Data. Cambridge, Massachusetts: The MIT Press.2001.
  15. 15. Macinko J, de Oliveira VB, Turci MA, Guanais FC, Bonolo PF, et al. The influence of primary care and hospital supply on ambulatory care-sensitive hospitalizations among adults in Brazil, 1999–2007. Am J Public Health. 2011 Oct;101(10):1963–1970. pmid:21330584
  16. 16. Macinko J, Dourado I, Aquino R, Bonolo Pde F, Lima-Costa MF, et al. Major expansion of primary care in Brazil linked to decline in unnecessary hospitalization. Health Aff (Millwood). 2010 Dec;29(12):2149–2160.
  17. 17. Arulampalam W, Bhalotra S. Persistence in infant mortality: Evidence for the Indian states. Discussion Paper 2488—IZA. 2006.
  18. 18. Menicucci TM. [The Unified National Health System, 20 years: assessment and perspectives]. Cad Saude Publica. 2009 Jul;25(7):1620–1625. pmid:19578584
  19. 19. Starfield B. Is primary care essential? Lancet. 1994 Oct 22;344(8930):1129–1133. pmid:7934497
  20. 20. Escorel S, Giovanella L, Magalhaes de Mendonca MH, de Castro Maia Senna M. [The Family Health Program and the construction of a new model for primary care in Brazil]. Rev Panam Salud Publica. 2007 Feb-Mar;21(2–3):164–176. pmid:17565803
  21. 21. Macinko J, Harris MJ. Brazil’s Family Health Strategy. N Engl J Med. 2015 Sep 24;373(13):1278. pmid:26398087
  22. 22. Starfield B, Shi L, Macinko J. Contribution of primary care to health systems and health. Milbank Q. 2005;83(3):457–502. pmid:16202000
  23. 23. Macinko J, Guanais FC, de Fatima M, de Souza M. Evaluation of the impact of the Family Health Program on infant mortality in Brazil, 1990–2002. J Epidemiol Community Health. 2006 Jan;60(1):13–19. pmid:16361449
  24. 24. Schell CO, Reilly M, Rosling H, Peterson S, Ekstrom AM. Socioeconomic determinants of infant mortality: a worldwide study of 152 low-, middle-, and high-income countries. Scand J Public Health. 2007;35(3):288–297. pmid:17530551
  25. 25. Santos ÁdS, Miranda SMRCD A enfermagem na gestão em atenção primária à saúde. Barueri, SP: Editora Manole 2007.
  26. 26. Saparolli ECL, Adami NP. Assessment of the quality of pediatric nursing consultation in Family Health Program. Acta Paulista de Enfermagem. 2007;20(1):7.
  27. 27. Aquino R, de Oliveira NF, Barreto ML. Impact of the family health program on infant mortality in Brazilian municipalities. Am J Public Health. 2009 Jan;99(1):87–93. pmid:19008516
  28. 28. Frias PG, Szwarcwald CL, Lira PI. [Evaluation of information systems on live births and mortality in Brazil in the 2000s]. Cad Saude Publica. 2014 Oct;30(10):2068–2280. pmid:25388310
  29. 29. Rasella D, Harhay MO, Pamponet ML, Aquino R, Barreto ML. Impact of primary health care on mortality from heart and cerebrovascular diseases in Brazil: a nationwide analysis of longitudinal data. BMJ. 2014 Jul 03;349:g4014. pmid:24994807
  30. 30. Andrade CL, Szwarcwald CL. [Socio-spatial inequalities in the adequacy of Ministry of Health data on births and deaths at the municipal level in Brazil, 2000–2002]. Cad Saude Publica. 2007 May;23(5):1207–1216. pmid:17486242
  31. 31. Arellano M, Bover O. Another look at the instrumental variables estimation of error-components models. Journal of Econometrics. 1995;68(1):30.
  32. 32. Blundell R, Bond S. Initial conditions and moment restrictions in dynamic panel data models. Journal of Econometrics. 1998;87(1):28.
  33. 33. Nickell S. Biases in dynamic models with fixed effects. Econometrica. 1981;49(6):9.
  34. 34. Arellano M, Bond S. Some Tests of Specification for Panel Data: Monte Carlo Evidence and an Application to Employment Equations. The Review of Economic Studies. 1991;58(2):20.
  35. 35. Roodman D. How to do xtabond2: An introduction to difference and system GMM in Stata. The Stata Journal 2009;9(1):50.
  36. 36. Roodman D. Practitioners’ Corner: A Note on the Theme of Too Many Instruments. Oxford Bulletin of Economics and Statistics. 2009;71(1):23.
  37. 37. Windmeijer F. A finite sample correction for the variance of linear efficient two-step GMM estimators. Journal of Econometrics. 2005;126(1):26.
  38. 38. Macinko J, Starfield B, Erinosho T. The impact of primary healthcare on population health in low- and middle-income countries. J Ambul Care Manage. 2009 Apr-Jun;32(2):150–171. pmid:19305227
  39. 39. Macinko J, Mde Marinho de Souza F, Guanais FC, da Silva Simoes CC. Going to scale with community-based primary care: an analysis of the family health program and infant mortality in Brazil, 1999–2004. Soc Sci Med. 2007 Nov;65(10):2070–2080. pmid:17689847
  40. 40. Rocha R, Soares RR. Evaluating the impact of community-based health interventions: evidence from Brazil’s Family Health Program. Health Econ. 2010 Sep;19 Suppl:126–158.
  41. 41. O’Hare B, Makuta I, Chiwaula L, Bar-Zeev N. Income and child mortality in developing countries: a systematic review and meta-analysis. J R Soc Med. 2013 Oct;106(10):408–414. pmid:23824332
  42. 42. Preston SH. The Changing Relation between Mortality and level of Economic Development. Population Studies 1975;29(2).
  43. 43. Kassouf AL. Acesso aos serviços de saúde nas áreas urbana e rural do Brasil. Rev Econ Sociol Rural. 2005;43(1).
  44. 44. ANS. National Regulatory Agency for Private Health Insurance and Plans 2018.
  45. 45. Cambota JN, Rocha FF. Determinantes das desigualdades na utilização de serviços de saúde: análise para o Brasil e regiões. Pesquisa e Planejamento Econômico. 2015;45:24.
  46. 46. Andrade MV, Maia AC. Diferenciais de Utilização do Cuidado de Saúde no Sistema Suplementar Brasileiro. Est econ. 2009;39(1):31.
  47. 47. Ezeh OK, Agho KE, Dibley MJ, Hall J, Page AN. The impact of water and sanitation on childhood mortality in Nigeria: evidence from demographic and health surveys, 2003–2013. Int J Environ Res Public Health. 2014 Sep 05;11(9):9256–9272. pmid:25198687
  48. 48. Fink G, Gunther I, Hill K. The effect of water and sanitation on child health: evidence from the demographic and health surveys 1986–2007. Int J Epidemiol. 2011 Oct;40(5):1196–1204. pmid:21724576
  49. 49. Barros FC, Matijasevich A, Requejo JH, Giugliani E, Maranhao AG, et al. Recent trends in maternal, newborn, and child health in Brazil: progress toward Millennium Development Goals 4 and 5. Am J Public Health. 2010 Oct;100(10):1877–1889. pmid:20724669
  50. 50. Guimarães MJ, Marques NM, Melo Filho DA, Szwarcwald C. [Living conditions and infant mortality: intra-urban differentials in Recife, Pernambuco State, Brazil]. Cad Saude Publica. 2003 Sep-Oct;19(5):1413–1424. pmid:14666223
  51. 51. Volpe FM, Abrantes MM, Capanema FD, Chaves JG. The impact of changing health indicators on infant mortality rates in Brazil, 2000 and 2005. Rev Panam Salud Publica. 2009 Dec;26(6):478–484. pmid:20107701
  52. 52. Van de Poel E, O’Donnell O, Van Doorslaer E. What explains the rural-urban gap in infant mortality: household or community characteristics? Demography. 2009 Nov;46(4):827–850. pmid:20084831
  53. 53. Gamper-Rabindran S, Khan S, Timmins C. The impact of piped water provision on infant mortality in Brazil: A quantile panel data approach. Journal of Development Economics. 2010;92(2):12.
  54. 54. UNICEF W. Diarrhoea: why children are still dying and what can be done. 2009.
  55. 55. Guanais FC. The combined effects of the expansion of primary health care and conditional cash transfers on infant mortality in Brazil, 1998–2010. Am J Public Health. 2013 Nov;103(11):2000–2006. pmid:24028257
  56. 56. Rasella D, Aquino R, Santos CA, Paes-Sousa R, Barreto ML. Effect of a conditional cash transfer programme on childhood mortality: a nationwide analysis of Brazilian municipalities. Lancet. 2013 Jul 6;382(9886):57–64. pmid:23683599
  57. 57. Rasella D, Aquino R, Barreto ML. Reducing childhood mortality from diarrhea and lower respiratory tract infections in Brazil. Pediatrics. 2010 Sep;126(3):e534–540. pmid:20679307
  58. 58. Vieira FS, Benevides RPdSe. Os impactos do novo regime fiscal para o financiamento do Sistema Único de Saúde e para a efetivação do direito à saúde no Brasil. IPEA—Instituto de Pesquisa Econômica Aplicada, Nota ténica 28. 2016.
  59. 59. Massuda A, Hone T, Leles FAG, de Castro MC, Atun R. The Brazilian health system at crossroads: progress, crisis and resilience. BMJ Glob Health. 2018;3(4):e000829. pmid:29997906
  60. 60. Rasella D, Basu S, Hone T, Paes-Sousa R, Ocke-Reis CO, Millett C. Child morbidity and mortality associated with alternative policy responses to the economic crisis in Brazil: A nationwide microsimulation study. PLoS Med. 2018 May;15(5):e1002570. pmid:29787574
  61. 61. Fontes LFC, Conceicao OC, Jacinto PA. Evaluating the impact of physicians’ provision on primary healthcare: Evidence from Brazil’s More Doctors Program. Health Econ. 2018 Aug;27(8):1284–1299. pmid:29770534
  62. 62. Feres BCJ. More Doctors, Better Health? Evidence from a Physician Distribution Policy. HEDG, WP 17/29 University of York. 2017.