Skip to main content

On the use of the outcome variable “small for gestational age” when gestational age is a potential mediator: a maternal asthma perspective

A Correction to this article was published on 29 January 2018

This article has been updated

Abstract

Background

The variable “small for gestational age,” frequently defined as birth weight below the 10th percentile in a gestational age and sex-normalized population, is nowadays generally perceived as a more adequate measure than birth weight or low birth weight (birth weight < 2500 g) to capture fetal growth. However, the use of small for gestational age rather than birth weight or low birth weight as an outcome (dependent) variable may have important impacts on the interpretation of analyses aimed at estimating the causal effect of an exposure of interest on infants. We hypothesized potential differences in both types of effects estimated (direct or total) and in ability to control for confounding bias.

Methods

We first examined the use of outcome variables birth weight and small for gestational age to get insights on modeling practices within the field of maternal asthma. Using directed acyclic graph simulations where gestational age was a potential mediator, we then compared estimated exposure effects in regression models for birth weight, low birth weight, and small for gestational age. Graphs with and without confounding were considered.

Results

Our simulations showed that the variable small for gestational age captures the direct effect of exposure on birth weight, but not the indirect effect of exposure on birth weight through gestational age. Interestingly, exposure effect estimates from small for gestational age models were found unbiased whenever exposure effect estimates from birth weight models were affected by collider bias due to conditioning on gestational age in the models.

Conclusions

The sole consideration of the outcome small for gestational age in a study may lead to suboptimal understanding and quantification of the underlying effect of an exposure on birth weight-related measures. Instead, our results suggest that both outcome variables (low) birth weight and small for gestational age should minimally be considered in studies investigating perinatal outcomes.

Peer Review reports

Background

Perinatal outcomes birth weight (BW) and low birth weight (LBW; BW < 2500 g) have a long history of use in public health and medical studies [1]. Since the past several decades, there has been an increased awareness and understanding of the limitation of these variables to convey notions of prematurity and fetal growth [1,2,3]. Although no numeric cut-offs were proposed at that time, the concept of being “small for gestational age” (SGA) can be traced back to the 1960s [1]. Nowadays, the variable SGA, frequently defined as BW below the 10th percentile in a gestational age (GA) and sex-normalized population [4, 5], is generally accepted as a more adequate measure than BW or LBW to characterize intrauterine growth [2, 6]. Indeed, (L)BW can be viewed as a heterogeneous variable influenced by two distinct processes, GA and fetal growth, which can complicate the interpretation of study results. However, while much of recent focus is put on the epidemiology of preterm birth and SGA (e.g., [1, 7]), the use of SGA rather than (L)BW itself as an outcome (dependent) variable may have analytical consequences which can pose difficulties if they are not well known or understood.

In this work we have investigated the practical implications, both in terms of types of estimated effects and potential biases, of using different birth weight-related parametrizations (that is either (L)BW or SGA) as outcome variables when assessing the effect of exposures in analyses which use these outcomes. As mentioned previously, the rationale for not exclusively using (L)BW in studies is now well understood. Although SGA is perceived as a more interpretable outcome than (L)BW since internally adjusted for GA, it remains that this variable is a mere statistical construct (e.g., [8]) and could also present some limitations.

Direct and indirect exposure effects are intuitive concepts which are at the heart of mediation analyses [9,10,11]. Fundamentally, one often desires to decompose the total effect of an exposure on an outcome in one effect that is mediated by an intermediate variable (indirect effect) and one effect that does not arise through that variable (direct effect). The potential for GA to lie in the causal pathway between an exposure of interest (e.g., inhaled corticosteroids (ICS) for treating asthmatic pregnant women) and BW is clearly evident. Indeed, any effect of the exposure on GA necessarily entails an indirect effect of the exposure on BW because of the strong causal association between GA and BW. Using mediation ideas, we have conceived a simulation study with the goal to shed further light on the advantages and disadvantages of considering either (L)BW or SGA as outcome variable of interest in exposure effect analyses. Our objective is to emphasize the statistical implications of using standard modelling approaches for variables that could formally be cast into mediation models.

Using directed acyclic graphs (DAGs) [12, 13] where GA is a potential mediator between the exposure to ICS and BW, we first examined different scenarios wherein the effect of ICS on the outcome was either direct, indirect or both. Because SGA is strictly a function of BW and GA, we then assessed the corresponding interpretations for this outcome. Using mediation DAGs with confounders, we also examined different scenarios which could bias similarly or dissimilarly the estimates of the exposure-(L)BW and exposure-SGA associations. To our knowledge, no comprehensive simulation-based study has yet brought evidence concerning the differential implications of using either parametrization in the presence of confounding biases.

The paper is divided as follows. We first offer a glance at practice in maternal asthma research with respect to the use of (L)BW, SGA and GA as outcome variables. While one can argue that these are relatively crude measures of fetal and subsequent infant health, these outcomes are widely used in this specific research area. Then we describe the two sets of simulated scenarios considered, the processes used to generate the data, and the analyses performed. The presentation of the results and a discussion conclude this work.

Method

A glance at practice in maternal asthma research

To get insights about the choice of (L)BW and SGA as outcome variables of interest, we selected two relatively recent articles which performed a meta-analysis or a systematic review on populations of asthmatic pregnant women. Our selection is not deemed exhaustive but rather insightful of preferred habits in the use of (L)BW and SGA as outcome variables in maternal asthma research.

The first article, authored by Murphy et al. [14], presents a meta-analysis of adverse perinatal outcomes in women with asthma. The meta-analysis is based on cohort studies published between 1975 and 2009 and for which the effect of asthma was assessed for at least one of the following outcomes: LBW, SGA, and preterm birth (PTB; GA < 37 weeks). Table 1 presents the selection of these three outcomes by included study; this table was constructed from the studies reported in relative risk Figures 1–3 in Murphy et al. [14].

Table 1 Meta-analysis by Murphy et al. [14]: a summary of selected outcome variables by included study

The second article, by Eltonsy et al. [15], presents a systematic review for the use of β2-agonists during pregnancy and their effects on perinatal outcomes. This systematic review searched for articles published before 2013 regarding the effects of β2-agonists on congenital malformations, SGA, mean and low BW, GA, and PTB. Table 2 presents the selection of the three types of outcomes investigated by included study; this table was constructed from the studies reported in Tables 4–8 in Eltonsy et al. [15].

Table 2 Systematic review by Eltonsy et al. [15]: summary of selected outcome variables by included study

We see that the GA-related variables (GA, PTB) most often appear in the studies listed in Tables 1 and 2. Thus, when only two variables were considered in a study, it usually included a GA-related variable. In such a situation, the other covariate selected was seen to split in a larger proportion for (L)BW as opposed to SGA (8/12 vs 4/12, respectively; studies counted once only). In the meta-analysis results (see Table 1), 5 out of 18 studies (27.8%) reported all three outcomes, and 4 out of 10 (40%) for the systematic review (see Table 2). Three studies, all in the meta-analysis from Murphy et al. [14], only reported either LBW or SGA (without GA).

No clear time trend with respect to the patterns of inclusion of the variables is seen in Tables 1 and 2. More recent articles also show differential preferences regarding the use of variables (L)BW, SGA, and GA/PTB. For instance, in [16] and [17], the authors used all these three variables as outcomes in analyses, while [18] and [19] only used SGA and GA, and LBW and GA, respectively. Although the use of (L)BW without SGA (and vice versa) is not widespread based on these articles and the selected meta-analysis and systematic review, it is frequent enough to establish the relevance of investigating (L)BW and SGA jointly.

In the next section, we introduce a series of DAGs to help interpret estimated causal effects of an exposure on (L)BW and SGA. At the same time, we also incorporate the GA dimension since the constructed variable SGA is a function of both BW and GA.

Design of primary simulations

We first present basic scenarios and corresponding data generation processes. These were used to gain insights on the interpretation of the associations between the exposure and the outcomes of interest ((L)BW, SGA). Then, to investigate the potential for confounding or collider bias [20, 21] when using (L)BW or SGA, we pursue with more complex scenarios that incorporate a single confounder of the relationships between the exposure, GA, and BW. Collider bias can be described as a spurious association between the exposure and the outcome which arises when studies, at the design or the analysis stage, stratify or adjust on a collider; in a DAG, a collider is a common effect of two variables lying on a path linking the exposure and the outcome [11]. The simulations were performed using R [22], version 3.0.2. We used the same initial seed to generate the data for all scenarios.

Basic scenarios

We considered four basic causal DAGs with three nodes: the exposure node ICS, the mediator node GA, and the outcome node BW. Figure 1 presents the DAGs showing the posited links between ICS, GA, and BW. Each of the DAGs corresponds to one scenario and is the basis for the generation of the corresponding variables of interest. First, a causal link between GA and BW is assumed in all four scenarios. Basic Scenario 1 depicts the case where ICS has no effect on BW, either direct or indirect. Basic Scenario 2 is a scenario in which the effect of ICS on BW is fully mediated by GA; in other words, the effect of ICS on BW exclusively occurs through a modification in GA. Basic Scenario 3 represents the situation where ICS only has a direct effect on BW; in this case, the effect of ICS on BW does not occur through a modification in GA. Finally, in Basic Scenario 4, ICS has a direct effect on BW, in addition to having an indirect effect mediated through GA. It can be noted that there are two scenarios which do not feature GA as a mediator of the relationship between ICS and BW (Basic Scenarios 1 and 3). These were considered as a benchmark for interpreting the results under scenarios in which GA is a mediator (Basic Scenarios 2 and 4).

Fig. 1
figure 1

Directed acyclic graphs for the Basic Scenarios 1–4. From top to bottom, Basic Scenario 1: Null effect of ICS on BW; Basic Scenario 2: Indirect effect of ICS on BW; Basic Scenario 3: Direct effect of ICS on BW; Basic Scenario 4: Direct and indirect effects of ICS on BW

Generation of outcomes (L)BW and SGA (basic scenarios)

For each basic scenario, we simulated data for 20,000 babies, where this process was repeated 1000 times to constitute 1000 samples of size 20,000. Exposures to ICS were generated from Bernoulli experiments with probability 0.5 (ICS = 1 if exposed and ICS = 0 if unexposed).

For each sample in the basic scenarios, the GAs for all babies were initially generated from a multinomial distribution with support between 26 and 43 weeks and probabilities corresponding to the estimated probability of delivering at each of these weeks according to Table 1 from Kramer et al. [4]. More precisely, the GAs were generated with theoretical probabilities coinciding with the empirical probabilities found in that table (that is, number of pregnancies observed for a given GA value divided by the total number of pregnancies). To generate the GAs for exposed babies in scenarios in which ICS had a direct effect on GA (Basic Scenarios 2 and 4), we subtracted 2 weeks from the GAs that were generated in the first place. The magnitude of this effect was chosen to well illustrate the concepts described herein and is likely larger than the real effect of ICS on GA. Indeed, low-to-moderate doses of ICS are generally regarded as safe regarding GA- and BW-related outcomes while some uncertainty remains regarding the effects of larger doses of ICS [23].

The BWs (in grams) were then generated independently according to a normal distribution with mean

$$ {\mu}_{BW}=-3703.3+183.25\ GA+{\beta}_{ICS}\ ICS $$
(1)

and standard deviation σ BW  = 333.82. The value of β ICS in Eq. (1) was set to −100 or 0, depending on whether ICS had a direct effect on BW or not. The values of the intercept and GA coefficients were defined on the basis of the data found in Table 1 from Kramer et al. [4]. Specifically, these two coefficients had been a priori calculated by fitting a linear regression model for BW versus GA on a large simulated sample (n = 100,000), where the GAs had also been generated according to tabulated empirical probabilities and the BWs generated from a normal distribution according to the GA-specific mean and standard deviation values found in that table. The value 333.82 for σ BW corresponds to the average of the standard deviation values found in the table and is substantially smaller than the residual standard error returned by the aforementioned large sample regression analysis (450.80). This value for σ BW was selected so that less variability is observed for the BW distributions conditional on the smallest GA values. Although the assumptions of a linear effect of GA on BW and of common variance of errors are not satisfied in the data summarized in Table 1 from Kramer et al. [4], these were made in the primary simulations for simplicity.

We created the SGA variables for all babies based on their values for GA and BW. The binary variable SGA was determined by comparing a baby’s BW to the 10th percentile of the normal BW distribution conditional on GA. For example, an unexposed baby born at 28 weeks of gestation was found small for his GA (SGA = 1) if his BW was less than 999.89 g and not small for his GA otherwise (SGA = 0), where 999.89 is the 10th percentile of a normal distribution with mean equal to −3703.3 + 183.25 × 28 = 1427.7 and standard deviation equal to 333.82. An exposed baby born at 28 weeks of gestation was also said small for his GA if his BW was less than 999.89 g. A common BW threshold was thus used to determine the SGA value of every baby born at the same GA. Finally, a baby was said having LBW if his BW was smaller than 2500 g.

Confounding scenarios

Next, we considered four additional DAGs with four nodes: a dichotomous confounder node V and the same three nodes as before (ICS, GA, BW). Figure 2 presents the DAGs depicting the causal links between V, ICS, GA, and BW in these confounding scenarios. All the DAGs feature a causal effect of ICS on BW fully mediated by GA. The DAGs differ by the posited relationships between V and the nodes ICS, GA, and BW. In Confounding Scenario 1, V is a confounder between ICS and GA, while V is a confounder between GA and BW in Confounding Scenario 2. In Confounding Scenario 3, V is a confounder between ICS and BW. Confounding Scenario 4 encompasses all previous scenarios as V is a common cause of ICS, GA, and BW simultaneously.

Fig. 2
figure 2

Directed acyclic graphs for the Confounding Scenarios 1–4. From top to bottom, Confounding Scenario 1: Indirect effect of ICS on BW with V confounder between ICS and GA; Confounding Scenario 2: Indirect effect of ICS on BW with V confounder between GA and BW; Confounding Scenario 3: Indirect effect of ICS on BW with V confounder between ICS and BW; Confounding Scenario 4: Indirect effect of ICS on BW with V common confounder between ICS, GA, and BW

Generation of outcomes (L)BW and SGA (confounding scenarios)

We also generated 1000 samples according to each of the four confounding scenarios. Each sample included 20,000 babies, among which half (10,000) had the value V = 1 and half had the value V = 0. Whenever there was an arrow from V to ICS in a DAG (Confounding Scenarios 1, 3, 4), the probability of being exposed to ICS was set to 0.7 for those with V = 1 and to 0.3 for those with V = 0. Otherwise (Confounding Scenario 2), the probability of being exposed to ICS was 0.5, independently of the value of V. Exposures to ICS were then generated according to these probabilities. Of note, the former values of 0.70 and 0.30 were selected to induce significant confounding arising through covariate V and exposure ICS. In general, in our simulations, we have allowed for strong relationships between variables to more comprehensively interpret the results.

For each sample, the GAs for the babies with V = 0 were generated from the multinomial distribution previously described. In Confounding Scenario 3, in which V had no direct effect on GA, the generated GA values for the babies with V = 1 were taken to be the same as those for the babies with V = 0. These values constituted the intermediate GA values for all babies. To create the GAs for the babies with V = 1 in scenarios in which V had a direct effect on GA, we subtracted 2 weeks from the GAs that were generated for the babies with V = 0. Lastly, to generate the final GA values for the exposed babies in scenarios in which ICS had a direct effect on GA, we subtracted 2 weeks from each exposed baby’s intermediate GA. For all other babies, their intermediate GA value was set to their final GA value.

The BWs (in grams) were generated independently according to a normal distribution with mean

$$ {\mu}_{BW}=-3703.3+183.25\ GA+{\beta}_{ICS}\ ICS+{\beta}_V\ V $$
(2)

and standard deviation σ BW  = 333.82. The value of β ICS (β V ) in Eq. (2) was equal to - 100 or 0, depending on whether ICS (V) had a direct effect on BW or not. The variables SGA and LBW were created as in the Basic Scenarios 1–4.

Primary analyses (all scenarios)

On each sample generated according to each of the four basic scenarios, we fitted four linear or logistic models (M), depending or SGA was modeled:

$$ {M}_1:\kern0.5em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS; $$
$$ {M}_2:\kern0.5em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_{GA}\ GA; $$
$$ {M}_3:\mathrm{logit}\left(P\left( LBW=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS; $$
$$ {M}_4:\mathrm{logit}\left(P\left( SGA=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS. $$

For each model under each basic scenario, we computed the mean of the 1000 estimated ICS coefficients, \( {\widehat{\beta}}_{ICS} \), where each coefficient represents either a mean difference or a log odds ratio. The standard error of the mean of \( {\widehat{\beta}}_{ICS} \) was calculated as the standard deviation of the 1000 estimated ICS coefficients divided by \( \sqrt{1000} \).We also assessed whether the Monte Carlo 95% normal confidence interval for β ICS included zero.

Similarly, we fitted the following eight models on each sample generated according to each of the four confounding scenarios:

$$ {M}_1:\kern1.25em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS; $$
$$ {M}_{1,V}:\kern0.75em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_V\ V; $$
$$ {M}_2:\kern1.5em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_{GA}\ GA; $$
$$ {M}_{2,V}:\kern0.75em E\left[ BW\right]={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_{GA}\ GA+{\beta}_V\ V; $$
$$ {M}_3:\kern1.25em \mathrm{logit}\left(P\left( LBW=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS; $$
$$ {M}_{3,V}:\kern0.75em \mathrm{logit}\left(P\left( LBW=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_V\ V; $$
$$ {M}_4:\kern1.25em \mathrm{logit}\left(P\left( SGA=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS; $$
$$ {M}_{4,V}:\kern0.75em \mathrm{logit}\left(P\left( SGA=1\right)\right)={\beta}_0+{\beta}_{ICS}\ ICS+{\beta}_V\ V. $$

We calculated the mean of the 1000 estimates \( {\widehat{\beta}}_{ICS} \) obtained for each model under each confounding scenario, with corresponding standard error. Again, we assessed whether the confidence interval for β ICS included zero. As an aid for interpretation, we also indicated whether the estimators of the ICS effect on (L)BW (\( {\widehat{\beta}}_{ICS} \) in models M1 − 3, M1 − 3, V) are biased according to the DAGs considered (e.g., see Pearl [24] for graphical causal rules). Whether or not the estimators of the ICS effect on SGA (\( {\widehat{\beta}}_{ICS} \) in models M4, and M4, V) are biased was also determined; in this case our insight relied on the estimates obtained from the simulations under the unconfounded basic scenarios.

Sensitivity analyses (all scenarios)

In the primary simulations, we made the assumptions of a linear effect of GA on BW and of common variance of BW errors. We subsequently decided to create additional simulation scenarios to examine the consequence of fitting the standard regression models previously described when these assumptions were not verified in the data.

We assessed the potential impact of the heterogeneity of BW errors on the basic and confounding scenarios’ results. Instead of simulating the BWs with constant errors (σ BW  = 333.82 for all GA), we simulated them according to the GA-specific standard deviation values found in our perinatal reference table (that is, Table 1 from Kramer et al. [4]). For instance, we used σ BW  = 241.00 for a baby born at 28 weeks and σ BW  = 447.00 for a baby born at 40 weeks. Modification to the way that the SGA variable was created was done accordingly. More precisely, the 10th percentile cut-off at a given GA was determined using the GA-specific standard deviation value used for simulating the BWs at that GA. No changes were made regarding the linearity of the effect of GA on BW in these new analyses; that is, Eq. (1) and (2) still apply for describing the mean BW formula used for generating this outcome.

We also assessed the potential impact of the nonlinearity of the effect of GA on BW on the basic and confounding scenarios’ results. To do that, we took an approach similar to what was done to calculate the intercept and GA coefficient values in Eq. (1). However, instead of fitting a simple linear model on the large sample with GA as a single explanatory variable, we considered a model with additional quadratic and cubic GA terms. The coefficients returned by this model were used to specify the equations used to generate the BW values for these sensitivity analyses. For both the basic scenarios and the confounding scenarios with nonlinearity of GA effect on BW, we thus had

$$ {\mu}_{BW}=29415-2922.4\ GA+94.46\ {GA}^2-0.9378\ {GA}^3+\varphi \left( ICS,V\right) $$

where φ(ICS, V) is the last term in Eq. (1) for the basic scenarios, and φ(ICS, V) is the sum of the last two terms in Eq. (2) for the confounding scenarios; the standard deviation of errors remained the same at σ BW  = 333.82. For these analyses, we also fitted either one or two additional models \( \left({M}_2^{\ast },{M}_{2,V}^{\ast}\right) \) where polynomial terms in GA up to degree 3 were included as regressors.

Results

Primary analysis results

The results for the basic and confounding scenarios for the primary simulations are presented in Tables 3 and 4, respectively.

Table 3 Basic scenarios without confounding
Table 4 Confounding scenarios for indirect ICS effect on BW

For Basic Scenario 1 (Table 3), all mean estimated ICS effects on (L)BW or SGA are close to zero, as is expected since ICS has no effect, either direct or indirect, on BW. For Basic Scenario 2, in which ICS has an indirect effect on BW, the ICS effect in the GA-unadjusted models M1 and M3 is significantly different from zero. In model M1, this mean estimate reflects the diminution of 2 weeks in GA for babies exposed to ICS, which in turn decreases the average BW by (−2) × 183.25 = −366.50 g (see Eq. (1)). Further, we observe that ICS is not associated with SGA (M4), nor is it with BW when we condition on GA (M2). For Basic Scenario 3, all mean estimated ICS effects on (L)BW or SGA are significantly different from zero. The mean estimated ICS effect on BW is close to −100 in both models M1 and M2, and reflects the direct effect of ICS on BW (see Eq. (1)). In that scenario, being exposed to ICS was found to increase the odds of having a LBW baby by exp (0.3283) = 1.3886 and the odds of having a SGA baby by exp (0.5611) = 1.7526. In Basic Scenarios 1 and 3, a reduction in variability is observed for the ICS estimates obtained from model M2 as opposed to model M1. In these scenarios, the mean estimated ICS effect is the same for both models M1 and M2, but the effect is more accurately estimated from the model including GA (M2). For Basic Scenario 4, the mean estimated ICS effect on BW in M1 and M2 are different, as in Basic Scenario 2. In this case, the estimated ICS effect on BW in M1 represents both the direct and indirect effects specified in simulation (−100–366.50 = −466.50 g), while the ICS effect on BW observed from M2 only reflects the direct effect. The estimated effect of ICS on SGA in Basic Scenario 4 is the same as in Basic Scenario 3. Together, the pairs of Basic Scenarios (1, 2) and Basic Scenarios (3, 4) thus reveal that the strength of the association between the exposure and SGA is driven by the size of the direct effect of the exposure on BW.

For interpreting the results from the confounding scenarios (Table 4), we take as reference the results obtained for Basic Scenario 2. In Confounding Scenario 1, in which V is a confounder between ICS and GA, only the estimated ICS effects from model M1 and M3 are biased. These results are interesting since they show that common causes of exposure and GA, and not only common causes of exposure and BW, create bias in the exposure-BW relationship when the (L)BW model do not adjust for them. In such a fully mediated scenario, this situation arises because of the presence of the open back-door path ICS ← V → GA → BW which creates a spurious association between ICS and BW; this path becomes closed when conditioning on GA. Still in Confounding Scenario 1, the mean estimate in model M3,V (1.3296) is somewhat different from the mean estimate for M3 in the Basic Scenario 2 (1.2676). This difference is attributed to the well-known non-collapsibility of the odds ratios, where, on this scale, marginal effects are generally different than conditional effects [25]. In Confounding Scenario 2, in which V is a confounder between GA and BW, the estimated ICS effects from model M2 are biased. This bias arises since adjusting for GA in the BW model opens the back-door path ICS → GA ← V → BW which creates a spurious association between ICS and BW; including V in the model closes the path and thus eliminates the bias (see corresponding result for model M2,V). In Confounding Scenario 3, in which V is a confounder between ICS and BW, all estimated ICS effects from non V-adjusted models (i.e., models M1, M2, M3,M4) are biased. Without surprise, all estimated ICS effects from non V-adjusted models are also biased in Confounding Scenario 4 (V confounder for all ICS, GA, and BW). Overall, there are thus two confounding scenarios that affect the estimates differently depending on whether we use (L)BW or SGA as outcome variables. Unlike (L)BW, our results indicate that SGA is not affected by confounders between the exposure and GA, nor by those between GA and BW. In particular, it appears that using this outcome variable prevents the collider bias problem that was seen under model M2 for the conditional effect of ICS on BW in Confounding Scenario 2.

Sensitivity analysis results

The results for the sensitivity analyses regarding the heterogeneity of errors in basic and confounding scenarios are found in Tables 5 and 6, respectively.

Table 5 Basic scenarios without confounding with heterogeneity of errors
Table 6 Confounding scenarios for indirect ICS effect on BW with heterogeneity of errors

In the basic scenarios with heterogeneity of errors (Table 5), the same interpretation as in the basic scenarios with homogeneity of errors can be done regarding the type of effect (total, direct) estimated in the different (L)BW and SGA models.

In the confounding scenarios with heterogeneity of errors (Table 6), the estimates obtained behaved similarly to those obtained under the confounding scenarios with homogeneity of errors from one scenario to the other. One notable difference is with regard to the non-null effect for some SGA models in Confounding Scenarios 2–4. Specifically, Models M4, M4,V in Confounding Scenarios 2 and M4,V in Confounding Scenarios 3 and 4 all featured a mean ICS estimate very close but significantly different from zero. For example, the mean estimate from M4,V in Confounding Scenarios 2 with heterogeneity of error was 0.0127 and the standard error of the mean was 0.0013, thus only slightly, but significantly, departing from the null. This is opposed to the null effect found in the corresponding M4,V result when assuming homogeneity of errors. A common feature of Confounding Scenarios 2–4, which could be the source of this very small discrepancy in these fully mediated scenarios, is the GA-dependent heterogeneity of effect of V on SGA. Indeed, in these scenarios, the direct effect of the V covariate on the mean BW is non-null and the same for all GA values (−100); however, when assuming heterogeneity of errors, this effect will translate differently across GA for SGA since the BW standard deviation, which intervenes in the definition of the SGA variable, varies across GA.

The results for the sensitivity analyses regarding the nonlinear GA effect on BW in basic and confounding scenarios are found in Tables 7 and 8.

Table 7 Basic scenarios without confounding with nonlinear GA effect on BW
Table 8 Confounding scenarios for indirect ICS effect on BW simulated with nonlinear GA effect on BW

With a nonlinear GA effect on BW, changes in results were observed for Basic Scenarios 2 and 4 which both feature an indirect effect of GA on BW (Table 7). For these scenarios, the effect obtained under model M2, which adjusts for a linear GA term only, did not unbiasedly represent the direct effect of exposure on BW. A positive residual bias of about 40 g was observed, but this bias vanished when additionally including the quadratic and cubic GA terms in the model (see results for \( {M}_2^{\ast } \) in these scenarios).

Compared to all previous scenarios, the interpretation of results for the confounding scenarios under nonlinearity of GA effect on BW is more complicated (Table 8). It is noted that even BW models adjusted for V (M1,V) did not estimate the total effect of exposure completely without bias when the confounder was associated with GA (Confounding Scenarios 1, 2, 4). However, model M1,V did unbiasedly estimate the total effect of exposure under Confounding Scenarios 3 (refer to M1,V in Basic Scenario 2, Table 7, for comparison), scenario in which the confounder is associated with exposure and outcome only. Nonetheless, adjusting for V in the (L)BW models M1,V and M3,V yielded estimates overall closer to the total effect of exposure, as expected. The unbiased estimation of the direct effect of exposure in the BW models was achieved under model \( {M}_{2,V}^{\ast } \), where both the correct functional form of GA and the confounder was accounted for. In these sensitivity analyses, no changes in results and interpretation were observed for SGA as compared with the Confounding Scenarios in the primary analyses (Table 4).

Discussion

Our study found that the sole consideration of SGA in a study may lead to suboptimal understanding and quantification of the underlying effect of an exposure on BW-related measures. Using DAGs where GA was a potential mediator between the exposure and BW, we have confirmed that SGA is an absorbing variable: the observed association between the exposure and SGA solely reflects the direct effect of the exposure on BW, effect which could be interpreted as a manifestation of intrauterine growth retardation. In the situation where the effect of exposure on BW was fully mediated by GA, the exposure and SGA were not causally linked. Therefore, an analyst may have concluded for no exposure effect on the BW-axis, but what should really be concluded is the absence of a direct effect of the exposure on BW. While one could argue this is precisely the purpose of using SGA, we believe that being able to assess the total (direct and indirect) effect of the exposure on BW is at least as valuable.

If the exposure has a direct beneficial effect but an overall detrimental effect, the direct effect has less importance. However, precise quantification of the total effect is obtained with the use of BW as outcome variable. From our perspective, it would thus not be advisable to only consider SGA and GA as outcome variables since the ability to precisely assess the indirect and total effects of exposure on BW would be diminished. Although the presence of an indirect effect of the exposure on BW could be deduced from results looking at the effect of the exposure on GA, our study suggests that considering BW as outcome variable is the most straightforward way to thoroughly investigate this issue. Indeed, recall that the strength of an indirect effect of an exposure on an outcome is a combination of two measures: 1) the strength of the association between the exposure and the mediator (GA) and 2) the strength of the association between the mediator and the outcome (BW) [26].

One advantage we found with respect to the use of SGA as outcome variable, as opposed to BW, is that is less prone to bias. Interestingly, exposure effect estimates from SGA models were found unbiased whenever exposure effect estimates from BW models were affected by collider bias due to conditioning on GA in the model. Therefore, accounting for GA internally (through the use of a GA-adjusted BW measure) or externally (by conditioning on GA in a model for BW) are two competing strategies that are not equally robust to bias. This finding is particularly important when confounders for the GA and BW association are unmeasured and cannot be adjusted for in the model for BW. In our study, only the SGA model, and not the BW model, yielded unbiased results interpretable as a direct effect of exposure on BW in the absence of such a confounder in the model. We also found that whenever one missing covariate was a common cause between exposure and BW, the exposure-BW association as well as the exposure-SGA association were biased. Therefore, one should not make the distinction between these outcome variables when selecting such confounders and adjusting for them in models. However, we found that in fully mediated sets-up, common causes between exposure and GA biased the exposure-(L)BW relationship unlike the exposure-SGA relationship when a model with the former variable (BW) did not adjust for GA. Distinguishing between common causes of exposure and GA and common causes of exposure and (L)BW could, however, be a rather difficult task in practice.

Conclusions

In light of our simulations and current wisdom, we recommend that, in addition to GA, both outcome variables (L)BW and SGA be considered in studies that rely on these perinatal outcomes. Alternatively to considering all three outcomes (GA, SGA, and (L)BW) in standard separate analyses, mediation analyses with GA as a mediator could be used to better understand the direct and indirect effects of an exposure on BW. When adopting a mediation strategy, the use of SGA could be omitted as one would be able to make the distinction between an exposure effect on (L)BW arising through a diminution in GA and one external to this mechanism. However, such a mediation model would nevertheless be inadequate to provide unbiased results in the presence of unmeasured mediator-outcome confounders [11] and could also be sensitive to linearity assumptions.

Globally, our study has highlighted the complexity of perinatal outcome modeling. Although our findings are directly relevant to the field of maternal asthma, we believe they are applicable to other research areas or specific types of studies where (L)BW and SGA, as measures of perinatal health, are the most useful and feasible. As a matter of fact, these variables have been recently considered in meta-analyses (e.g., [27]), large clinical trials in less-resourced countries (e.g., [28]), and large-scale perinatal studies based on administrative databases (e.g., [29]).

Change history

  • 29 January 2018

    Following publication of the original article [1], the authors reported that the following four references in Table 2 are incorrect:

Abbreviations

BW:

Birth weight

DAG:

Directed acyclic graph

GA:

Gestational age

ICS:

Inhaled corticosteroids

LBW:

Low birth weight

PTB:

Preterm birth

SGA:

Small for gestational age

References

  1. Hughes MM, Black RE, Katz J. 2500-g low birth weight cutoff: history and implications for future research and policy. Matern Child Health J. 2017;21(2):283–9.

    Article  PubMed  Google Scholar 

  2. Urquia ML, Ray JG. Seven caveats on the use of low birthweight and related indicators in health research. J Epidemiol Community Health. 2012;66(11):971–5.

    Article  PubMed  Google Scholar 

  3. Wilcox AJ. On the importance--and the unimportance--of birthweight. Int J Epidemiol. 2001;30(6):1233–41.

    Article  CAS  PubMed  Google Scholar 

  4. Kramer MS, Platt RW, Wen SW, Joseph KS, Allen A, Abrahamowicz M, et al. A new and improved population-based Canadian reference for birth weight for gestational age. Pediatrics. 2001;108(2):E35.

    Article  CAS  PubMed  Google Scholar 

  5. Campbell MK, Cartier S, Xie B, Kouniakis G, Huang W, Han V. Determinants of small for gestational age birth at term. Paediatr Perinat Epidemiol. 2012;26(6):525–33.

    Article  PubMed  Google Scholar 

  6. Savitz DA, Hertz-Picciotto I, Poole C, Olshan AF. Epidemiologic measures of the course and outcome of pregnancy. Epidemiol Rev. 2002;24(2):91–101.

    Article  PubMed  Google Scholar 

  7. Statistics Canada. Perinatal health indicators for Canada 2017. https://www.canada.ca/en/public-health/services/injury-prevention/health-surveillance-epidemiology-division/maternal-infant-health/perinatal-health-indicators-2017.html. Accessed 11 Oct 2017.

    Google Scholar 

  8. Kjøbli E, Bach R, Skogseth H, Jacobsen GW. The Scandinavian Small-for-Gestational Age (SGA) pregnancy and birth cohort – a source to continual insight into fetal growth restriction and long term physical and neurodevelopmental health in mother and offspring. Norw J Epidemiol. 2016;26(1–2):145–54.

    Google Scholar 

  9. MacKinnon DP, Fairchild AJ, Fritz MS. Mediation analysis. Annu Rev Psychol. 2007;58:593–614.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Kenny DA. MEDIATION. http://davidakenny.net/cm/mediate.htm Accessed 11 Oct 2017.

  11. Richiardi L, Bellocco R, Zugna D. Mediation analysis in epidemiology: methods, interpretation and bias. Int J Epidemiol. 2013;42(5):1511–9.

    Article  PubMed  Google Scholar 

  12. Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology. 1999;10(1):37–48.

    Article  CAS  PubMed  Google Scholar 

  13. VanderWeele TJ, Robins JM. Directed acyclic graphs, sufficient causes, and the properties of conditioning on a common effect. Am J Epidemiol. 2007;166(9):1096–104.

    Article  PubMed  Google Scholar 

  14. Murphy VE, Namazy JA, Powell H, Schatz M, Chambers C, Attia J, et al. A meta-analysis of adverse perinatal outcomes in women with asthma. BJOG. 2011;118(11):1314–23.

    Article  CAS  PubMed  Google Scholar 

  15. Eltonsy S, Kettani FZ, Blais L. Beta2-agonists use during pregnancy and perinatal outcomes: a systematic review. Respir Med. 2014;108(1):9–33.

    Article  PubMed  Google Scholar 

  16. Firoozi F, Lemiere C, Ducharme FM, Beauchesne MF, Perreault S, Berard A, et al. Effect of maternal moderate to severe asthma on perinatal outcomes. Respir Med. 2010;104(9):1278–87.

    Article  PubMed  Google Scholar 

  17. Rejno G, Lundholm C, Gong T, Larsson K, Saltvedt S, Almqvist C. Asthma during pregnancy in a population-based study--pregnancy complications and adverse perinatal outcomes. PLoS One. 2014;9(8):e104755.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Hodyl NA, Stark MJ, Scheil W, Grzeskowiak LE, Clifton VL. Perinatal outcomes following maternal asthma and cigarette smoking during pregnancy. Eur Respir J. 2014;43(3):704–16.

    Article  PubMed  Google Scholar 

  19. Mendola P, Laughon SK, Mannisto TI, Leishear K, Reddy UM, Chen Z, et al. Obstetric complications among US women with asthma. Am J Obstet Gynecol. 2013;208(2):127.e1–8.

    Article  Google Scholar 

  20. Elwert F, Winship C. Endogenous selection bias: the problem of conditioning on a collider variable. Annu Rev Sociol. 2014;40(1):31–53.

    Article  Google Scholar 

  21. Cole SR, Platt RW, Schisterman EF, Chu H, Westreich D, Richardson D, et al. Illustrating bias due to conditioning on a collider. Int J Epidemiol. 2010;39(2):417–20.

    Article  PubMed  Google Scholar 

  22. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. Available from: http://www.R-project.org.

    Google Scholar 

  23. Gregersen TL, Ulrik CS. Safety of bronchodilators and corticosteroids for asthma during pregnancy: what we know and what we need to do better. J Asthma Allergy. 2013;6:117–25.

    PubMed  PubMed Central  Google Scholar 

  24. Pearl J. Causality: models, reasoning, and inference. 2nd ed. Cambridge: Cambridge University Press; 2009.

    Book  Google Scholar 

  25. Hernán MA, Clayton D, Keiding N. The Simpson’s paradox unraveled. Int J Epidemiol. 2011;40(3):780–5.

    Article  PubMed  PubMed Central  Google Scholar 

  26. MacKinnon DP, Lockwood CM, Hoffman JM, West SG, Sheets V. A comparison of methods to test mediation and other intervening variable effects. Psychol Methods. 2002;7(1):83–104.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ruiz M, Goldblatt P, Morrison J, Kukla L, Svancara J, Riitta-Jarvelin M, et al. Mother’s education and the risk of preterm and small for gestational age birth: a DRIVERS meta-analysis of 12 European cohorts. J Epidemiol Community Health. 2015;69(9):826–33.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Muhihi A, Sudfeld CR, Smith ER, Noor RA, Mshamu S, Briegleb C, et al. Risk factors for small-for-gestational-age and preterm births among 19,269 Tanzanian newborns. BMC Pregnancy Childbirth. 2016;16:110.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Marzouk A, Filipovic-Pierucci A, Baud O, Tsatsaris V, Ego A, Charles MA, et al. Prenatal and post-natal cost of small for gestational age infants: a national study. BMC Health Serv Res. 2017;17(1):221.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Lao TT, Huengsburg M. Labour and delivery in mothers with asthma. Eur J Obstet Gynecol Reprod Biol. 1990;35(2–3):183–90.

    Article  CAS  PubMed  Google Scholar 

  31. Perlow JH, Montgomery D, Morgan MA, Towers CV, Porto M. Severity of asthma and perinatal outcome. Am J Obstet Gynecol. 1992;167(4 Pt 1):963–7.

    Article  CAS  PubMed  Google Scholar 

  32. Doucette JT, Bracken MB. Possible role of asthma in the risk of preterm labor and delivery. Epidemiology. 1993;4(2):143–50.

    Article  CAS  PubMed  Google Scholar 

  33. Jana N, Vasishta K, Saha SC, Khunnu B. Effect of bronchial asthma on the course of pregnancy, labour and perinatal outcome. J Obstet Gynaecol (Tokyo 1995). 1995;21(3):227–32.

    Article  CAS  Google Scholar 

  34. Schatz M, Zeiger RS, Hoffman CP, Harden K, Forsythe A, Chilingar L, et al. Perinatal outcomes in the pregnancies of asthmatic women: a prospective controlled analysis. Am J Respir Crit Care Med. 1995;151(4):1170–4.

    CAS  PubMed  Google Scholar 

  35. Stenius-Aarniala B, Riikonen S, Teramo K. Slow-release theophylline in pregnant asthmatics. Chest. 1995;107(3):642–7.

    Article  CAS  PubMed  Google Scholar 

  36. Demissie K, Breckenridge MB, Rhoads GG. Infant and maternal outcomes in the pregnancies of asthmatic women. Am J Respir Crit Care Med. 1998;158(4):1091–5.

    Article  CAS  PubMed  Google Scholar 

  37. Liu S, Wen SW, Demissie K, Marcoux S, Kramer MS. Maternal asthma and pregnancy outcomes: a retrospective cohort study. Am J Obstet Gynecol. 2001;184(2):90–6.

    Article  CAS  PubMed  Google Scholar 

  38. Bracken MB, Triche EW, Belanger K, Saftlas A, Beckett WS, Leaderer BP. Asthma symptoms, severity, and drug therapy: a prospective study of effects on 2205 pregnancies. Obstet Gynecol. 2003;102(4):739–52.

    PubMed  Google Scholar 

  39. Dombrowski MP, Schatz M, Wise R, Momirova V, Landon M, Mabie W, et al. Asthma during pregnancy. Obstet Gynecol. 2004;103(1):5–12.

    Article  PubMed  Google Scholar 

  40. Acs N, Puho E, Banhidy F, Czeizel AE. Association between bronchial asthma in pregnancy and shorter gestational age in a population-based study. J Matern Fetal Neonatal Med. 2005;18(2):107–12.

    Article  PubMed  Google Scholar 

  41. Bakhireva LN, Jones KL, Schatz M, Johnson D, Chambers CD. Asthma medication use in pregnancy and fetal growth. J Allergy Clin Immunol. 2005;116(3):503–9.

    Article  CAS  PubMed  Google Scholar 

  42. Sheiner E, Mazor M, Levy A, Wiznitzer A, Bashiri A. Pregnancy outcome of asthmatic patients: a population-based study. J Matern Fetal Neonatal Med. 2005;18(4):237–40.

    Article  PubMed  Google Scholar 

  43. Clark JM, Hulme E, Devendrakumar V, Turner MA, Baker PN, Sibley CP, et al. Effect of maternal asthma on birthweight and neonatal outcome in a British inner-city population. Paediatr Perinat Epidemiol. 2007;21(2):154–62.

    Article  PubMed  Google Scholar 

  44. Enriquez R, Griffin MR, Carroll KN, Wu P, Cooper WO, Gebretsadik T, et al. Effect of maternal asthma and asthma control on pregnancy and perinatal outcomes. J Allergy Clin Immunol. 2007;120(3):625–30.

    Article  PubMed  Google Scholar 

  45. Källén B, Otterblad Olausson P. Use of anti-asthmatic drugs during pregnancy. 2. Infant characteristics excluding congenital malformations. Eur J Clin Pharmacol. 2007;63(4):375–81.

    Article  PubMed  Google Scholar 

  46. Karimi M, Davar R, Mirzaei M, Mirzaei M. Pregnancy outcomes in asthmatic women. Iran J Allergy Asthma Immunol. 2008;7(2):105–6.

    PubMed  Google Scholar 

  47. Breton MC, Beauchesne MF, Lemiere C, Rey E, Forget A, Blais L. Risk of perinatal mortality associated with asthma during pregnancy. Thorax. 2009;64(2):101–6.

    Article  PubMed  Google Scholar 

  48. Schatz M, Zeiger RS, Harden KM, Hoffman CP, Forsythe AB, Chilingar LM, et al. The safety of inhaled beta-agonist bronchodilators during pregnancy. J Allergy Clin Immunol. 1988;82(4):686–95.

    Article  CAS  PubMed  Google Scholar 

  49. Schatz M, Zeiger RS, Harden K, Hoffman CC, Chilingar L, Petitti D. The safety of asthma and allergy medications during pregnancy. J Allergy Clin Immunol. 1997;100(3):301–6.

    Article  CAS  PubMed  Google Scholar 

  50. Alexander S, Dodds L, Armson BA. Perinatal outcomes in women with asthma during pregnancy. Obstet Gynecol. 1998;92(3):435–40.

    CAS  PubMed  Google Scholar 

  51. Olesen C, Thrane N, Nielsen GL, Sorensen HT, Olsen J. A population-based prescription study of asthma drugs during pregnancy: changing the intensity of asthma therapy and perinatal outcomes. Respiration. 2001;68(3):256–61.

    Article  CAS  PubMed  Google Scholar 

  52. Schatz M, Dombrowski MP, Wise R, Momirova V, Landon M, Mabie W, et al. The relationship of asthma medication use to perinatal outcomes. J Allergy Clin Immunol. 2004;113(6):1040–5.

    Article  PubMed  Google Scholar 

  53. Clifton VL, Rennie N, Murphy VE. Effect of inhaled glucocorticoid treatment on placental 11beta-hydroxysteroid dehydrogenase type 2 activity and neonatal birthweight in pregnancies complicated by asthma. Aust N Z J Obstet Gynaecol. 2006;46(2):136–40.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Fonds de recherche Québec-Santé. The sponsor had no role in study design, interpretation of results, in the writing of the article, and in the decision to submit the article for publication.

Availability of data and materials

All codes for simulations are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

GL conceived the study, designed the simulation and analyzed the results. MS and GL performed the literature review. GL wrote the manuscript. MS and GL edited the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Geneviève Lefebvre.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

The original version of this article was revised: The following four references in Table 2 have been corrected: [30], [38], [41] and [43].

A correction to this article is available online at https://doi.org/10.1186/s12874-018-0477-y.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lefebvre, G., Samoilenko, M. On the use of the outcome variable “small for gestational age” when gestational age is a potential mediator: a maternal asthma perspective. BMC Med Res Methodol 17, 165 (2017). https://doi.org/10.1186/s12874-017-0444-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12874-017-0444-z

Keywords