![]() |
|
|
,3

* Department of International Health,
Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD 21205-2103
3 To whom correspondence should be addressed. E-mail: jkatz{at}jhsph.edu.
| ABSTRACT |
|---|
|
|
|---|
KEY WORDS: birth weight micronutrients pregnancy infant mortality Nepal
Many studies use birth weight as an outcome of nutritional or other interventions in pregnancy because it is a measure of intrauterine growth retardation in the absence of gestational age data and also a strong predictor of early infant survival. The usual analytic approaches to the estimation of treatment effects involve calculating the mean difference in birth weight between control and treatment groups. This assumes that treatment shifts the distribution of birth weights by a constant amount. Alternately, researchers examine the relative risk of low or very low birth weight but this compares only the left tails of the 2 distributions, and may not adequately capture changes in the middle or upper end of the distribution.
Previously, we published the effects of daily antenatal supplementation with different combinations of micronutrients on birth weight (1). In that community-based, randomized trial, despite the apparently smaller adjusted birth weight increase with folic acid + iron alone (37 g) compared with a multiple micronutrient that included folic acid + iron (64 g), the percentage reduction in low birth weight (<2500 g) was very similar with folic acid + iron (16%) compared with the multiple micronutrients (14%), suggesting that the birth weight distributions of the treatment groups may have had different shapes from that of the control group, especially at the lower end of the distributions. Furthermore, the multiple micronutrient supplement increased the proportion of babies with birth weight
3300 g (7.7% compared with 5.3% in the control group), an increase not observed with folic acid + iron supplementation (6.0% compared with 5.3%) (1). Thus, one issue that may be relevant to randomized trials of nutritional interventions is whether treatments change the shape of the birth weight distribution in a way that may not be reflected when comparing means of treatment and control distributions.
In this paper, we explore these findings further and describe an analytic approach to addressing this issue in trials of nutritional interventions to improve birth weight. We use data from our community-randomized trial that examined the effect of different maternal micronutrient supplements on birth weight in rural Nepal to illustrate this approach.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The trial received ethical approval from the Committee on Human Research of the Johns Hopkins Bloomberg School of Public Health and the Nepal Health Research Council.
Pregnant women were interviewed at the time of enrollment and maternal height, weight, age, date of last menstrual period, parity, smoking history, and other characteristics were recorded. The main outcomes of the study were birth weight and infant survival. Because 95% of births occurred in the home, a female supplement distributor resident in the village reported the birth to a supervisor who dispatched an anthropometrist to the home to obtain "birth weight." The aim was to weigh the infant as soon after birth as possible. Weights were measured using a digital (Seca 727) infant weighing scale, accurate to 2 g. We enrolled a total of 4096 pregnancies that resulted in 4130 live-born infants;
80% of weights were obtained within 72 h.
Because weight is measured at variable times after birth, we wished to construct a statistical model that provided a better estimate of true birth weight using the age of the infant at the time of measurement. In addition,
7% of infants were missing weight measurements; 34% of these infants died very soon after birth (2.4% of all those missing birth weight), before the anthropometrist was able to reach the household. These infants are likely to have been smaller at birth than those who survived beyond the time they were weighed. We developed a measurement error model that allowed us to estimate the weight of the baby at birth for measurements made at variable times after birth and to impute missing measurements. Exploratory analysis of the relation between birth weight and time of the measurement indicated that birth weight dropped
50 g between 24 and 48 h, was constant between 48 and 72 h, and increased linearly after 72 h. We assumed the expected birth weight for baby
at time
had a main effect for the treatment and the vital status, was modeled as a natural cubic spline with knots placed at 24, 48, and 72 h for weights taken between birth and 72 h, varied linearly with time after 72 h and with number of cigarettes smoked, but varied nonlinearly with gestational age and maternal age, weight and height (modeled as natural splines).
Specifically, let
and
be the time and the corresponding birth weight measurement for infant
and let
be an indicator of vital status (alive or died) at time
. The following model was fit:
![]() |
where
is the vector of covariates for the
th woman (gestational age, maternal weight in the first trimester, maternal height, maternal age, and number of cigarettes smoked in the 7 d before the first trimester interview).
was the type of supplement the ith infant's mother received, and
was a cubic spline with
df (3 in this case). The gestational age of the infant was based on date of the last menstrual period from the first trimester interview, or if not available, the date of first positive pregnancy test was used. The imputation assumes that those with missing birth weight had weights similar to those with the same covariates who did have birth weights measured. Treatment group was one covariate in this model because treatment affected birth weight. Including treatment in the model would prevent the imputed values from being biased by the treatment effect (5).
Missing birth weights were imputed using a multiple imputation method (6). Specifically, let
be the predicted birth weight at time
conditioned on the vital status of the infant, supplement group, and maternal covariates. Let
be the estimated residual variance of the regression model. Imputed data sets (n = 50) were created by sampling
from a normal distribution with mean
and SD
for
. If the weight was taken after 72 h, then the birth weight was recalibrated by taking
A small percentage of infants (4%) had missing values for some of the mother's covariates. We imputed missing data in the mother's characteristics with the corresponding mean in the population. Figure 1 denotes the observed (circles), predicted (triangles) and imputed (squares) birth weights at time
for the babies that were measured after (A) and before (B) the 72 h [
].
|
Specifically, we assumed that the treatment effect was a smooth function of the percentiles of the birth weight distribution:
![]() |
where
and
are the quantile functions of the birth weight distribution for treatment and control groups, respectively, and s is a natural cubic spline of the percentile p with
df, which we have chosen to equal 4. To estimate
, we did the following: 1) calculated the percentiles
with
(the smallest number of infants across treatment groups); 2) calculated the corresponding empirical quantiles of the birth weights
under the treatment and control groups, respectively; and 3) smoothed the difference between quantiles
across the percentiles
.
Note that if we set P = 0.5, then our approach reduces to the common method of estimating a treatment effect by comparing the central values of the birth weight distribution (means or medians) between the treatment and control groups. The smoothing of quantile differences across percentiles to improve estimation of the mean difference between 2 outcomes was introduced by Dominici and Cope (7) for continuous and skewed data. This approach was implemented for estimating mean medical expenditures between diseased and nondiseased patients (8). In this paper, we modified this idea for a continuous and symmetric outcome to estimate percentile-specific treatment effects. For more details of this approach, see Dominici et al. (9).
To account for the uncertainty in the imputation of the missing values, we repeated steps 13 separately for 50 imputed data sets. We then calculated the percentile-specific treatment effect and its corresponding total statistical variance by using standard multiple imputation methods (6). More specifically, let
and
be the point estimate and the bootstrap variance of
for the j-imputed data set, respectively. For each
, the overall estimate of the treatment effect and its total variance, denoted by
, were obtained as follows:
![]() |
![]() |
where
![]() |
![]() |
where
![]() |
To test whether there was a constant treatment effect across the distribution, we performed a permutation test. Specifically, for
, we randomly reassigned the birth weights to the 2 treatment groups (each of the 4 active treatment groups against control) and calculated the test statistics
where
. The 1-sided P-value was calculated as the probability that
exceeded the observed test statistics
where
.
Because the correlation of birth weight within clusters that were randomized to different treatments was essentially zero, these analyses did not take into account the cluster randomization (1).
| RESULTS |
|---|
|
|
|---|
|
|
|
with the corresponding 95% bootstrap confidence bands for each treatment group (Fig. 3AD). The circles denote the difference in the empirical quantile functions between the 2 groups (
) as a function of the percentiles. As previously reported, there was a mean treatment effect for the folic acid + iron and multiple micronutrient groups (1). The 1-sided P-values of the permutation tests for the hypothesis of a constant treatment effect across percentiles of the birth weight distribution were 0.10 for folic acid + iron, 0.96 for multiple nutrients, 0.04 for folic acid + iron + zinc, and 0.60 for folic acid. Therefore, for folic acid + iron and for folic acid + iron + zinc, the treatment effect is significantly different across the distribution of birth weight at the 5 and 10% level, respectively.
|
2800 g with a decline in the treatment effect for the larger babies; 2) the multiple micronutrients increased birth weight across the entire distribution of weights; 3) folic acid + iron + zinc increased birth weight only for babies in the middle of the birth weight distribution (between 2400 and 2900 g approximately); and 4) folic acid alone did not significantly increase the birth weight across the entire distribution. | DISCUSSION |
|---|
|
|
|---|
Our data demonstrated that the mean treatment effect with iron folic acid was half that of the mean multiple micronutrient effect but in environments such as rural Nepal, it may be more important to affect the lower than the upper part of the birth weight distribution. In fact, affecting the upper part of the distribution may be harmful to the mother and infant. In our study, the multiple micronutrient supplement was associated with a slightly elevated risk of early infant mortality [relative risk (RR): 1.07, 95% CI: 0.75, 1.58], especially among term births (RR: 1.74, 95% CI: 1.00, 3.04) (2). This was despite the significant 14% reduction in low birth weight. The risk of birth asphyxia as a cause of neonatal mortality also appeared to be higher in the group administered the multiple micronutrient supplement. On the other hand, folic acid + iron was associated with an overall 21% reduction (95% CI: 48% reduction to increase of 20%), and
47% reduction in 3-mo mortality (95% CI: 8, 70%) among preterm births. Given an improvement in birth weight at the lower end of the distribution, this intervention may have produced improved survival overall, whereas the multiple micronutrient appeared to have no effect on survival because deaths averted in the smaller infants were negated by higher mortality at the upper end of the distribution. In a recent study examining a multiple micronutrient supplement relative to a control supplement of iron + folic acid in Janakpur, Nepal, there was a higher mean birth weight of 77 g (95% CI: 24, 130) in the multiple micronutrient group compared with the iron + folic acid group, but higher neonatal mortality in the multiple micronutrient group (RR: 1.53, 95% CI: 0.72, 3.2) (10). Pooled analyses of the Sarlahi and Janakpur studies found higher neonatal mortality (RR: 1.52 (1.03, 2.25) among those in the multiple micronutrient group relative to the folic acid + iron group (2,11).
When examined cross sectionally, our data also found evidence at a population level for a decline in weight of
50 g between 24 and 48 h after birth. This is comparable to the drop in weight seen in a population of newborns in Bangladesh (12). The weights rise thereafter, and from 72 h through 3 mo of age, weights can be modeled as a linear function of age in this population. For those who died before a weight measurement could be obtained, or who survived but were missing weight measurements due to migration, birth weights can be imputed using a variety of maternal and other covariates that are strong predictors of birth weight. Using recalibrated birth weights for those measured at varying times after birth, and imputed weights for those with missing weights, the estimated effect of maternal supplementation was slightly larger for all treatment groups than when the treatment effects were estimated without imputation (Fig. 3 and Table 2). The CI for the treatment effects using imputed values were wider, reflecting the uncertainty associated with the imputed birth weights. Although there is an assumption in imputation that the infants missing birth weight would have a similar weight to infants with similar covariates who were not missing weight, there are more infants missing weight because they died soon after birth, and such infants might have been smaller than those who survived long enough to be weighed. However, we do have birth weights for some infants who later died, and this imputation model does use those infants to help predict the birth weight of those who died before we could weigh them. The birth weight differential predicted from the model between those who died and those who survived is >300 g, although the babies we did weigh who later died may not adequately capture this difference in birth weight. It is also possible that babies who were large relative to the size of their mothers may have died so soon after birth that they were also not weighed. This may be true given a reverse J-shaped risk profile for mortality of these infants by birth weight (13). Thus, there may be some bias, but the direction of the bias is uncertain. However, only 2.4% of missing birth weights were imputed because the child had died before weight could be obtained; therefore, it is likely that the bias from imputing weights for these infants on overall birth weight is small. Given the constraints in obtaining birth weight in settings in which births occur at home, future studies or programs are likely to face the same situation with missing birth weight. Our measurement error model for imputing the missing birth weights and predicting the weights "at birth" for measurements made after the 72 h may be relevant for future studies that examine intervention effects on birth weight.
Because of the strong association between birth weight and survival, which may be an inverse J-shape in this population (13), understanding the effect of a treatment on different parts of the birth weight distribution may provide more insight into the effect of the treatment on both the surrogate (birth weight) and primary outcomes (survival) (14). Such analyses can help explain contradictions in treatment effects on surrogate and primary outcomes, and may help identify the extent to which a treatment effect acts directly on survival or through its effect on birth weight. Direct analyses of treatment effects on mortality by birth weight distribution are currently underway (9).
| FOOTNOTES |
|---|
2 Supported by USAID co-operative agreements HRN-A-00-97-00015-00 and GRA-A-00-03-00019-00, UNICEF Nepal country office, Bill and Melinda Gates Foundation. The study was a joint undertaking between the Johns Hopkins Bloomberg School of Public Health and the Nepal National Society for Comprehensive Eye Care, Kathmandu, Nepal. ![]()
Manuscript received 25 August 2005. Initial review completed 6 October 2005. Revision accepted 6 February 2006.
| LITERATURE CITED |
|---|
|
|
|---|
1. Christian P, Khatry SK, Katz J, Pradhan EK, LeClerq SC, Shrestha SR, Adhikari RK, West KP Jr. Effects of alternative maternal micronutrient supplements on low birth weight in rural Nepal. A double-masked randomized community trial. BMJ. 2003;326:5717.
2. Christian P, West KP Jr, Khatry SK, LeClerq SC, Pradhan EK, Katz J, Shrestha SR, Sommer A. Effects of maternal micronutrient supplementation on fetal loss and infant mortality: a cluster-randomized trial in Nepal. Am J Clin Nutr. 2003;78:1194202.
3. Christian P, Shrestha J, LeClerq SC, Khatry SK, Jiang T, Wagner T, Katz J, West KP Jr. Supplementation with micronutrients beyond iron and folic acid fails to further improve the hematologic status of pregnant women in rural Nepal. J Nutr. 2003;133:34928.
4. Freeman JV, Christian P, Khatry SK, Adhikari RK, LeClerq SC, Katz J, Darmstadt GL. Evaluation of neonatal verbal autopsy using physician review versus algorithm-based cause-of-death assignment in rural Nepal. Paediatr Perinat Epidemiol. 2005;19:32331.[Medline]
5. Raghunathan TE. What do we do with missing data? Some options for analysis of incomplete data. Annu Rev Public Health. 2004;25:99117.[Medline]
6. Rubin DB. Inference and missing data. Biometrika. 1976;63:58192.
7. Dominici F, Cope L, Naiman D, Zeger SL. Smooth quantile ratio estimation (SQUARE). Biometrika. 2005;92:54357.
8. Dominici F, Zeger SL. Smooth quantile ratio estimation with regression: estimating medical expenditures for smoking attributable diseases. Biostatistics. 2005;6:50519.
9. Dominici F, Zeger SL, Parmigiani G, Katz J, Christian P. Estimating percentile-specific effects: a case study of micronutrient supplementation, birth weight, and infant mortality. J R Stat Soc Ser C Appl Stat. 2006;50:120.
10. Osrin D, Vaidya A, Shrestha Y, Baniya RB, Manandhar DS, Adhikari RK, Filteau S, Tomkins A, Costello AM. Effects of antenatal micronutrient supplementation on birthweight and gestational duration in Nepal: double-blind randomized controlled trial. Lancet. 2005;365:95562.[Medline]
11. Christian P, Osrin D, Manandhar DS, Khatry SK, de L Costello AM, West KP Jr. Antenatal micronutrient supplements in Nepal. Lancet. 2005;366:7112.[Medline]
12. Arifeen SE. Birth weight, intrauterine growth retardation and prematurity: a prospective study of infant growth and survival in the slums of Dhaka, Bangladesh [dissertation]. Baltimore (MD): Johns Hopkins University School of Hygiene and Public Health; 1997.
13. Christian P. Effect of maternal micronutrient supplementation on newborn size and infant health and survival. In: Delange FM, West KP Jr, editors. Micronutrient deficiencies in the first months of life. Nestle Nutrition Workshop Series Pediatric Program, Vol 52, Nestle, 2003.
14. Wilcox AJ. On the importanceand the unimportanceof birthweight. Int J Epidemiol. 2001;30:123341.
This article has been cited by other articles:
![]() |
A. C. Lee, G. L. Darmstadt, S. K. Khatry, S. C. LeClerq, S. R. Shrestha, and P. Christian Maternal-Fetal Disproportion and Birth Asphyxia in Rural Sarlahi, Nepal Arch Pediatr Adolesc Med, July 1, 2009; 163(7): 616 - 623. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Roberfroid, L. Huybregts, H. Lanou, M.-C. Henry, N. Meda, J. Menten, P. Kolsteren, and for the MISAME Study Group Effects of maternal multiple micronutrient supplementation on fetal growth: a double-blind randomized controlled trial in rural Burkina Faso Am. J. Clinical Nutrition, November 1, 2008; 88(5): 1330 - 1340. [Abstract] [Full Text] [PDF] |
||||
![]() |
P Christian, G L Darmstadt, L Wu, S K Khatry, S C LeClerq, J Katz, K P West Jr, and R K Adhikari The effect of maternal micronutrient supplementation on early neonatal morbidity in rural Nepal: a randomised, controlled, community trial Arch. Dis. Child., August 1, 2008; 93(8): 660 - 664. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||