Abstract
Objective:
The purpose of this longitudinal study was to examine the effects of state-level socioeconomic status (SES), the density of fast food restaurants and walking to work on body mass index (BMI) among US adults aged 50 years and older. The study sought further to account for the interaction effects of three different hierarchical levels of variables, including time-varying variables and time-invariant variables such as individual and state-level variables. Confounding related to environmental/neighbourhood effects was controlled for through sample selection.
Methods:
Data were drawn from the 2000–2010 Health and Retirement Study (HRS) and the 2010 US Census. HRS survey participants aged 50 years and older in 2000, residing in the same census tract during 2000–2010 (N = 6,156) comprised the study sample. A three-level growth model was fitted to BMI trajectories.
Results:
BMI significantly increased over time showing both a linear and quadratic decelerating trajectory (p < .001) where BMI grew faster among sampled adults in their early 50s. Women living in the US states with a lower state-level SES had higher BMIs, while men with higher state-level SES had a higher BMI (b = −.242, p = .013). In states with a higher proportion of people who walked to work, Hispanics had higher BMIs at baseline and a lower BMI growth rate, compared with non-Hispanic whites (b = −.033, p < .001).
Conclusions:
Obesity interventions specific to adults aged 50 years and older require greater emphasis. Potential health inequalities regarding weight gain can be ameliorated through tailoring interventions based on sex, environmental and state-level SES influences.
Introduction
Obesity is a well-known public health issue both in the USA and worldwide (Flegal et al., 2012; Popkin et al., 2012). Middle-aged and older adults show the highest obesity rate compared with those of other age groups regardless of sex and race/ethnicity (Flegal et al., 2012). A strong body of evidence shows that clinical consequences of overweight and obesity, such as type 2 diabetes, are more serious and particularly problematic for older adults (Han et al., 2011). Obesity coupled with ageing increases healthcare service expenditure (Han et al., 2011), which increases the healthcare burden to society. Therefore, researchers continue to investigate a variety of obesity determinants such as built environment or obesogentic environment (e.g. access to fast food outlets and neighbourhood walkability), as well as micro- and macrolevel socioeconomic status (SES) to tackle this issue (Lakerveld et al., 2012), especially for middle-aged and older adults.
The built environment can be considered all buildings and space created, or modified to form the physical characteristics of a community. It may include a variety of items such as roads, homes, buildings, parks and recreation areas (Centers for Disease Control and Prevention, 2013). An obesogenic environment refers to an environment which is defined as the ‘sum of influences that the surroundings, opportunities, or conditions of life have on promoting obesity in individuals or populations’ (Swinburn et al., 1999). Although many obesogenic environment studies have been conducted, the association of obesogenic environments and obesity has not been fully studied due to the complexity of measurement options and potential confounding effects from the interaction between environmental and socioeconomic factors (Kirk et al., 2010), especially for middle-aged and older adults. Thus, a better understanding of obesity determinants for adults aged 50 years and older including obesogenic/socioeconomic factors is necessary to help reduce the high obesity prevalence and alleviate the overall burden of obesity-related diseases and associated healthcare costs.
Although there have been a few longitudinal studies concerning body mass index (BMI) trajectories among middle-aged and older adults, their results are inconclusive. One longitudinal study reports that BMI, on average, continues to increase and plateaus around 60–70 years of age (Villareal et al., 2005). Another longitudinal study reports that BMI increases up to 75 years of age (Botoseneanu and Liang, 2010). Furthermore, research findings are unclear concerning sex differences in weight changes among middle-aged and older adults (Botoseneanu and Liang, 2010) and on the associations between obesity rates and macrolevel environmental variables such as food outlets or walking environment in a geographic area (Ding and Gebel, 2012; Fan et al., 2011; Ruel et al., 2010; Winkleby and Cubbin, 2003).
Moreover, previous studies have revealed that obesity or BMI trends varied across geographic areas such as states and regions (Wang and Beydoun, 2007). States are unique entities with their own governance and authority to legislate state-specific laws and health policies. While some studies have examined the associations between weight gain and macrolevel environmental variables such as density of fast food restaurants and proportion of people who walk to work among younger adults, few studies have examined how these factors affect adults over the age of 50 (Mujahid et al., 2005; Pucher et al., 2010; Ruel et al., 2010). In addition, there are very limited numbers of studies regarding microlevel obesity factors (e.g. racial/ethnic segregation and sex) and their moderating effects on the association between macrolevel obesity determinants (e.g. state-level SES) and BMI (Mobley et al., 2006).
All these findings in the literature converge on the need to investigate the effect of state-level variables on the trends of BMI among a nationally representative sample of middle-aged and older adults. Given the aforementioned literature findings, albeit some are limited, our research questions were as follows:
How do longitudinal trajectory patterns of BMI differ by sex and increase over time in a sample of US adults aged 50 years or older?
How are state-level variables (i.e. state-level SES, density of fast food restaurants and proportion of people who walk to work) associated with longitudinal trajectory patterns of BMI?
Do the effects of the state-level variables on individuals’ BMI trajectories vary by demographics?
Does sex or race/ethnicity have a moderating effect on the association between state-level SES and obesity?
Methods
Data
The data used for this study were drawn from the US Health and Retirement Study (HRS), a nationwide longitudinal study of health, retirement and ageing among US adults aged over 50 years. It has been a leading source for information on the health and economic well-being of the aforementioned age group with funding support by the US National Institute on Aging (NIA U01AG009740) and the Social Security Administration (Health and Retirement Study, 2016). The HRS is an ongoing biennial research project that has been collecting multiple waves of data since 1990. Moreover, it takes 2 years to complete data collection for a specific wave of data in order to obtain intensive demographic and health-related data (e.g. age, sex, education, race/ethnicity, marital status, physical and mental health status, BMI). Therefore, data collected from 1990 to 1992 were considered the first wave of data. The HRS has collected 12 waves of data (i.e. from the first wave [1990–1992] to the 12th wave [2012–2014]) since the inception of the study.
The HRS survey response rates (2000–2010) ranged from 81.0% (2008–2010) to 88.9% (2004–2006) (Health and Retirement Study, 2017). The investigators linked the data (2000–2010) to the 2010 US Census data based on the census tract information of the HRS respondents. Although census data were available in 2000 and 2010, only the ‘restricted’ geographic information system (GIS) data regarding HRS subjects’ census tract information from 2010 Census were accessible to the researchers in the present study. Our research team had to develop an extensive and secure data protection plan to get the approval from the HRS to access the above GIS data. Therefore, 2010 Census data and 2010 HRS GIS information were linked for data analyses. The linkage allowed for investigation of the impact of statewide percentage of people who walk to work and density of limited-service restaurants, which is a proxy measure for density of fast food restaurants, on study participants’ BMI trajectories. Permission to use the participants’ census tract information was obtained from the HRS and the study protocols were approved by the Indiana University Institutional Review Board.
Study sample
Our goal was to select HRS respondents living in the same census tract from 2000 to 2010, which gave this study a good degree of control of each participant’s living environment and surrounding neighbourhood. This sample selection approach minimised a variety of confounding effects due to exposure to different environmental and neighbourhood factors (e.g. density of fast food restaurants and neighbourhood walkability). Figure 1 below illustrates each step of the sample selection process.

Flow chart of the study sample selection.
Based on the above selection criterion, the initial sample consisted of 7,664 participants. Since the HRS recruited study participants aged 50 years and older and their spouses of any age, some HRS respondents’ spouses were less than 50 years old. Thus, we excluded 225 primary participants’ spouses who were ineligible because they were younger than 50 as of Year 2000. Also excluded were the 622 individuals who died between 2000 and 2010. Another 556 respondents were excluded due to incomplete surveys resulting in missing data. Our analyses were limited to non-Hispanic whites, non-Hispanic blacks and Hispanics due to a very small percentage of people reporting their race/ethnicity as ‘others’; therefore, an additional 102 subjects were eliminated. Finally, three multivariate outliers determined by Mahalanobis distances were also excluded.
Analysis was conducted based on listwise deletion for cases with missing data which did not show systematic pattern according to Little’s MCAR (Missing Completely at Random) test (Dickson and Flynn, 2008; Little, 1988; Little and Rubin, 2002; Outhwaite and Turner, 2007; Seo and Li, 2009). In terms of quantity, less than 10% of the study participants had missing data. Thus, potential loss of statistical power due to the exclusion of missing data may not have biased estimation (Acock, 2015). In summary, the sample selection process resulted in a final study sample of 6,156 participants at baseline who were followed for 10 years.
Variables and measurement
The dependent variable in this study was BMI computed by weight (lb)/[height (in)]2 × 703 using self-reported weight and height information. There were six waves of BMI data from 2000 to 2010. The independent variables were categorised into three different hierarchical levels: time-varying individual changes (Level 1 variables), time-invariant individual variables (Level 2 variables), and time-invariant state-level variables (Level 3 variables). Year was used as a time variable, which was coded from 0 to 10. Level 2 variables included age in 2000, number of years of schooling completed, sex and race/ethnicity. Sex was coded as 0 for men and 1 for women. Race/ethnicity was coded as 0 for non-Hispanic whites, 1 for non-Hispanic blacks, and 2 for Hispanics. Level 3 variables included state-level SES, density of fast food restaurants (per 100,000 people), and percent of people who walked to work from the census data in 2010. All Level 3 variables were grand-mean centred to facilitate interpretation of multilevel models (Hox, 2002). The state-level rather than census tract level was used as the unit of analysis because (1) the states may have unique policies with their own governing power and influence and (2) there were sufficient sample sizes at the state level, not at the census tract level, to provide reliable estimates for the parameters in multilevel modelling (Hox, 2002).
The state-level SES was computed by taking the mean of SES indices of all census tracts in each state. The SES indices were based on median household income, educational attainment and poverty defined by the US Census Bureau (Gomez et al., 2007). Moreover, the state-level SES indices were reverse coded so that a high score indicated a higher level of SES. Consequently, the state-level SES represented a composite scale score reflecting SES for each state and ranges from a score of 3 (the worst) to 12 (the best).
In addition, limited-service restaurants from the census data have been used as a proxy measure for fast food restaurants (Lopez, 2007). Thus, in this study, state-level density of fast food restaurants was calculated from the number of limited-service restaurants divided by the total population in each state. Finally, state-level proportion of people who walk to work is a measure of built environment for physical activity (Pucher et al., 2010). In other words, this variable can be referred to as an obesogenic environment of physical inactivity (Waine, 2005). A higher proportion of people who walk to work indicates an environment with better walkability that encourages people to walk (Lovasi et al., 2009; Smith et al., 2008). State-level proportions of people who walked to work were computed by aggregating percentage of people who walked to work at the census tract level.
Statistical analysis
We fitted three-level hierarchical linear models using HLM 7.0 (Scientific Software International, Lincolnwood, IL, USA) to analyse longitudinal trajectories of BMI over 10 years with different levels of impact from temporal effects (the time-varying variable [i.e. BMI] at Level 1), individual sociodemographic characteristics (time-invariant variables [e.g. sex, race/ethnicity] at Level 2) and state-level characteristics such as state-level SES, density of fast food restaurants and proportion of people who walk to work (time-invariant variables [i.e. state-level SES, density of fast food restaurants and proportion of people who walk to work] at Level 3) where all 50 states had data for our analyses. HLM 7.0 statistical software was used as it provides great flexibility in capturing the changing patterns from the repeated measures within individuals and appropriately accounts for correlated errors within nested units of data collection (Raudenbush and Bryk, 2002).
Our modelling strategies were guided by the modelling steps in multilevel modelling literature (Raudenbush and Bryk, 2002). Full maximum likelihood estimation was used to describe the fit of both fixed effects and variance components. Log-likelihood ratio tests for model chi-square differences were used for model selection decisions. As the weighting criteria in the HRS were not based on the dependent variable (i.e. BMI), but rather the independent variables in the present study, sampling weights were not applied for analysis. According to Winship and Radbill (1994), ‘When sampling weights are only a function of independent variables included in the model being estimated, unweighted ordinary least squares (OLS) regression will be the appropriate course to take’.
Results
Table 1 outlines the baseline descriptive statistics of the study sample. Of the 6,156 participants, 61.6% were women; 79.0% were non-Hispanic whites. The minimum age was 50 years, the maximum age was 93 years at baseline (i.e. 2000), and the mean age was 64.5 years (SD = 7.9). The mean BMI was 27.5 (SD = 5.1) in 2000. The mean state-level density of fast food restaurants was 89 per 100,000 people (SD = 1.40) and the mean state-level percentage of people who walked to work was 3.45% (SD = 1.56). Figure 2 illustrates BMI trajectories by sex over the 10 years. Men consistently showed a higher level of BMI than women.
Baseline characteristics of the study sample (N = 6,156, age ≥50 years), Health and Retirement Study a .
SES: socioeconomic status.
The reported percentages are weighted percentages.
State-level SES is a composite scale score reflecting socioeconomic status for each state and ranges from 3 (the worst) to 12 (the best).

Trends in body mass index from 2000 to 2010: Results of six waves of data from Health and Retirement Study (HRS, 2016).
Table 2 illustrates the results from three-level hierarchical linear models. The first intraclass correlation (ICC) coefficient was 0.81 (25.27/[3.56 + 25.27 + 2.34]), which means that 81% of the total variation in BMI trajectory was attributable to differences between different individuals. The second ICC was 0.08 (2.34/[3.56 + 25.27 + 2.34]), which means that 8% of the total variation in BMI trajectory was attributable to differences between different states. Both ICCs justified the use of multilevel models. The unconditional growth model (Model 2) showed significant linear and quadratic decelerating growth in BMI. Women had a lower initial BMI (p < .001), but showed a faster linear growth rate of BMI (p < .001) than men in Models 3 and 4 (b = −.669 and b = −.668, respectively), supporting the first hypothesis that individual trajectory patterns differ by sex and tend to increase over time in those aged 50 years and older.
Estimates from the three-level multilevel models for body mass index, Health and Retirement Study (2000–2010).
SES: socioeconomic status.
State-level SES is a composite scale score reflecting socioeconomic status for each state and ranges from 3 (worst) to 12 (best).
Each model was compared to the immediate previous model.
p < .05, **p < .01, ***p < .001.
The second hypothesis that state-level variables are associated with individuals’ BMI trajectories was partially supported. When sex differences were not accounted for, BMIs at baseline and their trajectories did not vary by the state-level SES. Neither the state-level density of fast food restaurants nor the proportion of people who walked to work was significantly associated with individual BMI trajectories. However, the effect of state-level SES on individual BMI was significantly different by sex (b = −.242, p = .013) in Model 4. As shown in Figure 3, sex moderated the effect of state-level SES on BMI. There was a gradient in BMI across levels of state-level SES for both women and men, but in the opposite direction. Based on the study results, BMI was similar for men and women in states with low SES, but was increasingly divergent, where BMI was higher for men and lower for women, as state-level SES increased.

Cross-level interaction between sex and state-level SES for BMI.
The three-way interaction between year, Hispanics (vs. non-Hispanic whites), and state-level percentage of people who walked to work was significant (b = −.033, p < .001). To probe this interaction, the relations between the three variables were plotted in Figure 4. Hispanics living in states with a lower proportion of people who walked to work showed a higher growth in BMI over time than Hispanics living in states with a higher proportion of people who walked to work. On the other hand, the other three-way interactions involving Level-3 variables were not significant. Based on our findings, the third hypothesis that the effects of the three state-level variables on individuals’ BMI trajectories vary depending on demographic factors was partially supported.

Three-way cross-level interaction between year, Hispanic (vs. Non-Hispanic white), and the proportion of people who walk to work on BMI trajectories.
Discussion
The findings of this study contribute to existing literature on BMI-related health disparities by investigating cross-level interactions between microlevel and macrolevel obesity determinants in a nationally representative sample of US middle-aged and older adults. This study investigated when, where and for whom certain obesogenic environments are more influential to weight gain and obesity using a longitudinal three-level hierarchical linear modelling approach. It examined multiple contextual relationships between time, state-level SES and obesogenic environment, individuals’ demographic factors, and BMI changes over a 10-year period of time among US adults aged 50 years and older.
Ageing and BMI
Findings revealed evidence that individuals’ BMIs increase over time at a slightly decelerating rate among people aged 50 years and older. This is consistent with the results from Botoseneanu and Liang (2010), and further strengthens the evidence that suggests both adult men and women continue to gain weight, although more slowly, as they age. Knowing that BMI trajectory for older persons is growing at a decelerating rate, the obesity prevention effort should focus more on late middle-aged and young older. Proactive and aggressive actions should be taken to reduce obesity prevalence for this age group.
State-level SES and BMI
Findings also suggest that there is a sex difference in the influence of state-level SES on BMI trajectories, which supports previous research findings that women living in more socioeconomically disadvantaged areas have higher BMIs (Mujahid et al., 2005). A possible explanation is that women might perceive the influence from their surrounding environment (e.g. perceived safety to go for a walk) differently from men. Such differential perception has been reported to be associated with different health outcomes (Bengoechea et al., 2005). Another possible explanation is that women might be more dependent than men upon resources from the surrounding environment to maintain their effort to control weight (Stafford et al., 2005). Women are usually more connected than men to the local environment and spend more time on food acquisition/preparation and shopping; thus, women’s dependence on surrounding environment tends to be stronger than men (Matheson et al., 2008). Women living in socially disadvantaged states might be particularly vulnerable to unhealthy conditions (Winkleby and Cubbin, 2003). By contrast, men who lived in states with higher state-level SES had higher BMIs. This phenomenon may have to do with a lack of time for exercise, sedentary work environment or more opportunities to eat at restaurants for men compared with women (Chang and Lauderdale, 2005). These findings imply that obesity prevention interventions and initiatives targeted at middle-aged or older adults need to place greater emphasis on services for women in states with low SES and services for men in states with high SES. More research is warranted to investigate the mechanisms by which state-level SES affects weight gains differently between men and women.
Obesogenic environment and BMI
Our findings indicate that state-level density of fast food restaurants and proportions of people who walk to work might not be associated with individual’s BMI or increases in obesity among people aged 50 years and older. However, it deserves mention that there was a difference between racial/ethnic groups in the way that proportions of people who walked to work affected longitudinal BMI changes. As mentioned earlier, the percentage of people who walk to work represents a measure of obesogenic environment of physical inactivity or area walkability (Lovasi et al., 2009; Pucher et al., 2010).
More specifically, Hispanics living in states with a higher proportion of people who walked to work had a lower BMI growth rate. According to a recent American Community Survey Report (McKenzie, 2014), Hispanics comprise one of the top two racial/ethnic groups that have the highest percentages of people who walk to work. The relatively high percentages might be due to financial necessity, lower rates of car ownership or may be more likely to live near walkable urban areas (Besser and Dannenberg, 2005). Plaut (2004) also pointed out that people who walk to work usually have lower incomes and are less educated, indicators which are usually associated with a higher BMI. The finding of the present study indicates that higher area walkability may be beneficial to weight control at the population level. Given that Hispanics overall have higher obesity rates than other racial/ethnic groups and that the fastest growing segment of the US population is Hispanic, this is a significant piece of information that emphasises the potential for regular physical activity to promote weight control among Hispanics. Even without vigorous-intensity physical activity or lengthy workouts, people may be able to lose weight or decelerate weight gains with consistent physical activity (Haskell et al., 2007). Another implication of this finding is the importance of the built environment. In order to improve the built environment for physical activity, policies that encourage development of a more walkable environment need to be implemented and promoted. State governments need to consider investing in age-friendly walkability programmes and pedestrian-oriented commercial spaces. For instance, a greater integration of land use and street connectivity in community planning and housing development could bring commercial and residential areas closer to remove barriers to walking (Saelens and Handy, 2008).
With regard to the density of fast food restaurants, the hypothesis that such a state-level variable would be associated with longitudinal trajectory patterns of BMI was not supported. Actually, there are contradictory findings on the relationship between the density of fast food restaurants and an increase in body weight of people living in the respective geographical area. Some studies reported significant associations between density of fast food restaurants and obesity, but they were conducted either with cross-sectional designs (Fleischhacker et al., 2011) or in self-selected areas, which may have introduced selection bias (Ding and Gebel, 2012). By contrast, a number of other studies did not find a significant relationship between the two (Lopez, 2007). Given that no population-based studies to date have revealed the extent to which the density of fast food restaurants is linked to the frequency of individual visits to fast food restaurants or caloric intake among patrons in this 50-plus-year-old population, future research is warranted to investigate these immediate connections that may lead to weight gains as well as to corroborate the finding of the present study regarding the longitudinal effect of the density of fast food restaurants in this population (Pérez-Escamilla et al., 2012).
Limitations
Our findings should be interpreted in the light of several study limitations. First, BMIs for HRS respondents were calculated from self-reported weights and heights. Underestimation of weight and overestimation of height might have confounded our results (Gunnell et al., 2000). However, the longitudinal design in the present study may have compensated for the bias because comparing the same person’s weight and height over time is likely to minimise the errors from self-reported figures (Bowman and DeLucia, 1992). Second, over the study period, some state-level data may have varied but the present study treated state-level data as static. Thus, any possible differential influence of changes in the state-level data might have confounded our findings. Third, the death of some cohort participants during the study period, albeit a small proportion of the cohort, may have introduced bias.
Conclusion
Despite the above limitations, this study contributes to the literature by shedding light on the heterogeneity of BMI change trajectories during the transition to older adulthood. To the best of our knowledge, this study was one of the first attempts to investigate BMI trajectories in a US representative sample of adults aged 50 years and older using three-level growth models where cross-level interactions were probed. BMI significantly increased over time showing a quadratic decelerating trajectory (i.e. growing faster when people were at their younger age). Thus, more intensive obesity prevention interventions should be held for middle-aged adults because BMI growth rates peak at those ages. In addition, evidence suggests that health inequality may exist due to the variation in weight increasing patterns from different sex, environmental and socioeconomic factors. In order to decrease health inequality, public health professionals and policymakers should take into account the aforementioned factors carefully and create more effective obesity interventions such as providing better walking environments. Finally, our study expands existing research on macrolevel obesity determinants by examining the impact of obesogenic environments. Future studies may need to investigate macrolevel obesity factors using different measures of obesogenic environment and by attempting to track BMI for a longer time period.
Footnotes
Acknowledgements
The data use of the Health and Retirement Study (HRS) respondents’ census tract information was approved by the HRS.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
