 Research
 Open Access
 Open Peer Review
 Published:
Measuring and tracking obesity inequality in the United States: evidence from NHANES, 19712014
Population Health Metrics volume 14, Article number: 12 (2016)
Abstract
Background
Because people care about their weight relative to peers and society, obesity inequality plays a role in explaining obesity incidence and the impacts of being obese on subjective wellbeing. While the increase in obesity prevalence and mean body mass index (BMI) is well documented, the measurement of distributional changes and corresponding obesity inequality is yet to be fully explored.
Methods
The present study analyzed BMI data for adults aged 20 to 74 from the National Health and Nutritional Examination Survey (NHANES) I (19711974), II (19761980), III (19881994), and continuous NHANES (19992014). We applied tools developed to measure income inequality to analyze the intertemporal variation in the BMI distribution among US adults. Using stochastic dominance tests, we construct partial orderings on cumulative BMI distributions during the study period. Shapley decompositions and inequality indices are employed to quantify the source and extent of temporal variation and decompose the inequality into within and betweengroup components considering age, gender, and race.
Results
The BMI distribution of each NHANES study firstorder stochastically dominated the BMI distribution of the previous wave from 19711974 to 20032006, whereas more recent comparisons failed to reject the null hypothesis of nondominance. The Shapley decomposition analysis revealed that horizontal shifts of BMI distributions accounted for a majority of the increase in obesity prevalence since 19881991. Especially in recent years when the rate of obesity growth has slowed down, the contribution of the redistribution component dropped significantly and even became negative between 20072010 and 20112014. The inequality indexes consistently show a worsening of obesity inequality from the mid1970s to the mid2000s regardless of population subgroups, and this disproportionate shift of the BMI distribution is unlikely to be a result of a changing ethnic composition of the US population.
Conclusion
Our findings demonstrate that seemingly similar increases in obesity prevalence can be accompanied by very different patterns of distribution change. We find that the early phase of the obesity epidemic in the US was largely driven by increasing skewness, whereas more recent growth is a populationwide experience, regardless of demographic characteristics. Increasing morbid obesity certainly played an important role in the initial phase of the epidemic, but more recently the BMI distribution has largely horizontally shifted to the right.
Background
Many studies have documented a marked increase in obesity prevalence and mean body mass index (BMI) in the US over the last four decades [1–6]. This significant and consistent rise in bodyweight has been termed an "obesity epidemic," spreading across all gender, age, and ethnic groups. The use of such language evokes the idea of obesity being contagious, spreading from one person to another. For instance, gains in weight appear to spread through social ties, with friends and relatives apparently influencing others in their social network, in a way reminiscent of a contagious disease [7].
In conjunction with an overall rise in obesity prevalence, there has been an even more significant increase in the percentage of US adolescents [8, 9] and adults [4, 10, 11] who are morbidly obese. This is reflected in a rightward shift of the BMI distribution, more pronounced at its upper tail [12]. Such changes in the shape of the BMI distribution, and in particular the disproportionate growth in the distribution’s upper tail, have been explained by models that interact the effects of economic change (e.g., falling food prices) with social and physiological processes. In the social process, a person’s body weight standard depends on other people’s weight, and a relaxed standard can lead to weight increases [13–15]. In a society where one’s weight does not conform to the socially ideal weight, social pressure may exist, and result in a disutility cost to individuals [16].
In a study of 29 European countries, evidence suggests that overweight perceptions and dieting are influenced by a person’s relative BMI [17]. The authors suggest that, for a variety of reasons, it may be easier to be fat in a society that is fat, and provide empirical evidence that relative BMI influences subjective wellbeing. A more recent study demonstrates that the degree to which obesity is negatively associated with life satisfaction can be mitigated by the prevalence of obesity in a given geographic context [18].
Despite the importance of relative obesity in explaining the causes and consequences of the obesity epidemic, little is known about the evolution of BMI distributions over time. A notable exception is Contoyannis and Wildman in which relative distribution methods are employed to track changes in BMI using nonparametric methods [19]. Focusing on Canada and England, they found that the increase of obesity in England is characterized by more polarized growth towards the rightend of the BMI distribution, whereas the increase of obesity in Canada is driven primarily by an overall upward shift. While growing obesity inequality is believed to be a populationwide experience [12], a recent study of Americans (19992006) revealed a different pattern of polarization, with a more pronounced shift among ethnic minorities and the less educated [20].
Building upon this work, we expand on the focus, methods, and time horizon of earlier studies to present a longterm picture (19712014) of obesity inequality in the US incorporating a wide array of quantitative methods used in the study of income inequality [21, 22]. In addition to presenting an encompassing measurement of the transition of BMI inequality during the US obesity epidemic, the economic tools employed to analyze changes in obesity offer several new insights.
First, borrowing tools common in the study of poverty, we employ Stochastic Dominance (SD) tests. SD is very useful when making nonparametric comparisons between distributions of continuous variables such as income or, in our case, BMI. SD tests offer an ordinal comparison between distributions; a comparison that ranks the distributions but that does not estimate the magnitude of the differences between them. Because SD tests involve comparison over the BMI domain, they are independent of the choice of an obesity threshold. Moreover, as we explain below, the computation of statistics of dominance at multiple test points covering the whole range of the BMI distributions enables us to assess whether the shift of the BMI distribution was driven primarily by one part of the distribution or by the entire population [21, 22]. Second, we apply the Shapleyvalue based decomposition technique to decompose total change in obesity prevalence into a meangrowth effect and a redistribution component [21, 23, 24]. While the SD test primarily focuses on determining the ordinal dominance of BMI distributions, this second approach shows how much of the growth in obesity inequality is attributable to a horizontal shift of the distribution (an increase in average BMI), and how much is due to a change in the shape of the distribution (in particular, an increased skewness towards the right tail of the distribution). Third, we provide a singlevalue quantitative assessment of the degree of inequality measured by conventional inequality indices (Gini and Generalized Entropy). When analyzing univariate inequality measures, we pay particular attention to decomposing obesity inequality into inequalities within each segment of the population and inequalities between subgroups. This complementary approach allows us to reveal detailed aspects of the transition across subpopulations and examine whether a disproportionate growth of obesity is due to populationwide shifts of the BMI distribution or to a changing contribution by the different demographic groups.
Several studies point out that focusing on prevalence estimates is at best a crude approach to understanding the obesity epidemic since it ignores much of the available information, and the measurement of obesity rates depends heavily on somewhat arbitrary thresholds and does not correctly reflect the clinical implications of obesity around cutoff values [19, 21, 22, 25, 26]. Since the BMI distribution in the US shows that most individuals are centered around the overweight category (i.e., a BMI greater than or equal to 25), even a minor rightward shift of the BMI distribution would result in significantly higher prevalence estimates, which might overestimate the seriousness of obesity. If the recent growth of BMI is more pronounced for those at the right tail of the distribution, tracking only prevalence estimates over time does not correctly reflect accompanying mortality and morbidity risks. Complementing obesity prevalence estimates with distributionindependent techniques such as those proposed by this paper is thus critical to understanding the longterm pattern of obesity change in the US.
Methods
Data sources and study population
Baseline data were drawn from the National Health and Nutritional Examination Survey (NHANES) I (19711974), II (19761980), III (19881994), and continuous NHANES (19992014). The NHANES is a series of nationally representative crosssectional surveys of the US population, conducted by the National Center for Health Statistics at the Centers for Disease Control and Prevention. The NHANES data include clinical measurements of the respondents’ height and weight obtained using mobile examination centers and standardized procedures. This is an attractive feature of the data because selfreports of height and weight tend to be biased, leading to an underestimate of BMI [27–29]. This is particularly important when observations are clustered around the middle of the distribution since individuals whose BMI is close to the obesity threshold are more likely to underreport their weight [21].
For this study, we restricted the analysis to 2074 year olds because only persons aged 174 years were eligible to be interviewed for NHANES I and II, and different weight classification criteria are used for people under 20. Following NHANES analytic guidelines [30], respondents were classified into three age groups, 2039 years, 4059 years, and 6074 years, based on age at the examination date. For NHANES III and continuous NHANES, race and ethnicity were classified as nonHispanic white, nonHispanic black, Hispanic, and other. To facilitate comparability across waves, given the larger number of years covered by the initial waves, we aggregated two adjacent waves of continuous NHANES into a single survey. Since NHANES III covered a six year period (1988 to 1994) when the largest increase in obesity was experienced, we split NHANES III into two different periods, Phase I (19881991) and Phase II (19911994), according to NHANES III analytic guidelines [31]. The multistage sampling design of NHANES III selects 81 primary sampling units for a full six year survey, and then randomly assigns each primary sampling units to Phase I and Phase II, which makes each subsample representative of the US civilian noninstitutionalized population during the given period. In addition, we excluded female respondents who were pregnant at the time of the survey. The final sample for empirical analysis did not include observations with missing values or irregular responses in height and weight.
We adopted the clinical definition of obesity proposed by the International Obesity Task Force of the World Health Organization in 1997. Underweight is defined as BMI < 18.5, normal weight as BMI in the interval [18.5, 25.0), overweight as BMI ∈[25.0, 30.0), class I obesity as BMI ∈[30.0, 35.0), class II obesity as BMI ∈[35.0, 40.0), and class III obesity as BMI ≥40.0. The prevalence or increase in BMI of morbid or clinically severe obesity across waves was of special interest to this research since many direct medical costs associated with obesity are most pronounced for those at the higher spectrums of the obesity scale.
To account for the complex, stratified, multistage probability cluster sampling design of NHANES, we applied mobile examination centers sampling weights throughout the analysis. For NHANES III, a 3year sampling weight was applied to each Phase I and II, which benchmarked the 1990 and 1993 Current Population Survey (CPS), respectively [31, 32]. For the first combined sample of continuous NHANES, we used fouryear sample weights for 19992002, which was preadjusted by NHANES to account for the difference in the population base of the 19992000 and 20012002 surveys [30]. By rescaling a twoyear weight of adjacent surveys, this sample weight allowed us to make our sample representative of the population at the midpoint of the two surveys, even if different population bases were considered. For the subsequent fouryear datasets, we created a fouryear sample weight variable that assigned half of the twoyear weight for each period, as recommended by NHANES analytic guidelines [30]. We could then compare the distributions over time, since these weighting schemes were designed to ensure that the weighted sample was representative of the US civilian noninstitutionalized population; that is, it reflected the relative proportion of each demographic group to ensure equal selection probability of an individual given that some groups were oversampled.
Stochastic dominance test
The Stochastic Dominance test is an approach that allows an ordinal assessment of whether a cumulative distribution significantly differs from another without considering the shape of the distribution [33]. We applied the test here to determine the dominance of BMI distributions of US adults over time. Although commonly used by economists in poverty and economic inequality studies [33–35], it has only recently been applied to the study of obesity [21, 22].
Let \( {F}_{t_{n1}}(x) \) and \( {F}_{t_n}(x) \) denote two cumulative distribution functions (CDF) of BMI to be compared to each other, where t_{ n1 } and t_{ n } refer to time, i.e. to different NHANES waves, and\( {D}_t^1(x)={F}_t(x)={\displaystyle \underset{0}{\overset{x}{\int }}d{F}_t(y)} \)and\( {D}_t^s(x)={\displaystyle \underset{0}{\overset{x}{\int }}{D}_t^{s1}(y) dy} \) for any integer s ≥ 2. The distribution at time t_{ n }dominates the distribution t_{n − 1 }at order s if \( {D}_{t_{n1}}^s(x)\ge {D}_{t_n}^s(x) \) and strictly dominates if \( {D}_{t_{n1}}^s(x)>{D}_{t_n}^s(x) \), for all possible BMI values over the domain [21, 33].
Simple tstatistics were constructed to test the null hypothesis of nondominance (\( {H}_0:\kern1em {D}_{t_{n1}}^s(x){D}_{t_n}^s(x)=0 \)), for a series of test points up to the maximum BMI in the distribution. Unlike other studies testing dominance only within a range of interest, we tested the significance over the entire domain (we used 30 test points from the minimum to the maximum BMI) in order to investigate which part of the distribution changed most. Dominance of order s was declared if the null hypothesis was rejected for at least one test point at the 1 % significance level without any reversal in the signs of difference [21]. The stochastic dominance does not hold if, for instance, the difference is not significant or two cumulative distributions cross each other. In general, it has been shown that the stochastic dominance of one distribution over another can always be declared at a high enough order, provided that infinite comparisons between CDFs can be made [37]. The interpretation of higherorder comparisons is, however, less intuitive [38] and, in practice, comparisons are limited to thirdorder stochastic dominance [33]. We followed the convention of testing up to s = 3, i.e. thirdorder stochastic dominance, after which "no dominance" is declared [22, 35, 36].
Growthinequality decomposition
The growthinequality decomposition method allows researchers to decompose overall changes in a distribution into a meangrowth component and a redistribution component [39]. In the context of obesity, the meangrowth component captures the change in obesity prevalence attributable to a horizontal shift of the BMI distribution while holding the shape of the distribution constant at the reference year. The redistribution component represents the change in obesity as a result of a redistribution in the BMI curve while the mean BMI is kept constant. A third component, by definition, is the residual that cannot be exclusively attributed to the previous two elements. When applied to our study, the obesity rate at time t, Obs_{ t }, can be represented as
where Obs denotes obesity prevalence, T is the obesity threshold (30 for class I obesity), μ is the mean BMI, and L is the Lorenz curve representing the CDF of the empirical probability distribution of BMI. Letting t_{n − 1} be the base year, changes in obesity prevalence between two timeperiods can then be decomposed as
where G(⋅), R(⋅), and ε(⋅) represent the growth, redistribution, and residual components, respectively. Specifically, the growth and redistribution terms were defined as
\( G\equiv Obs\left(T\Big{\mu}_{t_n};{L}_{t_{n1}}\right)Obs\left(T\Big{\mu}_{t_{n1}};{L}_{t_{n1}}\right) \), and\( R\equiv Obs\left(T\Big{\mu}_{t_{n1}};{L}_{t_n}\right)Obs\left(T\Big{\mu}_{t_{n1}};{L}_{t_{n1}}\right) \). That is, G was the change in obesity driven by overall growth in population weight while holding relative position fixed, and R was the observed variation in relative position with no growth in mean BMI.
In empirical studies, a residual term controls for misspecified components in the decomposition analysis, which confound the interpretation of decomposition results, particularly when the residual term is relatively large [39]. A more desirable method would decompose changes in prevalence measures exactly into growth and redistribution factors without a residual term. In this context, a Shapleyvalue based decomposition approach which takes an equally weighted average of two decompositions, one at a reference point and the other at a later year, has been proposed [23, 24]:
where G^{s} and R^{s} represent the Shapley value of growth and distribution components of changes in obesity prevalence, and are given by:
Obesity inequality indices (Gini and Generalized Entropy)
While distributional dominance tests offer a partial ranking of BMI distributions, they do not measure cardinal differences between distributions and there are cases in which stochastic dominance cannot be determined. Decomposition analysis is also limited in that it relies heavily on the specific obesity threshold, T, to decompose the variation. In this section, we supplement our previous findings by summarizing obesity inequality into a univariate concentration index, the Gini coefficient, and track the cardinal growth of obesity inequality. Typically used to quantify income and wealth inequality, the Gini coefficient measures the statistical dispersion in a given distribution. It varies between 0, which reflects complete equality, and 1, which indicates complete inequality (one person has all the income or wealth, all others have none). Our approach closely followed Sahn where the Gini index was employed to track the obesity inequality in developing countries [22]. The Gini coefficient is computed as follows:
where N is the sample size, μ denotes the mean BMI, x_{ i } and r_{ i } represent individual BMI and corresponding rank of the i^{th} observation in ascending order. Considering the sample size in our study, we referred to a computationefficient formula, which approximates the Gini coefficient using a fast optimized algorithm [40].
In addition to the common Gini coefficient as a summary measure of inequality, generalized entropy (GE) inequality measures have been proposed. Compared to the Gini coefficient, which is more sensitive to variations around the mode of the distribution, the GE measures are more flexible allowing greater sensitivity away from the middle of the distribution [41]. This is an attractive feature given the focus of recent studies on the rise in morbid obesity at the upper tail of the BMI distribution and its contribution to the overall rise in obesity. The Generalized Entropy index can be expressed as:
where θ is a scaling parameter that represents the weight given to distances between individuals' BMI at different parts of the BMI distribution. For θ = 1 we obtained the Theil index, which treats differences between individuals' BMI levels at different points of the BMI distribution equally. The variation at the left tail of the distribution is given more weight with parameter values smaller than 1, whereas larger parameter values give more weight to the upper tail. We set θ equal to 0 and 2 for robustness and for comparison with a previous study examining this issue [41].
With distributional dominance tests and decomposition analysis, we were unable to split the analytic sample by demographic category and explore whether the growing obesity inequality is due to changing characteristics of particular population segments, or due to a populationwide shift of the BMI distribution. In the analysis of inequality over time, it could be the case that a growing inequality is influenced by greater disparities among different segments of the population, or by variation in the distribution of BMI within each subpopulation (provided the relative weight of different groups in the total population does not change). The GE class of inequality measures can be decomposed into within and betweengroup inequality such that
GE_{ t }(θ)=GE_{ t }(θ)_{ within } + GE_{ t }(θ)_{ between } [41, 42]. Specifically, GE_{ t }(θ)_{ within } = \( {\displaystyle \sum_j\left(\frac{BM{I}_{t,j}}{BM{I}_t}\right)G{E}_{t,j}} \) and GE_{ t }(θ)_{ between } = \( {\displaystyle \sum_j\left(\frac{BM{I}_{t,j}}{BM{I}_t}\right) \ln \left(\frac{BM{I}_{t,j}/BM{I}_t}{N_{t,j}/{N}_t}\right)} \)
where GE_{ t,j } and BMI_{ t,j } denote the GE index and BMI of subgroup j at time t, respectively; N_{ t,j } represents the number of respondents in subgroup j at time t; and BMI_{ t } represents the BMI of the total population at time t. That is, the first term indicates the weighted sum of inequalities within groups, whereas the second term captures the proportion attributable to the heterogeneity in inequality across the groups. If the contribution of betweengroup inequalities to total obesity inequality is negligible, and the evolution of withingroup inequality is comparable across groups, this indicates that worsening obesity inequality is not a result of changing demographic composition (even if the weight of different groups in the total population has indeed changed), but rather more of a populationwide experience.
Results and discussion
Studies reporting on the obesity epidemic have documented a dramatic increase in obesity prevalence over the last four decades [1–6]. In particular, the percentage increase in the higher obesity classes is substantial, suggesting that the population weight distribution has been disproportionately shifted rightward [11, 12, 43]. As illustrated in Fig. 1, which presents the kernel density and cumulative distribution of BMI over time, the increase in BMI was more pronounced between the late 1970s to early 2000s, whereas significant wavetowave differences have been small or not found in recent periods.
Stochastic dominance tests
To better assess the longterm transition of US population weight, we present the results of SD tests on BMI distributions over time (Table 1). Unlike ttests on prevalence estimates, stochastic dominance provides nonparametric pairwise comparisons of entire distributions (in our case at 30 points over the BMI domain) so that the comparisons of whether one CDF is greater in magnitude to the other can be made. Table 1 shows a clear pattern of BMI distribution dominance for all yeartoyear comparisons up until the mid2000s. Until the 20032006 survey, each BMI distribution firstorder dominated the previous distribution, which indicates a significant difference for at least one test point without a significant crossing of distributions. That is, the temporal increase in cumulative distributions during this period was greater than or equal to zero over the domain of the BMI distribution.
Although we observed significant temporal shifts in a statistical sense, the nature of the transition depends upon which part of the distribution most contributes to the dominance of one distribution over another. Table 2 shows tstatistics of firstorder dominance tests at 30 test points covering the whole range of the BMI distributions. For instance, the first column compares the CDF of the 19761980 survey to the previous period, and dominance is declared only at the second test point. Although we rejected the null hypothesis of nondominance at the firstorder of comparison, this difference did not seem meaningful in an economically significant sense. From 19761980 through 19992002, a significant increase was observed across most of the domain, indicating populationwide upward shifts of the BMI distribution. More specifically, from 19761980 to 19881991 and 19992002 to 20032006, a significant increase was found at the very upper tail of the distribution (at the 29^{th} and 30^{th} test points), indicating an even more disproportionate shift than in other periods. The shift of the BMI distribution across 19881991 and 19911994 and 19992002 and 20032006 was relatively more pronounced around the middle of distribution, although the distribution was becoming more skewed in the latter period. The null hypothesis of nondominance was not rejected during the 20032006 through 20072010 while the most recent comparison between 20072010 and 20112014 found a disproportionate downward shift across the very top end of the distribution.
Overall, results from the SD tests indicated that the distribution of BMI has disproportionately shifted upwards between 1971 and 2003, but this shift stalled in the mid2000s.
Growthinequality decomposition
In addition to dominance test results, Growth Incidence Curves graphically describe which part of the BMI distribution contributed more to the overall growth between two sampling periods (Fig. 2). They show the percentage change at each BMI percentile with reference to the horizontal line representing the rate of prevalence growth [44]. Consistent with the SD test results, Fig. 2 shows a moderate BMI increase from NHANES I (19711974) to II (19761980) caused by the shift of the lower tail of the distribution. Since then, we observe a clear pattern of a rapid rise in obesity prevalence between 19761980 and 19911994 that was most pronounced in higher obesity percentiles. Interestingly, over 19911994 and 19992002 when the highest increase was experienced, the rate of obesity growth was approximately the same across all BMI levels.
Table 3 presents the results of the Shapley decomposition, which decomposed the total change in the obesity rate between sampling periods into growth and redistribution components. Not surprisingly, our results supported the SD test results in Table 2. For instance, between 19761980 and 19881991 when we observed an upward shift at the right tail of distribution, approximately 26.6 % of the increase in obesity prevalence was explained by the redistribution component. Through 19881991 and 19992002, the BMI distribution of subsequent periods firstorder dominated the previous distribution mostly around the middle of the domain, and this was reflected in a smaller redistribution effect in the corresponding time frame. Similarly, the disproportionate shift from 19992002 to 20032006 being more pronounced at the righttail was reaffirmed by a sizable redistribution effect in Table 3. Throughout the decomposition analysis, two clear results stand out. First, except for the initial sample period comparison of 19711975 to 19761980, the increase in obesity in the US has been predominately due to the meangrowth effect. This implies that the recent rise in obesity in the US has not been a statistical artifact applicable only to a particular population group or due to the arbitrariness of the obesity threshold. However, while playing a lesser role, the redistribution component was positive and nontrivial up until the last sampling period. This indicates that between the 1970s and 2000s there was a continual increase in obesity inequality. Second, an interesting change in dynamics was observed between 2003 and 2012. During this time span, when the change in the obesity rate slowed substantially compared to the previous two decades, the contribution of the redistribution component shrank substantially and became negative between 20072010 and 20112014. If this trend continues, this could indicate that the rise in obesity inequality observed in the previous three decades could be stalling or even reversing despite continued increases in obesity prevalence.
Obesity inequality indices
Table 4 shows the historical trends for two different indices of BMI inequality: the Gini coefficient and Generalized Entropy when θ = 0 or 2. According to the Gini coefficient, there has been a steady and significant increase in obesity inequality in the US since the 1970s. Specifically, the degree of inequality measured by the Gini index increased most rapidly from the 19761980 to the 19881991 waves, followed by a relatively moderate but significant rise until 19992002. This pattern of transition is consistent with the SD test where firstorder dominance was observed at the upper tail of distribution between 19761980 and 19881991. Similarly, the increase in obesity inequality from 19992002 to 20032006 was slightly greater than that of more recent periods as the top end of distribution significantly shifted.
Mirroring the results of the Gini coefficient, the Generalized Entropy index also indicated that obesity inequality increased significantly, but suggested that the rate of growth was even greater (about twice that of the Gini index). The growing obesity inequality measured by GE(2) and GE(0) corresponded closely, indicating the robustness of our findings regardless of the relative importance of the lower or upper tails of the distribution. Overall, our analysis of obesity inequality suggests that the US adult population has experienced growing obesity inequality.
Breaking down the Gini coefficient based upon age, gender, and race categories indicates that the growth in obesity inequality has been a populationwide phenomenon across subpopulations (Fig. 3). For instance, both males and females have experienced a substantial increase in obesity inequality, with the rate of growth being nearly identical, although females started from a higher level of inequality. Consistent with Flegal and Troiano, we also found evidence of a more disproportionate shift of BMI distribution among younger adults [12]. We found no evidence that this disproportionate growth is due to particular ethnic groups, although the increase was less pronounced among the Hispanics.
Table 5 reports the obesity inequality decomposed into the within and betweengroup component. Given the little difference between GE(0) and GE(2), we present a decomposition of GE(1) (i.e., the Theil index) by age, gender, race and combinations of these categories. As expected, most obesity inequality was due to the inequality within groups, and this pattern did not vary significantly over time. The betweengroup inequality was very small. For instance, approximately 0.66 % (=0.00016/0.02429)  2.09 (=0.00056/0.02685) of total obesity inequality was attributable to race/ethnicity. Gender accounted for approximately 0.010.33 % of total obesity inequality, whereas age explained about 0.933.49 % of total obesity inequality. The combination of age, gender, and race accounted for only 3.384.76 % of inequality. More importantly, we did not find a systematic increase or decrease in the obesity inequality attributable to particular betweengroup components, although disparities among the ethnic groups appear to have increased slightly since 19992002. That is, this unequal shift of the BMI distribution has been more of a populationwide experience across the US.
Conclusions
Contributing to the growing literature focusing on different forms of inequality (e.g., income, wealth), this study quantifies through three alternative, complementary methods commonly employed in the economics literature, the trends in obesity inequality that have been experienced in the US during the past four decades. The methods range from ordinal comparisons of BMI distributions (through nonparametric tests of stochastic dominance between BMI distributions), to cardinal comparisons of inequality summarized in Gini and Generalized Entropy indices, including a decomposition of changes in the BMI distribution into mean growth, and a redistribution component driven by changes in the shape of the distribution. We used BMI data from the representative sample of adults in the NHANES (19712014), and found evidence consistent across the three methods that the rapid growth in obesity prevalence in the US has been accompanied by growing obesity inequality. Further, we found that a disproportionate shift of BMI distribution occurred when US obesity was increasing the most. This growth in inequality, which is not simply contained to an expansion in the righttail of the BMI distribution of extremely morbidly obese adults, is found across the distribution of BMI and population subgroups.
While the increase in obesity rates in the US is a clear economic and medical concern due to the direct linkages between obesity and chronic illnesses and increased medical costs, the increase in obesity inequality is a problem as well for two key reasons: First, there is evidence that life satisfaction is influenced by the prevalence of obesity surrounding an individual [18]. This suggests that as obesity inequality expands, thus driving a greater wedge between individuals lower versus higher on the BMI distribution of society, the negative influence of obesity on life satisfaction could increase. Thus, the growth in obesity inequality raises the specter that the wide array of negative social and psychological consequences that have been linked to individual obesity levels could be exacerbated in an increasingly obesityunequal nation. Identified consequences of obesity such as discrimination [45–48], depression [49–53], and increased social stigma [54], could potentially be intensified with growing obesity inequality. Second, because being obese may be easier in a fatter society as individuals judge their BMI relative to their peers, a more unequal BMI distribution may lead to further growth in obesity prevalence. A relaxation in relative weight standards may create a vicious circle, which imposes additional social costs on heavier societies. Overall, it is hoped that the presented measurements of obesity inequality will help spur further research to better understand not only the consequences of increased obesity prevalence but also how increased obesity inequality is affecting individuals across the BMI spectrum and society as a whole.
Abbreviations
 BMI:

body mass index
 NHANES:

National Health and Nutrition Examination Survey
 SD:

stochastic dominance
 CPS:

Current Population Survey
 CDF:

cumulative distribution function
References
 1.
Flegal KM, Carroll MD, Kuczmarski RJ, Johnson CL. Overweight and obesity in the United States: Prevalence and trends, 19601994. Int J Obes Relat Metab Disord. 1998;22:39–47.
 2.
Flegal KM, Carroll MD, Ogden CL, Curtin LR. Prevalence and trends in obesity among US adults, 19992008. JAMA. 2010;303:235–41.
 3.
Kuczmarski RJ, Flegal KM, Campbell SM, Johnson CL. Increasing prevalence of overweight among US adults: The National Health and Nutrition Examination Surveys, 1960 to 1991. JAMA. 1994;272:205–11.
 4.
Ogden CL, Carroll MD. Prevalence of overweight, obesity, and extreme obesity among adults: United States, trends 1960–1962 through 2007–2008. National Cent Health Stat. 2010;6:1–6.
 5.
Ogden CL, Carroll MD, Kit BK, Flegal KM. Prevalence of childhood and adult obesity in the United States, 20112012. JAMA. 2014;311:806–14.
 6.
Wang Y, Beydoun MA. The obesity epidemic in the United States  gender, age, socioeconomic, racial/ethnic, and geographic characteristics: A systematic review and metaregression analysis. Epidemiol Rev. 2007;29:6–28.
 7.
Christakis NA, Fowler JH. The spread of obesity in a large social network over 32 years. N Engl J Med. 2007;357:370–9.
 8.
Ogden CL, Carroll MD, Flegal KM. High body mass index for age among US children and adolescents, 20032006. JAMA. 2008;299:2401–5.
 9.
Skelton JA, Cook SR, Auinger P, Klein JD, Barlow SE. Prevalence and trends of severe obesity among US children and adolescents. Acad Pediatr. 2009;9:322–9.
 10.
Sturm R. Increases in clinically severe obesity in the United States, 19862000. Arch Intern Med. 2003;163:2146–8.
 11.
Sturm R. Increases in morbid obesity in the USA: 2000–2005. Public Health. 2007;121:492–6.
 12.
Flegal KM, Troiano RP. Changes in the distribution of body mass index of adults and children in the US population. Int J Obes Relat Metab Disord. 2000;24:807–18.
 13.
Burke MA, Heiland F. Social dynamics of obesity. Econ Inq. 2007;45:571–91.
 14.
Oswald AJ, Powdthavee N. Book review feature: Two reviews of the challenge of affulence: Self‐control and well‐being in the United States and Britain since 1950. Econ J. 2007;117:F441–54.
 15.
Etilé F. Social norms, ideal body weight and food attitudes. Health Econ. 2007;16:945–66.
 16.
Dragone D, Savorelli L. Thinness and obesity: A model of food consumption, health concerns, and social pressure. J Health Econ. 2012;31:243–56.
 17.
Blanchflower DG, Landeghem B, Oswald AJ. Imitative obesity and relative utility. J Eur Econ Assoc. 2009;7:528–38.
 18.
Wadsworth T, Pendergast PM. Obesity (sometimes) matters: The importance of context in the relationship between obesity and life satisfaction. J Health Soc Behav. 2014;55:196–214.
 19.
Contoyannis P, Wildman J. Using relative distributions to investigate the body mass index in England and Canada. Health Econ. 2007;16:929–44.
 20.
Houle BC. Measuring distributional inequality: Relative Body Mass Index distributions by gender, race/ethnicity, and education, United States (1999–2006). J Obes. 2010;959658.
 21.
Madden D. A profile of Obesity in Ireland, 2002–2007. J R Stat Soc Series A. 2012;175:893–914.
 22.
Sahn DE. Weights on the rise: Where and for whom? J Econ Inequal. 2009;7:351–70.
 23.
Kolenikov S, Shorrocks A. A decomposition analysis of regional poverty in Russia. Rev Dev Econ. 2005;9:25–46.
 24.
Shorrocks AF. Decomposition procedures for distributional analysis: A unified framework based on the Shapley value. J Econ Inequal. 2013;11:99–126.
 25.
Deurenberg P. Universal cutoff BMI points for obesity are not appropriate. Br J Nutr. 2001;85:135–6.
 26.
Jolliffe D. The income gradient and distributionsensitive measures of overweight in the US. National Poverty Center Working Paper Series 0727 [http://ageconsearch.umn.edu/bitstream/25677/1/cp060622.pdf]. Accessed 26 Nov 2013.
 27.
Burkhauser RV, Beyond CJ, BMI. The value of more accurate measures of fatness and obesity in social science research. J Health Econ. 2008;27:519–29.
 28.
Ezzati M, Martin H, Skjold S, Vander Hoorn S, Murray CJ. Trends in national and statelevel obesity in the USA after correction for selfreport bias: Analysis of health surveys. J R Soc Med. 2006;99:250–7.
 29.
Gorber SC, Tremblay M, Moher D, Gorber B. A comparison of direct vs. self‐report measures for assessing height, weight and body mass index: A systematic review. Obes Rev. 2007;8:307–26.
 30.
Johnson CL, PauloseRam R, Ogden CL, Carroll MD, KruszanMoran D, Dohrmann SM, et al. National health and nutrition examination survey. Analytic guidelines, 19992010. Vital Health Stat. 2013;2:161.
 31.
National Center for Health Statistics. The Third National Health and Nutrition Examination Survey, NHANES III (1988–1994). http://www.cdc.gov/nchs/nhanes/nh3data.htm. Accessed on Feb 01 2016.
 32.
Mohadjer L, Montaquila J, Waksberg J, Bell B, James P, FloresCervantes I, Montes M. National Health and Nutrition Examination Survey III, Weighting and Estimation Methodology. Hyattsville, MD: National Center for Health Statistics; 1996.
 33.
Davidson R, Duclos JY. Statistical inference for stochastic dominance and for the measurement of poverty and inequality. Econometrica. 2000;68:1435–64.
 34.
Anderson G. Nonparametric tests of stochastic dominance in income distributions. Econometrica. 1996;64:1183–93.
 35.
Sahn DE, Stifel DC. Poverty comparisons over time and across countries in Africa. World Dev. 2000;28:2123–55.
 36.
Sahn DE, Stifel DC. Robust comparisons of malnutrition in developing countries. Am J Agr Econ. 2002;84:716–35.
 37.
Foster JE, Shorrocks AF. Poverty orderings. Econometrica. 1988;56:173–7.
 38.
Sahn DE, Stifel DC, Younger SD. Intertemporal changes in welfare: Preliminary results from nine African countries. CFNPP Working Paper No. 94 [http://pdf.usaid.gov/pdf_docs/pnacm776.pdf]. Accessed 01 Dec 2013.
 39.
Datt G, Ravallion M. Growth and redistribution components of changes in poverty measures: A decomposition with applications to Brazil and India in the 1980s. J Dev Econ. 1992;38:275–95.
 40.
Karagiannis E, Kovacevic M. A method to calculate the Jackknife variance estimator for the Gini coefficient. Oxford B Econ Stat. 2000;62:119–22.
 41.
Sehili S, Elbasha EH, Moriarty DG, Zack MM. Inequalities in self‐reported physical health in the United States, 1993‐1999. Health Econ. 2005;14:377–89.
 42.
Shorrocks, AF. Inequality decomposition by population subgroups. Econometrica. 1984;52:1369–1385.
 43.
Ruhm CJ. Current and future prevalence of obesity and severe obesity in the United States. Forum Health Econ Pol. 2007;10:1–26.
 44.
Ravallion M, Chen S. Measuring propoor growth. Econ Lett. 2003;78:93–9.
 45.
Frieze IH, Olson JE, Good DC. Perceived and actual discrimination in the salaries of male and female managers. J Appl Psychol. 1990;20:46–67.
 46.
Pingitore R, Dugoni BL, Tindale RS, Spring B. Bias against overweight job applicants in a simulated employment interview. J Appl Psychol. 1994;79:909–17.
 47.
Puhl R, Brownell KD. Bias, discrimination, and obesity. Obes Res. 2001;9:788–805.
 48.
Roehling MV. Weight‐based discrimination in employment: Psychological and legal aspects. Pers Psychol. 1999;52:969–1016.
 49.
Carpenter KM, Hasin DS, Allison DB, Faith MS. Relationships between obesity and DSMIV major depressive disorder, suicide ideation, and suicide attempts: Results from a general population study. Am J Public Health. 2000;90:251–7.
 50.
Friedman MA, Brownell KD. Psychological correlates of obesity: Moving to the next research generation. Psychol Bull. 1995;117:3–20.
 51.
Istvan J, Zavela K, Weidner G. Body weight and psychological distress in NHANES I. Int J Obes Relat Metab Disord. 1992;16:999–1003.
 52.
Onyike CU, Crum RM, Lee HB, Lyketsos CG, Eaton WW. Is obesity associated with major depression? Results from the Third National Health and Nutrition Examination Survey. Am J Epidemiol. 2003;158:1139–47.
 53.
Radloff LS. The CESD scale a selfreport depression scale for research in the general population. Appl Psych Meas. 1977;1:385–401.
 54.
Wang SS, Brownell KD, Wadden TA. The influence of the stigma of obesity on overweight individuals. Int J Obes. 2004;28:1333–7.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
TYP conceptualized the study, drafted the manuscript, and takes responsibility for the integrity of data handling and accuracy of statistical analysis. SF and GC contributed to the study conception and design, interpretation of data, critical revision of the manuscript, and study supervision. All authors read and approved the final manuscript.
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.
About this article
Cite this article
Pak, T., Ferreira, S. & Colson, G. Measuring and tracking obesity inequality in the United States: evidence from NHANES, 19712014. Popul Health Metrics 14, 12 (2016). https://doi.org/10.1186/s1296301600815
Received:
Accepted:
Published:
Keywords
 BMI distribution
 Obesity inequality
 Stochastic dominance
 Shapley decomposition
 Gini coefficient