# Racial disparities in infant mortality: what has birth weight got to do with it and how large is it?

- Timothy B Gage
^{1, 2}Email author, - Fu Fang†
^{1}, - Erin K O'Neill†
^{1}and - A Gregory DiRienzo†
^{2}

**10**:86

**DOI: **10.1186/1471-2393-10-86

© Gage et al; licensee BioMed Central Ltd. 2010

**Received: **19 February 2010

**Accepted: **28 December 2010

**Published: **28 December 2010

## Abstract

### Background

It has been hypothesized that birth weight is not on the causal pathway to infant mortality, at least among "normal" births (i.e. those located in the central part of the birth weight distribution), and that US racial disparities (African American versus European American) may be underestimated. Here these hypotheses are tested by examining the role of birth weight on racial disparities in infant mortality.

### Methods

A two-component Covariate Density Defined mixture of logistic regressions model is used to decompose racial disparities, 1) into disparities due to "normal" versus "compromised" components of the birth cohort, and 2) further decompose these components into indirect effects, which are associated with birth weight, versus direct effects, which are independent of birth weight.

### Results

The results indicate that a direct effect is responsible for the racial disparity in mortality among "normal" births. No indirect effect of birth weight is observed despite significant disparities in birth weight. Among "compromised" births, an indirect effect is responsible for the disparity, which is consistent with disparities in birth weight. However, there is also a direct effect among "compromised" births that reduces the racial disparity in mortality. This direct effect is responsible for the "pediatric paradox" and maybe due to differential fetal loss. Model-based adjustment for this effect indicates that racial disparities corrected for fetal loss could be as high as 3 or 4 fold. This estimate is higher than the observed racial disparities in infant mortality (2.1 for both sexes).

### Conclusions

The results support the hypothesis that birth weight is not on the causal pathway to infant mortality among "normal" births, although birth weight could play a role among "compromised" births. The overall size of the US racial disparities in infant mortality maybe considerably underestimated in the observed data possibly due to racial disparities in fetal loss.

## Background

It has been argued that birth weight may not be on the "causal pathway" to infant mortality [1–3]. The best developed argument, originating with the Wilcox-Russell hypothesis [2, 4, 5], is supported by qualitative analyses using directed acyclic graphs [6]. Both of these approaches are based on simple graphical observations of the response of birth weight and birth weight specific infant mortality to exogenous stressors, such as smoking or altitude. The Wilcox-Russell hypothesis [2, 4, 5] suggests that in response to a stressor, the birth weight specific infant mortality curve and birth weight distribution appear to shift right or left together resulting in no change in total mortality. Consequently, there is no indirect effect of the stressor due to the shift in birth weight. Any changes in mortality due to the stressor are hypothesized to be due to the entire mortality curve shifting up or down independently of birth weight, i.e. direct effects of the stressor. Complementary analyses using directed acyclic graphs identify three plausible models which could account for the dynamics of birth weight specific mortality [6]. One model supports the Wilcox-Russell hypothesis [2, 4, 5], while the other two include birth weight on the "causal pathway" [6].

If birth weight is not on the "causal pathway", i.e. it does not mediate the effect of race on infant mortality (Figure 1c), then the US national policy of reducing infant mortality [7] in general, and racial disparities in particular, by reducing the low birth weight rate might not be effective. On the other hand, the Wilcox-Russell hypothesis [2, 4, 5] only applies to "normal" births (which they defined as those in the central part of the birth weight distribution) and not all births [5], so birth weight could still mediate the effect of race on infant mortality among the remaining births. An initial quantitative statistical test of the Wilcox-Russell [2, 4, 5] and Hernández-Diaz et al. [6] hypotheses using Covariate Density Defined mixture of logistic regressions (CDDmlr) with maternal age as a stressor supports the argument that birth weight is not on the causal pathway to infant mortality for either "normal" or the remaining "compromised" births [8].

Wilcox and others [5, 9–11] have also argued that racial disparities in infant mortality may be underestimated. This view is based on the simple graphical observation that lower birth weight African American births have better survival than their European American peers with similar birth weight despite much higher mortality overall, i.e. the racial "birth weight or pediatric paradox". The hypothesis is that unmeasured (and hence uncontrolled) heterogeneity between the racial groups might mask part of the true racial disparities. It has been shown that CDDmlr isolates the race "pediatric paradox" within the "compromised" subpopulation, allowing better control of this phenomenon [10].

The objective of this paper is to quantitatively document the role that birth weight plays in racial disparities in infant mortality using the 2001 United States non-Hispanic African and European American birth cohorts controlling for sex. In particular, we statistically test the hypothesis that birth weight is on the "causal pathway" to infant mortality and decompose racial disparities in infant mortality into effects, which are independent of birth weight (direct effects of race) and effects, which are due to the racial disparities in birth weight (indirect effects of race mediated by birth weight). A secondary aim is to estimate the magnitude of the racial disparities in infant mortality while controlling for the "pediatric paradox". We do not propose that "race" is the cause of these disparities, but simply a proxy for a collection of stressors (e.g. socio-economic status, education, and genetic etc, some of which may be unobserved), which are the underlying causes of these differences.

## Methods

### Data Source

Descriptive statistics for the 2001 sample populations

Birth Cohort | # Births | # Deaths | CDR | Birth Weight (grams) | |||
---|---|---|---|---|---|---|---|

Min | Mean | Median | Max | ||||

Eur. Am. F. | 1,023,583 | 3,558 | 3.48 | 500 | 3342 | 3365 | 6350 |

Eur. Am. M. | 1,076,814 | 4,880 | 4.53 | 500 | 3461 | 3487 | 7858 |

Af. Am. F. | 255,758 | 1,865 | 7.29 | 500 | 3092 | 3135 | 7002 |

Af. Am. M. | 264,130 | 2,545 | 9.64 | 500 | 3200 | 3260 | 7220 |

### Statistical Model - CDDmlr

Covariate Density Defined mixture of logistic regressions (CDDmlr), while a generally applicable statistical procedure, was specifically designed to test the Wilcox-Russell hypothesis [8]. It decomposes the birth weight distribution into a number of subpopulations, using standard mixture of Gaussian distributions, and simultaneously fits a separate birth weight specific mortality curve to each of the subpopulations identified by the birth weight density submodel [10]. A two-component CDDmlr model using Gaussian distributions (truncated at 500 grams) and logistic regressions (a 2^{nd} degree polynomial of birth weight) is the parsimonious model, that fits birth weight distributions [12] and birth weight specific mortality curves [9, 10] remarkably well. One subpopulation accounts for most births in the center of the birth weight distribution and appears to identify "normal" births, while the other accounts for most low and macrosomic births, and is hence called "compromised" births [9, 10, 12]. Clearly the "compromised" subpopulation represents a heterogeneous group, i.e. births "compromised" by a variety of potential factors. However, increasing the number of subpopulations does not resolve the "compromised" subpopulation into separate groups [13] and placing constraints on the fitting process [4] simply reduces the goodness of fit. The model represents the maximum likelihood division of the birth weight distribution given the assumption that the birth weight distribution is the sum of two Gaussian subpopulations. Furthermore, the "compromised" subpopulation differs slightly from Wilcox's "residual" subpopulation [4] given that it also accounts for births in the normal birth weight range, where as Wilcox's "residual" births [4] were restricted to the lower tail. However, a number of clinicians have argued that "compromised" births do occur in the normal birth weight range, but are not recognized as "compromised" when using the arbitrary low birth weight standards (i.e. <2500 grams) and are hence understudied [14, 15]. Given that the Reverse-J-shaped birth weight specific mortality curves fitted to each of the two subpopulations (i.e. "normal" and "compromised") is parsimonious [9, 10], we assume that the reverse-J shape is due to other unspecified covariates and not a "causal" effect of birth weight. This is consistent with Hernández-Diaz and her colleagues' assumption [6] that the reverse-J shape of the mortality curve is due to other unmeasured covariates, such as the theory of Basso and Wilcox [16, 17] that the reverse-J shape is due to confounding. Here, we use CDDmlr to statistically examine the Wilcox-Russell [2, 4, 5] and Hernández-Diaz et al. [6] hypotheses for both "normal" and "compromised" births. In addition, since the "pediatric paradox" is associated with the "compromised" subpopulation [10], CDDmlr can control for this phenomenon as well.

*θ*) in the birth weight density submodel and the six parameters (referred to collectively as

*β*or

*β**, representing the two 2

^{nd}degree polynomials of birth weight or standardized birth weight, respectively) in the mortality submodel of the basic CDDmlr model [10] as linear functions of a dummy variable (e.g. race). Thus this stratified model can quantify the differences in the birth weight distribution (i.e. the proportion of "compromised" births, and the means and standard deviations of both subpopulations) and the (standardized) birth weight specific mortality characteristics between African and European American birth cohorts. In this study, birth weight is standardized (Z-scored) for each subpopulation based on the subpopulation specific mean and variance. This step essentially breaks the association of race and birth weight so that we can estimate the birth weight independent effect (direct effect) and any remaining birth weight dependent effect [18, 19]. The latter may be potentially due to a direct effect of birth weight on infant mortality, or uncontrolled confounding between birth weight and infant mortality, or an interaction of race and birth weight on infant mortality. In particular, we investigate the effects of race on:

- (i)
the logit of minimum mortality (i.e. a vertical shift of the mortality curve by race, the direct effect of race);

- (ii)
the optimal standardized birth weight (i.e. a horizontal shift of the mortality curve by race, the indirect effect of race described by Wilcox-Russell [2, 4, 5]); and

- (iii)
the particular shape of the reverse-J-shaped standardized birth weight specific mortality curve (i.e. a second possible indirect effect of race, not considered by Wilcox-Russell [2, 4, 5] but equivalent to the interaction of the stressor and birth weight proposed by Hernández-Diaz et al. [6] as a possible alternative cause of the reverse-J shape of birth weight specific infant mortality).

This second indirect effect of race through birth weight (iii) occurs when a change in the variance of birth weight is not reflected in a compensatory change in the shape of the birth weight specific mortality curve. So that the standardized birth weight specific mortality curve changes by race. Finally, the mixing proportion may contribute to the overall observed racial disparities in infant mortality. This is an additional effect of race, which was not discussed in the Wilcox-Russell hypothesis [2, 4, 5] or its extension by Hernández-Diaz et al. [6]. However, it is similar to the concept of "confounding" in Basso and Wilcox [16, 17]. The mixing proportion does involve the birth weight density, nonetheless, the role of birth weight in this case is unclear, depending upon whether birth weight is the cause or the effect of being "compromised." In summary, the CDDmlr provides a reasonable statistical examination of the Wilcox-Russell hypothesis [2, 4, 5] concerning the potential effects of a stressor, i.e. its direct and/or indirect effects on infant mortality among "normal" as well as "compromised" births [8]. It can potentially distinguish between the "plausible" directed acyclic graphs identified by Hernández-Diaz et al. [6] (Figure 1).

*x*) only CDDmlr model (i.e. CDDmlr without any exogenous covariate) of infant mortality (

*y*) is formally defined as a product of the conditional mortality submodel

*f*

_{2}(

*y*|

*x*;

*θ, β*) and the birth weight density submodel

*f*

_{1}(

*x*;

*θ*):

*f*

_{1}(

*x*;

*θ*) is given by

*π*

_{ s }, the mixing proportion, is defined as the proportion of births belonging to the less numerous of the two subpopulations, that is, the secondary subpopulation (

*s*, "compromised" subpopulation) as opposed to the primary subpopulation (

*p*, "normal" subpopulation). The reparameterization of

*π*

_{ s }(Eq. 3) transforms the 0 and 1 bounds on

*π*

_{ s }to minus and plus infinity, respectively. For

*i*=

*s*and

*p*, $\tilde{N}\left(x;{\mu}_{i},{\sigma}_{i}^{2}\right)$ represents the Gaussian density, truncated at 500 grams, with mean

*μ*

_{ i }and variance ${\sigma}_{i}^{2}$. The conditional mortality submodel

*f*

_{ 2 }(

*y*|

*x*;

*θ*,

*β*) with two subpopulations is given by

*q*

_{ s }(

*x*;

*θ*) is the probability that an infant with birth weight

*x*belongs to the

*s*subpopulation. For

*i*=

*s*and

*p*,

*f*

_{1}(

*x*;

*θ*) (Eq. 2) determines that

Overall, there are 11 parameters, five defining the birth weight distribution, and six defining the subpopulation-specific mortalities.

*z*) on each of the 11 parameters in the basic CDDmlr model. Second, for

*i*=

*s*and

*p*, standardized birth weight (${x}_{i}^{*}$, i.e.

*x*is standardized according to the respective subpopulation mean and standard deviation) is used in the corresponding logistic regression function. Thus

This extended model includes 22 parameters, 11 representing the characteristics of European American births, and 11 representing the differences of African compared to European American birth outcomes, that is, the "race" effect. The 5 indicator variable terms in the density submodel (i.e. *η*_{1}, *μ*_{
i
} _{, 1}, and *σ*_{
i, 1
} for *i* = *s* and *p*) account for the effects of "race" on the birth weight distribution, while the 6 indicator variable terms in the mortality submodel (i.e. ${A}_{i,\mathit{1}}^{*}$, ${B}_{i,1}^{*}$, and ${C}_{i,\mathit{1}}^{*}$ for *i* = *s* and *p*) account for the racial difference in the standardized birth weight specific mortality curves.

### Model Fitting

^{nd}degree polynomial of standardized birth weight specific mortality curves are fitted in linear form, and then transformed to non-linear form after fitting. This significantly reduces the computational resources necessary to fit the model. The resulting parameter estimates are presented in Table 2. This model shows no evidence of lack of fit based on the Hosmer-Lemeshow statistic (with a p-value of 0.66 and 0.41 for females and males, respectively). Bias-adjusted 95% confidence intervals are estimated from 200 bootstrap samples of 200,000 births each, which are randomly generated from the entire birth cohort (as opposed to the more conventional procedure of re-sampling with replacement from the original sample a sample the same size as the original sample). The conventional procedure requires excessive computational resources. An independent study using maternal education as a binary exposure variable suggests that our bootstrap results are consistent with results from the conventional bootstrap method.

Parameter estimates for the 2001 sample populations

Birth Cohort | z = 0 | z = 1 | z = 0 | z = 1 |
---|---|---|---|---|

Eur. Am. F. | Af. Am. F. | Eur. Am. M. | Af. Am. M. | |

| -2.75 | -2.62 | ||

| 0.46 | 0.35 | ||

| 2678 | 2739 | ||

| -639 | -706 | ||

| 1098 | 1098 | ||

| 205 | 238 | ||

| 3380 | 3509 | ||

| -211 | -221 | ||

| 455 | 474 | ||

| 0 | -2 | ||

| -7.12 | -7.11 | ||

| 1.06 | 1.29 | ||

| 0.77 | 0.68 | ||

| -1.73 | -0.83 | ||

| 0.36 | 0.11 | ||

| 0.54 | 0.62 | ||

| -6.85 | -6.57 | ||

| 0.97 | 1.22 | ||

| 0.25 | 0.20 | ||

| 0.75 | 0.69 | ||

| -0.15 | -0.44 | ||

| 0.01 | 0.05 |

### Decomposition of the Racial Disparity

Decomposition of the racial disparity is carried out in two steps. First, the total absolute racial disparity in infant mortality is decomposed into deaths attributable to differences in the mixing proportion and rate effects for "normal" and "compromised" births using standard Kitagawa decomposition [21].

*i*subpopulation of European American births is the weighted average probability of infant death across all birth weights, that is

*i*subpopulation, the probability of death for an African American birth with is given by

*i*subpopulation of African American births is the weighted average probability of infant death across all birth weights, that is

*i*subpopulation is given by

*F*_{
i
} _{,1} is referred to as the direct factor of "race" in the *i* subpopulation. It is a constant, and independent of birth weight. *F*_{
i
} ,2 is referred to as the indirect factor of "race". It represents the combined effect of all birth weight related factors on the racial disparity in infant mortality of the *i* subpopulation. In particular, birth weight related factors include differences in the shape and the horizontal shift of the reverse-J-shaped standardized birth weight specific mortality curve, the non-linear transformation between the probability, and the logit of infant death at any standardized birth weight, as well as the difference in the truncating value of the standardized birth weight distributions between African and European American births.

## Results

### Characteristics of Race Specific Birth Weight Distributions and Infant Mortality

Model-estimated birth weight distribution and mortality characteristics for the 2001 sample populations

Birth Cohort | Eur. Am. F. | Af. Am. F. | Eur. Am. M. | Af. Am. M. | ||||
---|---|---|---|---|---|---|---|---|

estimate | LCI UCI | estimate | LCI UCI | estimate | LCI UCI | estimate | LCI UCI | |

"Normal" Subpopulation | ||||||||

Proportion (%) | 94.0 | (93.5; 94.3) | 90.7 | (90.5; 91.0) | 93.2 | (92.7; 93.6) | 90.6 | (90.5; 90.8) |

Mean (g) | 3380 | (3378; 3382) | 3169 | (3168; 3170) | 3509 | (3507; 3512) | 3288 | (3287; 3289) |

Standard Deviation (g) | 455 | (453; 457) | 455 | (454; 457) | 474 | (472; 476) | 472 | (471; 473) |

LBW Rate (%) | 2.8 | (2.7; 2.9) | 7.0 | (6.9; 7.0) | 1.8 | (1.7; 1.9) | 4.8 | (4.7; 4.9) |

Death Rate | 2.3 | (2.1; 2.6) | 4.4 | (4.2; 4.5) | 2.9 | (2.7; 3.2) | 5.2 | (5.0; 5.4) |

Percent of Total DR (%) | 63.1 | (58.2; 67.8) | 54.7 | (53.5; 56.1) | 60.5 | (55.6; 65.0) | 49.0 | (47.5; 50.5) |

"Compromised" Subpopulation | ||||||||

Proportion (%) | 6.0 | (5.7; 6.5) | 9.3 | (9.0; 9.5) | 6.8 | (6.4; 7.3) | 9.4 | (9.2; 9.5) |

Mean (g) | 2678 | (2630; 2730) | 2039 | (1995; 2072) | 2739 | (2684; 2785) | 2034 | (1995; 2072) |

Standard Deviation (g) | 1098 | (1067; 1126) | 1303 | (1289; 1320) | 1098 | (1076; 1121) | 1336 | (1323; 1350) |

LBW Rate (%) | 42.3 | (40.4; 43.8) | 58.5 | (57.7; 59.2) | 40.4 | (38.7; 42.4) | 58.3 | (57.5; 59.1) |

Death Rate | 21.2 | (18.3; 24.9) | 35.5 | (33.7; 37.1) | 26.3 | (23.0; 30.0) | 52.4 | (50.3; 55.0) |

Percent of Total DR (%) | 36.9 | (32.2; 41.8) | 45.2 | (43.9; 46.5) | 39.5 | (35.0; 44.4) | 51.0 | (49.5; 52.5) |

Total Population | ||||||||

LBW Rate (%) | 5.2 | (5.1; 5.3) | 11.8 | (11.7; 11.9) | 4.5 | (4.4; 4.6) | 9.8 | (9.8; 9.9) |

Death Rate | 3.5 | (3.2; 3.7) | 7.3 | (7.1; 7.4) | 4.5 | (4.3; 4.8) | 9.6 | (9.4; 9.9) |

### Racial Differences in Birth Weight Distributions

Race has substantial effects on the distribution of birth weight (Table 3 Figure 2a-2b, and 3a-3b). For both sexes, the proportion of "normal" births is approximately 3% smaller and the means of both subpopulations are significantly smaller in African American births compared to European American births. The standard deviation of "compromised" African American births is significantly larger compared to European American births. However, there is no difference in the standard deviation of the "normal" subpopulation between African and European American infants of the same sex. Collectively these differences account for the larger low birth weight rates generally observed in African American birth weight distributions (Table 3) [7, 22].

### Racial Differences in Infant Mortality

There are substantial racial differences in infant mortality as well (Table 3 Figures 2c-2e, and 3c-3e). The subpopulation specific results show that African American birth weight specific "normal" mortality is larger than European American mortality (Figures 2c and 3c), while African American birth weight specific "compromised" mortality is generally smaller than European birth weight specific mortality (Figures 2d and 3d). Birth weight specific total mortality shows the "pediatric paradox", that is significantly smaller African American mortality at lower birth weights but larger mortality in the larger birth weight range (Figures 2e and 3e). The lower mortality of African Americans at smaller birth weights is accounted for by the lower mortality of "compromised" African American births compared to European American births at the smaller birth weights where the "compromised" subpopulation predominates (Figures 2d and 3d). Similarly the excess mortality of African Americans in the normal birth weight range is accounted for by the larger mortality of African American "normal" births compared to European American "normal" births in the central part of the birth weight range where "normal" births predominate (Figures 2c and 3c).

Kitagawa decomposition of the observed racial disparities in infant mortality (death per 1000 births)

Decomposition | Females | Males | ||
---|---|---|---|---|

estimate | LCI UCI | estimate | LCI UCI | |

Mixing Proportion Effect | 0.81 | (0.67; 0.95) | 0.90 | (0.72; 1.08) |

Rate Effect | ||||

"Compromised" | 1.09 | (0.79; 1.37) | 2.11 | (1.75; 2.48) |

"Normal" | 1.90 | (1.67; 2.19) | 2.08 | (1.77; 2.39) |

Total Disparity | 3.81 | (3.53; 4.07) | 5.09 | (4.81; 5.40) |

### Birth Weight and the Racial Disparity

Subpopulation specific racial effect (relative risk) of infant mortality # decomposed into direct and indirect multiplicative factors

Racial Effect | Females | Males | ||
---|---|---|---|---|

estimate | LCI UCI | estimate | LCI UCI | |

"Normal" Subpopulation | ||||

Relative Risk | 1.89 | (1.71; 2.14) | 1.77 | (1.63; 1.98) |

Direct Factor | 2.11 | (1.87; 2.49) | 2.00 | (1.71; 2.68) |

Indirect Factor | 0.89 | (0.78; 1.01)* | 0.88 | (0.71; 1.00)* |

"Compromised" Subpopulation | ||||

Relative Risk | 5.19 | (4.29; 6.16) | 5.12 | (4.37; 6.08) |

Direct Factor | 0.18 | (0.01; 0.69) | 0.43 | (0.11; 1.55)* |

Indirect Factor | 29.38 | (13.83; 91.15) | 11.80 | (4.31; 33.46) |

Total Birth Cohort | ||||

Relative Risk | 3.10 | (2.83; 3.41) | 3.09 | (2.89; 3.37) |

Among "normal" births, there is a significant direct effect of being African American that contributes to excess mortality in African American births (Table 5). The indirect effect among "normal" births is marginally insignificant in both males and females and tends to **reduce** African American mortality! The direction of this association is surprising given that mean birth weight of "normal" African American births is significantly smaller than that of European Americans (Table 3, Figures 2a and 3a).

Among "compromised" births, on the other hand, the indirect effect is significant and contributes to the excess African American infant mortalities (Table 5). This excess infant mortality is consistent with expectation, that is higher mortality is associated with a significantly lower birth weight (Table 3). In addition, a direct effect among "compromised" African American births **reduces** infant mortality (Table 5). It is significant for females, but not for males. Since the direct and indirect effects tend to compensate for each other the true size of these effects may exceed the absolute effect predicted for each subpopulation.

Overall, a direct effect on the "normal" subpopulation is responsible for the higher infant mortality of African American births in the normal birth weight range (Figures 2e and 3e), while a direct effect on the "compromised" subpopulation is responsible for the lower infant mortality of African American births at lower birth weights (Figures 2e and 3e). As a result, the race "pediatric paradox" (i.e. African Americans have lower mortality at lower birth weights compared to their European American peers), is due to this beneficial direct effect of being an African American "compromised" birth (Figure 2f and 3f). Finally, a large indirect effect occurs in the "compromised" subpopulation (Table 5).

## Discussion

CDDmlr was designed to examine the Wilcox-Russell hypothesis [2, 4, 5], and its extensions, e.g. Hernández-Diaz et al. [6], and to provide quantitative estimates of the direct effects, which are independent of birth weight and the indirect effects that may operate through birth weight. As described above we have implemented the same assumptions as Wilcox-Russell [2, 4, 5] and Hernández-Diaz et al. [6]. Nevertheless, application of a quantitative model has some additional limitations over qualitative models, e.g. data quality and quantity, as well as the details of the implementation.

The analyses are based on the public use samples of the NCHS linked birth death files. These have very large sample sizes (Table 1) so there are unlikely to be issues with power. Birth weight is considered to be reliably measured. Mortality estimates may be slightly biased due to problems associated with linking birth and death certificates. However, these are the same data, with the same problems, that most representative analyses of the US are based upon. For our purposes the most troubling defects are that births at <500 grams and LMP gestational ages <20 weeks are not consistently reported by all states [23]. Following many analyses of these data we have truncated the data to avoid this problem. Consequently, we have used Gaussian distributions truncated at 500 grams to match the data.

One technical difficulty in models of this kind is estimating unbiased direct and indirect effects. The qualitative analysis in Hernández-Diaz et al. [6] is based on the assumptions of counterfactuals [24–26]. Here we take an alternative approach, developed from statistical decision theory [18]. In our case, we have modelled the birth weight density as the sum of two Gaussian distributions and the subpopulation specific mortality curves as a 2^{nd} degree polynomial of Z-scored birth weight standardized with respect to these Gaussians. This eliminates the main effects (associations) of race and birth weight and the logistic regressions can then estimate the direct effect of race on infant mortality versus potential interaction effects of race and birth weight on infant mortality. Direct and indirect effects can be estimated using procedures similar to direct standardization [18]. The result is called a "generated direct effect" by Geneletti [18], which is similar to Pearl's "natural direct effect" [24]. Since the "normal" and "compromised" subpopulations are defined as Gaussian distributions, the appropriate distribution is theoretically available for direct standardization. In this regard, truncation of the data at 500 grams creates a significant truncation difference in the standardized birth weight distributions between African and European American "compromised" births. Consequently, the results based on a common reference population (i.e. the European American distribution, Table 5) may be preferred. Identification issues concerning "generated direct effects" are discussed by Geneletti [18].

One advantage of the decision theory approach is that the assumptions concerning the existence of counterfactuals are not necessary. However, like counterfactual methods, the same strong unmeasured covariate assumptions are required. In particular; a) no unmeasured covariates which affect the stressor (race in this case) and the racial disparities in infant mortality, b) no unmeasured confounding of race and birth weight, and c) no unmeasured confounding of birth weight and infant mortality. Assumption a is necessary to estimate total racial disparities, all three are needed to estimate "generated direct effects" [18].

These assumptions may be less of a problem with race than with other variables such as smoking, which have more precise definitions. Race is typically considered to be socially constructed and defined as that collection of variables (some of which may be observable and some of which are currently unobservable) that are associated in some way with reported race. Given this view, all confounders of racial effects on birth weight or infant mortality, are integral parts of the definition of race. This is the assumption generally used when reporting total "racial disparities", such as those presented in Table 1. Of course it is possible to partial out the effects of measured confounders on racial disparities, e.g. the effects of maternal age, but what are left in this case are simply all the unmeasured and unknown effects of race. The results presented above are uncorrected for confounders, and consequently represent the sum total of all direct and indirect effects associated in some way with race. This should be considered when interpreting the results.

Based on the "pediatric paradox", Wilcox has argued that racial disparities may be underestimated due to unmeasured confounding [5, 9–11]. Gage has hypothesized that the lower birth weight specific mortality of African compared to European American "compromised" birth cohorts[10] is due to the heavier fetal loss and selection documented among African Americans [27, 28]. If this assumption is correct, then differential fetal loss is associated with the direct effect of being African American in the "compromised" subpopulation and with the "pediatric paradox". This interpretation is also consistent with Platt et al.'s finding [29] that the race "birth weight paradox" disappears when observable fetal deaths (total fetal loss is not observable) are included (as well as live births and infant deaths) in the analysis of racial disparities in infant mortality. Should this selection bias be included in the definition of "race" or should differential fetal loss be excluded from the definition of race? The answer depends upon the question, but CDDmlr potentially makes it possible to correct for this "unmeasured" source of confounding.

Model-based adjustment of this effect yields relative risks of 4.2 and 3.6 for African American female and male births, respectively. These are higher than the predicted total relative risks in Table 5, and much higher than the observed relative risk of 2.1 for both sexes derived from Table 1. This adjusted racial disparity needs to be considered with some caution, since it assumes that the direct effect in the "compromised" subpopulation is completely due to selection bias and can be reduced to zero while all other modelled effects remain the same. Nevertheless, it is possible that a substantial part of the racial disparity in infant mortality is hidden by differential fetal loss.

We assume that unmeasured confounding of birth weight and infant mortality (assumption c) is responsible for the reverse-J shape of the birth weight specific mortality curve [16, 17] and that the reverse-J shape is not a "causal" effect of birth weight. We have implemented the characteristic reverse-J shape of birth weight specific infant mortality using a second-degree polynomial to account for this unmeasured confounding. This could cause some error if it cannot adequately represent the shape determined by the unmeasured covariates assumed to be responsible for this phenomenon (Figure 1). A 2^{nd} degree polynomial, however, is a relatively flexible function, and is considered to provide an optimal fit to birth weight specific mortality in the homogeneous case [30].

Moreover, the CDDmlr model corrects for some unmeasured confounding of birth weight and infant mortality, referred to as "normal" versus "compromised" births. It is unlikely that dividing birth cohorts into two Gaussian subpopulations will account for all of the unmeasured confounding between birth weight and infant mortality. Nevertheless, the two subpopulations display significantly different mortality patterns indicating that the CDDmlr model accounts for some otherwise unmeasured heterogeneity [9, 10]. In particular, we have argued that the generally higher "normal" birth weight specific mortality compared to "compromised" birth weight specific mortality is due to greater fetal loss among "compromised" births, resulting in a highly selected "compromised" sample at live birth [9, 10] similar to the hypothesis concerning the "pediatric paradox". If correct, this effect would violate assumption c, unless the two subpopulations are examined separately, as they are here.

The statistical results presented above (Tables 2 and 5) are consistent with the Wilcox-Russell hypothesis [2, 4, 5], and its extensions [6] (Figure 1c) that suggest that birth weight is not on the "causal pathway" to infant mortality at least for "normal" births. The racial disparity in birth weight has no significant association with the racial disparity in infant mortality after controlling for the other paths in Figure 1c. There is no evidence of any residual difference in infant mortality between birth weight and infant mortality over and above the direct effect and the reverse-J shape of the standard population, European American births in this case. It is unlikely that this result is compromised by uncontrolled confounding of birth weight and infant mortality, since this would require that the sum total of associations generated by uncontrolled confounding equal zero. It is more likely that all of the effects of race on infant mortality in this subpopulation operate through pathways that do not include birth weight.

On the other hand, there is a substantial indirect effect, which disadvantages African American infant mortality among "compromised" births (Table 5). The results in Table 2 indicate that this association is largely due to a change in shape of the reverse-J-shaped birth weight specific mortality curve between the races. This could be due to an interaction of race and birth weight on infant mortality, or due to a violation of no unmeasured confounding assumptions b or c. It is also equivalent to the interaction [6] required by Figure 1a and also possible in Figure 1b, both of which require that birth weight be on the "causal pathway" to infant mortality. In any event an association between birth weight and infant mortality can not be excluded, and it remains possible that birth weight has a "causal" effect on infant mortality among these "compromised" births.

Overall, the findings suggest that interventions with respect to birth weight will not reduce racial disparities in mortality among "normal" births, but might reduce them among "compromised" births. Identification of the exact mechanisms and whether birth weight plays a "causal" role conditional on "compromised" birth will require additional analysis, i.e. control of potential confounding. The "compromised" subpopulation accounts for about 29-41% of the observed racial disparity for females and males respectively (Table 4).

If our hypothesis concerning the selection effects of fetal loss on observed racial disparities is correct, then the total racial disparity is higher than observed, and the proportion of the disparity due to the "compromised" subpopulation is larger than observed. The confounding, represented by the mixing proportion, accounts for an additional 17-21% of the observed racial disparity for males and females, respectively (Table 4). Nevertheless, completely eliminating the "compromised" subpopulation would a) reduce both the low and the macrosomic birth weight rates, which are generally associated with elevated infant mortality in both African and European American birth cohorts, b) reduce the size of the racial disparity if direct standardization based on the European American distribution are accepted, c) reduce the size of the disparity yet again if our hypothesis concerning the selection effects of fetal loss in the "compromised" subpopulation is correct and included as a potential bias, but d) still result in a population with a racial disparity of 1.9 and 1.8 for females and males, respectively (Table 5), about the level of the relative risk currently observed in the raw data (2.1 for both sexes, Table 1).

## Conclusions

Our results support the Wilcox-Russell [2, 4, 5] and Hernández-Diaz et al. [6] arguments that birth weight is not on the causal pathway to infant mortality at least among "normal" births. Improvements in birth weight may not necessarily impact infant mortality for these births! However, birth weight cannot be eliminated as a potential cause of infant mortality among a small subpopulation of "compromised" births, generally accounting for less than 10% of the birth cohort. Improvements in birth weight may reduce infant mortality among certain births.

The true racial disparity in infant mortality between African and European American birth cohorts may be obscured by unobserved heterogeneity. This heterogeneity may be due to differential fetal loss, which appears to account for the "pediatric paradox". The true racial disparities may also be obscured by lack of consistently reporting births at below 500 grams in the NCHS linked birth death files.

Part of the racial disparity is due to mixing proportion effects, i.e. a larger number of "compromised" births among African Americans than European Americans. Reducing the disparity in the size of "compromised" births will somewhat reduce racial disparities. If all "compromised" births could be eliminated (i.e. eliminating all possible statistically significant birth weight dependent effects), the racial disparities would decrease slightly (1.9 and 1.8 for females and males, respectively) from the currently observed level (2.1 for both sexes). Therefore, the complete elimination of racial disparities in infant mortality requires the elimination of birth weight independent (i.e. direct) effects, as well as any birth weight dependent (i.e. indirect) effects.

## Declarations

### Acknowledgements

This project is supported by NIEHS R01 HD037405. We'd like to thank Dr. Howard Stratton for his many useful discussions and extensive comments on an early draft. We'd also like to thank Kiersten Fussell for her work in preparing the final draft.

## Authors’ Affiliations

## References

- Mosely WH, Chen LC: An analytical framework for the study of child survival in developing countries. Population and Development Review Supplement to. Edited by: Mosely WH, Chen LC. 1984, Cambridge: Cambridge University Press, 10:
- Wilcox AJ: On the importance - and the unimportance - of birthweight. Int J Epidemiol. 2001, 30 (6): 1233-1241. 10.1093/ije/30.6.1233.View ArticlePubMed
- Wise PH: The anatomy of a disparity in infant mortality. Annu Rev Public Health. 2003, 24: 341-362. 10.1146/annurev.publhealth.24.100901.140816.View ArticlePubMed
- The Analysis of Birth Weight and Infant Mortality: an Alternative Hypothesis. [http://eb.niehs.nih.gov/bwt/index.htm]
- Wilcox AJ, Russell I: Why small black infants have a lower mortality-rate than small white infants - the case for population-specific standards for birth-weight. J Pediatr. 1990, 116 (1): 7-10. 10.1016/S0022-3476(05)81638-5.View ArticlePubMed
- Hernández-Diaz S WA, Schisterman EF, Hernán MA: From causal diagrams to birth weight-specific curves of infant mortality. 2008, 23 (3): 163-166.
- US-DOHHS: Healthy People 2010: Understanding and improving health. 2000, Washington, DC: US Department of Health and Human Services, 91-95.
- Gage TB, Fang F, O'Neill E, Stratton H: Maternal age and infant mortality: a test of the Wilcox-Russell hypothesis. Am J Epidemiol. 2009, 169 (3): 294-303. 10.1093/aje/kwn308.PubMed CentralView ArticlePubMed
- Gage TB: Birth-weight-specific infant and neonatal mortality: effects of heterogeneity in the birth cohort. Hum Biol. 2002, 74 (2): 165-184. 10.1353/hub.2002.0020.View ArticlePubMed
- Gage TB, Bauer MJ, Heffner N, Stratton H: Pediatric paradox: heterogeneity in the birth cohort. Hum Biol. 2004, 76 (3): 327-342. 10.1353/hub.2004.0045.View ArticlePubMed
- Wilcox AJ: Infant-mortality among blacks and whites. N Engl J Med. 1992, 327 (17): 1243-1243. 10.1056/NEJM199210223271714.View ArticlePubMed
- Gage TB, Therriault G: Variability of birth-weight distributions by sex and ethnicity: Analysis using mixture models. Hum Biol. 1998, 70 (3): 517-534.PubMed
- Gage TB: Classification of births by birth weight and gestational age: An application of multivariate mixture models. Ann Hum Bio. 2003, 30 (5): 589-604. 10.1080/03014460310001592678.View Article
- Cassady G, Strange M: The small-for-gestational-age (SGA) infant. 1987, Philadelphia: J. B. Lippincott
- Godfrey KM, Barker DJP: Fetal nutrition and adult disease. Am J Clin Nutr. 2000, 71 (5): 1344-1352.
- Basso O, Wilcox AJ, Weinberg CR: Birth weight and mortality: causality or confounding?. Am J Epidemiol. 2006, 164 (4): 303-311. 10.1093/aje/kwj237.View ArticlePubMed
- Basso O, Wilcox AJ: Intersecting birth weight-specific mortality curves: solving the riddle. Am J Epidemiol. 2009, 169 (7): 787-797. 10.1093/aje/kwp024.PubMed CentralView ArticlePubMed
- Geneletti S: Identifying direct and indirect effects in a non-counterfactual framework. J R Stat Soc Ser B. 2007, 69 (2): 199-215. 10.1111/j.1467-9868.2007.00584.x.View Article
- Schisterman EF, Whitcomb BW, Mumford SL, Platt RW: Z-scores and the birthweight paradox. Paediatr Perinat Epidemiol. 2009, 23 (5): 403-413. 10.1111/j.1365-3016.2009.01054.x.PubMed CentralView ArticlePubMed
- Bates DM, Chambers JM: Nonlinear models. Statistical Models in S. Edited by: Chambers JM, Hastie TJ. 1992, Pacific Grove, CA: Wadsworth and Brooks, 421-454.
- Gupta PD: A general method of decomposing a difference between two rates into several components. 1978, 15 (1): 99-112.
- Sappenfield WM, Buehler JW, Binkin NJ, Hogue CJR, Strauss LT, Smith JC: Differences in neonatal and postneonatal mortality by race, birth weight, and gestational age. Public Health Rep. 1987, 102 (2): 182-192.PubMed CentralPubMed
- US-DOHHS: State Definitions and Reporting Requirements for Live Births, Fetal Deaths, and Induced Terminations of Pregnancy. 1997, Washington, DC: US Department of Health and Human Services
- Pearl J: Causal inference in statistics: an overview. Stat Surv. 2009, 3: 96-146. 10.1214/09-SS057.View Article
- VanderWeele TJ: Marginal structural models for the estimation of direct and indirect effects. Epidemiology. 2009, 20 (1): 18-26. 10.1097/EDE.0b013e31818f69ce.View ArticlePubMed
- Robins JM, Hernán M, Brumback B: Marginal structural models and causal inference in epidemiology. Epidemiology. 2000, 11 (5): 550-560. 10.1097/00001648-200009000-00011.View ArticlePubMed
- Buck GM, Shelton JA, Mahoney MC, Michalek AM, Powell EJ: Racial variation in spontaneous fetal deaths at 20 weeks or older in upstate New York, 1980-86. Public Health Rep. 1995, 110 (5): 587-592.PubMed CentralPubMed
- Kallan JE: Race, intervening variables, and 2 components of low-birth-weight. Demography. 1993, 30 (3): 489-506. 10.2307/2061653.View ArticlePubMed
- Platt RW, Joseph KS, Ananth CV, Grondines J, Abrahamowicz M, Kramer MS: A proportional hazards model with time-dependent covariates and time-varying effects for analysis of fetal and infant death. Am J Epidemiol. 2004, 160 (3): 199-206. 10.1093/aje/kwh201.View ArticlePubMed
- Fryer JG, Hunt RG, Simons AM: Biostatistical considerations: the case for using models. Prevention of Perinatal Mortality and Morbidity. Edited by: Falkner F. 1984, Basel, Switzerland: Karger, 9-30.
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2393/10/86/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.