Abstract
Background: Previous research has suggested that obstetric outcomes are similar for deliveries by family physicians and obstetricians, but many of these studies were small, and none of them adjusted for unmeasured selection bias. We compared obstetric outcomes between these provider types using an econometric method designed to adjust for unobserved confounding.
Methods: We performed a retrospective population-based cohort study of all Canadian (except Quebec) hospital births with delivery by family physicians and obstetricians at more than 20 weeks gestational age, with birth weight greater than 500 g, between Apr. 1, 2006, and Mar. 31, 2009. The primary outcomes were the relative risks of in-hospital perinatal death and a composite of maternal mortality and major morbidity assessed with multivariable logistic regression and instrumental variable–adjusted multivariable regression.
Results: After exclusions, there were 3600 perinatal deaths and 14 394 cases of maternal morbidity among 799 823 infants and 793 053 mothers at 390 hospitals. For deliveries by family physicians v. obstetricians, the relative risk of perinatal mortality was 0.98 (95% confidence interval [CI] 0.85–1.14) and of maternal morbidity was 0.81 (95% CI 0.70–0.94) according to logistic regression. The respective relative risks were 0.97 (95% CI 0.58–1.64) and 1.13 (95% CI 0.65–1.95) according to instrumental variable methods.
Interpretation: After adjusting for both observed and unobserved confounders, we found a similar risk of perinatal mortality and adverse maternal outcome for obstetric deliveries by family physicians and obstetricians. Whether there are differences between these groups for other outcomes remains to be seen.
Over the past several decades in Canada, obstetric deliveries have increasingly been attended by specialist obstetricians rather than family physicians.1 Although specialized care is beneficial for high-risk mothers and their infants,2–4 there are concerns that it might increase risk for women whose deliveries could be safely managed without a specialized approach. Most prior studies have concluded that obstetric outcomes between family physicians and obstetricians are similar, but many of these studies were small, and none of them adjusted for unmeasured factors that might affect both the choice of delivery provider and outcomes.5–8
Obstetric risk is typically divided between providers, with family physicians and obstetricians sharing the lowest-risk patients, obstetricians caring for moderate-risk patients, and sub-specialized perinatologists caring for the highest-risk individuals. Although traditional statistical methods can be used to adjust for observed differences between these groups, they cannot be used to adjust for unobserved differences. For example, the presence of gestational diabetes mellitus is usually noted, but its severity is often not coded in administrative databases. Women with mild diabetes mellitus are usually eligible for delivery by family physicians, but those with severe diabetes (and the attendant increased risk of adverse outcome) are not. There are many prominent examples where traditional analyses of observational data produced results that were subsequently refuted by randomized trials, presumably because of these unmeasured or unknown factors that also affect treatment decisions or outcomes.4,9–11
The instrumental variable method from the field of econometrics is a technique designed to control for unmeasured covariates in regression analyses. Results from instrumental variable–controlled observational analyses of the effect of angiography after myocardial infarction9 and of long-acting bronchodilators on asthma control12 closely approximated those of randomized controlled trials, whereas analyses using traditional statistical methods differed substantially. Instrumental variable analyses of obstetric data have shown that traditional statistical approaches significantly underestimate the mortality benefit of high-volume hospitals for high-risk neonates.4 The objective of the current study was to compare perinatal mortality and maternal morbidity and mortality for deliveries by family physicians and obstetricians using instrumental variable methodology.
Methods
Study design, data sources and population
We collected maternal and neonatal data for all of Canada (except the province of Quebec) from the Canadian Institute for Health Information (CIHI) Discharge Abstracts Database for deliveries between Apr. 1, 2006, and Mar. 31, 2009. This database has been used for several surveillance reports13–15 and numerous studies16–22 of obstetric outcomes, and it captures clinically significant diagnoses with high sensitivity and specificity.23,24
We linked these records to Statistics Canada census-derived socioeconomic information using the maternal residential postal code.25 We also accessed records for discharged and readmitted infants. Infants with birth weight less than 500 g or gestational age less than 20 weeks at delivery were excluded, for consistency with other studies.2,26 We analyzed maternal data independently, whereas neonatal records were linked to the corresponding maternal record using a linkage variable provided by CIHI or by probabilistic linkage using other variables. We excluded infant records that could not be matched to a single mother. This research was approved by the provincial Health Research Ethics Authority (Newfoundland and Labrador).
Outcomes
Perinatal mortality and a composite of maternal morbidity and mortality were the primary outcome measures. The perinatal outcome was defined as in-hospital death after 20 weeks’ gestational age, and the maternal outcome was defined similarly to Joseph and associates,17 with some additions from other papers.21,27 Although most definitions of perinatal mortality exclude death after 7 days of age, we included such deaths for infants who were continuously hospitalized (including transfers), to protect against bias associated with technologically advanced hospitals and providers who had the capability to keep infants alive longer than that. The codes used to define the outcomes are listed in Appendices 1 and 2 (available at www.cmaj.ca/lookup/suppl/doi:10.1503/cmaj.141633/-/DC1).
Group assignment
The study database included fields recording the types of providers involved in delivering care during the admission and the role that each played. Mothers were enrolled in the family physician group if that provider was listed at any point as the most responsible provider. This approach appropriately assigns patients for whom a family physician delivery was planned but who experienced intrapartum complications requiring transfer to an obstetrician or other provider (e.g., for cesarean section). We believe that this is a conservative assignment of patients that may bias against family physicians in some models where care is shared between these providers and obstetricians. In these models, high-risk patients for whom an obstetrician delivery is planned are often admitted under the family physician.
Mothers attended by midwives were identified in a similar fashion but were excluded from analyses. The remaining patients were categorized into the obstetrician group if the delivery provider was an obstetrician, and any records remaining after these assignments were excluded. We conducted sensitivity analyses assigning patients solely to the practitioner coded as most responsible or to the practitioner coded as the delivery provider. We conducted an additional sensitivity analysis including midwife deliveries in the sample.
Instrumental variable
An instrumental variable is one that predicts the receipt of treatment but is not directly associated with outcomes, except through its effect on treatment. We used the proportion of women living within the catchment of the woman’s local hospital who were delivered by a family physician as an instrumental variable. We compared comorbidities across quintiles of the instrumental variable and calculated F statistics and partial correlations as measures of its suitability. We also confirmed that the F statistics exceeded the values defining an acceptable instrument.28 Additional information is provided in Appendix 3 (available at www.cmaj.ca/lookup/suppl/doi:10.1503/cmaj.141633/-/DC1).
Hospital and other covariates
We assigned women to hospital catchment areas using the Hospital Referral Region method,29 except that we did not adjust catchment areas for geographic contiguity. Briefly, this method assigns a postal code to a given hospital when a plurality of patients living within that postal code is admitted to that hospital for their acute inpatient care. All (not just obstetric) visits to acute care hospitals for the study period were used to assign postal codes in this fashion.
Each hospital in our data set was assigned a service level according to Canadian Pediatric Society guidelines.30 For additional details, see Appendix 3. We conducted sensitivity analyses including only a tertiary hospital term and excluding delivery hospital level entirely from the regression analyses. We also conducted sensitivity analyses including cesarean section or procedural delivery (forceps or vacuum delivery) as a covariate in our regression models.
Statistical analysis
Analyses were conducted using Stata software, version 13.1 (StataCorp LLP). We measured the effect of delivery provider on our outcomes initially with bivariable logistic regression, then with multivariable logistic regression controlling for the covariates listed in Appendix 3. We did not adjust logistic models for the instrumental variable. We estimated risk ratios from these logistic models as described previously.31 We also used 3 different instrumental variable methods because results sometimes vary with different approaches.32 For the first method, we divided the data into quintiles of the instrumental variable and estimated the relative risk of the outcome for each quintile, as described previously.31 We then analyzed the data using 2-stage least-squares regression, including a method that is robust to weak instruments,33 and the generalized method of moments (Stata syntax, including analytic first derivatives, as published previously34). All analyses were adjusted for clustering at the delivery hospital. For an expanded description of the instrumental variable analyses, please see Appendix 3.
Results
The cohorts and their exclusions are described in Figure 1. Of the 859 180 neonatal records accessed (before exclusions), 69.4% had delivery by obstetricians, 26.2% by family physicians and 3.7% by midwives. Across quintiles, midwives delivered 3.7%, 4.5%, 4.1%, 3.4% and 2.8% of the infants, respectively; the breakdown by quintile of deliveries by family physicians and obstetricians is provided in Table 1. For the remaining 0.7% of infants, either delivery was by another provider type or information on provider type was missing. We identified a total of 3600 (0.45%) perinatal deaths and 14 394 (1.82%) cases of maternal morbidity and mortality during the study period. Selected characteristics of the study population, delivery providers and hospitals are listed in Table 1. Notably, there was an apparent positive correlation between predicted risk of maternal morbidity and mortality and the proportion of the population who identified as Aboriginal, which was a significant predictor of perinatal mortality (data not shown).
Flow diagram for study cohort. FP = family physician, OB = obstetrician.
Selected characteristics of the study cohort, for records included in the final analysis
Strength of instrumental variable
Our instrumental variable predicted a wide range in the mean percentage of deliveries by a family physician (4.2% to 67.7% across quintiles). Whereas there were some differences in measured covariates across these quintiles (Table 1), there was essentially no correlation between mean perinatal mortality (r2 = 7.3 × 10−6) or maternal morbidity and mortality (r2 = 0.011) and the instrumental variable, a required characteristic to ensure unbiased results. The Kleibergen–Paap F statistics for our instrumental variable far exceeded the value necessary to define a strong instrument28 (see Appendix 3). The partial correlation coefficient between the delivery provider and the instrumental variable was 0.55 for both the neonatal and maternal analyses, indicating that 30% of the variation in the rate of delivery by family physicians was explained by the instrumental variable, which is also a marker of a strong instrument.
Effect of family physician as delivery provider
The risk ratios for the outcomes across the different family physician quintiles are presented in Table 2. Ad hoc risk ratios from the traditional 2-stage least-squares models were 0.95 for perinatal mortality and 1.01 for the maternal outcome, and the results were identical using methods that were robust to weak instruments.33 Results from logistic and generalized method of moments analyses are presented in Table 3. The sensitivity analyses described in the Methods section changed neither the direction of the estimated effects (risk ratio greater than or less than 1.0) nor the statistical significance of the association (data not shown).
Risk ratios for outcomes across FP delivery quintiles
Effect of delivery by family physicians on perinatal mortality and maternal morbidity and mortality using logistic and IV-adjusted regression
Interpretation
Using a statistical method that controls for both observed and unmeasured or unknown factors affecting obstetric outcomes, we found no difference in the risk of perinatal mortality or maternal morbidity and mortality between deliveries by family physicians and those by obstetricians. However, because of the limitations of this statistical method, the confidence intervals around the risk ratios are wide. These findings build on previous work that also supports the safety of obstetric delivery by family physicians.5–8 The current study is among the largest on this topic to date and strengthens these earlier findings by including an adjustment for unmeasured selection bias.
It is common to assume that more specialized or higher-volume medical care will result in improved outcomes. This has been most convincingly shown for highly technical, relatively uncommon procedures such as pancreatectomy, esophagectomy and elective repair of abdominal aortic aneurysm.35 Similarly, the obstetric literature has consistently shown that outcomes for high-risk newborns and mothers are best at higher-volume, more specialized hospitals;2–4 however, findings from the literature on low-risk deliveries are variable.26,36–39
Although there are exceptions, several studies have found differences between family physicians and obstetricians in the use of invasive procedures such as forceps, vacuum, episiotomy and cesarean section for the management of labour,6–8 and research suggests that this more invasive care may be harmful to mothers and infants.40,41 We chose not to adjust for procedures in our primary analyses because we wanted to capture differences in this covariate within the specialty variable. It is tempting to assume that the improved outcomes expected because of specialty training were offset in our study by an increased risk associated with a higher procedure rate. However, we observed a nonlinear trend in rates of cesarean section across the family physician quintiles (Table 1). We also conducted sensitivity analyses including cesarean section or procedural delivery as a covariate in our regression models, and these adjustments did not change our primary findings. We will further explore the effect of procedure use on obstetric outcomes in a forthcoming paper.
Instrumental variables analyses measure outcomes for the “marginal” population, which in our study consisted of those patients who would be delivered by a family physician in some jurisdictions but not in others. Although some clearly high-risk patients may on occasion be delivered by a family physician if labour proceeds precipitously and no obstetricians are available, most of these patients would have been delivered by a specialist, and our estimates of treatment effect do not apply to them. Because of this, instrumental variable–adjusted analyses do not easily translate into decisions that can be applied to an individual expectant mother, except if she fits the profile of the marginal population where she intends to deliver. Instead, instrumental variable–adjusted findings are more relevant for policy-makers who are making decisions at the health system level that will affect the choice of one treatment over another.
Some authors have found large differences between risk estimates with these 2 types of analyses;9,12 however, in our analyses, the results for instrumental variable–adjusted and comprehensively adjusted traditional models differed relatively little, particularly for the perinatal mortality outcome. This suggests that the covariates included in our analyses captured most of the variance associated with differences between our groups and/or that the indications for delivery by an obstetrician are relatively clear. Traditional statistical methods are more likely to produce unbiased estimates under one or both of these conditions.
Limitations
Several limitations of our study merit discussion. We excluded midwife deliveries from our analyses to obtain a “cleaner” comparison. Although midwives attend a substantial proportion of low-risk deliveries in some regions of the country, our statistical method adjusts for the between-group risk difference exacerbated by this exclusion. Including midwife deliveries in the sample in a sensitivity analysis did not alter our final conclusion. We have not presented separate analyses for the midwife group because they represent a small proportion of our sample overall, and there is much less variability in the proportion of deliveries by midwives across regions. Thus, the potential for bias, even with instrumental variable–adjusted analyses, is greater with this group. Furthermore, we did not have access to data for home deliveries, which make up a very small proportion of deliveries overall but often a substantial proportion of deliveries by midwives, and patients who choose a home delivery likely differ systematically from those who choose hospital delivery. Because home deliveries accounted for less than 1.5% of deliveries in Canada during our study period,42 we do not feel that this exclusion significantly biases the analyses presented here.
We did not attempt to measure perinatal outcomes other than death. Infants may benefit in other ways from delivery by obstetricians or family physicians, particularly the lower-risk infant population in whom the risk of death is very low. It is difficult to imagine important maternal complications that are not captured in our comprehensive definition of maternal morbidity. However, there is a small likelihood that we did not observe differences between groups in diagnoses that arose after hospital discharge, because we did not attempt to access maternal readmission records (as we did for the neonates).
Although our chosen variable met the commonly used thresholds for a strong instrument, there was some variability in measured covariates across quintiles of this variable (Table 1), which raises the possibility that unobserved covariates also vary across these levels. In particular, although the predicted risk of perinatal mortality was relatively consistent across quintiles, there was an apparent positive correlation between predicted risk of maternal morbidity and mortality and the instrument (although the r2 was very low). There was also a similar trend across quintiles in the proportion of the population who identified as Aboriginal, which was a significant predictor of perinatal mortality. Although we adjusted for these covariates, insofar as unobserved covariates vary in a similar pattern, the trend we observed would tend to bias results against family physicians. Additional unobserved covariates with a similar distribution may explain the trend toward worse maternal outcomes that we observed in the instrumental variable–adjusted analyses, despite the lower risk observed in logistic analyses (Table 3).
Finally we acknowledge that these data are now somewhat out of date. It is possible that there has been a cultural shift resulting in a change in the proportion of deliveries by family physicians or midwives in recent years. In addition, the increasing attention paid to procedural intervention in the obstetric literature (e.g., studies by Souza and colleagues40 and Villar and associates41) may have changed the management of deliveries by obstetricians.
Conclusion
In a large, population-based cohort of Canadian patients, we observed similar risks of perinatal mortality and adverse maternal outcome between obstetric deliveries by family physicians and those by obstetricians, using an econometric method designed to control for unmeasured bias. It remains to be seen whether there are differences between these groups in terms of other outcomes. Because of the analytical approach used, these findings apply only to mothers and infants who would be eligible for delivery by either family physicians or specialists in at least some jurisdictions of the country. These results do not apply to mothers who are consistently referred for delivery by obstetricians in all jurisdictions. Future research should explore the effect of different delivery providers on other outcomes and on health resource utilization.
Acknowledgement
The authors thank Dr. James Rourke for comments on an earlier version of this paper.
Footnotes
CMAJ Podcasts: author interview at https://soundcloud.com/cmajpodcasts/141633-res
Competing interests: None declared.
This article has been peer reviewed.
Contributors: Kris Aubrey-Bassler was the principal investigator and contributed to initial conception, funding application, study design, data release protocol, data retrieval, statistical analysis, manuscript writing and editing. Richard Cullen contributed to data linkage and cleaning, correspondence with the data holder, ethics approval, manuscript preparation and revision. Alvin Simms contributed to data linkage and cleaning and to manuscript revision, and contributed geographic expertise. Shabnam Asghari contributed epidemiologic and administrative data expertise and participated in statistical analysis and manuscript revision. Joan Crane contributed obstetric expertise and participated in initial conception, funding application, study design and manuscript revision. Peter Wang contributed statistical expertise and participated in study design and manuscript revision. Marshall Godwin contributed family medicine expertise and participated in the funding application and manuscript revision. All authors approved the final version for publication and agree to act as guarantors of the work.
Funding: Funding for this study was provided to the research team through the Dr. A.R. Cox Research Grant via the Medical Research Endowment Fund at Memorial University of Newfoundland. The research reported was conducted independently of the funder.
Disclaimer: Parts of this article are based on data and information provided by the Canadian Institute for Health Information (CIHI). However, the analyses, conclusions, opinions and statements expressed herein are those of the authors and not those of CIHI.
Data sharing: The Canadian Institute for Health Information (CIHI) does not allow data sharing. However, upon appropriate approvals, raw data are available from CIHI, and we are able to share statistical software syntax files.
- Accepted July 3, 2015.
Podcast