Effects of neighbourhood income on reported body mass index: an eight year longitudinal study of Canadian children
BMC Public Health volume 8, Article number: 16 (2008)
This study investigates the effects of neighbourhood income on children's Body Mass Index (BMI) from childhood (ages 2–3) to early adolescence (ages 10–11) using longitudinal data.
Five cycles of data from the Canadian National Longitudinal Survey of Children and Youth are analyzed for a sub-sample of children (n = 2152) aged 2–3 at baseline (1994) and assessed at two year intervals to 2002. Body mass index percentiles are based on height/weight estimates reported by proxy respondents (child's person most knowledgeable). Family and neighbourhood factors were assessed at baseline. The prevalence of neighbourhood low income was obtained from the 1996 Census and divided into three categories from 'most poor' to 'least poor'. Longitudinal modelling techniques were applied to the data.
After controlling for individual/family factors (age, sex, income, education, family structure) living in the 'most poor' neighbourhood was associated with increasing BMI percentile (1.46, 95% CI 0.16 to 2.75) over time compared to a 'middle' income neighbourhood. Living in an urban (vs. rural) neighbourhood was associated with a decreased BMI percentile (-3.57, 95% CI -6.38 to -0.76) across all time periods.
These findings provide evidence that effects of neighbourhood disadvantage on children's BMI occur between childhood and early adolescence and suggest that policies should target the conditions of childhood, including the neighbourhood environment.
The prevalence of childhood obesity and overweight has increased dramatically over the last three decades in most developed nations . The prevalence of overweight among Canadian children has doubled from 13% in 1977/8 to 26% in 2004 among 6-to 11-year olds . These increases are troubling because obesity in childhood persists into adulthood and may be linked to poor long-term health outcomes [3, 4]. Cross-sectional research from developed nations has demonstrated an association between neighbourhood disadvantage and an increased prevalence of obesity and overweight among children and youth [5–8]. Longitudinal studies are needed to establish causation and to examine the influence of neighbourhood disadvantage on children's body weights over time. Longitudinal studies of neighbourhood effects on children's health that include both individual and neighbourhood characteristics are rare .
Disparities in obesity may emerge between early childhood and adolescence . Neighbourhoods provide a social and physical context which structures opportunities for children to engage in behaviours that promote or inhibit weight gain [11, 12]. Residing in a disadvantaged neighbourhood may promote weight gain through access to a less healthy food supply, limited access to recreation facilities and increased safety concerns [6, 11, 13–15]. Neighbourhood factors may be more influential as children age and have more freedom to access the neighbourhood .
Several cross-sectional studies have found that relationships between neighbourhood deprivation and overweight are stronger for older children than younger children [17–19]. Collectively these findings suggest that disparities between overweight and neighbourhood deprivation may emerge between early childhood and adolescence.
The primary goal of this study was to assess the impact of the neighbourhood environment on children's BMI from early childhood to adolescence while controlling for family factors. It was hypothesised that from early childhood to adolescence there would be increasing disparity in body weight by neighbourhood income as neighbourhood factors become more influential and children are exposed to such environments over time.
Data used were from the Canadian National Longitudinal Survey of Children and Youth (NLSCY), a representative survey of Canadian children aged 0–11 beginning in 1994. Children are assessed bi-annually until age 25 and data used are from the first 5 cycles. The NLSCY uses a complex survey design based on the Statistics Canada Labour Force Survey and a full description of the data is available elsewhere . The initial cohort included 2,916 children aged 2–3 in 1994 (Cycle 1) and we use data for 2,229 respondents aged 10–11 in 2002 (Cycle 5) representing 76.4% of the original cohort.
In the NLSCY, children's heights and weights were reported by the person most knowledgeable (PMK) for ages 2–11 and self-report for older children. In Cycle 1 the PMK was the mother in 90% of cases and the father in 8% of cases. We limited our sample to a cohort of children aged 2–3 in Cycle 1. At Cycle 5 these children are aged 10–11 and heights and weights are reported by the PMK across all cycles. All individual and family variables were obtained from Cycle 1. Heights and weights were used to calculate Body Mass Index (BMI – kilograms (kg)/meters (m)2), which was then used to derive the outcome variable. We assess BMI percentile using age and sex specific values from the US Centre for Disease Control (CDC) which have been applied to Canadian data [22, 23]. This method assigns percentile values to age and sex adjusted BMI z-scores at one month intervals. Based on the CDC cut-offs we identify those with a BMI percentile ≥ 85th as overweight. The CDC categorization refers to the 85th percentile as 'at risk of overweight' and the 95th percentile as 'overweight'. In this study a BMI percentile equal to or greater than the 85th percentile is referred to as overweight. Individual characteristics included in the model are gender (coded 0 = male and 1 = female) and the child's age at Cycle1 (coded 0 = age 2 and 1 = aged 3).
A variable indicates if children are living in an 'intact family' with both biological parents. Previous research with the NLSCY has demonstrated that children living in an intact family have higher levels of health and well-being . Children living in an intact family were coded as 1 and others as 0.
Income adequacy is a categorical measure of 1994 income based on the PMK report of household income during the previous 12 months . The income categories are based on household income and household size. For example, a family of 2 would be classified in the highest group if making $60,000 CDN or more and a 3 person family would be classified in the highest group if making $80,000 CDN or more. This variable originally was created with 5 categories (lowest, low-middle, middle, high middle, and high). Lowest and low-middle were merged due to small sample sizes in the lowest group. Middle and high middle were combined to create three categories: low, middle (omitted reference group) and high.
Education of the PMK in 1994 was determined from a variable that groups education into 4 categories based on self-report: less than high school, high school, some postsecondary and postsecondary degree or diploma. Three analytic categories were created: less than high school, high school/some post secondary (omitted reference group), and post secondary degree/diploma.
Enumeration Areas, the smallest level of census geography in 1996 containing between 125 and 440 dwellings, were used as neighbourhood proxies. Neighbourhood low-income was assessed from the 1996 Census by calculating the proportion of the non-institutional population living below the low income cut-off. Proportions were divided into quintiles and the three middle quintiles were grouped resulting in three categories: 'least poor', 'middle' (omitted reference category) and 'most poor'.
Whether the neighbourhood was urban or rural (coded 0 = rural and 1 = urban) was determined by a variable indicating if it was in a Census Metropolitan Area (CMA). A CMA consists of one or more municipalities that form an urban core greater than 100,000 residents .
While data were available for 2,229 children at Cycle 5 some were excluded because of missing data for PMK education, family structure, or because a postal code could not be linked to neighbourhood data. There were 77 exclusions leaving a final sample of 2,152. Across all 5 time periods there are 8,915 valid BMI observations out of a total possible 10,760 (5 × 2152) so each child has an average of 4.14 out of 5 measurements. There was some variation in the number of valid BMI responses by family income: children from low income families have an average of 3.85 of 5 possible measures, middle income children have 4.19, and high income children have 4.30.
Individual growth modelling techniques were used to analyze longitudinal data and the 5 waves of data available is adequate for modelling linear change . A two-level model was specified in which observations at each time point (level 1) are clustered within individuals (level 2). The level-1 model describes within-individual change over time and the level-2 model assesses if predictor variables are related to inter-individual differences in change. An advantage of using this method is that children without complete BMI data for each survey cycle can be included .
A linear growth model was specified and the generic 2 level model used is:
BMIit = π 0i + π 1i (Timeit) + ε it
where BMIit refers to the outcome, BMI percentile for each child, i, at occasion t. Timeit represents each measurement occasion for the ith child at time t and the initial measurement (Cycle 1, 1994) was assigned a value of zero. As such, π 0i represents child i's true BMI percentile at Cycle 1 and π 1i represents the slope of the change trajectory over a 2 year period for the ith child. Predictor variables were added as β coefficients on the intercept and change trajectory. Random effects were specified on the intercept π 0i and slope π 1i of the individual change trajectory with other coefficients entered as fixed effects. Iterative Generalized Least Squares was used for statistical estimation and final models used an unstructured co-variance matrix . Sandwich estimators were used to calculate robust standard errors . Following Singer and Willett  we used a sequential model building process beginning with a null or empty model and subsequently adding individual, family and neighbourhood variables. Six models are presented.
Descriptive statistics were computed using SAS© and longitudinal models were computed using MLwin 2.02© . Variance estimates (not shown) of descriptive statistics were calculated using 1000 bootstrap weights and sampling weights were applied to longitudinal analysis .
Table 1 shows the weighted descriptive statistics for the sample. Figure 1 shows that the 'most poor' neighbourhood had the highest percent overweight across all time periods while the 'middle' and 'least poor' neighbourhoods had values similar to each other. Figure 2 shows that the 'most poor' category had a higher BMI percentile than the 'middle' and 'least poor' and this disparity increases as children age.
Tables 2 and 3 present the results of the initial and full models and the β's are the coefficient outcomes. Model A (Table 2) is an unconditional means model that has no predictors and is used to determine baseline estimates for comparison with subsequent models. The model goodness of fit statistic, the -2 log likelihood (88140.50), is compared with subsequent models to evaluate the model fit with the addition of predictor variables. Random effects are presented at the bottom of the table. Model B (Table 2) is an unconditional growth model which includes a single predictor of time. The fixed effects for this model include the intercept and change over time.
Model C (Table 3) presents the results of an initial status model that tests if predictor variables have an effect on children's initial status of BMI percentile and was used to gauge the effects of predictors without the addition of time. The estimate for age indicates that children aged 3 had a lower BMI percentile of -2.34 (95% CI -4.35 to -0.33) than children aged 2. Children in an intact family had an estimated differential in BMI percentile of -3.84 (95% CI -7.06 to -0.62) at initial status. The estimates for PMK education are in the expected direction but the confidence intervals include zero. Children living in the 'most poor' neighbourhood have an estimated differential in BMI percentile at initial status of 4.03 (95% CI 1.20 to 6.86) and children living in a CMA have an estimated BMI percentile difference of -3.51 (95% CI -5.55 to -1.46) at initial status. The -2 log likelihood for Model C was lower than Model B (88051.77 vs. 88102.66) indicating that the addition of predictor variables at initial status improves the model fit.
Model D (Table 3) includes effects of individual and family level predictors on both initial status and rate of change. Similar to the initial status model, being aged 3 and living in an intact family was associated with a decrease in BMI percentile at initial status. Effects over time show that being female was associated with a decreased rate of change.
Model E removes the family level variables and includes individual variables and neighbourhood income and urban/rural residence. This model shows the effect of neighbourhood income and urban/rural residence on BMI percentile change unadjusted for family factors. Living in a CMA decreased BMI percentile on initial status by -3.45 (95% CI -6.20 to -0.69) but the confidence interval for the estimate over time includes zero suggesting that there was no time varying effect. The initial estimate for the 'most poor' neighbourhood was 2.15 (95% CI -1.640 to 5.934) and the confidence interval included zero. The time varying estimate for the 'most poor' neighbourhood was 1.36 (95% CI 0.07 to 2.65) suggesting that neighbourhood income effects occur over time. When CMA residence was added as a time varying effect the initial status estimate decreased slightly to -3.45 (95% CI -6.20 to -0.69) and the confidence interval for the estimate over time included zero suggesting there was no time varying effect.
Model F includes both family and neighbourhood characteristics on initial status and rate of change. The -2 log likelihood was the lowest (88028.16) indicating this model best fits the data. While many of the family level variable estimates had confidence intervals that included zero they were retained for theoretical reasons and they contributed to the overall model fit . Living in an intact family and living in a CMA was associated with a reduced BMI percentile at initial status. Living in the 'most poor' neighbourhood increased BMI percentile by 1.46 (95% CI 0.161 to 2.75) over time.
The contributions of this study are twofold. First, it contributes to research on neighbourhood effects on children's health by adopting a longitudinal perspective. Second, it contributes to our understanding of young children's trajectories of BMI by examining associations with both family and neighbourhood characteristics.
The principal finding of this paper is that the early neighbourhood environment influences children's BMI percentiles. Children living in the 'most poor' neighbourhoods have an increased rate of change in BMI percentile relative to children living in a 'middle' income neighbourhood. Living in the 'least poor' neighbourhood did not confer benefits suggesting that it is the effect of neighbourhood poverty rather than affluence that may matter most. The final model showed that over time living in the 'most poor' neighbourhood increases BMI percentile which is consistent with our hypothesis that neighbourhood characteristics may have a greater influence as children age. It is possible that neighbourhoods may become more important as children age and have more freedom to access the neighbourhood. It is also possible that disparities emerge over time as children are exposed to such environments over a longer time period.
A strength of this study is the availability of data for a nationally representative longitudinal sample of Canadian children allowing an examination of the influence of the early neighbourhood environment on BMI. A limitation is the use of PMK reported heights and weights which have well-known limitations compared to direct anthropological measures . Evidence suggests parents underreport children's height leading to increased BMI values compared to direct measurements . The decrease in BMI percentile over time is likely because even small underreports of height can considerably overestimate BMI in young children who are shorter. Another limitation is that parental BMI was not assessed in the NLSCY. A statistical limitation is that the few cases per neighbourhood prevented specifying a three level model.
During the 8 year follow-up children may have experienced changes in family income, parental education, family structure (divorce, remarriage etc) or may have moved neighbourhoods. Because we did not account for children who move neighbourhoods the model examines the influence of neighbourhood income at ages 2–3 on trajectories of BMI. Future research should examine the effects of moving and if children's trajectories differ based on the number of moves, timing of moves and whether moves are associated with changes in socioeconomic status. It is possible that moving to a neighbourhood with a different income profile could influence children's BMI trajectories. Over the 8 year period the percent of children in the lowest income category decreased from 19% to 7% and children in the highest category increased from 15% to 34%. Percent of children living in an intact family decreased from 81% to 69%. Accounting for these changes would have introduced significant model complexity.
Existing studies examining the influence of neighbourhood deprivation on overweight have been cross-sectional or have not included both family and neighbourhood variables [5, 6, 17, 19]. These studies have found stronger effects for older children. Our study is the first we are aware of to examine longitudinal associations between the family and neighbourhood environment and children's body weight. A longitudinal study has found that socio-economic disparities in obesity did not increase as children aged from 11–16 . The period from childhood to early adolescence may represent a critical period in which disparities in overweight by neighbourhood income are established. We did not find statistically significant results related to family income or PMK education reported in other Canadian studies [2, 30].
The findings of this research suggest that obesity policies which focus on conditions of childhood including the places in which young children live may meet with the greatest success. Such policies may reduce the prevalence of obesity among all children and prevent the emergence of neighbourhood-based disparities in body weight as children age.
This study raises many questions for future research. This study found that children in rural areas have higher BMI percentiles than children in urban areas but the disparity, unlike neighbourhood income, does not increase with age. Future longitudinal studies examining urban-rural disparities in body weight are needed. This study did not evaluate the mechanisms and pathways through which low income neighbourhoods influence children's body weight over time. It is likely that differential access to food choices and opportunities for physical activity underlie this relationship. Cross-sectional studies have found that disadvantaged neighbourhoods have reduced access to healthy food options and parents report fewer safe parks and playgrounds [6, 31, 32]. However, longitudinal studies examining the influence of specific neighbourhood factors over time are needed to better understand how such factors influence children's bodyweight. At present we do not know if the influence of particular neighbourhood characteristics (e.g. parks and playgrounds, traffic) are greatest during certain periods in childhood and adolescence. This type of research may lead to the development of policies aimed at reducing neighbourhood disparities in overweight. Future research with a longer follow-up should investigate the effects of neighbourhood income on BMI from early adolescence to young adulthood.
In conclusion, cross-sectional studies have demonstrated that neighbourhood disadvantage is associated with increased risk of overweight among children and youth. This study is the first we are aware of that uses longitudinal data to demonstrate that the disparity in BMI by neighbourhood income emerges between childhood and early adolescence. As obesity in childhood often persists into adulthood the findings of this study suggest that policies to prevent neighbourhood disparities in overweight should focus on young children.
Ebbeling CB, Pawlak DB, Ludwig DS: Childhood obesity: public-health crisis, common sense cure. Lancet. 2002, 360: 473-82. 10.1016/S0140-6736(02)09678-2.
Shields M: Overweight and obesity among children and youth. Health Rep. 2006, 17: 27-42.
Dietz WH: Health consequences of obesity in youth: Childhood. Pediatrics. 1998, 101: 518-
Field AE: Epidemiology of the health and economic consequences of pediatric obesity. Handbook of pediatric obesity: etiology, pathophysiology, and prevention. Edited by: Goran MI, Sothern M. 2006, Boca Raton: CRC/Taylor & Francis Group, 1-19.
Janssen I, Boyce WF, Simpson K, Pickett W: Influence of individual- and area-level measures of socioeconomic status on obesity, unhealthy eating, and physical inactivity in Canadian adolescents. Am J Clin Nutr. 2006, 83: 139-45.
Oliver LN, Hayes MV: Neighbourhood socio-economic status and the prevalence of overweight Canadian children and youth. Can J Public Health. 2005, 96: 415-20.
Moffat T, Galloway T, Latham J: Stature and adiposity among children in contrasting neighborhoods in the city of Hamilton, Ontario, Canada. Am J Hum Biol. 2005, 17: 355-67. 10.1002/ajhb.20125.
Veugelers PJ, Fitzgerald AL: Prevalence of and risk factors for childhood overweight and obesity. CMAJ. 2005, 173: 607-13.
Schneiders J, Drukker M, van der Ende J, Verhulst FC, van Os J, Nicolson NA: Neighbourhood socioeconomic disadvantage and behavioural problems from late childhood into early adolescence. J Epidemiol Community Health. 2003, 57: 699-703. 10.1136/jech.57.9.699.
Strauss RS, Knight J: Influence of the home environment on the development of obesity in children. Pediatrics. 1999, 103: e85-10.1542/peds.103.6.e85.
Lake A, Townshend T: Obesogenic environments: exploring the built and food environments. The Journal of the Royal Society for the Promotion of Health. 2006, 126: 262-7. 10.1177/1466424006070487.
Macintyre S, Ellaway A, Cummins S: Place effects on health: how can we conceptualise, operationalise and measure them?. Social Science & Medicine. 2002, 55: 125-39. 10.1016/S0277-9536(01)00214-3.
Lee RE, Cubbin C: Neighborhood context and youth cardiovascular health behaviors. Am J Public Health. 2002, 92: 428-36.
Molnar BE, Gortmaker SL, Bull FC, Buka SL: Unsafe to play? Neighborhood disorder and lack of safety predict reduced physical activity among urban children and adolescents. Am J Health Promot. 2004, 18: 378-86.
Lumeng JC, Appugliese D, Cabral HJ, Bradley RH, Zuckerman B: Neighborhood safety and overweight status in children. Arch Pediatr Adolesc Med. 2006, 160: 25-31. 10.1001/archpedi.160.1.25.
Leventhal T, Brooks-Gunn J: The neighborhoods they live in: the effects of neighborhood residence on child and adolescent outcomes. Psychol Bull. 2000, 126: 309-37. 10.1037/0033-2909.126.2.309.
Kinra S, Nelder RP, Lewendon GJ: Deprivation and childhood obesity: a cross sectional study of 20,973 children in Plymouth, United Kingdom. J Epidemiol Community Health. 2000, 54: 456-60. 10.1136/jech.54.6.456.
Timperio A, Salmon J, Telford A, Crawford D: Perceptions of local neighbourhood environments and their relationship to childhood overweight and obesity. Int J Obes. 2005, 29: 170-5. 10.1038/sj.ijo.0802865.
Wardle J, Brodersen NH, Cole TJ, Jarvis MJ, Boniface DR: Development of adiposity in adolescence: five year longitudinal study of an ethnically and socioeconomically diverse sample of young people in Britain. BMJ. 2006, 332: 1130-5. 10.1136/bmj.38807.594792.AE.
Human Resources Development Canada and Statistics Canada: User's Handbook and Microdata Guide: National Longitudinal Survey of Children and Youth, Cycle 1. 1997, Ottawa: Special Surveys Division
Kuczmarski RJ, Ogden C, Guo S: 2000 CDC growth charts for the United States: Methods and development. National Center for Health Statistics. Vital Health Stat. 2001, 246: 1-190.
Janssen I, Katzmarzyk PT, Srinivasan SR, Chen W, Malina RM, Bouchard C, et al: Utility of childhood BMI in the prediction of adulthood disease: comparison of national and international references. Obes Res. 2005, 13: 1106-15.
Phipps SA, Lethbridge L, Burton P: Long-run consequences of parental paid work hours for child overweight status in Canada. Soc Sci Med. 2006, 62: 977-86. 10.1016/j.socscimed.2005.06.054.
Kerr D: Family Transformations and The Well-being of Children: Recent Evidence from Canadian Longitudinal Data. Journal of Comparative Family Studies. 2004, 35: 73-90.
Statistics Canada: 1996 Census Dictionary – Final Edition. 1999, Ottawa: Statistics Canada
Singer J, Willett J: Applied longitudinal data analysis: modeling change and event occurrence. 2003, New York: Oxford University Press
Rasbash J, Browne W, Goldstein D, Yang M, Plewis I, Healy M, Woodhouse G, Draper D, Langford I, Lewis T: A User's Guide to MLWIN. 2000, London: University of London
Power C, Lake JK, Cole TJ: Review: Measurement and long-term health risks of child and adolescent fatness. Int J Obes. 1997, 21: 507-10.1038/sj.ijo.0800454.
Human Resources Development Canada and Statistics Canada: Microdata User Guide: National Longitudinal Survey of Children and Youth, Cycle 5. 2003, Ottawa: Special Surveys Division
Willms JD, Tremblay MS, Katzmarzyk PT: Geographic and demographic variation in the prevalence of overweight Canadian children. Obes Res. 2003, 11: 668-73.
Pearce J, Blakely T, Witten K, Bartie P: Neighborhood deprivation and access to fast-food retailing: a national study. Am J Prev Med. 2007, 32: 375-82. 10.1016/j.amepre.2007.01.009.
Reidpath D, Burns C, Garrard J, Mahoney M, Townsend M: An ecological study of the relationship between social and environmental determinants of obesity. Health & Place. 2002, 8: 141-5. 10.1016/S1353-8292(01)00028-4.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/8/16/prepub
We would like to thank the four reviewers for their comments and suggestions. Lisa Oliver was supported through the Statistics Canada Tom Symons Doctoral Fellowship Program.
The author(s) declare that they have no competing interests.
LO conceived of the idea for this project, conducted the analysis and prepared the manuscript. MH assisted with conceptualizing the paper and preparing the manuscript.
About this article
Cite this article
Oliver, L.N., Hayes, M.V. Effects of neighbourhood income on reported body mass index: an eight year longitudinal study of Canadian children. BMC Public Health 8, 16 (2008). https://doi.org/10.1186/1471-2458-8-16