Abstract

This study explores whether hospitals with higher increases in obesity levels have higher cesarean section (CS) rates and the consequential effects on maternal and newborn health in Mexico for 2008–2015. It models how changes in the obesity level of hospitals’ patient pools may affect the quantity and quality of care by focusing on the use of CS and the potential returns to specialization. And it creates a measure of hospital-level obesity, based on the fraction of obesity-related discharges for women of childbearing age. Exploiting temporal and hospital variation of this measure, results show that higher hospital-level obesity increases a woman’s probability of having a CS. Also, delivery-related birth outcomes improve: maternal mortality, birth injuries, and birth trauma decrease. The evidence is consistent with hospital-level specialization in CS leading to better birth outcomes.

1. Introduction

Cesarean section (CS) procedures increased worldwide from 7 percent to 21 percent between 1990 and 2018, with the largest increases in low- and middle-income countries (Betran et al. 2021). These rates are beyond the 10–15 percent rates recommended by the World Health Organization (WHO), as CS may have adverse effects on maternal and child health (Jachetta 2014; Costa-Ramón et al. 2018). Parallel to CS increasing trends, the prevalence of being overweight and obesity—a common risk factor for CS and adverse pregnancy outcomes (Kominiarek et al. 2010; Fyfe et al. 2011)—has also increased worldwide, becoming a major health issue in low- and middle-income countries (Popkin, Corvalan, and Grummer-Strawn 2020). Despite these trends, the literature on the determinants and consequences of CS rates pertains to high-income (Currie and MacLeod 2008; Kozhimannil, Arcaya, and Subramanian 2014; Costa-Ramón et al. 2018; Card, Fenizia, and Silver 2023) rather than developing countries, where the inefficiency and inequality of the health systems are more salient, and where maternal and infant mortality are high and remain targets of the Sustainable Development Goals (SDGs). As CS are surgical procedures that can save the lives of mothers and children, which factors play a role in maximizing their benefits in contexts where health systems are resource constrained is a pressing question.

Empirical evidence—mostly from the United States—has shown that substantial variation of CS rates across hospitals may be uncorrelated with medical needs, and that supply rather than demand factors might explain unnecessary CS procedures (Kozhimannil et al. 2014). On the one side, while there is some evidence of benefits for riskier patients (Doyle et al. 2015), the excessive use of unnecessary intensive treatments may be wasteful in terms of resources. On the other side, if underutilized, increasing their use efficiently may lead to important gains for population health. Thus, understanding the observed differences in the use of intensive treatments across hospitals which may have a limited impact on patients’ health is a key policy concern. This issue is especially important in developing countries with increasing CS rates and scarce resources. While in Sub-Saharan Africa, an average CS rate of 5 percent might indicate under utilization and unmet needs for CS, in Latin America and the Caribbean a 42 percent CS rate suggests overuse leading to the health consequences of unnecessary procedures, and waste of human and financial resources (Betran et al. 2021).

Recent evidence from the United States shows that the large variations in the use of CS across health facilities may be explained by hospital delivery practices and that CS hospital specialization may affect maternal and health outcomes (Card, Fenizia, and Silver 2023). However, even in high-income countries, the relationship between hospital-level patient characteristics that may increase medically indicated CS (e.g. obesity prevalence among mothers) and CS specialization is unclear. Similarly, the effects of hospital specialization on health outcomes of different types of patients are ambiguous.

This paper aims to fill these two gaps: the relative lack of evidence of the determinants and consequences of CS in low- and middle-income settings and the effects of changes in the patient pool on procedure use. It studies whether hospitals with a higher fraction of obesity-related discharges have higher CS rates and are more likely to become specialized in performing CS. In addition, it analyzes the effects of specialization on newborn and maternal health outcomes at birth. In order to help predict the effects of increases in obesity and hospital specialization, it develops a model of CS choice, specialization, and health outcomes based on Bozzoli, Deaton, and Quintana-Domeque (2009). CS choice is modeled as a function of the health (obesity) burden of a hospital’s population. A useful feature of the model is that its predictions are presented in terms of the mean hospital-level health (i.e. they are independent of an individual mother’s specific weight). Thus, it is straightforward to link the empirical findings to the model’s predictions and to interpret the results in terms of specialization.

The study focuses on Mexico, where CS rates doubled from 23 percent in the 1990s to 46 percent in 2016 and where at least 72 percent of the adult population is overweight or obese (OWOB) (Brenes-Monge et al. 2019). As supplementary online appendix fig. S1.1 shows, there is a positive correlation between these variables.1 It is well documented that obesity is a risk factor for CS (Kominiarek et al. 2010; Fyfe et al. 2011; Cnattingius et al. 2013). In Mexico, mother’s obesity is explicitly listed as a risk factor during pregnancy in terms of risk of preterm delivery and is highlighted as a condition under which VBAC (vaginal birth after CS) is not recommended by the clinical guidelines of the Mexican Health Ministry (IMSS 2013, IMSS 2017). Furthermore, anecdotal evidence from medical personnel suggests that higher CS rates are concentrated in larger hospitals that specialize in riskier pregnancies related to factors such as a high pregestational BMI and hypertension.2 Hence, obesity could be a potential factor in explaining the high CS rate increase in Mexico. In addition, studying the consequences of CS on maternal and neonatal health is particularly relevant in the Mexican case. Although the maternal mortality rate decreased substantially from 88 maternal deaths per 100,000 live births in 1990 to 38 in 2016,3 this rate is far above Mexico’s SDGs target of 22. Moreover, women who are poor, less educated, and indigenous are at higher risk of maternal mortality (Rodríguez-Aguilar 2018).

Time Variation in Hospital-Level CS and Obesity Rates
Figure 1.

Time Variation in Hospital-Level CS and Obesity Rates

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: Obesity index (left) and C-Section (CS) rate (right) averaged within group according to the hospital type specified in the legend, using year 2008 as the baseline year. Hospitals were classified as low/medium/high variable change according to their observed individual linear time trend for each variable, and sorted into groups of size according to it: low change (25 percent), medium (50 percent), high (25 percent). A hospital may be in one category for obesity changes and in another when it comes to CS changes (e.g. high obesity change but medium CS change).

Geographic Variation in CS and Obesity Time Trends
Figure 2.

Geographic Variation in CS and Obesity Time Trends

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: Hospitals’ linear trend for obesity index (left) or C-Section (CS) rate (right) over the study period (2008–2015). These are “Hexbin” maps which use a density 2d technique plotted on top of a map. Each figure shows the distribution of a variable (obesity and C-section) on a map, splitting the map into a set of hexagons. The bandwidth chosen in the figures is 0.6, and only the hexagons with a positive count of points (hospitals) are shown on the map.

The analysis merges the panel of all public hospitals with the census of birth records for 2008–2015, which contains information on birth-type delivery and several birth outcomes, as well as maternal characteristics. While the administrative health data contain detailed information on inpatient diagnoses, it lacks information on whether patients themselves are obese. To overcome this data limitation, this study constructs a novel obesity measure at the hospital level, using patient diagnoses to identify all women of reproductive age admitted at a health facility who were diagnosed with a condition highly related to obesity based on the medical literature (Kotchen 2010; Cnattingius et al. 2013; Kearns et al. 2014). The hospital-level obesity measure is strongly correlated with state-level female OWOB and obesity prevalence at the state level. This is an alternative measure of obesity levels that could be used in other settings when anthropometric data is unavailable.

The empirical strategy uses variation in the health-facility-level obesity index—both across hospitals and over time—to identify the effects of obesity prevalence at the hospital level on CS and on delivery-related health outcomes for newborns and mothers. The sample includes all mothers with a low risk of having a CS. The hospital-level obesity index is the main explanatory variable, conditioning on mother, pregnancy, and hospital-level characteristics, as well as time and hospital fixed effects. The main threats to identification come from endogenous selection and differential changes in patient composition at high- versus low-obesity hospitals. To alleviate these concerns, the analysis considers municipalities with only one clinic, as this limits patient choice in terms of place of delivery. In addition, as shown, the obesity index does not predict predetermined maternal characteristics such as age, marital status, or education. Moreover, the results are robust to alternative specifications of the hospital-level obesity index and fixed-effects models.

The study has two main findings. First, the individual probability of delivering via CS increases as the proportion of female obesity-related discharges in public hospitals increases. Overall, the positive effect of the hospital-level obesity measure on an individual woman’s likelihood of CS is consistent with a hospital-level CS specialization argument. If higher-obesity hospitals perform more CS due to medical reasons—as obesity is a known risk factor for CS—there may be productivity spillovers from knowledge, experience, and learning by physicians (as in Chandra and Staiger (2007)) that lead to specialization in CS.

Alternative mechanisms, such as women’s preferences over the type of birth delivery, do not seem to drive the results: the effects across mothers’ education levels are similar. Furthermore, doctors’ medical or monetary incentives or changes in hospital congestion are also unlikely to explain the results.

Second, increases in the hospital-level obesity index improve birth delivery-related outcomes: there is a significant reduction in maternal mortality and newborn birth injury and trauma for all mothers with low risk of CS. It is worth noting that the hospital-obesity measure does not affect delivery-unrelated outcomes, such as low-birth weight (LBW), as these should not be affected by CS practices. Moreover, the hospital-obesity measure increases mother survival and decreases the probability of birth injury and trauma for preterm newborns, who are almost always delivered via CS and are thus unaffected by the potential effects of obesity on delivery-type choice. Overall, this evidence supports the hospital-level CS specialization hypothesis.

This paper contributes to the emerging literature studying the determinants and consequences of increasing CS rates in low- and middle-income countries.4 Despite the high CS prevalence (Boerma et al. 2018), the evidence for Mexico is limited.5 This study is closest to Card, Fenizia, and Silver (2023) who show, for the United States, that proximity to hospitals with high CS rates leads to more cesarean deliveries and higher average Apgar scores among low-risk first births. Nevertheless, the evidence on the effects of CS on child health outcomes is not clear-cut. While Jensen and Wust (2015) show that breech births benefit from CS deliveries, Jachetta (2014) finds that CS leads to higher incidence of asthma and Costa-Ramón et al. (2018) show that non-medically indicated CS delivery leads to lower Apgar scores. This paper highlights that CS specialization can reduce maternal mortality. This is a relevant finding for low- and middle-income countries with high maternal mortality rates.

Finally, this study relates to the literature on the effects of medical practices on health outcomes (Chandra and Staiger 2007; Doyle et al. 2015; Card, Fenizia, and Silver 2023). The results extend beyond the effects of specialization on delivery mode and are applicable to less developed contexts. They also highlight how mothers’ obesity prevalence may affect treatment choice in a context where the health system is more resource constrained (Lagarde and Blaauw 2022).

2. Model

This section presents a simple model of procedure (CS) choice, specialization, and health outcomes. This paper builds on Bozzoli, Deaton, and Quintana-Domeque (2009) and adapts it to model CS choice as a function of the health burden of the population at a given hospital. It also proposes an extension where doctors may become specialized.

Suppose that each pregnancy carries a certain risk which is negatively correlated with the underlying health of the newborn, denoted by hi. Assume that health is distributed in the population with distribution function F(h) and considers populations to be defined at the hospital level. The decision to perform a CS depends on the pregnancy risk and a CS will be performed when pregnancy risk is higher than a given cutoff, i.e. when a newborn’s health is below a cutoff z. Hence, the CS rate is given by F(z).

Then consider an additional “health burden,” which reduces newborns’ health endowment and increases pregnancy risk. This additional burden is denoted by vt and, in the context of this study, is captured by obesity. The value of vt can vary year to year. Larger values of vt indicate riskier pregnancies. This increase in pregnancy risk leads to an increase in the CS rate, where a delivery will occur via CS if the health of the newborn “net” of the obesity burden is less than z, or if hivtz and the CS rate, given the burden of obesity, vt, is given by |$\%\mathrm{CS}_t = F(z+v_t)$|⁠. It is worth noting that z may vary either over time or across hospitals, as specialization, physician skills, or resources available may lead to different thresholds for using CS. These scenarios are considered in the next subsection.

To include the potential “stunting” (or “scarring” as in (Bozzoli, Deaton, and Quintana-Domeque 2009)) effect of obesity on newborns’ health, assume that a fraction 0 ≤ θ ≤ 1 of vt will be subtracted from a newborn’s health stock, hi. A newborn’s (observed) health will therefore be given by |$\tilde{h}_{it}=h_i-\theta v_t$|⁠.

Given the effective cutoff at which CS are performed, the average health of newborns born by procedure T ∈ {V, CS}, where V stands for vaginal and CS stands for C-section, is given by

(1)
(2)

The first term for both equations reflects how average health for a given procedure changes as the threshold varies, which may be interpreted as a selection effect. This pure selection effect could be observed if θ = 0, i.e. if no stunting was associated with obesity. It is worth noting that for θ = 0,

for both T ∈ {V, C}. Intuitively, as the threshold for performing a CS moves to the right, marginally healthier babies are delivered via CS (left of the threshold)—increasing the average health of those born via CS—while the least healthy babies are “removed” from those born via vaginal delivery (right of the threshold)—thus also increasing the average health of those born vaginally. This is referred to as a pure selection effect.

Given that it is likely that θ > 0, the effect of vt on average newborn health can be rewritten by differentiating equations (1) and (2), respectively:

(3)
(4)

For both equations (3) and (4), the first term on the right-hand side is always positive (at the limit, for CS the healthiest individual has hi = z + vt − θvt, while for vaginal births the least healthy individual has hi = z + vt − θvt). The final sign of the equation may be either positive or negative, as it can change sign over the range of vt. In this context, it implies that an increase in obesity may lead to more or less healthy newborns, and may actually have a non-monotonic relationship depending on the parameters (θ). The overall health of newborns is given by

Note that while the average health of both types of births may be increasing or decreasing with vt, the overall health is always decreasing in vt, as long as θ > 0.

2.1. Specialization

Specialization in CS (intensive treatment) is introduced in the model via an additional term, st, that modifies the threshold under which CS are performed (effectively moving it to the right). In particular, a delivery will be a CS if hivtz + st and the CS rate, given the burden of obesity, vt, and the specialization, st, is now given by |$\%\mathrm{CS}_t = F(z+v_t+s_t)$|⁠.

Consider also that specialization in CS will be associated with higher physician skill in performing the intensive treatment and potentially in decreased skill in the competing treatment (as in Chandra and Staiger (2007), Currie and MacLeod (2008), Currie and MacLeod (2017)). This feature is introduced by adding the term δst for CS deliveries and subtracting the term γst for vaginal deliveries to health outcomes that are delivery related. For outcomes that are delivery unrelated, such as birth weight, assume that δ = γ = 0. The value of st depends linearly on the obesity burden faced by the hospital, vt, so that st(vt) = πvt, with π ≥ 0, as doctors treating a pool of riskier pregnancies are more likely to become specialized. With this setup, the mean (delivery-related) health of newborns born via procedure T, |$\bar{h}^{T,s}_t(v_t,s(v_t))$|⁠, is

Note that for a given |$\bar{v}$|⁠, the mean health of newborns will be larger when specialization is allowed, as |$\bar{h}_t^{\mathrm{CS},s}(\bar{v},s(\bar{v}))\gt \bar{h}_t^{\mathrm{CS}}(\bar{v})$|⁠. However, the relationship between |$\bar{h}_t^{V,s}(\bar{v},s_t)$| and |$\bar{h}_t^{V}(\bar{v})$| will depend on the specific parameters. More formally, for a given |$\bar{v}$|⁠, it is possible to calculate how average health (overall and by type of delivery) varies as the degree of obesity-driven specialization, captured by the parameter π, increases (note that the case of no specialization is captured by π = 0). In particular,

(5)
(6)

As in equations (3) and (4), the first terms on the right-hand sides of equations (5) and (6) will always be positive, as the mean health of those born via CS/vaginally will always be lower/higher than that of the individual at the threshold, given by hi = z + (1 + π)vt − (θ − πδ)vt [hi = z + (1 + π)vt − (θ + πγ)vt]. Given that in equation (5) the second term is also always positive, for a given level of obesity (⁠|$\bar{v}$|⁠), more specialization unambiguously increases the average health of those born through CS. This is not the case for vaginal births, where the last term in equation (6) is negative and could thus lead to worse health outcomes.

Additionally, when considering the health of the whole population, independent of the procedure, increases in obesity without specialization always lead to worse health outcomes, as

while with specialization, the effect on overall health is ambiguous and depends on the value of the specific parameters and on the level of obesity, vt:

Finally, for a given |$\bar{v}$|⁠, the effect of obesity-driven specialization on overall health is given by

suggesting that overall health may improve with specialization (for a given obesity burden), but only for the case of delivery-related health outcomes. Note that for delivery-unrelated health outcomes, where δ = γ = 0, specialization should have no effect (as compared to no specialization, or π = 0).

2.2. Testable Predictions

By construction, the first prediction that arises from the model is that increases in the obesity burden lead to increases in the CS rate, with or without specialization (the increase would be larger with specialization). In addition, the model captures the different pathways through which an effect of hospital-level obesity on birth outcomes may operate. In particular, the direct effect of obesity on newborn health can be seen by comparing |$\bar{h}^{T}(v,0)$| to |$\bar{h}^{T}(0,0)$|⁠, while the model’s predictions regarding how |$\bar{h}^{T}(v,0)$| and |$\bar{h}^{T}(v,\pi v)$| compare, allow us to test whether there is evidence of specialization. It is worth noting that the model yields predictions on the average health of newborns (overall and by type) without requiring information on individual mothers’ weight. The model’s main predictions, for both delivery-related and delivery-unrelated health outcomes, are summarized in tables 1 and 2.

Table 1.

Delivery-Related Outcomes

Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|??
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|??
Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|??
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|??

Source: Authors’ calculations.

Table 1.

Delivery-Related Outcomes

Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|??
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|??
Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|??
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|??

Source: Authors’ calculations.

Table 2.

Delivery-Unrelated Outcomes

Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|?=
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|?x > B
Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|?=
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|?x > B

Source: Authors’ calculations.

Table 2.

Delivery-Unrelated Outcomes

Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|?=
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|?x > B
Baseline (B)
|$\bar{h}^{\mathrm{all}}(0,0)$||$\bar{h}^{\mathrm{CS}}(0,0)$||$\bar{h}^{V}(0,0)$||$\bar{h}^{\mathrm{all}}(v,0)$||$\bar{h}^{\mathrm{CS}}(v,0)$||$\bar{h}^{V}(v,0)$|
(x)(1)(2)(3)(4)(5)(6)
A.|$\bar{h}^{\mathrm{all}}(v,0)$|x < B=
B.|$\bar{h}^{\mathrm{CS}}(v,0)$|?=
C.|$\bar{h}^{V}(v,0)$|?=
D.|$\bar{h}^{\mathrm{all}}(v,\pi v)$|?=
E.|$\bar{h}^{\mathrm{CS}}(v,\pi v)$|?x > B
F.|$\bar{h}^{V}(v,\pi v)$|?x > B

Source: Authors’ calculations.

Considering delivery-related outcomes, by comparing cells A1 and D4, as obesity increases, a positive effect on overall newborn (delivery related) health may only be observed if there is obesity-driven specialization, for instance.

Regarding delivery-unrelated outcomes, cell D4 would suggest a null effect from specialization for delivery-unrelated health outcomes for overall newborn health with respect to the no-specialization case.

3. Data

This study uses two sources of data. The first data source comprises the administrative health data from the universe of public health hospitals managed by the Ministry of Health (SSA) in Mexico. These data include information on inpatient hospital discharges and diagnoses recorded using ICD-10 codes for all SSA hospitals between 2008 and 2015,6 as well as patients’ age, sex, and municipality of residence. The analysis considers all SSA public health facilities for which at least eight births were observed in every month between 2008 and 2015.7,8 In addition, to alleviate concerns of endogenous selection, the study focuses on a balanced panel of 226 hospitals which are the only public health facility in the municipality (as explained in the next Section).

The second data source is the Mexican Birth Certificate Data, a census of all the registries of live births.9 While this census contains information on all births, our analysis focuses only on those public hospitals administered by the Ministry of Health. The birth certificate contains information related to the type of birth delivery (vaginal or cesarean), the newborn’s sex, birth weight, gestational age (in weeks), and Apgar score.10

It also includes information on whether the birth was a multiple birth, birth order, and some basic socioeconomic characteristics of the mother (age, education, civil status, and the municipality of birth and residence). As birth certificates include a variable for any congenital anomaly, illness, or lesion of the newborn (also using ICD-10 codes), dummies for whether there was any birth injury or trauma are also included.11 Finally, birth certificates include a variable indicating whether the mother survived after delivery, which we use as an outcome.

The study considers as women of reproductive age (WRA) those 15–45 years of age. Following Card, Fenizia, and Silver (2023), the analysis focuses on low-risk mothers given the variables available in our data. Although women who previously had a CS are generally more likely to undergo a CS in subsequent births, the sample is not restricted to first-time mothers. As there are no differences in the data on the likelihood of CS by birth order, the analysis considers all low-risk mothers.12 The low-risk sample excludes women under 18 and over 35, multiple births, births under 37 weeks of gestation, and women with more than 20 prenatal visits. Selecting mothers on these criteria allows us to exclude the women at a higher risk for CS (high-risk births are defined as those that do not satisfy any one of the aforementioned conditions), although it does not necessarily allow us to focus on pregnancies that are at a very low risk for a CS.13

Only live births are considered, as the dataset available for fetal deaths does not include a hospital identifier, so it cannot be matched to the data. Excluding fetal deaths from the analysis may imply a survival bias for our estimates, although the direction of the bias is a priori unclear. In particular, a direct (negative) effect of obesity on health leading to an increase in overall fetal deaths may lead to overstated health effects, while (positive) effects of specialization leading to a higher number of surviving newborns would likely bias our estimates downwards.

Since the administrative data on hospital inpatient discharges and from birth certificates do not include information on patients’ individual weight and height, this study creates a measure of hospital-level obesity. It uses ICD-10 patient diagnoses to identify obesity-related hospital discharges based on the medical literature (Kotchen 2010; Cnattingius et al. 2013; Kearns et al. 2014). Then it calculates the fraction of obesity-related discharges for WRA per month. The hospital obesity measure considers patient diagnoses for the following obesity-related diseases: pregestational hypertension (I10–I15, O10, and O11), pregestational diabetes (E10–E14 and 0240–O243), preeclampsia (O14 and O15), gestational diabetes (O244), hypertension, diabetes (E11.9), heart disease (I15.9, I20–I25, I25, I25.1), emphysema (J43.9), elevated blood pressure (R03.0), knee osteoarthritis (M17), gallbladder disease, hyperlipidemia (E78.5), chronic bronchitis (J41.0), stroke (I63.9), and asthma (J44.9, J69.8, J82). It also considers pregnancy complications associated with obesity such as preeclampsia, gestational diabetes, pregestational diabetes, and pregestational hypertension—though these may not be linked to specific births. Arguably, using diagnoses information from discharge data is an innovative way to obtain a useful and alternative measure of hospital-level obesity that may be valuable in other settings when patient-level weight data are not available.

The information on obesity prevalence is only available at the state level from national health surveys,14 while the obesity index data come from the administrative inpatient hospital records. As a validation exercise, supplementary online appendix table S2.1 shows that the aggregated index at the state level has a positive and statistically significant correlation with the state-level obesity and overweight rates, by regressing state obesity on the obesity index.15 If the estimates were to be in levels, a coefficient closer to one would reveal that the index perfectly captures obesity rates. However, the empirical strategy uses changes in the index not levels. The raw correlation of the obesity index with state obesity (overweight and obesity) is 0.49 (0.35).

Table 3.

Summary Statistics by Hospital-Level Obesity Index

High-obesity hospitalsLow-obesity hospitals
MeanStd. dev.MeanStd. dev.Difft-test
Panel A: Mother characteristics
Age24.260.4024.450.420.19***(3.47)
Schooling8.730.758.380.77−0.35***(−3.44)
Married0.350.140.350.140.00(0.10)
Panel B: Birth outcomes
Delivery-related outcomes
CS0.320.100.280.10−0.04***(−3.07)
Apgar8.910.118.920.130.01(0.82)
Apgar <70.010.000.010.010.00(0.58)
Birth injury5.918.706.288.870.37(0.31)
Birth trauma2.273.672.735.430.46(0.75)
Mother survival996.804.65997.075.510.27(0.40)
No. prenatal visits6.320.696.270.74−0.05(−0.54)
Non-delivery-related outcomes
Deliveries124.40101.1589.6975.38−34.71***(−2.92)
Gestational weeks39.370.1439.400.150.03*(1.80)
Low birth weight0.010.000.010.00−0.00(−0.42)
Weight3,268.3484.543,248.5483.78−19.80*(−1.77)
Obesity index0.600.65−0.470.31−1.07***(−15.99)
No. hospitals113113226
High-obesity hospitalsLow-obesity hospitals
MeanStd. dev.MeanStd. dev.Difft-test
Panel A: Mother characteristics
Age24.260.4024.450.420.19***(3.47)
Schooling8.730.758.380.77−0.35***(−3.44)
Married0.350.140.350.140.00(0.10)
Panel B: Birth outcomes
Delivery-related outcomes
CS0.320.100.280.10−0.04***(−3.07)
Apgar8.910.118.920.130.01(0.82)
Apgar <70.010.000.010.010.00(0.58)
Birth injury5.918.706.288.870.37(0.31)
Birth trauma2.273.672.735.430.46(0.75)
Mother survival996.804.65997.075.510.27(0.40)
No. prenatal visits6.320.696.270.74−0.05(−0.54)
Non-delivery-related outcomes
Deliveries124.40101.1589.6975.38−34.71***(−2.92)
Gestational weeks39.370.1439.400.150.03*(1.80)
Low birth weight0.010.000.010.00−0.00(−0.42)
Weight3,268.3484.543,248.5483.78−19.80*(−1.77)
Obesity index0.600.65−0.470.31−1.07***(−15.99)
No. hospitals113113226

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: High-obesity hospitals are hospitals with an average standardized obesity index higher than the sample median. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. Schooling is defined as number of years of education. Married is an indicator equal to one if the mother is currently married. Apgar is based on a total score between 0 and 10, where a higher score reflects a better health outcome. Apgar <7 is a dummy variable for whether the Apgar is less than 7. Birth injury is a dummy for the ICD-10 codes P10–15, P209–P221, P219–P221, P228, P229, P240, P285, P011, P032, P368, P369, P399, P545, and P914 multiplied by 1,000. Birth trauma is a dummy for the ICD-10 codes P10–15 multiplied by 1,000. Mother survival is a dummy for mother survival multiplied by 1,000. Low birth weight is a dummy variable for having low birth weight (lower than 2,500 g). Obesity index: Standardized monthly obesity index lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 3.

Summary Statistics by Hospital-Level Obesity Index

High-obesity hospitalsLow-obesity hospitals
MeanStd. dev.MeanStd. dev.Difft-test
Panel A: Mother characteristics
Age24.260.4024.450.420.19***(3.47)
Schooling8.730.758.380.77−0.35***(−3.44)
Married0.350.140.350.140.00(0.10)
Panel B: Birth outcomes
Delivery-related outcomes
CS0.320.100.280.10−0.04***(−3.07)
Apgar8.910.118.920.130.01(0.82)
Apgar <70.010.000.010.010.00(0.58)
Birth injury5.918.706.288.870.37(0.31)
Birth trauma2.273.672.735.430.46(0.75)
Mother survival996.804.65997.075.510.27(0.40)
No. prenatal visits6.320.696.270.74−0.05(−0.54)
Non-delivery-related outcomes
Deliveries124.40101.1589.6975.38−34.71***(−2.92)
Gestational weeks39.370.1439.400.150.03*(1.80)
Low birth weight0.010.000.010.00−0.00(−0.42)
Weight3,268.3484.543,248.5483.78−19.80*(−1.77)
Obesity index0.600.65−0.470.31−1.07***(−15.99)
No. hospitals113113226
High-obesity hospitalsLow-obesity hospitals
MeanStd. dev.MeanStd. dev.Difft-test
Panel A: Mother characteristics
Age24.260.4024.450.420.19***(3.47)
Schooling8.730.758.380.77−0.35***(−3.44)
Married0.350.140.350.140.00(0.10)
Panel B: Birth outcomes
Delivery-related outcomes
CS0.320.100.280.10−0.04***(−3.07)
Apgar8.910.118.920.130.01(0.82)
Apgar <70.010.000.010.010.00(0.58)
Birth injury5.918.706.288.870.37(0.31)
Birth trauma2.273.672.735.430.46(0.75)
Mother survival996.804.65997.075.510.27(0.40)
No. prenatal visits6.320.696.270.74−0.05(−0.54)
Non-delivery-related outcomes
Deliveries124.40101.1589.6975.38−34.71***(−2.92)
Gestational weeks39.370.1439.400.150.03*(1.80)
Low birth weight0.010.000.010.00−0.00(−0.42)
Weight3,268.3484.543,248.5483.78−19.80*(−1.77)
Obesity index0.600.65−0.470.31−1.07***(−15.99)
No. hospitals113113226

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: High-obesity hospitals are hospitals with an average standardized obesity index higher than the sample median. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. Schooling is defined as number of years of education. Married is an indicator equal to one if the mother is currently married. Apgar is based on a total score between 0 and 10, where a higher score reflects a better health outcome. Apgar <7 is a dummy variable for whether the Apgar is less than 7. Birth injury is a dummy for the ICD-10 codes P10–15, P209–P221, P219–P221, P228, P229, P240, P285, P011, P032, P368, P369, P399, P545, and P914 multiplied by 1,000. Birth trauma is a dummy for the ICD-10 codes P10–15 multiplied by 1,000. Mother survival is a dummy for mother survival multiplied by 1,000. Low birth weight is a dummy variable for having low birth weight (lower than 2,500 g). Obesity index: Standardized monthly obesity index lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 4.

Reduced-Form Evidence: The Effect of Hospital-Level Obesity Index on Mothers’ Characteristics

AgeAge groupSchoolingMarried
(1)(2)(3)(4)
Panel A: All births
Obesity index0.0020.0000.010−0.001
(STD lagged one year)(0.009)(0.001)(0.009)(0.001)
Mean dep. variable (2008)24.381.588.170.38
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−0.015−0.0020.002−0.001
(STD lagged one year)(0.017)(0.002)(0.013)(0.002)
Mean dep. variable (2008)24.611.618.530.40
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.0080.0010.010−0.001
(STD lagged one year)(0.011)(0.002)(0.010)(0.001)
Mean dep. variable (2008)24.281.568.020.37
No. mothers1,383,1481,383,1481,383,1481,383,148
AgeAge groupSchoolingMarried
(1)(2)(3)(4)
Panel A: All births
Obesity index0.0020.0000.010−0.001
(STD lagged one year)(0.009)(0.001)(0.009)(0.001)
Mean dep. variable (2008)24.381.588.170.38
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−0.015−0.0020.002−0.001
(STD lagged one year)(0.017)(0.002)(0.013)(0.002)
Mean dep. variable (2008)24.611.618.530.40
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.0080.0010.010−0.001
(STD lagged one year)(0.011)(0.002)(0.010)(0.001)
Mean dep. variable (2008)24.281.568.020.37
No. mothers1,383,1481,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for month and hospital fixed effects. Standard errors clustered at the hospital-month level. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. The age group variable is defined as 1 for mothers 18–24, as 2 for mothers 25–30, and as 3 for mothers 30–35. Schooling is defined as the number of years of education. Married is an indicator equal to 1 if the mother is currently married. Obesity index: standardized monthly lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 4.

Reduced-Form Evidence: The Effect of Hospital-Level Obesity Index on Mothers’ Characteristics

AgeAge groupSchoolingMarried
(1)(2)(3)(4)
Panel A: All births
Obesity index0.0020.0000.010−0.001
(STD lagged one year)(0.009)(0.001)(0.009)(0.001)
Mean dep. variable (2008)24.381.588.170.38
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−0.015−0.0020.002−0.001
(STD lagged one year)(0.017)(0.002)(0.013)(0.002)
Mean dep. variable (2008)24.611.618.530.40
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.0080.0010.010−0.001
(STD lagged one year)(0.011)(0.002)(0.010)(0.001)
Mean dep. variable (2008)24.281.568.020.37
No. mothers1,383,1481,383,1481,383,1481,383,148
AgeAge groupSchoolingMarried
(1)(2)(3)(4)
Panel A: All births
Obesity index0.0020.0000.010−0.001
(STD lagged one year)(0.009)(0.001)(0.009)(0.001)
Mean dep. variable (2008)24.381.588.170.38
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−0.015−0.0020.002−0.001
(STD lagged one year)(0.017)(0.002)(0.013)(0.002)
Mean dep. variable (2008)24.611.618.530.40
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.0080.0010.010−0.001
(STD lagged one year)(0.011)(0.002)(0.010)(0.001)
Mean dep. variable (2008)24.281.568.020.37
No. mothers1,383,1481,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for month and hospital fixed effects. Standard errors clustered at the hospital-month level. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. The age group variable is defined as 1 for mothers 18–24, as 2 for mothers 25–30, and as 3 for mothers 30–35. Schooling is defined as the number of years of education. Married is an indicator equal to 1 if the mother is currently married. Obesity index: standardized monthly lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 5.

Higher Hospital-Level Obesity Index Increases the Probability of CS

All mothersAll mothersFTMFTM
(low risk)(low risk)
(1)(2)(3)(4)
Obesity index0.003**0.003**0.005**0.004**
(STD one year lagged)(0.001)(0.001)(0.002)(0.002)
Mean dep. variable (2008)0.290.310.340.33
No. mothers1,988,8712,805,181641,2901,029,699
All mothersAll mothersFTMFTM
(low risk)(low risk)
(1)(2)(3)(4)
Obesity index0.003**0.003**0.005**0.004**
(STD one year lagged)(0.001)(0.001)(0.002)(0.002)
Mean dep. variable (2008)0.290.310.340.33
No. mothers1,988,8712,805,181641,2901,029,699

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for a woman’s age, education, marital status, gestational weeks, number of prenatal visits, number of births, and month and hospital fixed effects. Standard errors clustered at the hospital-month level. Low-risk mothers are defined as mothers between 18 and 35 years old excluding multiple births, 37 or fewer gestation weeks and with fewer than 20 prenatal visits. FTM denotes first-time mothers. Obesity index: standardized monthly obesity index lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 5.

Higher Hospital-Level Obesity Index Increases the Probability of CS

All mothersAll mothersFTMFTM
(low risk)(low risk)
(1)(2)(3)(4)
Obesity index0.003**0.003**0.005**0.004**
(STD one year lagged)(0.001)(0.001)(0.002)(0.002)
Mean dep. variable (2008)0.290.310.340.33
No. mothers1,988,8712,805,181641,2901,029,699
All mothersAll mothersFTMFTM
(low risk)(low risk)
(1)(2)(3)(4)
Obesity index0.003**0.003**0.005**0.004**
(STD one year lagged)(0.001)(0.001)(0.002)(0.002)
Mean dep. variable (2008)0.290.310.340.33
No. mothers1,988,8712,805,181641,2901,029,699

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for a woman’s age, education, marital status, gestational weeks, number of prenatal visits, number of births, and month and hospital fixed effects. Standard errors clustered at the hospital-month level. Low-risk mothers are defined as mothers between 18 and 35 years old excluding multiple births, 37 or fewer gestation weeks and with fewer than 20 prenatal visits. FTM denotes first-time mothers. Obesity index: standardized monthly obesity index lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 6.

Reduced-Form Evidence: The Effect of Hospital-Level Obesity Index on Delivery-Related Birth Outcomes

Mother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)
Panel A: All births
Obesity index1.239***−0.533*−0.130−0.000
(STD one year lagged)(0.207)(0.279)(0.141)(0.000)
Mean dep. variable (2008)991.754.181.480.01
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index1.798***−1.066**−0.144−0.001*
(STD one year lagged)(0.315)(0.427)(0.141)(0.000)
Mean dep. variable (2008)989.275.451.180.01
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.995***−0.240−0.118−0.000
(STD one year lagged)(0.183)(0.285)(0.177)(0.000)
Mean dep. variable (2008)992.773.651.600.01
No. mothers1,383,1481,383,1481,383,1481,383,148
Mother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)
Panel A: All births
Obesity index1.239***−0.533*−0.130−0.000
(STD one year lagged)(0.207)(0.279)(0.141)(0.000)
Mean dep. variable (2008)991.754.181.480.01
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index1.798***−1.066**−0.144−0.001*
(STD one year lagged)(0.315)(0.427)(0.141)(0.000)
Mean dep. variable (2008)989.275.451.180.01
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.995***−0.240−0.118−0.000
(STD one year lagged)(0.183)(0.285)(0.177)(0.000)
Mean dep. variable (2008)992.773.651.600.01
No. mothers1,383,1481,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, birth weight, gestational weeks, number of births, number of prenatal consults, and month and hospital fixed effects. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. Standard errors clustered at the hospital-month level. Apgar is based on a total score between 0 and 10, where a higher score reflects a better health outcome. Apgar <7 is a dummy variable for whether the Apgar is less than 7. Birth injury is a dummy for the ICD-10 codes P10–15, P209–P221, P219–P221, P228, P229, P240, P285, P011, P032, P368, P369, P399, P545, and P914 multiplied by 1,000. Birth trauma is a dummy for the ICD-10 codes P10–15 multiplied by 1,000. Mother survival is a dummy for mother survival multiplied by 1,000. The (monthly) obesity index is standardized and lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 6.

Reduced-Form Evidence: The Effect of Hospital-Level Obesity Index on Delivery-Related Birth Outcomes

Mother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)
Panel A: All births
Obesity index1.239***−0.533*−0.130−0.000
(STD one year lagged)(0.207)(0.279)(0.141)(0.000)
Mean dep. variable (2008)991.754.181.480.01
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index1.798***−1.066**−0.144−0.001*
(STD one year lagged)(0.315)(0.427)(0.141)(0.000)
Mean dep. variable (2008)989.275.451.180.01
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.995***−0.240−0.118−0.000
(STD one year lagged)(0.183)(0.285)(0.177)(0.000)
Mean dep. variable (2008)992.773.651.600.01
No. mothers1,383,1481,383,1481,383,1481,383,148
Mother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)
Panel A: All births
Obesity index1.239***−0.533*−0.130−0.000
(STD one year lagged)(0.207)(0.279)(0.141)(0.000)
Mean dep. variable (2008)991.754.181.480.01
No. mothers1,988,8711,988,8711,988,8711,988,871
Panel B: CS births
Obesity index1.798***−1.066**−0.144−0.001*
(STD one year lagged)(0.315)(0.427)(0.141)(0.000)
Mean dep. variable (2008)989.275.451.180.01
No. mothers605,723605,723605,723605,723
Panel C: Vaginal births
Obesity index0.995***−0.240−0.118−0.000
(STD one year lagged)(0.183)(0.285)(0.177)(0.000)
Mean dep. variable (2008)992.773.651.600.01
No. mothers1,383,1481,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, birth weight, gestational weeks, number of births, number of prenatal consults, and month and hospital fixed effects. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only considers municipalities with one clinic. Standard errors clustered at the hospital-month level. Apgar is based on a total score between 0 and 10, where a higher score reflects a better health outcome. Apgar <7 is a dummy variable for whether the Apgar is less than 7. Birth injury is a dummy for the ICD-10 codes P10–15, P209–P221, P219–P221, P228, P229, P240, P285, P011, P032, P368, P369, P399, P545, and P914 multiplied by 1,000. Birth trauma is a dummy for the ICD-10 codes P10–15 multiplied by 1,000. Mother survival is a dummy for mother survival multiplied by 1,000. The (monthly) obesity index is standardized and lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 7.

The Effect of Hospital-Level Obesity Index on Delivery-Unrelated Birth Outcomes

Birth weightLow birth weightMacrosomia
(1)(2)(3)
Panel A: All births
Obesity index−0.728−0.309*0.223***
(STD one year lagged)(0.890)(0.165)(0.063)
Mean dep. variable (2008)3,263.919.332.12
No. mothers1,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−2.159−0.3500.323**
(STD one year lagged)(1.430)(0.298)(0.140)
Mean dep. variable (2008)3,316.228.743.90
No. mothers605,723605,723605,723
Panel C: Vaginal births
Obesity index−0.896−0.2620.176***
(STD one year lagged)(1.006)(0.192)(0.058)
Mean dep. variable (2008)3,242.399.581.38
No. mothers1,383,1481,383,1481,383,148
Birth weightLow birth weightMacrosomia
(1)(2)(3)
Panel A: All births
Obesity index−0.728−0.309*0.223***
(STD one year lagged)(0.890)(0.165)(0.063)
Mean dep. variable (2008)3,263.919.332.12
No. mothers1,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−2.159−0.3500.323**
(STD one year lagged)(1.430)(0.298)(0.140)
Mean dep. variable (2008)3,316.228.743.90
No. mothers605,723605,723605,723
Panel C: Vaginal births
Obesity index−0.896−0.2620.176***
(STD one year lagged)(1.006)(0.192)(0.058)
Mean dep. variable (2008)3,242.399.581.38
No. mothers1,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, number of births, number of prenatal consults, gestational weeks, and month and hospital fixed effects. Standard errors clustered at the hospital-month level. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only using municipalities with one clinic. Low birth weight is a dummy variable for having low birth weight (lower than 2,500 g). Macrosomia is a dummy variable for having fetal macrosomia (birth weight higher than 4,500 g). LBW and macrosomia are multiplied × 1, 000. Obesity index: standardized monthly lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 7.

The Effect of Hospital-Level Obesity Index on Delivery-Unrelated Birth Outcomes

Birth weightLow birth weightMacrosomia
(1)(2)(3)
Panel A: All births
Obesity index−0.728−0.309*0.223***
(STD one year lagged)(0.890)(0.165)(0.063)
Mean dep. variable (2008)3,263.919.332.12
No. mothers1,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−2.159−0.3500.323**
(STD one year lagged)(1.430)(0.298)(0.140)
Mean dep. variable (2008)3,316.228.743.90
No. mothers605,723605,723605,723
Panel C: Vaginal births
Obesity index−0.896−0.2620.176***
(STD one year lagged)(1.006)(0.192)(0.058)
Mean dep. variable (2008)3,242.399.581.38
No. mothers1,383,1481,383,1481,383,148
Birth weightLow birth weightMacrosomia
(1)(2)(3)
Panel A: All births
Obesity index−0.728−0.309*0.223***
(STD one year lagged)(0.890)(0.165)(0.063)
Mean dep. variable (2008)3,263.919.332.12
No. mothers1,988,8711,988,8711,988,871
Panel B: CS births
Obesity index−2.159−0.3500.323**
(STD one year lagged)(1.430)(0.298)(0.140)
Mean dep. variable (2008)3,316.228.743.90
No. mothers605,723605,723605,723
Panel C: Vaginal births
Obesity index−0.896−0.2620.176***
(STD one year lagged)(1.006)(0.192)(0.058)
Mean dep. variable (2008)3,242.399.581.38
No. mothers1,383,1481,383,1481,383,148

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, number of births, number of prenatal consults, gestational weeks, and month and hospital fixed effects. Standard errors clustered at the hospital-month level. The sample consists of all low-risk mothers (i.e. mothers 18–35 years of age excluding multiple births, 37 or fewer gestational weeks, and with fewer than 20 prenatal visits) and only using municipalities with one clinic. Low birth weight is a dummy variable for having low birth weight (lower than 2,500 g). Macrosomia is a dummy variable for having fetal macrosomia (birth weight higher than 4,500 g). LBW and macrosomia are multiplied × 1, 000. Obesity index: standardized monthly lagged one year. *p < 0.1, **p < 0.05, ***p < 0.01.

Table 8.

Robustness Tests: The Effect of Hospital-Level Obesity Index on Delivery-Related Birth Outcomes in Premature Births

CSMother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)(5)
Panel A: Premature birth def. 1
Obesity index0.0001.241−1.885−3.907**−0.002
(STD one year lagged)(0.017)(2.643)(5.894)(1.817)(0.004)
Mean dep. variable (2008)0.66986.0125.172.800.02
No. mothers5,9985,9985,9985,9985,998
Panel B: Premature birth def. 2
Obesity index0.0013.216***−1.068−0.1960.001
(STD one year lagged)(0.005)(0.842)(3.547)(0.341)(0.002)
Mean dep. variable (2008)0.57988.5038.791.390.05
No. mothers66,94066,94066,94066,94066,940
CSMother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)(5)
Panel A: Premature birth def. 1
Obesity index0.0001.241−1.885−3.907**−0.002
(STD one year lagged)(0.017)(2.643)(5.894)(1.817)(0.004)
Mean dep. variable (2008)0.66986.0125.172.800.02
No. mothers5,9985,9985,9985,9985,998
Panel B: Premature birth def. 2
Obesity index0.0013.216***−1.068−0.1960.001
(STD one year lagged)(0.005)(0.842)(3.547)(0.341)(0.002)
Mean dep. variable (2008)0.57988.5038.791.390.05
No. mothers66,94066,94066,94066,94066,940

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, gestational weeks, number of births, number of prenatal visits, and month and hospital fixed effects. For birth outcomes, birth weight is included as a control. Two definitions of premature births are used. Premature birth definition 1: first time mothers (age >30) gestational age between 30 and 37 weeks. Premature birth definition 2: all mothers, gestational age 30–37, and only low-birth-weight births. Significance levels: *p < 0.1, **p < 0.05, p < 0.01.

Table 8.

Robustness Tests: The Effect of Hospital-Level Obesity Index on Delivery-Related Birth Outcomes in Premature Births

CSMother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)(5)
Panel A: Premature birth def. 1
Obesity index0.0001.241−1.885−3.907**−0.002
(STD one year lagged)(0.017)(2.643)(5.894)(1.817)(0.004)
Mean dep. variable (2008)0.66986.0125.172.800.02
No. mothers5,9985,9985,9985,9985,998
Panel B: Premature birth def. 2
Obesity index0.0013.216***−1.068−0.1960.001
(STD one year lagged)(0.005)(0.842)(3.547)(0.341)(0.002)
Mean dep. variable (2008)0.57988.5038.791.390.05
No. mothers66,94066,94066,94066,94066,940
CSMother survivalBirth injuryBirth traumaApgar <7
(1)(2)(3)(4)(5)
Panel A: Premature birth def. 1
Obesity index0.0001.241−1.885−3.907**−0.002
(STD one year lagged)(0.017)(2.643)(5.894)(1.817)(0.004)
Mean dep. variable (2008)0.66986.0125.172.800.02
No. mothers5,9985,9985,9985,9985,998
Panel B: Premature birth def. 2
Obesity index0.0013.216***−1.068−0.1960.001
(STD one year lagged)(0.005)(0.842)(3.547)(0.341)(0.002)
Mean dep. variable (2008)0.57988.5038.791.390.05
No. mothers66,94066,94066,94066,94066,940

Source: Authors’ analysis based on data from birth/hospitalization registries from Mexico’s Ministry of Health (Secretaria de Salud, SSA).

Note: All models control for woman’s age, schooling, marital status, gestational weeks, number of births, number of prenatal visits, and month and hospital fixed effects. For birth outcomes, birth weight is included as a control. Two definitions of premature births are used. Premature birth definition 1: first time mothers (age >30) gestational age between 30 and 37 weeks. Premature birth definition 2: all mothers, gestational age 30–37, and only low-birth-weight births. Significance levels: *p < 0.1, **p < 0.05, p < 0.01.

Table 3 shows descriptive statistics of the main variables used in the empirical models in 2008, for high- and low-obesity hospitals, as measured by the obesity index. A hospital is classified as high obesity if its average obesity index in 2008 is higher than the median of the sample. As mentioned earlier, the estimation sample includes all mothers who are at low risk of having a CS birth as earlier defined. Regarding mothers’ characteristics, there are small albeit statistically significant differences between high- and low-obesity hospitals in women’s age and schooling, but not in marital status. Mothers are slightly more educated in high-obesity hospitals but, this difference is relatively small in magnitude (0.3 of a year of schooling). In any case, any selection concerns in the patient pool are addressed in our empirical strategy and the models control for maternal characteristics.

Hospitals with high obesity prevalence tend to have a higher number of monthly CS deliveries than hospitals with low obesity prevalence in 2008, the baseline year for our analysis period. The fraction of birth deliveries by CS is 32 percent in high-obesity hospitals and 28 percent in low-obesity hospitals. There are no statistically significant differences between high- and low-obesity hospitals in delivery-related birth outcomes, such as average Apgar scores as well as a dummy variable for Apgar scores lower than 7, frequencies of birth injury, birth trauma, and mother survival. Similarly, differences in delivery-unrelated outcomes, such as low birth weight (LBW), are statistically indistinguishable between high- and low-obesity hospitals. Nevertheless, newborns are heavier in high-obesity hospitals, which is consistent with results linking heavier mothers to heavier newborns and fetal macrosomia (Khashan and Kenny 2009). Finally, by construction, the obesity index is higher for hospitals with high obesity prevalence. As the empirical strategy exploits changes in the obesity index (and not levels), and time-invariant hospital characteristics are absorbed by hospital fixed effects, concerns about differences in both types of hospitals are limited.

During the study period, there is considerable variation in CS rates and obesity prevalence, both over time and across hospitals. Using the sample of hospitals, fig. 1 illustrates how CS rates and the hospital-level obesity measure varies over time, with respect to 2008 and shows that substantial variation existed in both variables. Furthermore, fig. 2 presents the geographic variation in trends for both CS rates and obesity levels and shows that increases in neither CS nor obesity index were concentrated in a single geographic area. Together, these figures indicate that there is enough time and spatial variation to identify the effects of the hospital-level obesity measure on a woman’s probability of having a CS and her newborn’s delivery-related birth outcomes.

4. Empirical Strategy

The analysis starts by considering the effect of hospital-level obesity prevalence on the likelihood of a woman having a cesarean birth. The identification comes from variation across time and hospitals in the obesity index, conditioning on mother, pregnancy, and newborn characteristics, using time and hospital fixed effects. In order to mitigate concerns regarding endogenous selection into hospitals, the main analysis focuses on municipalities where there is only one SSA hospital. The following model is estimated:

(7)

where Yitc is a dummy variable that captures whether birth i at time t (month) in hospital c was via CS. The variable Obesityindext − 1, c is the standardized obesity index for clinic c at time t − 1 (i.e. lagged one year).16,17 The variable Xi is a set of time-varying control variables associated with birth i at time t in hospital c, including mother’s age, education, and marital status, number of prenatal visits, weeks of gestation, and total number of births in hospital c in month t. The variables δt are time fixed effects (year-month) that capture countrywide time trends and the seasonality of births during the period of analysis. The variables θc are hospital fixed effects that control for all observed and unobserved time-invariant differences between hospitals. Standard errors are clustered at the hospital-month level.

As mentioned earlier, the main sample includes all low-risk mothers, defined as the sample of mothers aged 18–35 years, excluding multiple births, 37 or fewer gestation weeks, and with fewer than 20 prenatal visits. The analysis focuses on low-risk births, as they are less likely ex ante to require a medically recommended CS. And it is not limited to first-time mothers. The results of the obesity index on CS are presented for all births and first-time mothers, as well as low-risk first-time mothers.

The main threats to identification are endogenous patients’ selection into hospitals and differential changes in patient composition at high- versus low-obesity hospitals. Arguably, if women were to endogenously sort into high/low specialization hospitals according to their underlying health risks, such endogenous selection would be likely to bias the results against finding a positive effect of obesity on health outcomes. That is because women more at risk would choose to deliver in highly specialized CS hospitals (and if our hypothesis holds, in high-obesity hospitals), thus leading to worse health outcomes in general. Nevertheless, in order to mitigate concerns regarding endogenous selection, the analysis excludes hospitals in municipalities that have more than one public hospital and where maternal choice in terms of delivery place may be more concerning. Arguably, endogenous hospital choice is less likely to be present in municipalities with only one hospital. Additionally, these concerns are also alleviated as the study setting is limited to women delivering in public hospitals where women are assigned the clinic they must attend and therefore have limited choice. Hence, the analysis considers the panel of 226 hospitals for which births are observed for every time period and that are in municipalities served only by one hospital. The results are also robust to including the unbalanced panel sample of these hospitals across the period of interest.

Another threat to identification is related to hospitals’ patient composition. If the overall health or individual characteristics of a hospital’s patient pool have different trends for hospitals with different obesity, the estimates from equation (7) could be biased. To present evidence of the validity of the design, the main models are estimated using a set of predetermined maternal characteristics as the dependent variable. In particular, changes in the hospital’s measure of obesity are regressed against maternal age, marital status, and education as dependent variables. The estimates presented in table 4 show that changes in the obesity index do not appear to predict mothers’ age, marital status, and education. Even without correcting for multiple hypotheses testing, the coefficients for these dependent variables are not significant at conventional levels. These findings mitigate potential concerns that the results might be driven by changes in the population within a hospital catchment area. Additionally, in Section 6, Robustness Checks, the preferred specification is estimated including state-level time trends to further alleviate concerns that the results might be driven by secular trends in the health of patients in the hospital catchment area.

Using the same specification as in equation (7), the study estimates the reduced-form effects of hospital-level obesity prevalence on newborn health outcomes that could be affected by CS specialization during delivery: (a) whether the Apgar is less than 7 (Apgar <7)18; (b) maternal survival; (c) whether there was birth trauma; and (d) whether there was birth injury. Any positive effects of high-obesity hospitals on birth outcomes are likely to be lower bounds of “true” effects, as in general, increases in mothers’ weight may be associated with worse health outcomes in newborns (Khashan and Kenny 2009; Gallardo et al. 2015) and macrosomic newborns—more frequently born to overweight or obese mothers—are more likely to suffer from birth trauma or injury (Wollschlaeger et al. 1999; McDonald et al. 2010). Such bias would work against finding any positive effects of hospital-level obesity on newborns’ health.

The results on birth outcomes are presented for all births and separately for CS and vaginal births. Although presenting results by type of delivery is conditioned on an outcome, showing results separately may help depict a more complete picture of the mechanisms driving the results. Additionally, the model gives a theoretical framework with predictions regarding the average health of CS and vaginal births at the hospital level (with and without specialization). As predicted in the theoretical model, while specialization should not affect birth outcomes for newborns delivered vaginally directly, selection may lead to observing better health outcomes for vaginal births as marginally risky mothers who would have otherwise delivered vaginally are induced to deliver via CS, and the average risk of women delivering vaginally falls. This would be consistent with an improvement in newborn health outcomes for newborns born vaginally—for both delivery-related and delivery-unrelated health outcomes—due to selection. A similar (and positive) selection effect would be observed for CS newborns as the marginally risky mother switching from vaginal birth to CS would be relatively less risky than the average mother having a CS in the absence of specialization, thus decreasing average risk. However, specialization may also lead to worse delivery-related birth outcomes for vaginal births (the competing treatment) if physicians become less skilled at this procedure as they specialize in CS or mothers’ obesity is negatively correlated with newborn health. Hence, the expected direction of specialization on delivery-related birth outcomes for vaginal births is a priori ambiguous, as predicted by the theoretical model.

The study also explores the effects of obesity on delivery-unrelated outcomes. The theoretical model predicts that these outcomes should not be affected by specialization, but could be affected by changes in mothers’ health burden (e.g. obesity). For this, birth weight and dummies for LBW and macrosomia are used. It is worth noting that the literature suggests that there is no evidence of LBW being affected by mothers’ obesity when controlling for gestational age (McDonald et al. 2010) (in terms of the model, θ = 0) while there is evidence of obesity being associated with a higher probability of macrosomia (Wollschlaeger et al. 1999). When estimating the models using these birth outcomes that are arguably not determined at time of delivery, one would expect to see no (positive) effects from specialization if the specialization hypothesis is correct (although negative direct effects from increases in average mothers’ obesity in the case of macrosomia may be observed).

5. Results

Table 5 reports the effects of the hospital-level obesity measure on women’s likelihood of having a CS. It shows the estimates for all low-risk mothers (column 1), all mothers (column 2), first-time low-risk mothers (column 3), and all first-time mothers (column 4). For all four samples, there is a positive and statistically significant relationship between changes in obesity at the hospital level and changes in the individual probability that a woman delivers via CS. Considering all low-risk mothers, an increase of one standard deviation in hospital-level obesity (with respect to 2008, the baseline year) leads to a 0.3 percentage points (pp) increase in the probability of having a CS. Given the 2008 CS rates of 29 percent, this corresponds to an increase in CS of approximately 1 percent. This estimate is relevant as it corresponds to the sample of women who ex ante have a lower risk of having a CS. For this reason, the main analysis considers the sample of all low-risk mothers, as this may present a complete picture of the potential effects of CS specialization.

Moreover, supplementary online appendix table S2.2 explores whether the results vary according to whether deliveries happened during the “daytime” versus “nighttime and weekends” CS birth delivery, as birth certificates do not specify whether a CS was scheduled or not. Thus, following Costa-Ramón et al. (2018), who find that the proportion of women that deliver via an unplanned CS is highest during the nighttime, one may think of CS as unplanned if it occurred during night hours (i.e. 8pm–8am) or weekends. The 8pm–8am range is used, as the data show a discontinuity in CS rates at these times (see fig. S1.2), that most likely is related to changes in medical staff shifts. There is no evidence showing that CS specialization might be driven by doctors working during the “daytime.” Although the coefficient for the “daytime CS” is slightly larger than the effect for all births (column 1), it is only statistically significant at the 10 percent level and not statistically different from the coefficient for “nighttime and weekend,” reported in column 3.

Next, the reduced-form effects of our hospital-level obesity measure on newborn health for delivery-related and delivery-unrelated health outcomes are considered. The estimates in table 6, panel A, highlight a significant increase in mother survival for all births, as well as a decrease in the probability of experiencing birth injury or birth trauma, although the latter reduction is not statistically significant at conventional levels. No effects on Apgar scores were found. These results are consistent with higher physicians’ specialization in CS due to an increase in risky pregnancies at the hospital level from a heavier population of women of reproductive age. As doctors are more likely to perform CS, they become more skilled at this procedure and develop a preference for it. Delivery-related outcomes such as the likelihood of birth trauma or injury thus decrease with specialization. In particular, the results on the reduction of maternal mortality suggest remarkable returns of specialization in the Mexican context where maternal mortality is still a prevalent issue.

Panels B and C of table 6 show the effects of the obesity measure on health outcomes by delivery type: CS births and vaginal births. As these analyses are conditioned on an outcome, they should be interpreted with caution and taking into account that selection into CS is not exogenous: as discussed above, increases in obesity are likely to increase the likelihood of a woman receiving a CS. Nevertheless, these estimates are presented by delivery types to get a clearer picture of the potential mechanisms driving the health effects on all births. Panel B, column 1 shows that the increase in the probability of mother survival after delivery is mostly driven by CS. Similarly, panel B, column 2 indicates that the lower probability of birth injuries is driven by CS births. These effects are consistent with the hypothesis of specialization. Furthermore, panel B, column 4 shows a decrease in the probability of having an Apgar score <7 for CS births. Low Apgar scores are often the result of prolonged labor (Altman et al. 2015; Card, Fenizia, and Silver 2023), which is less likely for CS birth as they tend to be quicker (Costa-Ramón et al. 2018). Again, this finding supports the specialization hypothesis: a lower probability of a low Apgar in CS births can be explained by doctors becoming specialized in CS.

Furthermore, panel C shows improvements in mother survival for vaginal births and although we find reductions in birth injury and birth trauma, these are not statistically significant. These findings could be interpreted as a selection effect: to the extent that specialization leads marginally riskier women to deliver via CS, the average risk of women delivering vaginally falls—thus improving delivery-related health for newborns born vaginally as well.

In addition, supplementary online appendix table S2.3 shows the reduced form for birth delivery outcomes for all CS births (panel A), and separately for CS during the day (panel B) and CS at nighttime and on weekends (panel C). Arguably, the latter two categories may proxy for scheduled and unscheduled CS, respectively. A similar caveat to the one for table 6 applies: although the time of delivery is plausibly exogenous, the estimates are conditional on having a CS birth, thus should be interpreted with caution. As shown in column 1, CS births increase the probability of mother survival regardless of the time of day. The estimates in column 2 indicate a decrease in the probability of birth injury after delivery is stronger for “daytime CS,” and, even though it is not statistically significant for “nighttime and weekend CS,” both coefficients have the same sign and are not statistically different. This finding is also consistent with a specialization hypothesis.

Finally, table 7 shows the effects of our hospital-level obesity measure on delivery-unrelated health outcomes, namely on a set of transformations of the birth weight variable. Columns 1 and 2 show that the obesity index does not have a statistically significant effect on birth weight and low birth weight. Column 3 reports a positive effect of our hospital-level obesity index on the probability of a newborn being macrosomic (i.e. having a birth weight higher than 4,500 g). These results provide a falsification test for the specialization hypothesis, adding validity to the interpretation of the results as evidence of CS specialization. While one would expect an association between hospital-level obesity and macrosomia as mothers’ obesity is a risk factor for fetal macrosomia (Khashan and Kenny 2009), one would not expect birth weight or the probability of low birth weight (conditional on gestational age) to be affected by our hospital-obesity index (McDonald et al. 2010).

5.1. Mechanisms

The findings support the hypothesis that an increase in hospital-level obesity prevalence increases a woman’s probability of CS and affects newborns’ delivery-related outcomes through the specialization channel. However, it is possible that other mechanisms might be at play, such as preferences over the type of birth delivery, supply driven changes, or hospital congestion.

Supplementary online appendix table S2.4, columns 2–4 analyzes differential effects by mothers’ age. Variation in the measure of facility-level obesity on CS is concentrated among women under 30—although there is also a positive effect, albeit not statistically significant, for women 30–35. This result suggests that the effect of CS specialization is probably concentrated in marginal pregnancies, where the underlying risk of having a CS is not too high. In addition, supplementary online appendix table S2.4, columns 5–6 shows that the effect of hospital-level obesity measure on CS is both larger and statistically stronger for women with lower-than-average schooling (i.e. fewer than eight years of education). Furthermore, supplementary online appendix table S2.4, panel B shows that the returns of CS specialization in maternal mortality reduction are found across maternal age cohorts and levels of education. Overall, estimates from supplementary online appendix table S2.4 suggest that the results are unlikely to be driven by higher-educated women’s preferences over the type of birth delivery or by selection.

It is also possible that changes in CS are supply driven. In particular, if doctors in hospitals with higher increases in obesity face differential monetary incentives to perform CS, the estimates could be driven by the supply side. However, this is unlikely in this context: physicians working at public hospitals are salaried workers who earn a monthly wage. They do not receive any additional benefits based on the number or type of deliveries they perform.

An alternative explanation could be that the results are driven by greater incentives to perform CS due to space limitations, if changes in obesity at the hospital level are correlated with changes in hospital congestion. For instance, CS deliveries take less time than vaginal ones (Costa-Ramón et al. 2018). The specifications however, control for changes in the contemporaneous monthly volume of births at the hospital level, which would capture changes in congestion. Hence, the findings are not likely to be driven by congestion.

6. Robustness Checks

The analysis first considers the effects of the hospital-level obesity measure on delivery-related birth outcomes for a subgroup of newborns that are most frequently (almost always) born via CS: preterm newborns. In these cases, the use of CS is more frequent and arguably would be less likely to be affected by the physician’s choice. The models are estimated using two definitions of the sample preterm newborns. The first definition includes all infants born between 30 and 37 gestational weeks to first-time mothers 30 years old and older. The second definition includes babies with the same gestational age but conditions to low-birth-weight newborns from all mothers. Table 8 shows that, regardless of the definition used, the hospital-level obesity index is uncorrelated to the probability of a CS for this subset of births (and that the CS rate is in fact higher; i.e. 57–66 percent depending on the sample definition). However, when looking at delivery-related birth outcomes, table 8 shows evidence consistent with a specialization hypothesis: there is a positive association between the obesity index and delivery-related birth outcomes for preterm newborns—where no effect on the probability of CS is observed. A higher obesity index at the hospital level is associated with a lower probability of birth trauma and maternal mortality for preterm newborns. These results would (only) be expected if higher obesity improves physicians’ CS skills through CS specialization.

Next, supplementary online appendix table S2.5 shows that the results are robust to a series of alternative specifications. First, one might be concerned that the results are exclusively driven by the short-run fluctuation of the obesity index. The main analysis uses the previous year as the relevant measure for learning at time t as it allows for a long period for doctors to be specialized (and frequent enough to have sufficient variation in the data).19 Nevertheless, to mitigate concerns that the time window choice is not driving the results, supplementary online appendix table S2.5, panel A uses a moving average of the obesity index over 12 months instead of the lagged index. Estimates show that the main effects of the obesity index on CS, mother survival, birth injury, and birth trauma are overall robust to the measure used.

Second, the estimation sample is limited to the public hospitals that have at least eight births per month in every year of our study, and as such, it is a balanced panel of facilities. Panel B estimates the models using all public health facilities included in the SSA data by not restricting our estimation to hospitals that had at least eight monthly births during every year over our study period, i.e. the unbalanced panel of hospitals. Supplementary online appendix table S2.5, panel B shows that the results remain robust to this specification.

Third, panel C estimates model specifications that control for specific linear time trends to alleviate concerns that the results are driven by secular trends in the health of patients in the hospital catchment areas. As the estimates in supplementary online appendix table S2.5, panel C show, the effects of the hospital-level obesity measure on the probability of having a CS are robust to adding state-specific linear time trends. Therefore, there are no concerns that seasonal or omitted trends at the state level might be driving the results.

Fourth, supplementary online appendix table S2.5 panel D explores the effects of the obesity index on high-risk births. In line with the model’s predictions, one would expect to see a weak or null effect of obesity on the probability of CS for riskier pregnancies, whereas positive effects on this subgroup are likely to be evidence of specialization. Estimates show that there is only a marginally significant effect on CS for high-risk births, but a large and statistically significant improvement in maternal mortality. This result can be interpreted as further evidence of specialization.

As discussed previously, the exact relevant window for learning at time t for doctors is unknown. However, one would expect doctors to specialize in CS over time if they are exposed to riskier pregnancies (in this context, to more obese mothers). The empirical strategy captures both an increase in CS from treating more obese patients and a decrease in CS from treating fewer obese patients. In order to explore whether the results are being driven by increases or decreases in CS specialization (and in order to alleviate potential concerns regarding the definition of the “learning period”) the analysis presents the results from estimating a regression of the hospital-level CS time trend (over the entire study period) on its obesity index time trend. Intuitively, this analysis shows whether, in the longer term, trends in obesity are correlated with trends in CS. Furthermore, this regression can be estimated separately for those clinics where the estimated linear time trend for obesity was positive (overall increase in obesity) or negative (overall decrease in obesity). By doing this, it is possible to present some evidence on whether the results are being driven by doctors specializing in CS when they treated a riskier patient pool or whether they are driven by doctors performing fewer CS as obesity falls. Arguably, if it is hard to de-specialize once a doctor becomes specialized, one should see a positive effect for hospitals with increasing obesity, and a smaller or null effect for hospitals with a decrease in obesity. The estimates from this analysis are reported in supplementary online appendix table S2.6. Results show that a hospital’s obesity time trend is correlated with its CS time trend: increases in obesity are positively correlated with increases in CS. Furthermore, columns 2 and 3 show that this association is positive and statistically significant only for hospitals with an upward obesity trend, while it is indistinguishable from zero for those hospitals in which obesity trended downwards. These findings suggest that treating more obesity cases is likely to lead to performing more CS overall.

Finally, to further mitigate concerns regarding patient selection, the analysis explores whether selection into private hospitals could affect the results. Using the birth certificates data, one can calculate the share of births in a municipality that happened in private hospitals. As the analysis focuses on municipalities with only one public hospital, the share of births occurring at a private clinic is regressed on the obesity index of the (only) public hospital in a municipality (and municipality and year/month fixed effects are included to mimic the preferred specification). Supplementary online appendix table S2.7 shows that changes in obesity in the public hospital do not predict changes in private hospital share, suggesting that selection to/from the private sector is unlikely to be driving our results. Columns 1 and 2 include all municipalities in our sample (with column 2 including the total number of births as a control), while columns 3 and 4 focus on municipalities where the share of births in private facilities is nonzero. Coefficients are not significant across all specifications.

7. Conclusion

Over the past decades, CS rates have increased worldwide—particularly in middle- and low-income countries. Similarly, obesity prevalence—a risk factor for CS—has also dramatically increased. The literature on the determinants and consequences of increasing CS rates mostly focuses on high-income countries. Hence, the evidence for low- and middle-income countries, where maternal and neonatal mortality are still high and health systems are resource constrained, is lacking.

This paper studies how variation in hospital practices induced by high obesity levels might affect hospital CS rates and whether that has an impact on maternal and neonatal delivery-related outcomes in Mexico, a country where CS rates are the second largest in Latin America, 70 percent of the population is obese, and maternal and infant mortalities are still relatively high.

The analysis uses a novel measure of hospital-level obesity by exploiting data from patient diagnoses. It leverages the obesity-index variation across time and hospitals to identify the effect of obesity on a woman’s likelihood of having a CS.

The study finds that hospital-level obesity increases a woman’s probability of CS among a sample of mothers with low risk of CS. The findings suggest a potential story of CS specialization in which doctors develop a preference for CS procedures in hospitals with high female-obesity prevalence and become more skilled at performing CS.

While other mechanisms might be at play, the findings are not consistent with demand-driven channels (e.g. higher-educated women’s preferences over the type of birth delivery), with differential monetary incentives to perform CS in higher-obesity hospitals, or with greater incentives to perform CS due to hospital congestion. The results thus speak to how the use of intensive treatments may affect health outcomes (as in Doyle et al. (2015)), which is a key policy concern given the ample variation in treatment choice (i.e. CS) observed in Mexico and particularly relevant for resource-constrained settings.

Furthermore, CS specialization leads to better birth outcomes for all mothers with low risk of CS, consistent with evidence on productivity spillovers in healthcare (Chandra and Staiger 2007). Delivery-related birth outcomes improve for CS: there is a reduction in birth injuries, birth trauma, and an increase in mother survival. There is also a reduction in maternal mortality among mothers who deliver vaginally. The results might inform other low- and middle-income countries on the path of increasing CS and obesity rates.

To the extent that there are positive returns to specialization in CS, a policy recommendation derived from the results is that the need for CS increases as population health changes (e.g. obesity increases). Thus, having hospitals that are more/less specialized in performing CS and sorting women efficiently into these according to their pregnancy risk, could be welfare enhancing, particularly in resource-constrained settings.

While this study analyzes the immediate effects of CS on newborn and maternal health, it does not speak to the literature on the long-term effects of CS on mothers’ (e.g. Halla et al. (2020)) or children’s (e.g. Costa-Ramón et al. (2022)) health. While CS specialization may improve health outcomes in the very short term, improvements in some health outcomes, such as Apgar score, should be weighed against potential long-term costs associated with increased CS use. However, the benefits of increased mother survival from CS specialization are unlikely to be outweighed by other long-term costs associated with CS use.

Data Availability Statement

The data underlying this article are available in Zenodo, at https://zenodo.org/record/8270682. The datasets were derived from sources in the public domain. The publicly available data sources used in the paper are: 1. Egresos Hospitalarios, Secretaría de Salud, México (Dirección General de Información en Salud) 2008—2016. http://www.dgis.salud.gob.mx/contenidos/basesdedatos/da_egresoshosp_gobmx.html 2. Nacimientos, Subsistema de Información sobre Nacimientos (SINAC). Secretaría de Salud, Dirección General de Información en Salud. 2008—2016 http://www.dgis.salud.gob.mx/contenidos/basesdedatos/da_nacimientos_gobmx.html 3. Componente de Nutrición, Instituto de Salud Pública, México (Encuesta Nacional de Salud y Nutrición) 2012. https://ensanut.insp.mx/encuestas/ensanut2012/descargas.php 4. Componente de Nutrición, Instituto de Salud Pública, México (Encuesta Nacional de Salud y Nutrición) 2016. https://ensanut.insp.mx/encuestas/ensanut2016/descargas.php.

Footnotes

1

The 2012 and 2016 female OWOB prevalence are the latest estimates at the state and national levels overlapping with the study period.

2

Executive Commission of Attention to Victims, “Diagnosis on Victimization Due to Obstetric Violence in Mexico,” last retrieved at www.gob.mx/cms/uploads/attachment/file/194701/Diagno_stico_VO_port.pdf.

3

Dirección General de Información en Salud, Ministry of Health.

4

Recent studies have analyzed the effects of policy changes in the health sector on CS rates and newborn health in Brazil (De Oliveira, Lee, and Quintana-Domeque 2022; Melo and Menezes-Filho 2023) and Chile (de Elejalde and Giolito 2021).

5

Guendelman et al. (2017), using data from 2014 birth certificates, show that the type of insurance coverage and type of facilities (public versus private) are factors associated with women’s CS deliveries. Nevertheless, there is no evidence on the effects of hospital practices related to obesity levels on CS rates.

6

ICD-10 is the 10th revision of the International Statistical Classification of Diseases and Related Health Problems. It is a medical classification list by the World Health Organization (WHO). The data contain codes for diseases, signs and symptoms, abnormal findings, complaints, social circumstances, and external causes of injury or diseases.

7

Following Heckman (1981) and Greene (2001)’s rule of thumb.

8

While the vast majority of these are “second- and third-level healthcare facilities”—general and specialized hospitals—some are smaller clinics where births sometimes take place. Healthcare facilities are referred to as hospitals hereafter.

9

This census comes from the National Health Information System Birth Certificates (SINAC), collected by the Ministry of Health.

10

The Apgar score is a widely used diagnostic indicator of newborn health, ranging from 0 to 10 (Card, Fenizia, and Silver 2023). A higher score reflects a better health outcome.

11

Birth trauma is defined as a dummy variable for any of the ICD-10 codes from P10 to P15, following the WHO definition. Birth injury is defined as a dummy variable for any of the following ICD-10 codes: P10–P15, P209–P211, P219–P221, P228, P229, P240, P284, P285, P011, P032, P368, P369, P399, P545, and P914. The dummy variables for birth injury, birth trauma, and mother survival after birth delivery are multiplied by 1,000.

12

There is no evidence in the data that higher-order births are more likely to be CS (as compared to first births). In fact, there is evidence supporting the opposite: CS rates seem to be decreasing in birth order. In 2008 the CS rate was 33.3 percent for first-time mothers, 32.7 percent for second births, 32.4 percent for third births, and 31.1 percent for higher-order births.

13

Due to data limitations, it is not possible to identify vertex first births, mothers with BMI above the 90th percentile, and mothers with eclampsia, preeclampsia, or growth restrictions, all of which are also excluded from the analysis in Card, Fenizia, and Silver (2023).

14

The 2012 Encuesta Nacional de Salud y Nutricin (ENSANUT) is the latest available round of this survey overlapping with the study period. ENSANUT is representative at the state level.

15

Different specifications of the obesity index are tested by including different subsets of the obesity-related hospital diagnoses and the index that includes all diagnoses discussed above has the strongest correlation with obesity at the state level; see supplementary online appendix table S2.1.

16

The obesity index is standardized by subtracting the 2008 mean level across hospitals and dividing by the 2008 standard deviation. Thus, all estimates can be interpreted in terms of an increase equivalent to one standard deviation at baseline.

17

The lagged index is calculated as the average of the lagged calendar year (i.e. if the birth occurred on February 2010, the average index for 2009 was assigned).

18

The outcomes considered are a dummy variable for whether the Apgar is less than 7 (Apgar <7), as the American Academy of Pediatrics defines a five-minute Apgar score of 7–10 as reassuring, a score of 4–6 as moderately abnormal, and a score of 0–3 as low in the term infant and late-preterm infant (American Academy of Pediatrics 2015).

19

Results remain mostly unchanged when using the previous quarter.

Author Biography

Catalina Herrera-Almanza (corresponding author) is an assistant professor at the Department of Agricultural and Consumer Economics, University of Illinois at Urbana-Champaign; her email address is [email protected]. Fernanda Marquez-Padilla is an assistant professor at the Centro de Estudios Economicos, El Colegio de Mexico, Mexico. Her email address is [email protected]. Silvia Prina is an associate professor at the Department of Economics, Northeastern University. Her email address is [email protected]. The authors thank Ana Costa-Ramón, Janet Currie, Susan Parker, Heather Royer, Felix Matthys; conference participants at the University of Illinois, Urbana-Champaign, AAEA, ASHECON, PAA, Colmex Seminar Series, and at the ITAM Alumni Conference; the editor; and three anonymous referees for helpful comments. We thank Ana Paola Paredes Rodriguez for outstanding research assistance. Catalina Herrera-Almanza thanks the International Seed Grant through the College of ACES Office of International Programs at University of Illinois, Urbana-Champaign. All errors are our own. A supplementary online appendix is available with this article at The World Bank Economic Review website.

References

A Kotchen
 
T.
 
2010
. “
Obesity-Related Hypertension: Epidemiology, Pathophysiology, and Clinical Management
American Journal of Hypertension
.
23
(
11
):
1170
78
.

Altman
 
M.
,
Sandström
 
A.
,
Petersson
 
G.
,
Frisell
 
T.
,
Cnattingius
 
S.
,
Stephansson
 
O.
.
2015
. “
Prolonged Second Stage of Labor Is Associated with Low Apgar Score
European Journal of Epidemiology
.
30
(
11
):
1209
15
.

American Academy of Pediatrics
.
2015
. “
ACOG Committee Opinion. Number 644, October 2015 (Replaces No. 333, May 2006): The Apgar Score
.”.
Committee on Fetus and Newborn, ACOG and Committee on Obstetric Practice and American Academy of Pediatrics and others
.

Betran
 
A. P.
,
Ye
 
J.
,
Moller
 
A.-B.
,
Souza
 
J. P.
,
Zhang
 
J.
.
2021
. “
Trends and Projections of Caesarean Section Rates: Global and Regional Estimates
BMJ Global Health
.
6
(
6
):
e005671
.

Boerma
 
T.
,
Ronsmans
 
C.
,
Melesse
 
D. Y.
,
Barros
 
A. J.
,
Barros
 
F. C.
,
Juan
 
L.
,
Temmerman
 
M.
.
2018
. “
Global Epidemiology of Use of and Disparities in Caesarean Sections
Lancet
.
392
(
10155
):
1341
48
.

Bozzoli
 
C.
,
Deaton
 
A.
,
Quintana-Domeque
 
C.
.
2009
. “
Adult Height and Childhood Disease
Demography
.
46
(
4
):
647
69
.

Brenes-Monge
 
A.
,
Saavedra-Avendaño
 
B.
,
Alcalde-Rabanal
 
J.
,
Darney
 
B. G.
.
2019
. “
Are Overweight and Obesity Associated with Increased Risk of Cesarean Delivery in Mexico? A Cross-Sectional Study from the National Survey of Health and Nutrition
BMC Pregnancy and Childbirth
.
19
(
1
):
239
.

Card
 
D.
,
Fenizia
 
A.
,
Silver.
 
D.
 
2023
. “
The Health Impacts of Hospital Delivery Practices
.”.
American Economic Journal: Economic Policy.
 
15
(
2
):
42
81
.

Chandra
 
A.
,
Staiger
 
D. O.
.
2007
. “
Productivity Spillovers in Health Care: Evidence from the Treatment of Heart Attacks
Journal of Political Economy
.
115
(
1
):
103
40
.

Cnattingius
 
S.
,
Villamor
 
E.
,
Johansson
 
S.
,
Bonamy
 
A. K. E.
,
Persson
 
M.
,
Wikström
 
A. K.
,
Granath
 
F.
.
2013
. “
Maternal Obesity and Risk of Preterm Delivery
Lancet
.
309
(
22
):
2362
70
.

Costa-Ramón
 
A. M.
,
Rodríguez-González
 
A.
,
Serra-Burriel
 
M.
,
Campillo-Artero
 
C.
.
2018
. “
It’s about Time: Cesarean Sections and Neonatal Health
.”.
Journal of Health Economics
.
59
:
46
59
.

Costa-Ramón
 
A.
,
Kortelainen
 
M.
,
Rodríguez-González
 
A.
,
Sääksvuori
 
L.
.
2022
. “
The Long-Run Effects of Cesarean Sections
Journal of Human Resources
.
57
(
6
):
2048
85
.

Currie
 
J.
,
MacLeod
 
W. B.
.
2008
. “
First Do No Harm? Tort Reform and Birth Outcomes
Quarterly Journal of Economics
.
123
(
2
):
795
830
.

Currie
 
J.
,
MacLeod
 
W. B.
.
2017
. “
Diagnosing Expertise: Human Capital, Decision Making, and Performance among Physicians
Journal of Labor Economics
.
35
(
1
):
1
43
.

de Elejalde
 
R.
,
Giolito
 
E.
.
2021
. “
A Demand-Smoothing Incentive for Cesarean Deliveries
.”.
Journal of Health Economics
.
75
:
102411
.

De Oliveira
 
V. H.
,
Lee
 
I.
,
Quintana-Domeque
 
C.
.
2022
. “
The Effect of Increasing Women’s Autonomy on Primary and Repeated Caesarean Sections in Brazil
.”.
Health Economics
.
31
(
8
):
1800
4
.

Doyle
 
J. J.
,
Graves
 
J. A.
,
Gruber
 
J.
,
Kleiner
 
S. A.
.
2015
. “
Measuring Returns to Hospital Care: Evidence from Ambulance Referral Patterns
Journal of Political Economy
.
123
(
1
):
170
214
.

Fyfe
 
E. M.
,
Anderson
 
N. H.
,
North
 
R. A.
,
Chan
 
E. H.
,
Taylor
 
R. S.
,
Dekker
 
G. A.
,
Screening for Pregnancy Endpoints (SCOPE)Consortium
.
2011
. “
Risk of First-Stage and Second-Stage Cesarean Delivery by Maternal Body Mass Index among Nulliparous Women in Labor at Term
Obstetrics & Gynecology
.
117
(
6
):
1315
22
.

Gallardo
 
J. M.
,
Gómez-López
 
J.
,
Medina-Bravo
 
P.
,
Juárez-Sánchez
 
F.
,
Contreras-Ramos
 
A.
,
Galicia-Esquivel
 
M.
,
Sánchez-Urbina
 
R.
 et al.  
2015
. “
Maternal Obesity Increases Oxidative Stress in the Newborn
Obesity
.
23
(
8
):
1650
54
.

Greene
 
W.
 
2001
.
“Estimating Econometric Models with Fixed Effects.”
.
Technical report
,
New York University
,
Leonard N. Stern School of Business
.

Guendelman
 
S.
,
Gemmill
 
A.
,
Thornton
 
D.
,
Walker
 
D.
,
Harvey
 
M.
,
Walsh
 
J.
,
Perez-Cuevas
 
R.
.
2017
. “
Prevalence, Disparities, and Determinants of Primary Cesarean Births among First-Time Mothers in Mexico
Health Affairs
.
36
(
4
):
714
22
.

Halla
 
M.
,
Mayr
 
H.
,
Pruckner
 
G. J.
,
García-Gómez
 
P.
.
2020
. “
Cutting Fertility? Effects of Cesarean Deliveries on Subsequent Fertility and Maternal Labor Supply
.”.
Journal of Health Economics
.
72
:
102325
.

Heckman
 
J. J.
 
1981
.
The Incidental Parameters Problem and the Problem of Initial Conditions in Estimating a Discrete Time-Discrete Data Stochastic Process
.
Cambridge
:
MIT Press
.

IMSS
 
2013
. “
Parto después de una cesárea: Evidencias y recomendaciones
.”.
Catálogo Maestro de Guías de Práctica Clínica Mexico

IMSS
 
2017
. “
Control prenatal con atención centrada en la paciente.” Guía de Evidencias y Recomendaciones: Guía de Práctica Clínica
.
CENETEC
 
Mexico
.

Jachetta
 
C.
 
2014
. “
Cesarean Sections and Later Child Health Outcomes
.”.

Jensen
 
V. M.
,
Wust
 
M.
.
2015
. “
Can Cesarean Section Improve Child and Maternal Health? The Case of Breech Babies
.”.
Journal of Health Economics
.
39
:
289
302
.

Kearns
 
K.
,
Dee
 
A.
,
Fitzgerald
 
A. P.
,
Doherty
 
E.
,
Perry
 
I. J.
.
2014
. “
Chronic Disease Burden Associated with Overweight and Obesity in Ireland: The Effects of a Small BMI Reduction at Population Level
BMC Public Health
.
14
(
1
):
143
.

Khashan
 
A. S.
,
Kenny
 
L. C.
.
2009
. “
The Effects of Maternal Body Mass Index on Pregnancy Outcome
European Journal of Epidemiology
.
24
(
11
):
697
.

Kominiarek
 
M. A.
,
VanVeldhuisen
 
P.
,
Hibbard
 
J.
,
Landy
 
H.
,
Haberman
 
S.
,
Learman
 
L.
,
Wilkins
 
I.
,
al
 
et
.
2010
. “
The Maternal Body Mass Index: A Strong Association with Delivery Route
American Journal of Obstetrics and Gynecology
.
203
(
3
):
264
e1
.

Kozhimannil
 
K. B.
,
Arcaya
 
M. C.
,
Subramanian
 
S. V.
.
2014
. “
Maternal Clinical Diagnoses and Hospital Variation in the Risk of Cesarean Delivery: Analyses of a National US Hospital Discharge Database
PLOS Medicine
.
11
(
10
):
e1001745
.

Lagarde
 
M.
,
Blaauw
 
D.
.
2022
. “
Overtreatment and Benevolent Provider Moral Hazard: Evidence from South African Doctors
.”.
Journal of Development Economics
.
158
:
102917
.

McDonald
 
S. D.
,
Han
 
Z.
,
Mulla
 
S.
,
Beyene
 
J.
.
2010
. “
Overweight and Obesity in Mothers and Risk of Preterm Birth and Low Birth Weight Infants: Systematic Review and Meta-Analyses
.”.
BMJ
.
341
.

Melo
 
C.
,
Menezes-Filho
 
N.
.
2023
. “
The Effects of a National Policy to Reduce C-Sections in Brazil
Health Economics
.
32
(
2
):
501
17
.

Popkin
 
B. M.
,
Corvalan
 
C.
,
Grummer-Strawn
 
L. M.
.
2020
. “
Dynamics of the Double Burden of Malnutrition and the Changing Nutrition Reality
Lancet
.
395
(
10217
):
65
74
.

Rodríguez-Aguilar
 
R
.
2018
. “
Maternal Mortality in Mexico, beyond Millennial Development Objectives: An Age-Period-Cohort Model
PloS one
.
13
(
3
):
e0194607
.

Wollschlaeger
 
K.
,
Nieder
 
J.
,
Köppe
 
I.
,
Härtlein
 
K.
.
1999
. “
A Study of Fetal Macrosomia
Archives of Gynecology and Obstetrics
.
263
(
1
):
51
55
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]