Research article | Open | Published:
Mediation analysis of gestational age, congenital heart defects, and infant birth-weight
BMC Research Notesvolume 7, Article number: 926 (2014)
In this study we assessed the mediation role of the gestational age on the effect of the infant’s congenital heart defects (CHD) on birth-weight.
We used secondary data from the Baltimore-Washington Infant Study (1981–1989). Mediation analysis was employed to investigate whether gestational age acted as a mediator of the association between CHD and reduced birth-weight. We estimated the mediated effect, the mediation proportion, and their corresponding 95% confidence intervals (CI) using several methods.
There were 3362 CHD cases and 3564 controls in the dataset with mean birth-weight of 3071 (SD = 729) and 3353 (SD = 603) grams, respectively; the mean gestational age was 38.9 (SD = 2.7) and 39.6 (SD = 2.2) weeks, respectively. After adjusting for covariates, the estimated mediated effect by gestational age was 113.5 grams (95% CI, 92.4–134.2) and the mediation proportion was 40.7% (95% CI, 34.7%–46.6%), using the bootstrap approach.
Gestational age may account for about 41% of the overall effect of heart defects on reduced infant birth-weight. Improved prenatal care and other public health efforts that promote full term delivery, particularly targeting high-risk families and mothers known to be carrying a fetus with CHD, may therefore be expected to improve the birth-weight of these infants and their long term health.
Congenital heart defects (CHD), also known as congenital cardiovascular malformations, are the most common type of structural birth defects and have substantial effects on pediatric mortality and morbidity [1, 2]. The prevalence of CHD in the live-born population of the United States is estimated to be 1 out of 125 births [3, 4], and approximately 40,000 infants are born annually with CHD . In addition, the 2009 annual summary of vital statistics in the United States showed that CHDs were the leading cause of all birth defect-related infant deaths .
CHDs originate in the early weeks of pregnancy when the heart is forming and are usually diagnosed either prenatally or soon after birth, although some defects are not diagnosed until later in childhood or even in adulthood. Among the most consistently observed findings in newborns with CHD are reduced birth weight [3, 7–11] and increased incidence of clinically low and very low birth weight infants, i.e. <2,500 g and <1,500 g, respectively [3, 10, 12, 13]. An additional highly consistent finding is that infants with CHD are often born prematurely, suggesting that the association of CHD with low birth weight may be attributed, at least in part, to prematurity or reduced gestational age [8, 14–18]. Whether CHD, low birth weight, and prematurity are all co-outcomes of a common, underlying teratogenic influence such as maternal environmental exposure is plausible but as yet unknown [19–21].
The aim of our study was to investigate the relationship between the CHD status of infants and their birth weight by considering gestational age as a mediating variable (mediator), i.e., a variable that “mediates” the relation between the predictor (an infant’s CHD status) and the outcome (birth weight). Using statistical mediation analysis, we considered gestational age as the mediator, we assessed its indirect effect on the relationship of interest, and we estimated the proportion of the effect of CHD status on infant birth weight accounted for by this mediator, taking into account infant and maternal covariates known to affect birth weight and CHD status.
The Baltimore-Washington Infant Study was a population-based case–control study of the potential risk factors of CHD, conducted from 1981 to 1989 in Maryland, the District of Columbia, and Northern Virginia. The study enrolled all live-born infants with CHD and a randomly selected sample of controls that was representative of all regional live births. Controls were randomly sampled at birth hospitals in the region during the same time period. Parents of cases and controls were interviewed at home by trained personnel within one year of the infant’s birth and asked questions concerning medical and demographic/socioeconomic factors, as well as questions regarding cigarette smoking and alcohol consumption. Mothers reported the infant’s birth weight in pounds and ounces, which we converted to grams for this analysis. Gestational age was calculated from the mother’s report of the infant’s date of birth in relation to her reported expected delivery date [1, 13]. The diagnosis of each case was made during the first year of life and confirmed by pediatric cardiologists using echocardiography, cardiac catheterization, surgery, or autopsy. Detailed information about the Baltimore-Washington Infant Study dataset as well as the study in general has been reported previously [1, 22, 23]. The study was approved by the Institutional Review Board of the University of Maryland, Baltimore.
From the extensive Baltimore-Washington Infant Study dataset, we selected variables that are known from prior reports to influence both infants’ CHD status and birth weight, as potential confounders for our study. These included an infant’s sex [24–26] and race [27–29], as well as maternal age [27, 30, 31], education , pre-existing maternal diabetes [32, 33], parity [30, 34, 35], and cigarette smoking during the first trimester [19, 31, 32, 35]. We excluded twins and other multiple births due to their much smaller birth weights relative to singletons and their high prevalence of premature births.
In mediation analysis, the main interest is to assess how much the effect of an exposure on a response variable is mediated through an intermediate variable, oftentimes called the mediator . Thus, the intermediate variable lies on the causal pathway between the exposure (independent variable) and the response (dependent variable). Based on the statistical mediation approach proposed by Baron and Kenny , the mediation analysis is carried out by partitioning the total effect of the independent variable on the dependent variable into its mediated (or indirect) effect and its direct effect . The mediated effect is the effect of the independent variable on the dependent variable accounted for by the mediator, whereas the direct effect is the effect of the independent variable on the dependent variable holding the mediator constant. Imai et al.  and MacKinnon  operationalized the procedure using a regression-based approach by proposing the implementation of three separate, but related, regression models. In the first model (Model 1), the dependent variable is regressed on the independent variable; in the second model (Model 2), the mediator is regressed on the independent variable; and in the third model (Model 3), the dependent variable is regressed on both the mediator and the independent variable. An attractive feature of this regression-based approach is that additional potential confounders and other terms (e.g., interactions between covariates) can be included.
After obtaining parameter estimates from these three regression models, we estimated the mediated effect as the product of the estimate of the effect of the independent variable on the mediator from Model 2 with the estimate of the effect of the mediator on the dependent variable from Model 3 [37, 40–43]. The 95% confidence intervals for the mediation effect were calculated and compared using three different approaches: (1) the Wald statistic approach, based on the Sobel test ; (2) the likelihood ratio test approach , which is based directly on the likelihood function rather than on the quadratic approximation of the log-likelihood function used by the Wald statistic approach; and (3) the bootstrap approach [43, 45–47].
From an epidemiological point of view, in order to explore the mechanism of exposure-disease relations, it is important to assess the mediated (indirect) effect. A relevant way of expressing mediated effects is by considering the “proportion mediated”, which is the proportion of the total effect explained by a particular mediator [38–40, 48]. The mediation proportion was estimated by taking the ratio of the estimate of the mediated effect to the estimate of the total effect (i.e., the effect from Model 1). We used and compared three different approaches  to calculate the 95% confidence intervals for the mediation proportion: (1) the Delta method  to compute the standard error of the mediation proportion so as to construct the confidence intervals for the mediation proportion; (2) the Fieller’s method ; and (3) the bootstrap approach [45, 51–53]. The results from the bootstrap approach were compared, for the purpose of justification, with the results from the other two approaches. During the bootstrap approach, we used 500 draws (i.e., bootstrap resamples) for point estimation as well as construction of confidence intervals. The mediation analyses were performed using SAS 9.3 (SAS Institute Inc., Cary, NC) and R (R Foundation for Statistical Computing, Vienna, Austria).
The parents of all participants signed an informed consent agreement to participate in the original Baltimore-Washington Infant Study protocol, which was approved by the Instututional Review Board of the University of Maryland, Baltimore. The present analysis was performed on data abstracted without identifiers.
The Baltimore-Washington Infant Study enrolled a total of 3,377 cases and 3,572 controls from 1981 to 1989. Fifteen cases and 8 controls were dropped from the analyses due to missing values. Characteristics of the study population are shown in Table 1. The average birth weights of the cases and controls were 3071 (Standard deviation = 729) and 3353 (Standard deviation = 603) grams, respectively (p < 0.0001), confirming the expected deficit among the infants with CHD. The average gestational ages were found to be 38.9 (standard deviation = 2.7) and 39.6 (standard deviation = 2.2) weeks, respectively, and the difference was statistically significant (p < 0.0001). Infant sex, maternal education, and maternal cigarette smoking were distributed similarly in cases and controls (p > 0.05), whereas the case–control differences in infant race, maternal age, previous pregnancies, pre-pregnancy diabetes, and method of delivery were each statistically significant. Based on these results we explored these variables as potential covariates in the subsequent statistical models.
The parameter estimates, standard errors, and p- values from the three regression models (Model 1, Model 2, and Model 3) are shown in Table 2. During the model building steps, interaction terms between the variables were considered, but none of them were found to be statistically significant (p > 0.10; where p = 0.10 was used as statistical cut-point for the inclusion of interaction terms into the final models). Race of the infant, maternal age, education and parity were not statistically significant in any of the models considered and were not included in the final models. Two covariates – sex of the infant and maternal smoking status – were associated with infant birth weight, with or without gestational age in the model. In Model 1, the CHD status had a statistically significant detrimental effect on birth weight (p < 0.0001, Table 2(a)). Adjusted for sex of infants and maternal smoking status, on average, infants with CHD had a birth weight 279 g lower (95% confidence interval, 248–310) than that of infants without CHD. The CHD status also had a small but statistically significant association with gestational age in Model 2 (p < 0.0001, Table 2(b)). When CHD status and gestational age of infants were both considered as predictors of the birth weight of the infants (Model 3), both variables were found to be statistically significant (p < 0.0001 each, Table 2(c)). Adjusted for gender, maternal smoking, and gestational age, the difference between the birth weight of infants with CHD and infants without CHD decreased to 165.9 g (95% confidence interval, 141–190).
Analysis of the mediated effect and mediation proportion
From these three regression models, we observed that the total effect of CHD status on infant birth weight was 279 g (Table 2(a)). The direct effect of the CHD status on birth weight was 166 g (Table 2(c)), and the mediated effect was 164.13 × 0.69 = 113 g (95% confidence interval, 93.6–132.6 based on the Wald statistic approach, and 99.9–117.8 based on the likelihood-based approach; Table 3(a)). The point estimate of the mediated effect found from the bootstrap approach was 113.5 g, which is closely comparable to the estimates found from Wald statistic and likelihood-based approaches (95% confidence interval, 92.4–134.2, Table 3(a)). The 95% confidence interval from the likelihood-based approach was narrower compared to the Wald statistic and bootstrap approaches.
The estimated mediation proportion was 40.5% using the Wald statistic (Delta) and Fieller’s methods (95% confidence intervals, 32.2%–48.8% and 35.2%–45.9%, respectively), and 40.7% (95% confidence interval, 34.7%–46.6%) using the bootstrap method (Table 3(b)). In other words, adjusted for gender and maternal smoking, approximately 41% of the effect of the CHD status on birth weight was mediated by the gestational age of the infant, and the remaining 59% was a direct effect of CHD status on birth weight.
We also investigated the mediation proportion for the two sexes separately. The mediation proportion was 45.2% (95% Fieller’s confidence interval, 36.6%–54.0%) for males and 36.7% (95% Fieller’s confidence interval, 29.8%–43.5%) for females; the difference in these two proportions was not statistically significant (p = 0.3342).
In studying the association between an independent and dependent variable, there are situations in which the independent variable affects the outcome variable both directly and indirectly, by influencing a third variable (the mediator) which in turn influences the dependent variable. Thus, the relationships between variables are often complex and cannot be well modeled by a single regression equation, while consideration of a mediator allows researchers to more fully explore the separate and combined effects of interest. The problem of CHD as a potential cause of reduced birth weight in affected infants is amenable to this approach because birth weight itself is highly dependent on gestational age, which in turn may be influenced by CHD to some extent. Consistent with this view, several studies have previously documented the association of CHD in newborns with their decreased birth weight and gestational age, relative to unaffected infants, but to our knowledge, the mediation effect has not been directly assessed in prior published research [1, 13]. Additionally, we were able to account for the effects of additional covariates of birth weight in our models: males had consistently higher birth weights than females on average in this study, while maternal smoking resulted in decreased birth weight in a dose-dependent manner. Both effects were independent of gestational age. Caesarian section, a possible indicator for a planned inducement of early delivery for a fetus known or suspected to have birth defects, was more common among cases than controls. While it reduced the mean gestational age of cases by 0.5 weeks (p = 0.002) and had a lesser affect on controls (−0.3 weeks, p = 0.99), the use of caesarian section did not significantly affect the mean birth weights of infants in this study.
This study has several strengths and limitations. The data is from the Baltimore-Washington Infant Study, which had achieved complete ascertainment of cases in the region and a population-based random sample of controls; therefore selection bias is extremely unlikely as an explanation for our findings. On the other hand, we cannot rule out the possibility of recall bias of infant birth weight, nor random (non-systematic) errors in recalling the weight or the estimated date of delivery. There are also statistical issues that should be considered in evaluating our results. We confined our analysis to the linear effects of variables on birth weight, but future work should consider possible non-linear effects as well. While our use of regression models permitted us to assess these mediation effects, the technique can become unwieldy for more complex situations [40, 44, 54]. For example, the original Baltimore-Washington Infant Study dataset contains a large number of variables that could be explored jointly as mediators and/or covariates. Furthermore, it would be interesting to consider a common teratogen model, in which a factor such as maternal smoking was modeled as affecting both the occurrence of the malformation and contributing simultaneously to premature birth and decreased weight at birth.
In summary, we used data from a population-based epidemiological study to evaluate the mediation effect of gestational age in explaining the relationship between CHD and decreased birth weight of infants. The results suggest that this mediation accounts for a substantial proportion (about 41%) of the latter association, which was not apparent in previous research in this field due to the limitations of conventional statistical association methods. Further research is needed to confirm these findings in other populations and to explore even more complex methods such as structural equation modeling [55, 56].
Ferencz C, Rubin JD, Loffredo CA, Magee CM: The Epidemiology of Congenital Heart Disease: The Baltimore-Washington Infant Study (1981–1989). 1993, Mount Kisco, NY: Futura Publishing Company
Nieminen HP, Jokinen EV, Sairanen HI: Causes of late deaths after pediatric cardiac surgery: a population-based study. J Am Coll Cardiol. 2007, 50: 1263-1271. 10.1016/j.jacc.2007.05.040.
Archer JM, Yeager SB, Kenny MJ, Soll RF, Horbar JD: Distribution of and mortality from serious congenital heart disease in very low birth weight infants. Pediatrics. 2011, 127: 293-299. 10.1542/peds.2010-0418.
Tanner K, Sabrine N, Wren C: Cardiovascular malformations among preterm infants. Pediatrics. 2005, 116: e833-e838. 10.1542/peds.2005-0397.
Centers for disease control and prevention: five facts about congenital heart defects. [http://www.cdc.gov/ncbddd/heartdefects/specificdefects.html]
Kochanek KD, Kirmeyer SE, Martin JA, Strobino DM, Guyer B: Annual summary of vital statistics: 2009. Pediatrics. 2012, 129: 338-348. 10.1542/peds.2011-3435.
Oberg S, Cnattingius S, Sandin S, Lichtenstein P, Iliadou AN: Birth weight predicts risk of cardiovascular disease within dizygotic but not monozygotic twin pairs: a large population-based co-twin-control study. Circulation. 2011, 123: 2792-2798. 10.1161/CIRCULATIONAHA.110.987339.
Wren C, Irving CA, Griffiths JA, O’Sullivan JJ, Chaudhari MP, Haynes SR, Smith JH, Hamilton JR, Hasan A: Mortality in infants with cardiovascular malformations. Eur J Pediatr. 2012, 171: 281-287. 10.1007/s00431-011-1525-3.
Khalil A, Aggarwal R, Thirupuram S, Arora R: Incidence of congenital heart disease among hospital live births in India. Indian Pediatr. 1994, 31: 519-527.
Kramer MS: Birthweight and infant mortality: perceptions and pitfalls. Paediatr Perinat Epidemiol. 1990, 4: 381-390. 10.1111/j.1365-3016.1990.tb00661.x.
Levin DL, Stanger P, Kitterman JA, Heymann MA: Congenital heart disease in low birth weight infants. Circulation. 1975, 52: 500-503. 10.1161/01.CIR.52.3.500.
Godfrey M, Schimmel MS, Hammerman C, Farber B, Glaser J, Nir A: The incidence of congenital heart defects in very low birth weight and extremely low birth weight infants. Isr Med Assoc J. 2010, 12: 36-38.
Rosenthal GL, Wilson PD, Permutt T, Boughman JA, Ferencz C: Birth weight and cardiovascular malformations: a population-based study. the Baltimore-Washington Infant study. Am J Epidemiol. 1991, 133: 1273-1281.
Honein MA, Kirby RS, Meyer RE, Xing J, Skerrette NI, Yuskiv N, Marengo L, Petrini JR, Davidoff MJ, Mai CT, Druschel CM, Viner-Brown S, Sever LE, National Birth Defects Prevention Network: The association between major birth defects and preterm birth. Matern Child Health J. 2009, 13: 164-175. 10.1007/s10995-008-0348-y.
Malik S, Cleves MA, Zhao W, Correa A, Hobbs CA, National Birth Defects Prevention Study: Association between congenital heart defects and small for gestational age. Pediatrics. 2007, 119: e976-e982. 10.1542/peds.2006-2742.
Laas E, Lelong N, Thieulin AC, Houyel L, Bonnet D, Ancel PY, Kayem G, Goffinet F, Khoshnood B, EPICARD Study Group: Preterm birth and congenital heart defects: a population-based study. Pediatrics. 2012, 130: e829-e837. 10.1542/peds.2011-3279.
Cedergren MI, Kallen BA: Obstetric outcome of 6346 pregnancies with infants affected by congenital heart defects. Eur J Obstet Gynecol Reprod Biol. 2006, 125: 211-216. 10.1016/j.ejogrb.2005.07.006.
Wren C: Prematurity, low birth weight, and cardiovascular malformations. Pediatrics. 2011, 127: 385-386. 10.1542/peds.2010-3345.
Alverson CJ, Strickland MJ, Gilboa SM, Correa A: Maternal smoking and congenital heart defects in the Baltimore-Washington Infant Study. Pediatrics. 2011, 127: e647-e653. 10.1542/peds.2010-1399.
Malik S, Cleves MA, Honein MA, Romitti PA, Botto LD, Yang S, Hobbs CA, National Birth Defects Prevention Study: Maternal smoking and congenital heart defects. Pediatrics. 2008, 121: e810-e816. 10.1542/peds.2007-1519.
Kallen K: Maternal smoking and congenital heart defects. Eur J Epidemiol. 1999, 15: 731-737. 10.1023/A:1007671631188.
Loffredo CA: Epidemiology of cardiovascular malformations: prevalence and risk factors. Am J Med Genet. 2000, 97: 319-325. 10.1002/1096-8628(200024)97:4<319::AID-AJMG1283>3.0.CO;2-E.
Wilson PD, Loffredo CA, Correa-Villasenor A, Ferencz C: Attributable fraction for cardiac malformations. Am J Epidemiol. 1998, 148: 414-423.
Stevenson DK, Verter J, Fanaroff AA, Oh W, Ehrenkranz RA, Shankaran S, Donovan E, Wright L, Lemons J, Tyson J, Korones S, Bauer C, Stoll B, Papile L: Sex differences in outcomes of very low birthweight infants: the newborn male disadvantage. Arch Dis Child Fetal Neonatal Ed. 2000, 83: F182-F185. 10.1136/fn.83.3.F182.
Storms MR, Van Howe RS: Birthweight by gestational age and sex at a rural referral center. J Perinatol. 2004, 24: 236-240. 10.1038/sj.jp.7211065.
Brothwood M, Wolke D, Gamsu H, Benson J, Cooper D: Prognosis of the very low birthweight baby in relation to gender. Arch Dis Child. 1986, 61: 559-564. 10.1136/adc.61.6.559.
Swamy GK, Edwards S, Gelfand A, James SA, Miranda ML: Maternal age, birth order, and race: differential effects on birthweight. J Epidemiol Community Health. 2012, 66: 136-142. 10.1136/jech.2009.088567.
Mathews TJ, Minino AM, Osterman MJ, Strobino DM, Guyer B: Annual summary of vital statistics: 2008. Pediatrics. 2011, 127: 146-157. 10.1542/peds.2010-3175.
Conley D, Bennett NG: Race and the inheritance of low birth weight. Soc Biol. 2000, 47: 77-93.
Vohr BR, Tyson JE, Wright LL, Perritt RL, Li L, Poole WK, NICHD Neonatal Research Network: Maternal age, multiple birth, and extremely low birth weight infants. J Pediatr. 2009, 154: 498-503. 10.1016/j.jpeds.2008.10.044.
Nobile CG, Raffaele G, Altomare C, Pavia M: Influence of maternal and social factors as predictors of low birth weight in Italy. BMC Public Health. 2007, 7: 192-10.1186/1471-2458-7-192.
Bener A, Salameh KM, Yousafzai MT, Saleh NM: Pattern of maternal complications and low birth weight: associated risk factors among highly Endogamous Women. ISRN Obstet Gynecol. 2012, 540495:
Seghieri G, Anichini R, De Bellis A, Alviggi L, Franconi F, Breschi MC: Relationship between gestational diabetes mellitus and low maternal birth weight. Diabetes Care. 2002, 25: 1761-1765. 10.2337/diacare.25.10.1761.
Shah PS, Knowledge Synthesis Group on Determinants of LBW/PT births: Parity and low birth weight and preterm birth: a systematic review and meta-analyses. Acta Obstet Gynecol Scand. 2010, 89: 862-875. 10.3109/00016349.2010.486827.
Phung H, Bauman A, Nguyen TV, Young L, Tran M, Hillman K: Risk factors for low birth weight in a socio-economically disadvantaged population: parity, marital status, ethnicity and cigarette smoking. Eur J Epidemiol. 2003, 18: 235-243.
Preacher KJ, Kelley K: Effect size measures for mediation models: quantitative strategies for communicating indirect effects. Psychol Methods. 2011, 16: 93-115.
Baron RM, Kenny DA: The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol. 1986, 51: 1173-1182.
Hafeman DM: “Proportion explained”: a causal interpretation for standard measures of indirect effect?. Am J Epidemiol. 2009, 170: 1443-1448. 10.1093/aje/kwp283.
Imai K, Keele L, Tingley D: A general approach to causal mediation analysis. Psychol Methods. 2010, 15: 309-334.
MacKinnon DP: Introduction to Statistical Mediation Analysis. 2008, New York: Taylor & Francis
Freedman LS, Graubard BI, Schatzkin A: Statistical validation of intermediate endpoints for chronic diseases. Stat Med. 1992, 11: 167-178. 10.1002/sim.4780110204.
Buyse M, Molenberghs G: Criteria for the validation of surrogate endpoints in randomized experiments. Biometrics. 1998, 54: 1014-1029. 10.2307/2533853.
MacKinnon DP, Fairchild AJ, Fritz MS: Mediation analysis. Annu Rev Psychol. 2007, 58: 593-614. 10.1146/annurev.psych.58.110405.085542.
Sobel ME: Asymptotic confidence intervals for indirect effects in structural equation models. Sociological Methodology. Edited by: Leinhart S. 1982, San Francisco: Jossey-Bass, 290-312.
Preacher KJ, Hayes AF: SPSS and SAS procedures for estimating indirect effects in simple mediation models. Behav Res Methods Instrum Comput. 2004, 36: 717-731. 10.3758/BF03206553.
Zhang Z, Wang L: Methods for evaluating mediation effects: Rationale and comparison. New Trends in Psychometrics. Edited by: Shigemasu K, Okada A, Imaizumi T, Hoshino T. 2008, Tokyo: Universal Academy Press, 595-604.
Shrout PE, Bolger N: Mediation in experimental and nonexperimental studies: new procedures and recommendations. Psychol Methods. 2002, 7: 422-445.
Freedman LS: Confidence intervals and statistical power of the ‘Validation’ ratio for surrogate or intermediate endpoints. J Stat Plan Inference. 2001, 96: 143-153. 10.1016/S0378-3758(00)00330-X.
Casella G, Berger RL: Statistical Inference. 2002, Belmont, CA: Duxbury Press, 240-245. 2
Beyene J, Moineddin R: Methods for confidence interval estimation of a ratio parameter with application to location quotients. BMC Med Res Methodol. 2005, 5: 32-10.1186/1471-2288-5-32.
Efron B: Better bootstrap confidence intervals. J Am Stat Assoc. 1987, 82: 171-185. 10.1080/01621459.1987.10478410.
Carpenter J, Bithell J: Bootstrap confidence intervals: when, which, what? a practical guide for medical statisticians. Stat Med. 2000, 19: 1141-1164. 10.1002/(SICI)1097-0258(20000515)19:9<1141::AID-SIM479>3.0.CO;2-F.
Haukoos JS, Lewis RJ: Advanced statistics: bootstrapping confidence intervals for statistics with “difficult” distributions. Acad Emerg Med. 2005, 12: 360-365. 10.1111/j.1553-2712.2005.tb01958.x.
Brown RL: Assessing specific sediational effects in complex theoretical models. Struct Equ Model. 1997, 4: 142-156. 10.1080/10705519709540067.
Iacobucci D, Saldanha N, Deng JX: A meditation on mediation: evidence that structural equations models perform better than regressions. J Consum Psychol. 2007, 17: 140-154.
Little TD, Card NA, Bovaird JA, Bovaird KJ, Crandall CS: Structural equation modeling of mediation and moderation with contextual factors. Modeling Contextual Effects in Longitudinal Studies. Edited by: Little TD, Bovaird JA, Card NA. 2007, Mahwah, NJ: Lawrence Erlbaum Associates, 207-230.
The authors declare that they have no competing interests.
All of the authors participated in writing the manuscript. C.L. designed the study. Data analysis was performed by AFW, IB, and GL. All authors read and approved the final manuscript.