# Variation between hospital caesarean delivery rates when Robson's classification is considered: An observational study from a French perinatal network

Thibault Quibel, Patrick Rozenberg, Camille Bouyer, Jean Bouyer

► To cite this version:

Thibault Quibel, Patrick Rozenberg, Camille Bouyer, Jean Bouyer. Variation between hospital cae-sarean delivery rates when Robson's classification is considered: An observational study from a French perinatal network. PLoS ONE, 2021, 16 (8), e0251141. doi:10.1371/journal.pone.0251141

HAL Id: hal-03401827

https://hal.science/hal-03401827

Submitted on 31 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Variation between hospital caesarean delivery rates when Robson's classification is considered: An observational study from a French perinatal network

Thibaud Quibel<sup>1,2*</sup>, Patrick Rozenberg<sup>1,2</sup>, Camille Bouyer<sup>1,3</sup>, Jean Bouyer<sup>1,3</sup>

1 Department of Gynecology and Obstetrics, Intercommunal Hospital Centre of Poissy-Saint-Germain-en-Laye, Poissy, France, 2 EA 7285, Research Unit "Risk and Safety in Clinical Medicine for Women and Perinatal Health", Versailles-Saint-Quentin University (UVSQ), Montigny-le-Bretonneux, France, 3 Réseau de Périnatalité Maternité en Yvelines et Périnatalité Active (MYPA), Saint-Germain-en-Laye, France,

4 Université Paris-Saclay, UVSQ, Inserm, CESP, Villejuif, France

* thibaud.quibel@ght-yvelinesnord.fr

<figure><img src="image_2.png" /></figure>

## Abstract

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Quibel T, Rozenberg P, Bouyer C, Bouyer J (2021) Variation between hospital caesarean delivery rates when Robson's classification is considered: An observational study from a French perinatal network. PLoS ONE 16(8): e0251141.

https://doi.org/10.1371/journal.pone.0251141

**Editor:** Alireza Abdollah Shamshirsaz, Baylor

College of Medicine, UNITED STATES

**Received:** November 1, 2020

**Accepted:** April 20, 2021

**Published:** August 20, 2021

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0251141

**Copyright:** © 2021 Quibel et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its S1 File.

**Funding:** The authors received no specific funding for this work.

## Introduction

WHO has recommended using Robson's Ten Group Classification System (TGCS) to monitor and analyze CD rates. Its failure to take some maternal and organizational factors into account, however, could limit the interpretation of CD rate comparisons, because it may contribute to variations in hospital CD rates.

## Objective

To study the contribution of maternal socioeconomic and clinical characteristics and hospital organizational factors to the variation in CD rates when using Robson's ten-group classification system for CD rate comparisons.

## Methods

This prospective, observational, population-based study included all deliveries at a gestational age > 24 weeks at the 10 hospitals of the French MYPA perinatal network in the Paris area. CD rates were calculated for each TGCS group in each hospital. Interhospital variations in these rates were investigated with hierarchical logistic regression models to quantify the variation explained by differences in patient and hospital characteristics when the TGCS is considered. Variations in CD rates between hospitals were estimated with median odds ratios (MOR) to express interhospital variance on the standard odds ratio scale. The percentage of variation explained by TGCS and maternal and hospital characteristics was also calculated.
**Competing interests:** The authors have declared that no competing interests exist.

## Results

The global CD rate was 24.0% (interhospital range: 17–32%). CD rates within each TGCS group differed significantly between hospitals (*P*<0.001). CD was significantly associated with maternal age (>40 years), severe preeclampsia, and two organizational factors: hospital status (private maternities) and the deliveries per staff member per 24 hours. The MOR in the empty model was 1.27 and did not change after taking the TGCS into account. Adding maternal characteristics and hospital organizational factors lowered the MOR to 1.14 and reduced the variation between hospital CD rates by 70%.

## Conclusion

Maternal characteristics and hospital factors are needed to address variation in CD rates among the TGCS groups. Therefore, comparisons of these rates that do not consider these factors should be interpreted carefully.

## Introduction

The cesarean delivery (CD) has become the most commonly performed operation around the world [1]. Despite the 1985 recommendation of the World Health Organization (WHO) that "[t]here is no justification for any region to have a cesarean delivery rate higher than 10–15%" [2], CD rates have continued to increase worldwide and have led to major and controversial public health problems [3, 4]. Although CDs remain essential for some obstetric conditions (for example, placenta previa), many studies have highlighted harmful consequences of its rise, especially the notably higher risks for mothers and babies in subsequent pregnancies, such as abnormally invasive placenta, placental abruption, and stillbirth [5–8]. Consistent with the lack of clear evidence about the optimal mode of delivery, CD rates vary considerably between institutions, regions, and countries [9–12]. These differences may be due to many factors including variations in patient characteristics or preferences, access to care, clinician behavior, and hospital culture or policy [13–17].

Several classification systems have been devised to analyze and try to understand discrepancies in CD rates between hospitals [18–20]. In 2011, a systematic review of 27 CD classification systems identified the 10-group classification system (TGCS) proposed by Robson in 2001 as the most appropriate for classifying women according to their CD risk [21]. The TGCS classifies all deliveries into one of 10 groups based on 6 variables: parity, history of CD, mode of onset of labor, fetal presentation, multiple pregnancy, and gestational age. These 10 categories are mutually exclusive, totally inclusive and can be applied prospectively, since each woman admitted for delivery can be classified immediately according to a few routinely recorded variables [22]. In 2015, WHO recommended that healthcare facilities use the TGCS as a common starting point for comparing and analyzing CD rates (between hospitals), in the expectation that this classification would help healthcare providers to develop interventions to optimize CD rates. Indeed, its use implicitly suggests that women in each group have the same risk profile for CDs. However, the downside of the simplicity of this classification is that it does not take into account many other risk factors related to maternal characteristics, including but not limited to maternal age, obesity, diseases such as hypertensive disorders, and hospital organizational factors, all of which can also contribute to variations in CD rates. Therefore, the use of the single TGCS classification may lead to inappropriate conclusions in the comparison of CD rates between healthcare facilities.
The aim of this study was to examine the contribution of maternal socioeconomic and clinical characteristics and hospital organizational factors to the variation in CD rates when considering Robson's ten-group classification system.

## Material and methods

This study used data from a population-based cohort of births which occurred between the first January 2014 and the 31st December 2014 in the MYPA (Maternité en Yvelines et Péri-natalité Active) perinatal network, which covers the district of Yvelines (west of Paris). This network, which manages around 18,000 deliveries per year, is composed of 10 hospitals (referred to as hospitals A to J), including one hospital (A) affiliated with a university (Université Versailles-Saint Quentin). Half of the hospitals were public (A, B, C, E, H) and half private (D, F, G, I, J). Only two hospitals (A and B), both of which include a neonatal intensive care unit (NICU), manage all types of pregnancies. The other eight (C to J) are level-1 departments intended to provide services only for low-risk pregnancies. The number of deliveries/years differed notably between hospitals: units A, B, and C each had more than 2,000 deliveries/year; units D, E, F, G, and H had between 1,000 and 2,000; and units I and J had fewer than 1,000.

Data are extracted from the CoNaissance 78 program, which was created in 2008 to monitor maternal and perinatal morbidity and mortality in the perinatal network. This dataset contains all births in the district with fully anonymized demographic characteristics and medical information about each pregnancy and delivery, as well as about maternal, fetal, and neonatal health. Therefore, any of the authors had access to data with identifying information prior to de-identification in the database. Since its creation, data are continuously recorded from two health certificates completed for each birth in the network at the hospital of delivery:

* the "first health certificate" of infants born in network hospitals, which is completed during a medical examination that is compulsory nationwide, to be performed within eight days after birth, usually in the maternity ward;

* an additional local health certificate, specifically developed within the MYPA network and reporting additional anonymized data including severe maternal morbidity, social characteristics (e.g., lack of health insurance, single parenthood), and data for all fetal deaths and medical terminations of pregnancy at and after 22 weeks of gestation.

These certificates are completed by midwives and physicians in each hospital and then collected and recorded. Two research midwives provide assistance in collecting the missing data for the CoNaissance 78 program database and in controlling its quality to improve the completeness of the certificates. The proportion of missing data is < 3% for all variables in this study [23]. The National Committee for Data Protection (Commission Nationale de l'Informaticque et des Libertés, registration number 1295794) approved the study, which was conducted in accordance with French legislation. Under French law, this study is exempt from informed consent requirements because patients received standard care and because the dataset contains no information that enables patient identification. Similarly, ethics committee approval was not required because the study used an anonymized database and did not modify patient care.

The maternal risk profile included the following variables: maternal age, recorded and coded in two categories (<40 years yes/no), educational level (coded in three categories: university/high school/primary and middle school), occupational activity (yes/no), and severe preeclampsia, defined according to the ACOG criteria in the perinatal network [24]. Hospital organizational factors included hospital status (public or private), the presence of a neonatal intensive care unit (yes/no), and the ratio of deliveries to staff members (midwives and obstetricians) per 24 hours. This last variable was developed to create a quantitative measure for
appropriate staffing, in view of the lack of any standard for obstetric care in labor wards between countries, given that deliveries are managed by either midwives or obstetricians. In our analysis, we took the median of deliveries per staff member per 24 hours among MYPA hospitals for a reference.

For each TGCS group, we calculated its specific cesarean rate (that is, the number of women undergoing a CD divided by the total number of women in the group), and the relative size of the group (as a percentage of the total population of women within each hospital). We also calculated each group's contribution to the overall CD rate, also as a percentage, that is, the number of women with a CD in the group divided by the total number of CDs.

The variation in hospital CD rates was analyzed with a multilevel logistic regression model with a random intercept to take into account that women were nested within hospitals [25–27]. This model enabled us to estimate the area level variance to quantify the variation in CD rates between hospitals [28]. We then computed the median odds ratio (MOR) to express hospital-level variance according to the standard odds ratio (OR) scale, which can be interpreted consistently and intuitively. The MOR is defined as the median value of the odds ratio between the hospital at highest risk of CD and the hospital at lowest risk, with random selection of two hospitals. MOR can be conceptualized as the increased median risk that a patient would have if moved to another hospital with a higher risk. In this study, the MOR showed the extent to which the individual probability of a CD was determined by the hospital's level of care and was therefore appropriate for quantifying contextual phenomena [25]. If the MOR is equal to one, the probability of a CD does not differ between the hospitals.

Models were fitted with the three-stage approach described by Merlo et al. [25]. First, a null (or empty) model with hospital random intercepts only was fit to estimate the unadjusted hospital-level variation in CD rates. Then we added the TGCS group as a variable (that is, the group in which each woman is classified) in the model to investigate the extent to which the variation in CD rates between hospitals was modified when the TGCS classification was considered. Finally, we added the variables describing the women's clinical, social, and demographic characteristics as well as the hospitals' characteristics, to determine the amount of variation associated with these variables and to test the significance of the residual between-hospital variance. Maternal and hospital variables whose crude associations with CD occurrence had a p-value < .2 were selected for the multilevel model.

As more than 25% of deliveries took place in a single hospital (A, 4,530/17,433), we performed a sensitivity analysis by excluding it. Variation in CD rates among the nine remaining hospitals was again analyzed with the same process.

All analyses were performed with R Studio version 1.0.136. The analysis was conducted by applying a significance level of α = 0.05.

## Results

During 2014, 17,511 deliveries took place in the 10 maternity units of the perinatal network. After excluding 78 women (0.4%) due to missing data for one of the 6 TGCS variables, 17,433 women remained in the final sample. Two thirds of the deliveries occurred in a public hospital (66.0%, 11,512/17,433), and 40.0% (6,999/17,433) in the hospitals with a NICU. There were 4,182 CDs during this year (overall rate: 24.0%, interhospital range: 17.2%-32.6%). Fig 1 shows the CD rate for each hospital, and Table 1 summarizes the characteristics of the women and the hospitals by TGCS group.

Most women were considered at low risk, i.e., classified in TGCS groups 1 and 3 (56.2%; 9,804/17,433), including nulliparous and parous women with no previous CD, with a singleton in cephalic presentation ≥ 37 weeks, with spontaneous labor. Groups 2 and 4 (4,003/17,433;
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Cesarean delivery (CD) in each hospital of the MYPA perinatal network. The CD rates were significantly different between each hospital (p<0.001). Hospital A and B are second level maternities.</figcaption></figure>

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

nulliparas and paras with no previous CD, with a singleton in cephalic presentation ≥ 37 weeks, with induction of labor or a CD before labor) contained 22.9% of the women, while 10.3% were in group 5 (1,795/17,433; previous CD, singleton in cephalic presentation ≥ 37 weeks). Around 10% of the population was assigned to groups 6 to 10. (Table 2) The CD rates for each group differed significantly between the hospitals (P < .001). Fig 2 illustrates the variations in CD rates between hospitals within each TGCS group.

The crude associations of CD occurrence with maternal characteristics and hospital organizational factors are reported in Table 3. CD was significantly associated with a maternal age > 40 years (OR 1.77; 95% CI 1.54–2.01) and with severe preeclampsia (OR 4.2, 95% CI

<table><caption>**Table 1.** Maternal characteristics and hospital factors for each group in the Ten Group Classification System (TGCS).</caption><thead><tr><th>Maternal characteristics</th><th>Overall</th><th>Group 1</th><th>Group 2</th><th>Group 3</th><th>Group 4</th><th>Group 5</th><th>Group 6</th><th>Group 7</th><th>Group 8</th><th>Group 9</th><th>Group 10</th></tr></thead><tbody><tr><td>Maternal age (years): median (IQR)</td><td>31.4 (28–35)</td><td>28.8 (25.9–31.9)</td><td>29.4 (26.1–32.9)</td><td>32.5 (29.4–35.5)</td><td>33.4 (30.3–36.7)</td><td>33.4 (30.2–36.4)</td><td>30.2 (27.4–33.1)</td><td>34.0 (31.0–36.9)</td><td>31.8 (29.0–35.8)</td><td>35.1 (29.3–38.0)</td><td>31.3 (27.8–35.1)</td></tr><tr><td>age >40 years</td><td>3.9%</td><td>1.6%</td><td>3.4%</td><td>5.8%</td><td>8.6%</td><td>7.3%</td><td>3%</td><td>9.2%</td><td>10.4%</td><td>14.5%</td><td>6.7%</td></tr><tr><td>Woman's level of education</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>University (%)</td><td>61.4%</td><td>65.5%</td><td>62.7%</td><td>58.6%</td><td>57.2%</td><td>56.2%</td><td>72.6%</td><td>51.3%</td><td>66.0%</td><td>72.9%</td><td>53.2%</td></tr><tr><td>High school (%)</td><td>23.2%</td><td>22.5%</td><td>24.4%</td><td>23.3%</td><td>24.5%</td><td>24.3%</td><td>17.6%</td><td>27.0%</td><td>20.6%</td><td>18.7%</td><td>24.7%</td></tr><tr><td>Primary and middle school (%)</td><td>15.4%</td><td>12%</td><td>12.9%</td><td>18.1%</td><td>18.3%</td><td>19.5%</td><td>9.8%</td><td>21.7%</td><td>13.4%</td><td>8.3%</td><td>22.1%</td></tr><tr><td>Mother's occupational activity</td><td>69.7%</td><td>75.6%</td><td>73.4%</td><td>65.4%</td><td>65.2%</td><td>65.3%</td><td>81.8%</td><td>68.3%</td><td>65.4%</td><td>75.0%</td><td>68.2%</td></tr><tr><td>Severe preeclampsia</td><td>1.1%</td><td>0.2%</td><td>1.5%</td><td>0.5%</td><td>2.1%</td><td>1.5%</td><td>2.7%</td><td>0.7%</td><td>2.2%</td><td>2.0%</td><td>3.7%</td></tr><tr><td>Hospital factors</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Proportion of deliveries in private hospital (%)</td><td>34%</td><td>33%</td><td>38%</td><td>33%</td><td>39%</td><td>39%</td><td>41%</td><td>41%</td><td>21%</td><td>48%</td><td>24%</td></tr><tr><td>Proportion of deliveries in hospitals with a NICU (%)</td><td>40%</td><td>39.5%</td><td>34.4%</td><td>39.2%</td><td>36.7%</td><td>38.3%</td><td>33.1%</td><td>35.2%</td><td>65.3%</td><td>42.6%</td><td>59%</td></tr></tbody></table>

NICU: Neonatal intensive care unit.

Qualitative data are given with %.

Quantitative data are given with median and interquartile (IQR).

https://doi.org/10.1371/journal.pone.0251141.t001
<table><caption>Table 2. The relative size, the cesarean delivery rate, and the contribution of each group in Robson's classification to this rate in the MYPA network.</caption><thead><tr><th></th><th>Number of CDs</th><th>Number of patients</th><th>Group size (%)<sup>1</sup></th><th>Group CD rate (%)<sup>2</sup></th><th>Relative group Contribution to overall CD rate (%)<sup>3</sup></th></tr></thead><tbody><tr><td>Group 1: Nulliparous, single cephalic, ≥37 weeks, in spontaneous labor</td><td>482</td><td>4300</td><td>24.7%</td><td>11.2%</td><td>11.5%</td></tr><tr><td>Group 2: Nulliparous, singleton cephalic, ≥37 weeks, with elective delivery</td><td>752</td><td>1880</td><td>10.8%</td><td>40.0%</td><td>18.0%</td></tr><tr><td>Group 3: Multiparous, singleton cephalic, ≥37 weeks, in spontaneous labor</td><td>173</td><td>5504</td><td>31.6%</td><td>3.1%</td><td>4.1%</td></tr><tr><td>Group 4: Multiparous, singleton cephalic, ≥37 weeks, with elective delivery</td><td>453</td><td>2123</td><td>12.2%</td><td>21.3%</td><td>10.8%</td></tr><tr><td>Group 5: Multiparous, singleton cephalic, ≥37 weeks, with history of CD</td><td>1203</td><td>1795</td><td>10.3%</td><td>67.0%</td><td>28.8%</td></tr><tr><td>Group 6: Nulliparous, singleton, breech,</td><td>281</td><td>293</td><td>1.7%</td><td>95.9%</td><td>6.7%</td></tr><tr><td>Group 7: Multiparous, single breech,</td><td>157</td><td>196</td><td>1.1%</td><td>81.1%</td><td>3.8%</td></tr><tr><td>Group 8: Multiple pregnancies</td><td>313</td><td>467</td><td>2.7%</td><td>67%</td><td>7.5%</td></tr><tr><td>Group 9: Transverse presentation</td><td>46</td><td>49</td><td>0.2%</td><td>93.8%</td><td>1.1%</td></tr><tr><td>Group 10: Singleton, cephalic presentation &lt;37 weeks</td><td>322</td><td>825</td><td>4.7%</td><td>39.0%</td><td>7.7%</td></tr></tbody></table>

1. % = n of women in the group/total N women delivered in the setting x 100.

2. % = n of CS in the group/total N of women in the group x 100.

3. % = n of CS in the group/total N of CS in the setting x 100.

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

3.01–5.7). Hospital status and the ratio of deliveries per staff member per 24 hours were also associated with CD rates.

Table 4 presents the three stages of the analysis with a multilevel logistic regression model. In the empty model, the hospital-level variance was 0.065 and the MOR 1.27. When the TGCS was included in the model, the variance and the MOR remained similar (variance 0.077, MOR 1.33). However, in the final model, which also included maternal and hospital factors, the interhospital variance decreased substantially by more than 70% (variance 0.019), and the MOR declined to 1.14. In this final model, CD rates were significantly associated with maternal age >40 years (ORa 1.05; 95% CI 1.03–1.06), educational level (ORa 1.20; 95% CI 1.13–1.27), severe preeclampsia (ORa 2.20, 95% CI 2.07–2.33), and private maternity units (ORa 1.72; 95% CI 1.35–2.18). The lack of an NICU, however, was associated with a significantly lower CD rate (ORa 0.74; 95% CI 0.56–0.99).

After exclusion of births from hospital A, the analysis of the remaining 12,933 births produced similar results and thus confirmed that most of the variation in CD rates among hospitals was related to maternal characteristics and hospital factors.

# Discussion

## Principal findings

In this study, we considered a population in which CD rates varied between hospitals. We developed regression models to quantify the extent to which variations in CD rates between maternity units remained when TGCS was taken into account and how the consideration of women's socioeconomic and clinical characteristics and hospital organization factors might modify this variation. Our analyses showed that patient and hospital factors reduced the inter-hospital variance in CD rates by more than 70%, thereby demonstrating the need to integrate case-mix models when we compare CD rates within groups of the TGCS between hospitals.
<figure><img src="image_5.png" /><figcaption>**Fig 2. Cesarean delivery (CD) rate of each hospital.** For each histogram, the first column represents the mean CD rate of the Mypa network, followed by each CD rate of each hospital (A to J). CD rates the Group 9 (fetus in a transverse presentation) are not shown, as the relative size of Group 9 is 0,2% of overall population.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0251141.g002

## Implications of the work

As the difficulties in understanding the rise in CD rates worldwide become more generally recognized, the TGCS has emerged as the most appropriate—and the most frequently used—classification for assessing CD rates, as it builds clinically meaningful subgroups of the population of pregnant women and thus makes it possible to analyze variations in CD rates across institutions, countries, and time. Many classification systems are based on indications for CD, such as fetal distress, dystocia, failure to progress, cephalopelvic disproportion, failed induction, macrosomia, or failed trial of labor, but these classifications are limited by ambiguous terminology [21]. Over the last decade, therefore, the TGCS has been adopted for monitoring (and even comparing) CD rates, with the use of specific indicators such as the CD rate within each TGCS group and the relative size of the group (percentage of the total population of women).
<table><caption>Table 3. Crude odds ratios for cesarean deliveries by maternal socioeconomic and clinical characteristics and hospital organizational factors.</caption><thead><tr><th></th><th>OR</th><th>95% CI</th><th>P-value</th></tr></thead><tbody><tr><td><strong>Maternal characteristics</strong></td><td></td><td></td><td></td></tr><tr><td>Maternal age > 40 years</td><td>1.77</td><td>1.54-2.01</td><td>&lt;0.001</td></tr><tr><td><strong>Educational level</strong></td><td></td><td></td><td>0.11</td></tr><tr><td>University</td><td>ref</td><td>ref</td><td></td></tr><tr><td>High school</td><td>1.01</td><td>0.99-1.08</td><td></td></tr><tr><td>Primary and middle school</td><td>1.04</td><td>1.01-1.07</td><td></td></tr><tr><td><strong>Occupational activity</strong></td><td>0.97</td><td>0.90-1.04</td><td>0.47</td></tr><tr><td><strong>Severe preeclampsia</strong></td><td>4.16</td><td>3.01-5.7</td><td>&lt;0.001</td></tr><tr><td><strong>Hospital factor</strong></td><td></td><td></td><td></td></tr><tr><td>Private status</td><td>1.37</td><td>1.27-1.47</td><td>&lt;0.001</td></tr><tr><td>No NICU unit</td><td>0.93</td><td>0.87-1.04</td><td>0.06</td></tr><tr><td>&gt; 0.92 deliveries per staff member per 24h*</td><td>0.89</td><td>0.83-0.95</td><td>0.002</td></tr></tbody></table>

NICU: Neonatal intensive care unit.

Data are expressed in numbers, percentages, and odds ratios with their 95% confidence intervals

*0.92 deliveries/staff member/24 hours corresponds to the median of deliveries/staff member among MYPA hospitals.

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

This classification illustrates the variation in CD risk from one group to another and helps to identify groups of women at low and at high risk of CD.

However, while the TCGS is a starting point for describing CD practices in relevant obstetric groups, its utility in accounting for the wide variance of these practices between hospitals is limited. Moreover, use of the TGCS as an audit tool can lead to misinterpretations of CD rate variations between hospitals. Most studies evaluating the TGCS do not describe maternal characteristics, or the aspects of hospital organization that may influence CD rates [29-31]. For

<table><caption>Table 4. Measures of associations between individual/hospital characteristics and cesarean delivery (CD) and measures of variability of cesarean delivery rates in the MYPA network, France, 2014, obtained from multilevel logistic models.</caption><thead><tr><th></th><th>Empty model</th><th>Model using TGCS</th><th>Model using TGCS plus maternal characteristics and hospital factors</th></tr></thead><tbody><tr><td><strong>Measures of association with CD occurrence (OR, 95% CI)</strong></td><td></td><td></td><td></td></tr><tr><td>Maternal age >40 years</td><td></td><td>1.05(1.03-1.06)</td><td></td></tr><tr><td>Educational level</td><td></td><td>1.20 (1.13-1.27)</td><td></td></tr><tr><td>Severe preeclampsia</td><td></td><td>2.20 (2.07-2.33)</td><td></td></tr><tr><td>Private maternity</td><td></td><td>1.72 (1.35-2.18)</td><td></td></tr><tr><td>No NICU unit</td><td></td><td>0.74 (0.56-0.99)</td><td></td></tr><tr><td>&gt;0.92 deliveries/staff member/24h*</td><td></td><td>0.86 (0.56-1.16)</td><td></td></tr><tr><td><strong>Measures of variability of CD rates</strong></td><td></td><td></td><td></td></tr><tr><td>Area level variance</td><td>0.065</td><td>0.077</td><td>0.019</td></tr><tr><td>Proportional change in variance</td><td></td><td>-18%</td><td>71%</td></tr><tr><td>MOR</td><td>1.27</td><td>1.30</td><td>1.14</td></tr></tbody></table>

NICU: Neonatal intensive care unit, OR: Odds ratio, 95%CI: 95% confidence interval.

ICC, intraclass correlation; MOR, median odds ratio.

0.92 deliveries/staff member/24 hours corresponds to the median of deliveries/staff member among MYPA hospitals.

https://doi.org/10.1371/journal.pone.0251141.t004
example, we observed that the CD rate was significantly lower in hospitals without NICUs. Similarly, Le Ray et al. analyzed the influence of the maternity unit's perinatal care level on the rate of intrapartum cesarean delivery among women with low-risk pregnancies and reported that maternity units with NICUs that manage high-risk pregnancies have higher rates of cesareans during labor for their population of nulliparas at low risk than do facilities that deal mainly with low-risk pregnancies. They identified the maternity unit level as the only structural factor that was a significant risk factor for cesarean delivery during labor. Le Ray et al. suggested several hypotheses to explain the differences in these cesarean rates as a function of level of perinatal care: Perhaps, in high-risk facilities, physicians are more likely to expect problems and may encourage the use of cesarean deliveries while in low-risk facilities, midwives are less likely to expect problems and may thus encourage vaginal delivery, which is important because support for women during labor is essential for successful vaginal delivery. Another hypothesis centers on the organization of care: if the availability of personnel makes it easier to perform a cesarean delivery in a level-3 unit, this might affect obstetric decision-making [32].

Some authors have suggested that the association of the TGCS with maternal factors and prenatal risk factors may provide the most reliable model for exploring this variation. Studying CD rates in seven international tertiary referral hospitals (as well as a regional referral hospital and the Iceland 2005 birth cohort), Brennan et al. have suggested that the variation of CD rates among nulliparous and parous women in spontaneous labor with singletons in cephalic presentation (TGCS groups 1 and 3) might be related to maternal clinical, social, and demographic risk factors, which the TGCS does not consider [12]. The findings of Colais et al. support this hypothesis, in demonstrating that the residual variability in CDs is explained by clinical, social, and demographic variables and that the TGCS classification does not suffice to eliminate case-mix differences [33]. Pasko et al. also confirmed this point recently, in their estimates of the contributions of the characteristics of patients, healthcare providers, and hospitals to the variation in the frequency of nulliparous, term, singleton, vertex CDs [34]. In a cohort of 115,502 women, they observed that patient characteristics accounted for 24% of the variation in this group, while healthcare provider-hospital characteristics were not significantly associated with CD frequency. These findings, reinforced by ours, indicate that understanding CD rates requires recording not only clinical/obstetric factors, but also factors related to maternity organizational factors (ratio of deliveries/staff, ratio of deliveries/delivery room, number per 24 hours of inhouse attending obstetricians, anesthesiologists, and pediatricians). The third paper of a series on optimizing CD use in the *Lancet* illustrated the components affecting CD frequency in a schematic diagram in which the three outer rings represent the layers of complexity of these factors [35]. We agree that a classification requiring such data is time-consuming and complex and that in the end we still need to compare CD rates. The TGCS is an appropriate classification for highlighting differences in these rates, but we have demonstrated why comparisons of CD rates, even in TGCS groups, must be interpreted carefully.

## Strengths and limitations

The main strength of our study lies in its large population-based sample size, multicenter nature, and prospective design, including the collection of numerous data items that make it possible to adjust for a large set of confounding factors. The maternal clinical, social, and demographic data were prospectively recorded in the database by a research midwife and enabled an exclusion rate for missing data of less than 1%.

Our study used case-mix models, which provide a relevant statistical approach for comparing CD rates across hospitals and are suggested for assessing CD rates as institutional quality
indicators. These models help us to understand the determinants of the variability of health-
care quality indicators such as CD rates and thus to understand their crude rates and risk ratios
better. Precisely because outcomes may depend on preexisting patient characteristics, simply
measuring CD rates may not provide useful insight into quality of care. These models have
been widely employed in such clinical disciplines as cardiothoracic surgery, but not yet for
evaluating obstetric outcomes [36, 37]. Nonetheless Main et al. have suggested that controlling
for maternal age reduces CD rate variability in nulliparous term cephalic singleton pregnancies
[38].

We computed the median odds ratio (MOR), which is an epidemiologically more suitable option for obtaining measures of variance in logistic regression. It is not statistically dependent on the prevalence of the outcome and allows area-level variance to be expressed on the well-known OR scale. Therefore, it allows comparison of area-level variations (here the hospital level) with the impact of specific factors [28].

Limitations of this study should be noted. Even though this multilevel model may be gener-
alizable to other studies, the variation of CD rates associated with maternal and hospital factors
in our study should be interpreted with caution, as these factors might vary across populations.
Furthermore, more than 25% of deliveries occurred in a single hospital, which was the only
university hospital of the perinatal network. However, findings remained unchanged after its
exclusion in a sensitivity analysis.

We had very little information on other maternal clinical characteristics, specifically status
for diabetes mellitus or fetal growth restriction (as we did not have access to prenatal ultra-
sound findings), which are also important variables to consider in interpreting TGCS results.
The CoNaissance program does not collect the ultrasound findings, and although we had the
birth weight, we were not able to determine if the fetus was considered either small-for-gesta-
tional age or a growth-restricted fetus. The variable for diabetes was missing in more than 50%
of cases, and we were not able to differentiate gestational diabetes mellitus treated with or with-
out insulin from one another, or from preexisting (non-pregnancy-related) diabetes [39].
Moreover, there is no evidence that gestational diabetes is associated with a higher risk of CD
[40]. Finally, this study does not consider maternal body mass index (BMI) in the risk adjust-
ment, as unfortunately it was not recorded in birth certificates. The well-known association of
obesity with higher CD risk means that it might well be a factor of variability in CD rates
between hospitals [41]. We suppose that adding these maternal and obstetric risk factors to the
final model (which uses maternal characteristics, hospital factors, and TGCS) would probably
decrease the area-level variance of the model and therefore decrease a little bit more the
amount of variation explained by each hospital.

Conclusions

The Ten Group Classification System is an easy tool to implement for auditing CD and for highlighting differences in CD rates among maternity units. However, maternal characteristics and hospital organizational factors are needed before comparing CD rates within groups of the TGCS.

**Supporting information**

S1 File.
(XLTX)
## Author Contributions

**Conceptualization:** Thibaud Quibel.

**Data curation:** Camille Bouyer.

**Formal analysis:** Thibaud Quibel.

**Investigation:** Thibaud Quibel.

**Methodology:** Thibaud Quibel, Patrick Rozenberg, Jean Bouyer.

**Writing – review & editing:** Patrick Rozenberg, Jean Bouyer.

## References

1. Masdenet al., Acta Obstet Gynecol Scand 2013; 92:256-63.

2. WHO (1985) Appropriate technology for birth. Lancet 2: 436-437. PMID: 2863457

3. Gibbons L, Belizan JM, Lauer J, Betran AP, Merialdi M, Althabe F. (2010) The global numbers and costs of additionally needed and unnecessary caesarean sections performed per year: overuse as a barrier to universal coverage. World Health Report. Geneva, Switzerland: World Health Organization. Background Paper, 30.

4. Betran AP, Torloni MR, Zhang JJ, Gülmezoglu AM. WHO Working Group on Caesarean Section. WHO Statement on Caesarean Section Rates. BJOG. 2016; 123:667-70. https://doi.org/10.1111/1471-0528.13526 PMID: 26681211

5. Keag OE, Norman JE, Stock JE. Long-term risks and benefits associated with cesarean delivery for mother, baby, and subsequent pregnancies: Systematic review and meta-analysis. PLoS Med. 2018; 15: e1002494 https://doi.org/10.1371/journal.pmed.1002494 PMID: 29360829

6. Moraitis AA, Oliver-Williams C, Wood AM, Fleming M, Pell JP, Smith GCS. Previous caesarean delivery and the risk of unexplained stillbirth: retrospective cohort study and meta-analysis. BJOG. 2015; 122:1467-74 https://doi.org/10.1111/1471-0528.13461 PMID: 26033155

7. Esteves-Pereira AP, Deneux-Tharaux C, Nakamura-Pereira M, Saucedo M, Bouvier-Colle MH, Leal Mdo C. Caesarean Delivery and Postpartum Maternal Mortality: A Population-Based Case Control Study in Brazil. PLoS One. 2016; 11: e0153396 https://doi.org/10.1371/journal.pone.0153396 PMID: 27073870

8. Downes KL HS, Sjaarda LA, et al. Previous prelabor or intrapartum cesarean delivery and risk of placenta previa. Am J Obstet Gynecol. 2015; 212:669 e1-6.

9. Organization for Economic Cooperation and Development. OECD (2018), Caesarean sections (indicator). https://doi.org/10.1787/adc3c39f-en (Accessed on 06 May 2018)

10. Vogel JP, Betrán AP, Vindevoghel N, Souza JP, Torloni MR, Zhang J, et al.; WHO Multi-Country Survey on Maternal and Newborn Health Research Network. Use of the Robson classification to assess caesarean section trends in 21 countries: a secondary analysis of two WHO multicountry surveys. Lancet Glob Health. 2015; 3:e260-70. https://doi.org/10.1016/S2214-109X(15)70094-X PMID: 25866355

11. Betran AP, Ye J, Moller AB, Zhang J, Gulmezoglu AM, Torloni MR. The increasing trend in caesarean section rates: Global, regional and national estimates: 1990±2014. PLoS ONE. 2016; 11(2) https://doi.org/10.1371/journal.pone.0148343 PMID: 26849801

12. Brennan DJ, Robson MS, Murphy M, O'Herlihy C. Comparative analysis of international cesarean delivery rates using 10-group classification identifies significant variation in spontaneous labor. Am J Obstet Gynecol 2009; 201:308.e1-8. https://doi.org/10.1016/j.ajog.2009.06.021 PMID: 19733283

13. Plough AC, Galvin G, Li Z, Lipsitz SR, Alidina S, Henrich NJ, et al. Relationship Between Labor and Delivery Unit Management Practices and Maternal Outcomes. Obstet Gynecol. 2017; 130:358-365. https://doi.org/10.1097/AOG.0000000000002128 PMID: 28697107

14. Coulm B, Le Ray C, Lelong N, Drewniak N, Zeitlin J. Obstetric interventions for low-risk pregnant women in France: do maternity unit characteristics make a difference? Birth 2012; 39:183-91. https://doi.org/10.1111/j.1523-536X.2012.00547.x PMID: 23281900

15. Paranjothy S, Frost C, Thomas J. How much variation in CS rates can be explained by case mix differences? BJOG 2005; 112:658-66. https://doi.org/10.1111/j.1471-0528.2005.00501.x PMID: 15842293

16. Nakamura-Pereira M, do Carmo Leal M, Esteves-Pereira AP, Domingues RM, Torres JA, Dias MA, et al. Use of Robson classification to assess cesarean section rate in Brazil: the role of source of payment for childbirth. Reprod Health. 2016; 13:128. https://doi.org/10.1186/s12978-016-0228-7 PMID: 27766941
17. Mazzoni A, Althabe F, Liu NH, Bonotti AM, Gibbons L, Sánchez AJ, et al. Women's preference for caesarean section: a systematic review and meta-analysis of observational studies. BJOG. 2011; 118:391-9. https://doi.org/10.1111/j.1471-0528.2010.02793.x PMID: 21134103

18. Thomas J, Callwood A, Brocklehurst P, Walker J. The National Sentinel Caesarean Section Audit. BJOG 2000; 107:579-80. https://doi.org/10.1111/j.1471-0528.2000.tb13296.x PMID: 10826569

19. Notzon FC, Cnattingius S, Bergsjo P, Cole S, Taffel S, Irgens L, et al. Cesarean section delivery in the 1980s: international comparison by indication. Am J Obstet Gynecol 1994; 170:495-504. https://doi.org/10.1016/s0002-9378(94)70217-9 PMID: 8116703

20. Wilkes PT, Wolf DM, Kronbach DW, Kunzel M, Gibbs RS. Risk factors for cesarean delivery at presentation of nulliparous patients in labor. Obstet Gynecol 2003; 102:1352-7. https://doi.org/10.1016/j.obstetgynecol.2003.08.006 PMID: 14662226

21. Torloni MR, Betran AP, Souza JP, Widmer M, Allen T, Gulmezoglu M, et al. Classifications for cesarean section: a systematic review. PLoS One. 2011; 6:e14566. https://doi.org/10.1371/journal.pone.0014566 PMID: 21283801

22. Robson M. Classification of caesarean sections. Fetal Matern Med Rev 2001; 12:23-39.

23. Lebreton E, Rozenberg P, Chalavoux K, Huart F, Cotte B, Pineau C, et al. Evaluation of a perinatal network using the first certificates of health. Journal de Gynécologie Obstétrique et Biologie de la Reproduction 2014; 43:342-350 https://doi.org/10.1016/j.jgyn.2013.02.009 PMID: 23608346

24. ACOG Committee on Practice Bulletins—Obstetrics. ACOG practice bulletin. Diagnosis and management of preeclampsia and eclampsia. Number 33, January 2002. Obstet Gynecol. 2002; 99:159-67. https://doi.org/10.1016/s0029-7844(01)01747-1 PMID: 16175681

25. Merlo J. Multilevel analytical approaches in social epidemiology: measures of health variation compared with traditional measures of association. J Epidemiol Community Health 2003; 57:550-2. https://doi.org/10.1136/jech.57.8.550 PMID: 12883048

26. Larsen K, Merlo J. Appropriate assessment of neighborhood effects on individual health—integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol 2005; 161:81-8. https://doi.org/10.1093/aje/kwi017 PMID: 15615918

27. Goldstein H, Browne W, Rasbash J. Partitioning variation in generalised linear multilevel models. Understanding Statistics 2002; 1:223-32.

28. Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, et al. A brief conceptual tutorial of multi-level analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006; 60:290-7. https://doi.org/10.1136/jech.2004.029454 PMID: 16537344

29. Farine D, Shepherd D, Robson M. Classification of caesarean sections in Canada: the modified robson criteria. J Obstet Gynaecol Can. 2012; 34:1133. PMID: 25162083

30. Le Ray C, Blondel B, Prunet C, Khireddine I, Deneux-Tharaux C, Goffinet F. Stabilising the caesarean rate: which target population? BJOG. 2015; 122:690-9. https://doi.org/10.1111/1471-0528.13199 PMID: 25412695

31. Tapia V, Betran AP, Gonzales GF. Caesarean Section in Peru: Analysis of Trends Using the Robson Classification System. PLoS One. 2016; 11:e0148138. https://doi.org/10.1371/journal.pone.0148138 PMID: 26840693

32. Le Ray C, Carayol M, Zeitlin J, Bréart G, Goffinet F; PREMODA Study Group. Level of perinatal care of the maternity unit and rate of cesarean in low-risk nulliparas. Obstet Gynecol. 2006; 107:1269-77. https://doi.org/10.1097/01.AOG.0000218098.70942.a2 PMID: 16738151

33. Colais P, Fantini MP, Fusco D, Carretta E, Stivanello E, Lenzi J, et al. Risk adjustment models for inter-hospital comparison of CS rates using Robson's ten group classification system and other socio-demographic and clinical variables. BMC Pregnancy Childbirth. 2012; 12:54. https://doi.org/10.1186/1471-2393-12-54 PMID: 22720844

34. Pasko DN, McGee P, Grobman WA, Bailit JL, Reddy UM, Wapner RJ, et al. Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) Maternal-Fetal Medicine Units (MFMU) Network. Variation in the Nulliparous, Term, Singleton, Vertex Cesarean Delivery Rate. Obstet Gynecol. 2018; 131:1039-1048. https://doi.org/10.1097/AOG.0000000000002636 PMID: 29742665

35. Betrán AP, Temmerman M, Kingdon C, Mohiddin A, Opiyo N, Torloni MR, et al. Interventions to reduce unnecessary caesarean sections in healthy women and babies. Lancet. 2018 Oct 13; 392 (10155):1358-1368 https://doi.org/10.1016/S0140-6736(18)31927-5 PMID: 30322586

36. Grobman WA, Bailit JL, Rice MM, Wapner RJ, Varner MW, Thorp JM Jr, et al. Eunice Kennedy Shriver National Institute of Child Health and Human Development Maternal-Fetal Medicine Units Network. Can differences in obstetric outcomes be explained by differences in the care provided? The MFMU Network APEX study. Am J Obstet Gynecol. 2014; 211:147.e1-147.e16.
37. Bailit JL, Grobman WA, Rice MM, Spong CY, Wapner RJ, Varner MW, et al. Eunice Kennedy Shriver National Institute of Child Health and Human Development Maternal–Fetal Medicine Units Network. Risk-adjusted models for adverse obstetric outcomes and variation in risk-adjusted outcomes across hospitals. Am J Obstet Gynecol. 2013; 209:446.e1–446.e30.

38. Main EK, Moore D, Farrell B, Schimmel LD, Altman RJ, Abraham C, et al. Is there a useful cesarean birth measure? Assessment of the nulliparous term singleton vertex cesarean birth rate as a tool for obstetric quality improvement. Am J Obstet Gynecol. 2006; 194: 1655–1651

39. Martis R, Crowther CA, Shepherd E, Alsweiler J, Downie MR, Brown J. Treatments for women with gestational diabetes mellitus; an overview of Cochrane systematic reviews. Cochrane Database Syst Rev. 2018 Aug 14; 8(8):CD012327.

40. Crowther CA, Hiller JE, Moss JR, McPhee AJ, Jeffries WS, Robinson JS; Australian Carbohydrate Intolerance Study in Pregnant Women (ACHOIS) Trial Group. Effect of treatment of gestational diabetes mellitus on pregnancy outcomes. N Engl J Med. 2005 Jun 16; 352(24):2477–86. https://doi.org/10.1056/NEJMoa042973 PMID: 15951574

41. Rogers AJG, Harper LM, Mari G. A conceptual framework for the impact of obesity on risk of cesarean delivery. Am J Obstet Gynecol. 2018; 219:356–363. https://doi.org/10.1016/j.ajog.2018.06.006 PMID: 29902446