 Technical advance
 Open Access
 Published:
Thinking outside the curve, part I: modeling birthweight distribution
BMC Pregnancy and Childbirth volume 10, Article number: 37 (2010)
Abstract
Background
Greater epidemiologic understanding of the relationships among fetalinfant mortality and its prognostic factors, including birthweight, could have vast public health implications. A key step toward that understanding is a realistic and tractable framework for analyzing birthweight distributions and fetalinfant mortality. The present paper is the first of a twopart series that introduces such a framework.
Methods
We propose describing a birthweight distribution via a normal mixture model in which the number of components is determined from the data using a model selection criterion rather than fixed a priori.
Results
We address a number of methodological issues, including how the number of components selected depends on the sample size, how the choice of model selection criterion influences the results, and how estimates of mixture model parameters based on multiple samples from the same population can be combined to produce confidence intervals. As an illustration, we find that a 4component normal mixture model reasonably describes the birthweight distribution for a population of white singleton infants born to heavily smoking mothers. We also compare this 4component normal mixture model to two competitors from the existing literature: a contaminated normal model and a 2component normal mixture model. In a second illustration, we discover that a 6component normal mixture model may be more appropriate than a 4component normal mixture model for a general population of black singletons.
Conclusions
The framework developed in this paper avoids assuming the existence of an interval of birthweights over which there are no compromised pregnancies and does not constrain birthweights within compromised pregnancies to be normally distributed. Thus, the present framework can reveal heterogeneity in birthweight that is undetectable via a contaminated normal model or a 2component normal mixture model.
Background
The impact of birthweight on perinatal mortality and morbidity has been debated for decades [1–11]. Although advances in maternal and perinatal care have reduced overall mortality, infants with very low birthweights (10001500 g; VLBW) and extremely low birthweights (<1000 g; ELBW) remain at high risk. These infants require more intensive utilization of health resources, at increased costs relative to normal birthweight (NBW; 25004000 g) infants [12–14]. Even infants of moderately low birthweight (15002500 g; MLBW) and high birthweight (>4000 g; HBW) have elevated mortality and morbidity [15, 16]. Greater epidemiologic understanding of the relationships among fetalinfant mortality and its prognostic factors, including birthweight, could have vast public health implications. A key step toward that understanding is a realistic yet tractable framework for analyzing birthweight distribution and fetalinfant mortality.
Simple bell curves are inadequate characterizations of birthweight distributions [17, 11, 18–20]. Wilcox and Russell proposed a contaminated normal model, in which a predominant normal distribution accounts for most birthweights while a contaminating residual distribution yields most VLBW and ELBW cases [21]. The residual distribution does not have a specific structure and, in particular, is not normal. The contaminated normal model was later extended by Umbach and Wilcox to accommodate two residual distributions, one yielding excess births in the left tail and the other in the right tail [22].
Gage and Therriault took a different approach, employing a 2component normal mixture model [23]. A primary normal distribution accounts for most birthweights, while a secondary normal distribution is linked not only to most VLBW and ELBW cases but also to many HBW cases. The 2component normal mixture (resp., contaminated normal model) dichotomizes birthweights: those arising from the primary distribution (resp., predominant distribution) are conceptualized as reflecting ordinary fetal development, while the rest are considered to signify compromised fetal development [24]. Gage also formulated a parametric mixtures of logistic regressions (PMLR) technique to evaluate heterogeneity in mortality associated with this dichotomy [24].
While the aforementioned works demonstrate great insight, their statistical models have some limitations. In particular, the number of constituent distributions (predominant, residual, primary, secondary) is fixed a priori. If a constituent distribution can signify compromised fetal development [24], perhaps different biological mechanisms for compromised fetal development warrant a model with more than two or three constituent distributions. Likewise, perhaps more than two or three birthweightspecific mortality curves are needed to describe heterogeneity in mortality.
The present paper is the first in a twopart series that introduces a new framework for modeling birthweight distribution and fetalinfant mortality. We propose a normal mixture model for birthweight distribution in which the number of components is not fixed a priori but rather determined from the data using the Flexible Information Criterion (FLIC) (Pilla and Charnigo, Consistent estimation and model selection in semiparametric mixtures, submitted) or another model selection technique [25, 26]. In the companion paper, we show how to estimate birthweightspecific mortality within each component using a generalization of PMLR [24] and how to compare mortality across components within a single population or across populations within a single component. In both papers, we seek statistical models that provide an empirically reasonable fit to the data. However, the goal is not to find good fitting models for their own sake. Rather, such models may lead to better assessments of mortality.
Results
1. Pragmatics for mixture modeling
a. Finite normal mixture models
Many phenomena cannot be accurately described via a normal distribution. When no other commonly used probability distribution seems appropriate, a finite normal mixture model is often reasonable. We now briefly describe the model. Readers interested in theoretical developments may consult references [27–30] and works cited therein.
Let f(x;μ,σ) denote the probability density for the normal distribution with mean μ and standard deviation σ. A finite normal mixture model with k components has probability density
A common way to interpret Equation (1) is to imagine that the full population consists of k subpopulations. The proportion of individuals in the full population belonging to subpopulation j is p _{ j } . In subpopulation j, measurements are normally distributed with mean μ _{ j } and standard deviation σ _{ j } .
The mixture components may or may not represent subpopulations with obvious biological definitions outside the statistical model. For example, in a 2component normal mixture describing birthweights for white singletons in the United States, there is not an obvious biological characterization for the two components: we may say that the component with the smaller mean reflects compromised pregnancies, but we cannot immediately attribute the compromised pregnancies to a specific biological mechanism.
Ideally, modeling with finite normal mixtures may lead to discoveries of subpopulations with biological definitions that were not immediately obvious, although the mixture components themselves may still only be approximations to such subpopulations.
b. Order selection and the flexible information criterion
Equation (1) may be an imperfect description of real data regardless of k, but with k sufficiently large the description may be adequate to address a problem of scientific interest. Conversely, if k is too large, the model may become unwieldy. Hence, a researcher with real data must confront the problem of "order selection" (i.e., choosing an appropriate number of components).
Let M denote the maximum number of components that a researcher is willing to accept. For 1 ≤ m ≤ M, let L _{ m } denote the maximum value of the likelihood attainable by an mcomponent normal mixture. The Akaike Information Criterion (AIC) [25], Bayesian Information Criterion (BIC) [26], and Flexible Information Criterion (FLIC) (Pilla and Charnigo, Consistent estimation and model selection in semiparametric mixtures, submitted) are
Above, (3m  1) is the number of free parameters in an mcomponent normal mixture. Also, n denotes the sample size, δ the average fraction of withincomponent variability to total variability over the M normal mixtures fitted by maximum likelihood, and B(n,δ) a bivariate function taking values between 0 and 1 (Pilla and Charnigo, Consistent estimation and model selection in semiparametric mixtures, submitted). The criteria balance fidelity to the observed data against model complexity; models are preferred for which the criteria are smaller. Note that m indexes normal mixtures being judged by the criteria, while k pertains to a normal mixture that has been adopted for data analysis.
The FLIC is distinguished from the AIC and BIC in that its penalty term is determined not only by the sample size but also by the configuration of data points: a configuration suggesting greater heterogeneity allows a model with more components to be selected. The penalty term of the FLIC also depends on M, so that a researcher must specify M. In analyzing birthweight data, we fix M = 7 since having too many components would impede inference about mortality risk. The FLIC and AIC perform well for small samples, while the FLIC and BIC are better for large samples, so we prefer to rely on the FLIC (Pilla and Charnigo, Consistent estimation and model selection in semiparametric mixtures, submitted).
c. Computational procedures
To employ the FLIC, we must obtain maximum likelihood estimates of the proportions, means, and standard deviations in all finite normal mixture models under consideration. For models with more than one component, numerical optimization procedures must be used. We apply the expectation maximization (EM) algorithm to obtain preliminary estimates [31], followed by the optimization (optim) procedure in version 2.3.1 of R (R Foundation for Statistical Computing, Vienna, Austria, 2006) to acquire final estimates. Our R code is available upon written request to the corresponding author. See Section I of [Additional file 1] for details on using EM and optim, including initial value specification.
2. Analyzing birthweight data with the FLIC
a. A FLICselected model and competitors
To exemplify use of the FLIC, we draw a random sample of size 50,000 from the 202,849 white singletons who were born (or experienced fetal death) from 2000 to 2002 and whose mothers smoked heavily (at least twenty cigarettes per day). Since records with birthweights less than 500 grams or gestational ages less than 22 weeks were not consistently documented [32], we require infants in our sample to have known gestational ages of at least 22 weeks and birthweights between 500 and 5500 grams. The data source is the National Center for Health Statistics (NCHS) PublicUse Perinatal Mortality Data Files.
The FLIC selects a 4component model (Figure 1a),
Component 3 is loosely analogous to the predominant distribution in the contaminated normal model [22] and the primary distribution in the 2component model [23]. Component 1 in the 4component model includes ELBW and VLBW cases, component 2 contains mostly MLBW and NBW cases but also some VLBW and HBW cases, and component 4 comprises NBW and HBW cases.
Next we fit the contaminated normal and 2component models to the same data set. For the contaminated normal model, we take the bin width to be 200 grams and use the BIC to select the number of contaminated bins [22]. Approximately 2.5% of cases are assigned to the lower residual distribution (threshold: 1700 grams), 97.5% to the predominant distribution (estimated mean and standard deviation, 3168 and 488 grams), and less than 1 in 8700 to the upper residual distribution (threshold: 5300 grams). Regarding the 2component model, approximately 88.0% of cases are assigned to the primary distribution (estimated mean and standard deviation, 3186 and 458 grams) and 12.0% to the secondary distribution (estimated mean and standard deviation, 2617 and 951 grams).
The fitted contaminated normal, 2component, and 4component models are compared in Figure 2. The contaminated normal model fits the ELBW and VLBW data nicely but exhibits artifacts at the thresholds of 1700 and 5300 grams; the contaminated normal model also understates the HBW data. The 2component model provides a good fit at most birthweights but severely understates the ELBW data. The 4component model avoids these weaknesses but has an exaggerated peak near the component 1 mean.
b. Reproducibility of order selection
In the preceding example, the selection of a 4component model was based on a specific sample of 50,000 white singletons whose mothers smoked heavily. If we draw another sample of size 50,000, will the FLIC express the same preference?
We can address this question by drawing N _{ rep } samples of size 50,000 with replacement and applying the FLIC to each sample. Here "with replacement" means that an infant can appear in more than one sample, not that an infant can appear twice in the same sample. The frequency with which the FLIC prefers a 4component model indicates the reproducibility of order selection.
Table 1 shows the verdicts of the FLIC and other criteria for N _{ rep } = 25 samples of size 50,000. The FLIC prefers a 4component model for 22 out of 25 samples; for the other three samples, the FLIC narrowly prefers a 6component model. The verdicts of the BIC match those of the FLIC. The AIC is equivocal between 6component and 7component models. Table 1 also identifies the preferences of the FLIC for sample sizes smaller than 50,000. The tendency to favor simpler models at smaller sample sizes can be understood by analogy to a hypothesis test. Imagine testing a null hypothesis that there are two components against an alternative hypothesis that there are more than two components: as the sample size decreases, the power to reject a false null hypothesis also decreases.
c. Uncertainty in parameter estimation
Although we may be comfortable using a 4component model for the birthweights of white singletons whose mothers smoked heavily, Equation (5) does not convey the uncertainty in the parameter estimates for that model.
To assess uncertainty in parameter estimation, we fit kcomponent models using each of N _{ rep } samples of equal size; in our example, k = 4 and there are N _{ rep } = 25 samples of size 50,000. Let θ represent a parameter of interest, such as μ _{3}, and let represent estimates of θ from the N _{ rep } samples. With denoting the "metasample" mean of and serving as an overall estimate of θ, and with denoting the corresponding standard deviation, we can define a confidence interval via
If were normally distributed with expected value θ, then for 95% confidence we should choose C as the upper .025 quantile of the standard normal distribution (or of a T distribution); in the absence of normality, to be conservative we could choose based on Chebychev's inequality [33]. However, not even C = 5.0 yields a coverage probability of 95% (see Section 3b of Results). There are two problems.
First, mixture model parameter estimates may have nonnegligible bias; the expected value of may not be close to θ. Second, when each of the N _{ rep } samples constitutes a large fraction of the underlying population, are not independent due to the large overlaps among the N _{ rep } samples.
The first problem can be addressed by modifying Equation (6) to
where denotes the estimated absolute value of the bias [34]. Our approach to acquiring is simulationbased. We simulate a birthweight data set from , where are the overall estimates of their respective parameters, and then compare to its own estimate arising from the simulated data set: the "drift" from to should mirror the drift from θ to . However, since relying on a single simulated data set seems precarious, we define as the average value of over five simulated data sets.
The second problem can be resolved by choosing the value of C according to the fraction of the underlying population that each of the N _{ rep } samples constitutes. Let C _{0} denote the value of C that would be chosen if this fraction were negligibly small, and let C _{ φ } denote the value that would be chosen if this fraction were equal to φ, a positive number less than 1. In Section II of [Additional file 1], we show that
Section II of [Additional file 1] also explains why we sample with replacement, why we sample instead of using the full population, and how to compare parameters within and between populations.
Table 2 lists overall estimates and confidence intervals for parameters in a 4component model for the birthweights of white singletons born to heavilysmoking mothers, using Equations (7) and (8) with the same N _{ rep } = 25 samples of size 50,000 in Table 1, C _{0} = 2.5 (see Section 3b of Results), and φ = .2465 = 50,000/202,849. Figure 1b displays the mixture model implied by the overall estimates in Table 2. Section III of [Additional file 1] examines how the overall estimates and confidence intervals change when the sample size is less than 50,000.
3. Further illustrations
a. Simulation study on model selection
For our first simulation study we generated 25 nonoverlapping data sets of size 5000 from designs A through E in Table 3; see also Figure 3. Designs A through E represent the fitted 2 through 6component models derived from the 25 samples of size 50,000 in Table 1. Values in the data sets less than 500 or greater than 5500 were discarded since the 2 through 6component models were meant to mimic a birthweight distribution; new values were drawn as needed to complete the data sets. We assessed how often the FLIC, BIC, and AIC recovered the correct number of components. This was repeated for data sets of different sizes up to 100,000.
As shown in Table 4, the FLIC and BIC consistently returned the correct answer with the 2component model at a sample size of 5000, the 3component model at a sample size of 10,000, and the 4component model at a sample size of 25,000. The FLIC and BIC did not consistently return the correct answer for the 5component or 6component model at any sample size, although they occasionally detected components 5 and 6 at a sample size of 100,000. The AIC was erratic.
At larger sample sizes, the FLIC and BIC routinely claimed a third (nonexistent) component for the 2component model. We attribute this to the removal of values less than 500 or greater than 5500, after which the 2component model was, strictly speaking, no longer a normal mixture but rather a truncated normal mixture.
b. Simulation study on calibrating confidence intervals
For our second simulation study we generated 25 overlapping data sets of size 50,000 from design C in Table 3, the degree of overlap consistent with a population of 200,000. For each of various C between 2.0 and 5.0, we used Equation (7) to form confidence intervals for the mixture parameters p _{1}, p _{2}, p _{3}, p _{4}, μ _{1}, μ _{2}, μ _{3}, μ _{4}, σ _{1}, σ _{2}, σ _{3}, σ _{4}. We recorded how many of the mixture parameters were contained in their respective confidence intervals. This was repeated nine more times, and we tabulated how many of the 120 = 12 × 10 confidence intervals contained their targets. Confidence intervals were also formed using Equation (6) for comparative purposes. The above steps were repeated with overlapping data sets consistent with a population of 1,000,000 and with nonoverlapping data sets consistent with an effectively infinite population.
The results are summarized in Table 5. With an effectively infinite population, only 81.7% of the confidence intervals formed using Equation (6) contained their targets at C = 5.0. The confidence intervals formed using Equation (7) contained their targets 95.0% of the time at C = 2.5. The adjustment suggested by Equation (8) appears reasonable: φ = .05 = 50,000/1,000,000 and N _{ rep } = 25 yield C _{ φ } = 1.315 C _{0}, which accords with the 95.8% capture of mixture parameters at C = 3.5 ≈ 1.315 × 2.5 with a population of 1,000,000.
c. Another example with real data
We also drew 25 samples of size 50,000 from the 1,749,827 black singletons who were born (or experienced fetal death) from 2000 to 2002, regardless of maternal smoking status. Table 6 records the frequencies with which the FLIC selected the 2 through 7component models as well as the overall estimates of component proportions, means, and standard deviations for each of these models. The 6component model was overwhelmingly preferred by the FLIC. Figure 4 juxtaposes the fitted 4component and 6component models implied by the overall estimates. The four components in the 4component model are loosely analogous to the second through fifth components in the 6component model, so that the main rationale for adding two more components appears to be providing a more elaborate description of the far left and right tails of the birthweight distribution.
Discussion
Our approach to modeling birthweight distribution is distinguished from previous proposals in that the data determine the number of components in the normal mixture model. We have seen that data sets of size 50,000 for white singletons born to heavilysmoking mothers typically warrant 4 components, while data sets of size 50,000 for black singletons usually demand 6 components. These results underscore the idea that a one size fits all paradigm  whether that be a 2component normal mixture model or even the across the board use of a 4component normal mixture model  may lead to unreasonable representations of birthweight distribution for some populations. Our approach, on the other hand, allows birthweight distribution to be described differently for different populations. We also note here that, although results have not been presented in this paper for a full spectrum of populations, our experience has been that data sets of size 50,000 usually call for between 3 and 6 components.
The second paper in our twopart series will elucidate the main advantage of our approach over the contaminated normal model [21, 22] and the 2component model [23], namely its greater potential to expose heterogeneity in mortality risk. By this we mean that, even at a fixed birthweight, some infants may be at higher risk than others. While such heterogeneity seems plausible, if not altogether obvious, it may not be adequately expressed by either the contaminated normal model or the 2component model. Hence, allowing a model to have more than 2 components is not an intellectual exercise or fitting the data for the sake of fitting the data but rather a way to improve assessment of mortality.
Since gestational age is sometimes considered in tandem with birthweight [19, 20], we now comment on its relation to the methodology in this twopart series.
Our approach to modeling birthweight distribution does not explicitly consider gestational age. However, our experience is that the first component typically captures most very preterm births. For instance, the birthweight distribution for white singletons with gestational ages > 37 weeks is well approximated by a 3component model whose components resemble the second through fourth components of a 4component model for white singletons in general.
Even so, one may be interested in extending our methodology to explicitly consider gestational age and/or other covariates. We envisage at least two possible extensions. The first would generalize the work of Fang, Stratton, and Gage [19] in which the number of components had been constrained a priori to two, while the second would be novel.
The first extension would be to model the joint probability density of birthweight and gestational age as a bivariate normal mixture, with the number of components determined from the data using the FLIC rather than being constrained a priori to two. Then, instead of estimating the mortality risk within each component as a function of birthweight only, one could estimate the mortality risk within each component as a function of both birthweight and gestational age.
The second extension would be to retain the univariate normal mixture model for birthweight distribution but create auxiliary models to relate covariates, such as gestational age, to mixture components. The appeal of this extension is that it could allow some mixture components to be placed in approximate correspondence with identifiable subpopulations.
Conclusions
The present paper, the first in a twopart series, develops a new and flexible approach to modeling a birthweight distribution using a normal mixture model with the number of components determined from the data rather than fixed a priori. This approach allows the detection of heterogeneity in birthweight that cannot be found with a contaminated normal model or a 2component normal mixture model. Unlike a contaminated normal model, our approach does not assume the existence of an interval of birthweights over which there are no compromised pregnancies. Unlike a 2component normal mixture model, our approach does not constrain birthweights within compromised pregnancies to be normally distributed. Yet, better modeling of birthweight distribution is a means to an end, namely a greater understanding of fetalinfant mortality. The second paper in our twopart series reveals that, when coupled with methodology for estimating birthweightspecific mortality curves within each component, this paper's approach to describing a birthweight distribution can also reveal heterogeneity in mortality.
Methods
[Additional file 1] presents technical details on our methodology and its implementation.
Abbreviations
 AIC:

Akaike Information Criterion
 BIC:

Bayesian Information Criterion
 ELBW:

extremely low birthweight
 EM:

expectation maximization
 FLIC:

Flexible Information Criterion
 HBW:

high birthweight
 MLBW:

moderately low birthweight
 NBW:

normal birthweight
 NCHS:

National Center for Health Statistics
 PMLR:

parametric mixtures of logistic regressions
 VLBW:

very low birthweight
References
 1.
Brimblecombe F, Ashford J, Fryer J: Significance of Low Birth Weight in Perinatal Mortality: A Study of Variations within England and Wales. Br J Prev Soc Med. 1968, 22: 2735.
 2.
Rooth G: Low birthweight revised. Lancet. 1980, 1: 639641. 10.1016/S01406736(80)911307.
 3.
Goldstein H: Factors related to Birth Weight and Perinatal Mortality. Br Med Bull. 1981, 37: 259264.
 4.
Fryer J, Hunt R, Simons A: Biostatistical Considerations: The Case for Using Models. Child Health. 1984, 3: 930.
 5.
Kleinman JC: Methodological Issues in the Analysis of Vital Statistics. Reproductive and Perinatal Epidemiology. Edited by: Kiely M. 1991, Boca Raton: CRC Press, 453462.
 6.
Kiely JL, Kleinman JC: BirthWeightAdjusted Infant Mortality in Evaluations of Perinatal Care: Towards a Useful Summary Measure. Stat Med. 1993, 12: 377392. 10.1002/sim.4780120319.
 7.
Cogswell M, Yip R: The Influence of Fetal and Maternal Factors on the Distribution of Birthweight. Semin Perinatol. 1995, 19: 222240. 10.1016/S01460005(05)80028X.
 8.
Klebanoff MA, Schoendorf KC: What's So Bad about Curves Crossing Anyway? (Invited Commentary). Am J Epidemiol. 2004, 160: 211212. 10.1093/aje/kwh203.
 9.
Basso O, Wilcox A, Weinberg C: Birthweight and Mortality: Causality or Confounding?. Am J Epidemiol. 2006, 164: 303311. 10.1093/aje/kwj237.
 10.
Basso O: Birthweight is Forever. Epidemiology. 2008, 19: 204205. 10.1097/EDE.0b013e31816379d9.
 11.
Bjørstad AR, IrgensHansen K, Daltveit AK, Irgens LM: Macrosomia: Mode of Delivery and Pregnancy Outcome. Acta Obstet Gynecol Scand. 2010, 89: 664669. 10.3109/00016341003686099.
 12.
MacDonald H: Perinatal Care at the Threshold of Viability. Pediatrics. 2002, 110: 10241027. 10.1542/peds.110.5.1024.
 13.
Blackmon L, Batton DG, Bell EF, Denson SE, Engle WA, Kanto WP, Martin GI, Stark AR: Levels of Neonatal Care. Pediatrics. 2004, 114: 13411347. 10.1542/peds.114.1.229.
 14.
Russell RB, Green NS, Steiner CA, Meikle S, Howse JL, Poschman K, Dias T, Potetz L, Davidoff MJ, Damus K, Petrini JR: Cost of hospitalization for preterm and low birth weight infants in the United States. Pediatrics. 2007, 120: 19. 10.1542/peds.20062386.
 15.
Wilcox AJ, Russell IT: Birthweight and Perinatal Mortality: II. On WeightSpecific Mortality. Int J Epidemiol. 1983, 12: 319325. 10.1093/ije/12.3.319.
 16.
Escobar GJ, McCormick MC, Zupancic JAF, ColemanPhox K, Armstrong MA, Greene JD, Eichenwald EC, Richardson DK: Unstudied Infants: Outcomes of Moderately Premature Infants in the Neonatal Intensive Care Unit. Arch Dis Child Fetal Neonatal Ed. 2006, 91: 238244. 10.1136/adc.2005.087031.
 17.
Wilcox AJ: On the importanceand the unimportanceof birthweight. Int J Epidemiol. 2001, 30: 12331241. 10.1093/ije/30.6.1233.
 18.
Gage T, Bauer M, Heffner N, Stratton H: Pediatric Paradox: Heterogeneity in the Birth Cohort. Hum Biol. 2004, 76: 327342. 10.1353/hub.2004.0045.
 19.
Fang F, Stratton H, Gage T: Multiple mortality optima due to heterogeneity in the birth cohort: a continuous model of birthweight by gestational agespecific infant mortality. American Journal of Human Biology. 2007, 19: 475486. 10.1002/ajhb.20607.
 20.
Schwartz SL, Gelfand AE, Miranda ML: Joint Bayesian Analysis of Birthweight and Censored Gestational Age Using Finite Mixture Models. Stat Med. 2010, 29: 17101723. 10.1002/sim.3900.
 21.
Wilcox AJ, Russell IT: Birthweight and Perinatal Mortality: 1. On the Frequency Distribution of Birthweight. Int J Epidemiol. 1983, 12: 314319. 10.1093/ije/12.3.314.
 22.
Umbach D, Wilcox AJ: A Technique for Measuring Epidemiologically Useful Features of Birthweight Distributions. Stat Med. 1996, 15: 13331348. 10.1002/(SICI)10970258(19960715)15:13<1333::AIDSIM271>3.0.CO;2R.
 23.
Gage T, Therriault G: Variability of BirthWeight Distributions by Sex and Ethnicity: Analysis Using Mixture Models. Hum Biol. 1998, 70: 517534.
 24.
Gage T: BirthWeightSpecific Infant and Neonatal Mortality: Effects of Heterogeneity in the Birth Cohort. Hum Biol. 2002, 74: 165184. 10.1353/hub.2002.0020.
 25.
Akaike H: Information theory and an extension of the maximum likelihood principle. Second International Symposium on Information Theory. Edited by: Petrov BN, Csaki F. 1973, Akademiai Kiado, Budapest
 26.
Schwarz G: Estimating the dimension of a model. Annals of Statistics. 1978, 6: 461464. 10.1214/aos/1176344136.
 27.
Titterington D, Smith AFM, Makov U: Statistical Analysis of Finite Mixture Distributions. 1985, Wiley, New York
 28.
Lindsay BG: Mixture Models: Theory, Geometry and Applications. 1995, IMS NSFCBMS Regional Conference Series, Hayward
 29.
McLachlan G, Peel D: Finite Mixture Models. 2000, Wiley, New York
 30.
Charnigo R, Sun J: Testing homogeneity in a mixture distribution via the L^{2} distance between competing models. Journal of the American Statistical Association. 2004, 99: 488498. 10.1198/016214504000000494.
 31.
Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithm. J R Stat Soc. 1977, 39: 122.
 32.
Martin J, Hoyert D: The National Fetal Death File. Semin Perinatol. 2002, 26: 311. 10.1053/sper:2002.29834.
 33.
Casella G, Berger R: Statistical Inference. 2002, Duxbury, Pacific Grove, 2
 34.
Loader C: Local Regression and Likelihood. 1999, Springer, New York
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712393/10/37/prepub
Acknowledgements
The authors thank Vicki Flenady, Gerald Hoff, and an anonymous Associate Editor for feedback that led to improvement of this manuscript.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
RC  Concept and design, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, statistical analysis, read and approved final manuscript. LWC  Concept and design, acquisition of data, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, read and approved final manuscript. TL  Analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, read and approved final manuscript. RSK  Analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, read and approved final manuscript.
Electronic supplementary material
Technical Appendix
Additional file 1: . Additional file 1 presents technical details on our methodology and its implementation. (DOC 122 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Charnigo, R., Chesnut, L.W., LoBianco, T. et al. Thinking outside the curve, part I: modeling birthweight distribution. BMC Pregnancy Childbirth 10, 37 (2010). https://doi.org/10.1186/147123931037
Received:
Accepted:
Published:
Keywords
 Bayesian Information Criterion
 Normal Model
 Normal Mixture
 Residual Distribution
 Order Selection