Article Text

Download PDFPDF

Original research
Understanding neighbourhood retail food environmental mechanisms influencing BMI in the Caribbean: a multilevel analysis from the Jamaica Health and Lifestyle Survey: a cross-sectional study
  1. Colette Andrea Cunningham-Myrie1,
  2. Novie O Younger2,
  3. Katherine P Theall3,
  4. Lisa-Gaye Greene4,
  5. Parris Lyew-Ayee4,
  6. Rainford Wilks2
  1. 1Department of Community Health and Psychiatry, University of the West Indies, Mona, Jamaica
  2. 2Caribbean Institute for Health Research, University of the West Indies, Mona, Jamaica
  3. 3Global Community Health and Behavioral Sciences, Tulane University, New Orleans, Louisiana, USA
  4. 4Mona GeoInformatics Institute, University of the West Indies, Mona, Jamaica
  1. Correspondence to Dr Colette Andrea Cunningham-Myrie; colette.cunninghammyrie{at}uwimona.edu.jm

Abstract

Objective To derive estimates of the associations between measures of the retail food environments and mean body mass index (BMI) in Jamaica, a middle-income country with increasing prevalence of obesity.

Design Cross-sectional study.

Setting Data from the Jamaica Health and Lifestyle Survey 2008 (JHLS II), a nationally representative population-based survey that recruited persons at their homes over a 4-month period from all 14 parishes and 113 neighbourhoods defined as enumeration districts.

Participants A subsample of 2529 participants aged 18–74 years from the JHLS II who completed interviewer-administered surveys, provided anthropometric measurements and whose addresses were geocoded.

Primary outcome measure Mean BMI, calculated as weight divided by height squared (kg/m2).

Results There was significant clustering across neighbourhoods for mean BMI (intraclass correlation coefficients=4.16%). Fully adjusted models revealed higher mean BMI among women, with further distance away from supermarkets (β=0.12; 95% CI 8.20×10−3, 0.24; p=0.036) and the absence of supermarkets within a 1 km buffer zone (β=1.36; 95% CI 0.20 to 2.52; p=0.022). A 10 km increase in the distance from a supermarket was associated with a 1.7 kg/m2 higher mean BMI (95% CI 0.03 to 0.32; p=0.020) in the middle class. No associations were detected with fast-food outlets or interaction by urbanicity.

Conclusions Higher mean BMI in Jamaicans may be partially explained by the presence of supermarkets and markets and differ by sex and social class. National efforts to curtail obesity in middle-income countries should consider interventions focused at the neighbourhood level that target the location and density of supermarkets and markets and consider sex and social class-specific factors that may be influencing the associations.

  • epidemiology
  • nutrition & dietetics
  • public health
  • diabetes & endocrinology
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

View Full Text

Statistics from Altmetric.com

Strengths and limitations of this study

  • This study is the first in a Caribbean island to demonstrate the influence of the retail food environment on body mass index (BMI) using geocoded data and multilevel modelling.

  • This study provided a large sample size representative of Jamaicans 15–74 years.

  • Individual geocoded addresses from a nationally representative survey were linked with specific objective geographic information system-based retail food environmental measures.

  • Enumeration districts were used to define Jamaican neighbourhoods which are quite heterogeneous geographically in size, composition and context, and may not fully represent exposure to the food obesogenic environments.

  • The reliability and validity of the area-level environmental variables used were not ascertained for the local context and therefore they may not be the most effective in explaining any variance in BMI.

Introduction

Obesity has been increasing in the Caribbean1 2 and in Jamaica is now a major public health problem.3 Over the past five decades a rapid increase in obesity has been reported with women having consistently higher rates than men.4–6 The Jamaica Health and Lifestyle Survey 2008 (JHLS II)3 has documented increased prevalence over the earlier 2001 survey (JHLS I),7 in obesity as well as the comorbid chronic non-communicable disease (NCD) conditions of diabetes mellitus and hypertension. Approximately 99% of Jamaicans consumed below the daily recommended portions of fruits and vegetables and 30% of obese persons preferred fried protein in their diets.3

The presence of supermarkets/markets has been thought to indicate better access to and intake of healthier foods, given its association with higher intake of fruits and vegetables8 and inversely associated with obesity.9–11 For example, among Canadian children residing in Toronto, those who lived in close proximity to a supermarket had decreased odds of being overweight or obese.12

The presence of supermarkets has been shown to be inversely associated with neighbourhood socioeconomic status (SES) in the USA, whereby data have revealed greater poverty being associated with a decreased presence.13 14 With regard to sex differences, research by Wang et al15 among adults aged 25–74 years in California revealed that closer proximity to a supermarket and higher neighbourhood density of small grocery stores were associated with higher body mass index (BMI) among women.

Fast-food outlets (FFO) have increased in many countries and thought to be associated with the global rise in obesity. While there is no universally accepted agreement on what the definition of fast food is, most research includes foods sold that are low cost, energy dense with high fat and/or sugar content and low nutrient content. Studies have found that frequent consumption of fast foods in areas with a high density of FFOs has been found to increase body weight16 17 and a positive association of proximity to FFOs with measures of adiposity.18 19

The interaction between SES and the density of FFOs has also been investigated. In Australia, it was found that persons with poor SES (based on median weekly income) had 2.5 times exposure to FFOs than persons in the wealthiest SES category.20 Similar associations have been reported in the USA21 22 and Great Britain.23 In Europe, eating at restaurants (which included eating at FFOs) was positively associated with BMI among men.24 On the other hand, in the USA eating at FFOs was positively associated cross-sectionally with BMI among low-income women.25

The limited body of research on environment influences on the chronic NCDs in Jamaica and the developing world, as well as the apparent lack of lifestyle changes despite health promotion programmes targeting individual-level prevention, suggests that barriers to these changes may yet be unrecognised and accounted for in the traditional modelling of risk factors. The studies on obesity previously referenced3 4 6 7 have only assessed geographical variations according to the dichotomised classification of urban/rural area of residence. Within the Caribbean and Latin American regions, there are a limited number of studies assessing geographical variations in obesity or other measures of adiposity such as mean BMI, using multilevel modelling (MLM) statistical techniques and/or geographic information systems (GIS) to determine whether there are associations with the built environment in a middle-income country (MIC) or small island-developing state (SIDS) context.26–28

The aim of this study was to provide a unique and important opportunity to address these gaps in understanding the retail food environmental mechanisms influencing mean BMI in Jamaica, a small island MIC. Our objective was to derive estimates of the associations between measures of the retail food environments and mean BMI, using a combination of MLM and GIS-based methods for contextualising the national survey data and calculating objective community exposures.

We hypothesised that: (A) there was variability in the mean BMI across Jamaican neighbourhoods, (B) the pathway between greater presence/closer proximity to supermarkets/markets, and lower mean BMI would be stronger for those of higher SES, (C) the pathway between greater presence/closer proximity to FFOs and greater mean BMI would be stronger for those of low SES and residing in urban areas, and (D) there would be sex differences in these associations.

Materials and methods

Study design and sample

The JHLS II was a cross-sectional, interviewer-administered, island-wide survey over a 4-month period between 2007 and 2008. The sample of 2848 participants aged 15–74 years represented approximately 70% of the predominantly (94%) Black Jamaican population.3 A stratified random sample of clusters known as enumeration districts (ED) was selected using a probability proportionate to the size of population of the parishes in the island in order to yield a nationally representative sample. Trained interviewers administered a structured questionnaire on diseases and lifestyle behaviours and performed anthropometry. Further details on the sampling technique are provided elsewhere.3

A total of 2529 (or 89% of) participants from the JHLS II data set were geocoded out of the original 2848 participants. Kreft and De Leeuw29 suggest a ‘30/30 rule’ so that researchers should strive for a sample of at least 30 groups with 30 individuals per group. For this study, each of the 101 EDs (sampling units) had an average of 28 individuals, providing sufficient power for the proposed secondary multilevel analyses. Furthermore, we also calculated power to detect a difference in BMI from 2 to 10 units based on differences in food environment exposure, at alpha=0.05 and power of 80%, with a design effect employed ranging from an intraclass correlation coefficient (ICC) of 2% to 10% and in all scenarios our sample and number of groups were sufficient to detect this difference in BMI.

Patient and public involvement

No patients were involved in the study. The study participants were community residents and were not involved in the design, recruitment or the conduct of the study. The study findings will be disseminated to the Ministry of Health, Jamaica and general public, including the study participants.

Measures

Individual-level measures

The primary outcome was mean BMI, calculated as weight divided by height squared (kg/m2). Weight was measured using calibrated electronic scales (Tanita models HD 314 or 2204) to 0.1 kg precision and height measured using a portable stadiometer (Seca) to 0.1 cm precision.

Additional covariates included age, sex, educational attainment, occupation, urbanicity and perceived community safety, and were examined as potential confounders; sex and urbanicity were also examined as effect modifiers, based on a priori theory. Named jobs were first categorised using the Jamaica Standardised Occupational Classification codes for 199130 which comprised 16 categories. These were collapsed into four groups: (A) highly skilled/professional (legislators, senior officials and managers, professionals/technicians and associate professionals, clerks, service workers and shop and market sales workers), (B) skilled (skilled agricultural and fishery workers, craft and related trade workers, plant and machine operators and assemblers), (C) unskilled (elementary occupations), and (D) unemployed/other (armed forces, retired, unemployed, housewife, self-employed, student and unclassified).

Perception of community safety was determined by asking each participant how safe he or she felt to walk in the community.

Household-level measures

The number of possessions owned (including but not exclusive to owning a radio, telephone, refrigerator, television, computer or car ownership) was used as a proxy for SES3 and classification based on the following tertiles: first tertile=lower class=≤6 items, second tertile=middle class=7–9 items, third tertile=upper class=10–16 items. The tertile categorisation was based on the distribution of ownership of these items. SES, using this definition, was examined as a potential confounder or effect modifier, based on a priori theory.

Environment-level measures

Each observation was linked, through a geocoded residential address, to neighbourhood-level proximity and density measures for supermarkets, markets and FFOs. Neighbourhood was defined as the ED. The final choices of environment-level measures for investigation were based on a combination of previously derived GIS-based measures,31 32 documented associations seen with the outcome of interest31 33 and data availability.

The locations of supermarkets/markets (figure 1A) and FFOs (figure 1B) were identified from Mona GeoInformatics Institute’s (MonaGIS) proprietary JAMNAV database, and were collected in 2009. The proximities and densities of supermarkets/markets and FFOs were estimated using application of Spatial Analyst tool in ArcGIS to data from MonaGIS proprietary JAMNAV database. Proximity of supermarkets and markets, combined to represent good sources of fresh fruits and vegetables, was determined as the straight-line distance (km) from each geocoded address to the closest supermarket or market. Two density variables were created. The first was supermarkets/markets per km2, and the second density variable, the number of supermarkets/markets per 1000 persons in the corresponding ED according to the 2011 census from the Statistical Institute of Jamaica.34 Buffer zones were set at 1 km to reflect a short walking distance that could be completed in about 10–15 min.35

Figure 1

Spatial distribution of supermarkets/markets (A) and fast-food outlets in Jamaica (B).

FFOs were defined as places where high-caloric food could be obtained relatively quickly and excluded traditional cook shops, snack shops and sit-down restaurants. Proximity and density measures were created in a similar way as done for supermarkets/markets.

The above-mentioned retail food environments excluded informal food spaces (eg, street vendors).

Zero-inflated variables

The absence of the environmental-level measure based on the participant’s geocoded address was indicated by a large proportion of zero values for most density measures as shown in the online supplementary table. New indicator variables (dummy variables) were subsequently created and the specific dummy variable included in regression models alongside the original quantitative forms of the respective retail food environmental-level explanatory variable. These dummy variables are also referenced as the zero-inflated form of the density measures.

Missing variables

Addresses for 11% of the JHLS II study participants could not be geocoded and contributed to missing data in subsequent analyses. The age/sex population of the data set used for this JHLS II secondary analysis subsample was compared with that for the non-geocoded data and no key deviations were observed.

Geocoded and non-geocoded (missing) participants were compared with respect to age, sex, SES categories and the key outcome variable of mean BMI, using mixed effects models, regression models accounting for survey design and regression models that ignored survey design. No associations were detected between the geocoded status (present or missing) and these other variables. This data analysis was done only on geocoded data, on the assumption that these participants were representative of the target population.

Statistical analyses

A complex database was created that combined individual-level JHLS II data with contextual environment-level data.

Descriptive data analysis estimated sex-specific and total survey-weighted means, proportions and 95% CIs for the outcome, explanatory and confounding variables as well as age-adjusted mean BMI and prevalence of poor fruit and vegetable intake, a key cardiovascular disease (CVD) risk factor. Age adjustment used direct standardisation across the strata identified as 10 years age band with weights being survey-weighted population proportions of the respective 10 years age group, as estimated using the JHLS II data.

The adjusted Wald test and the Pearson’s χ2 test corrected for survey design were used to determine whether, respectively, the age-adjusted and unadjusted estimates differed with respect to sex.

ICCs from hierarchical models quantified the proportion of variation in mean BMI potentially explainable at the ED level.

To determine and account for the effect of clustering at the neighbourhood level, subsequent analyses used multilevel models based on EDs nested within parish and examined the stratum-specific estimates of the effect of the environment variables on mean BMI with and without adjustment for covariates. Strata were defined using the urbanicity, sex and SES variables and stratum-specific multilevel models estimated if terms for interaction between the environment and strata variables were statistically significant. The Akaike information criterion statistic was used to determine the final best models. Collinearity assessed using Goodman and Kruskal’s gamma coefficient γ36 was the basis for selection of model covariates. To assess the chance of false positive errors, p values from these models were compared with the Bonferroni corrected significance level.

All analyses were conducted using STATA V.12 and V.14 (StataCorp, College Station, Texas).

Results

Sample characteristics

The weighted total and sex-specific summary statistics are shown in table 1. Women had higher mean BMI than men and higher proportions in the highly skilled/professional, unskilled and unemployed categories and in both first and second SES tertiles based on number of possessions owned (p<0.001). However, a greater proportion of men perceived their communities as unsafe (males=86.29%; 95% CI 82.64% to 89.94% vs females=80.38%; 95% CI 75.63% to 85.13%; p=0.012). There were no sex differences in urbanicity nor among those who had not completed high school.

Table 1

Total and sex-specific weighted sample characteristics (95% CI) for Jamaicans (JHLS II, 2008)

There was clustering in mean BMI across neighbourhoods in Jamaica, with an ICC of 4.16%, the proportion of the variance in mean BMI that can be accounted for by the neighbourhood level.29 37 No associations were found between the retail food environment variables and mean BMI in unadjusted regression models. There was also no effect modification by urbanicity.

Sex-specific regression models

There was interaction between sex and the following variables in their relationship with mean BMI: supermarket proximity (p=0.023), absence of supermarkets within a 1 km buffer zone (p=0.008) and FFO proximity (p=0.031). Figure 2 reveals that for women, in fully adjusted models, a 10 km increase in distance from supermarkets (or further proximity) was associated with a 1.20 kg/m2 higher mean BMI (95% CI 8.20×10–3, 0.24; p=0.036); the absence of supermarkets within a 1 km buffer zone was associated with a 1.36 kg/m2 higher mean BMI (95% CI 0.20 to 2.52; p=0.022). Proximity to FFOs was not associated with mean BMI in any sex.

Figure 2

Sex-specific unadjusted and adjusted β coefficients for the association of retail food environments with mean body mass index. ǂDummy variable for zero-inflated predictor. Model 1—unadjusted. Model 2—age adjusted. Model 3—adjusted for age and number of possessions. Model 4—adjusted for age, number of possessions, urban, occupation, education, perception of unsafe community. *p<0.05.

SES tertile-specific regression models

There was interaction between SES of a participant and a few retail food environment variables in their relationship with mean BMI. These included supermarket proximity (p=0.015), absence of supermarkets/1000 people/ED (p=0.033) and FFO proximity (p=0.045). Figure 3 reveals that a kilometre increase in the distance from a supermarket was consistently associated with higher mean BMI for all models for persons within the middle class, with a 0.17 kg/m2 higher mean BMI (95% CI 0.03 to 0.32; p=0.020) in the final model. Among persons in the upper class, the absence of supermarkets/1000 people/ED was associated with a 2.00 kg/m2 higher mean BMI (95% CI 0.08 to 3.92; p=0.041) only in age-adjusted models. Proximity to FFOs was not associated with mean BMI in any of the SES classes

Figure 3

Socioeconomic status (SES)-specific unadjusted and adjusted β coefficients for the association of retail food environments with mean body mass index. ǂDummy variable for zero-inflated predictor. Model 1—unadjusted. Model 2—age adjusted. Model 3—adjusted for age and sex. Model 4—adjusted for age, sex, urban, occupation, education, perception of unsafe community. *p<0.05; **p<0.01. ED, enumeration district.

Discussion

This study is the first to examine the impact of the retail food environment on obesity-related outcomes in a small Caribbean island. While we observed no significant associations between the retail food environment variables and mean BMI in unadjusted regression models, results revealed significant sex differences in the impact of the food environment, particularly for supermarkets. The further distance away from supermarkets and markets, and their absence within a 1 km buffer zone from residences, were associated with higher mean BMI in women; and further proximity to supermarkets associated with higher mean BMI for the middle class. There was no association with proximity to or density of FFOs, nor urban-rural differences.

We also observed clustering of BMI at the community level, with approximately 4.0% of the variance in mean BMI potentially explainable by environmental influences outside of the individual or at the neighbourhood level. This is similar to those reported by Harrington and Elliott38 in a Canadian sample and by Masood and Reidpath39 for many of the countries that participated in the 2003 World Health Survey. It is also similar to ICC by neighbourhood seen for obesity-related outcomes among adolescents in the USA as part of the National Health and Nutrition Examination Survey (NHANES).40

Our study findings corroborate previous research conducted in developed countries that indicate closer proximity to supermarkets/markets12 19 20 and increased density of supermarkets/markets9–11 16 41 42 are associated with less obesity-related outcomes in adults and children. Of note however, Wang and colleagues found higher neighbourhood density of small grocery stores associated with higher BMI among a sample of US women,15 while our study found higher BMI for women being associated with further distance from supermarkets/markets (inclusive of small grocery stores). There are potential explanations for this sex-specific difference as some studies have found that residential environments have a greater effect on women’s health.43 For example, women as primary food providers may depend more on neighbourhood sources for food than men do, which may in turn have an effect on their health outcomes, such as BMI. However, the influence of gender-specific roles in the provision of food supplies for families was not the focus of our study and so additional studies would be needed to better understand the differences observed. Although the geographic variables are somewhat dissimilar (a density vs proximity measure), the difference in direction of the associations raises the question as to whether presence of grocery stores, supermarkets and markets across different geographical contexts, including small island states like Jamaica, represents the same degree of accessibility, availability and ultimate consumption of healthy foods. For example, we found no urban-rural differences as hypothesised, perhaps due to physical activity levels (PAL), diet or some other unknown confounder masking the association. For example, it is quite possible that the presence or absence of supplementary food sources, used in rural communities and not captured in the environmental-level variables used in our study, for example, small produce plots or seasonal vegetables/fruits that supplement diets, may have played a role in the lack of urban-rural differences seen.

Living a greater distance away from supermarkets/markets was positively associated with mean BMI in the middle class in our study, however many other studies found this association in the lower class.44 We are unclear as to the reasons for the association observed among the middle-class participants. We surmise this could possibly be due to the influence of shopping preferences in terms of location. For example, a study conducted in low-income neighbourhoods in Philadelphia, USA, found that many residents did most of their shopping outside of the neighbourhood.45 Other possible influences could be cost, types of food purchased or lower PALs in this group, the latter perhaps being influenced by car ownership and/or greater use of motorised transport. Inagami et al46 found in a study conducted in the Los Angeles area of California, USA, that BMI was higher (A) where persons frequented grocery stores located in more disadvantaged neighbourhoods, which usually have higher availability of relatively inexpensive, energy-dense foods and (B) among those who owned cars and travelled farther to their grocery stores. However, other studies reveal inconsistencies, with mixed findings across various neighbourhood food retail contexts which may be due to the heterogeneity in defining neighbourhoods using SES47 48 and methodological limitations in measuring the inter-relatedness of neighbourhood residence, determinants of purchasing choices within and outside of the residential neighbourhood and issues such as dietary preferences.49 In our study, PAL, motorised transport and diet were not adjusted for in the regression models. These omissions may have suppressed associations.

It is also unclear why associations between FFO proximity and mean BMI were neither seen in the sex-specific nor SES-specific models. This may be an indication that demonstration of associations may depend on which obesogenic retail food environment variable or outcome measure of adiposity is chosen for the regression models. For example, Menke et al,50 in analysing the third US NHANES conducted, found that WC maintained a stronger association with CVD risk factors than the other measures of adiposity. Further work is needed on assessing the quality and utility of the measures used in this study as well as the development of new ones for the Jamaican and developing world context.

The major strength of this study is that it represents pioneering work in a SIDS and MIC context. We linked individually geocoded addresses from a nationally representative survey, with specific objective GIS-based retail food environment measures, and provided empirical evidence using MLM to explore the association between objective neighbourhood-level retail food environment measures and BMI. Additionally, the sample characteristics that differed significantly with sex were also included as covariates in the final models in order to minimise overestimation or underestimation of the true strength of the associations detected between the neighbourhood food retail environment measures and BMI.

Despite strengths, there are limitations that deserve mention, including the inability to make causal inferences given the cross-sectional design of the JHLS II and our definition of neighbourhood. Using EDs to represent a Jamaican neighbourhood could be deemed inadequate as they (1) are quite heterogeneous geographically in size, composition and context, and (2) may not fully represent exposure to the obesogenic environments being investigated. Additionally, although the outcome was objectively assessed, the risk factors included self-reported data on a single individual representing a household, which may introduce information bias. Furthermore, there was temporal mismatch of the data collected from individual JHLS II participants with that for the retail food environment-level variables, most of which were collected by MonaGIS in and after the year 2009, subsequent to the end of data collection for the JHLS II in 2008. This may have biased the results as individual exposures may have varied after the survey period, although food consumption behaviours are believed to be relatively stable over time. Lastly, the reliability and validity of the area-level environmental variables were untested locally and therefore uncertainty remains as to whether they were most effective in explaining any variance in the obesity-related outcomes. For example, a study conducted in a low-income neighbourhood in Spain found a mismatch between GIS-based measures of the food environment and resident’s perceptions of the environment.51 Another challenge faced in Jamaica with accurately characterising the neighbourhood retail food environment relates to the presence and patronage of the foods provided by many itinerant vendors which are often energy-dense foods. The use of mixed methods approaches in future research within SIDS similar to Jamaica may improve understanding on the associations with observed health-related behaviours such as food purchases (cost and locations) and dietary choices with BMI.

In conclusion, we found that further distance away from supermarkets and markets and their absence within a 1 km buffer zone from residences were associated with higher mean BMI, with important sex and social class differences. There has been an increase in the prevalence of obesity in Jamaica3 despite the implementation of policies and programmes to ameliorate its impact on the continuum towards NCDs.52 Higher mean BMI in Jamaicans may be partially explained by the presence of supermarkets and markets and differ by sex and social class. National efforts to curtail obesity in SIDS, like Jamaica, should consider the inclusion of interventions and future studies focused at the neighbourhood level that target the location and density of supermarkets and markets and those that consider sex and social class-specific factors that may be influencing the associations.

References

View Abstract

Footnotes

  • Contributors CACM and KPT conceived the study. LGG and PLA geocoded the data. CACM, NOY and KPT analysed and interpreted the data. RW supervised the research. CACM wrote the manuscript. PLA, KPT and RW edited the manuscript. All authors provided critical intellectual contributions and read and approved the final manuscript.

  • Funding This secondary analysis research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

  • Map disclaimer The depiction of boundaries on this map does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. This map is provided without any warranty of any kind, either express or implied.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Patient consent for publication Not required.

  • Ethics approval The University of the West Indies/University Hospital of the West Indies Ethics Research Committee (ECP 169, 14/15).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement All data relevant to the study are included in the article or uploaded as supplementary information. No additional data are available.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.