Background: State-issued identification cards are a promising data source for neighborhood-level obesity estimates. Methods: We used information from three million Oregon state-issued identification cards to compute age-adjusted estimates of average adult body mass index (BMI) for each census tract in the state. We used multivariate linear regression to identify associations between weight status and population characteristics, food access, commuting behavior, and geography. Results: Together, home values, education, race, ethnicity, car commuting, and rural-urban commuting area (RUCA) explained 86% of the variation in BMI among tracts. BMI was lower in areas with higher home values and greater educational attainment, and higher in areas with more workers commuting by car. Discussion: Our findings are consistent with other research on socioeconomic disparities in obesity. This demonstrates state-issued identification cards are a promising data source for BMI surveillance and may offer new insight into the association between weight status and economic and environmental factors. Public health agencies should explore options for developing their own obesity estimates from identification card data.
Obesity prevention is a top public health priority; in 2010, one-third of adults and 17% of children and adolescents in the United States were overweight by an unhealthy amount [
Morris
Because DMV records are a novel data source for public health tracking, research is needed to determine whether estimates based on these data can be trusted. If BMI estimates from DMV records are reliably biased, patterns should match up with demographic and community factors known to be associated with obesity.
Like other public health threats, obesity prevalence is greater among groups with less education and lower incomes [
Weight loss or gain is largely a function of the foods people eat and the energy they expend, but genetic and environmental factors have a strong influence as well. Millions of Americans live in food deserts, where access to fresh, healthy, and affordable food is severely limited [
ESRI ArcGIS 10.0 was used to geocode DMV records and SPSS 19.0 was used for analysis. Data for this study came from 3 million Oregon driver licenses and identification cards issued to adults ages 18–84 years between 2005 and 2012. The Driver and Motor Vehicle Services Division (DMV) of the Oregon Department of Transportation provided data for this study. Obtaining the data was not difficult, because Oregon’s DMV recognized that the state public health agency had a legitimate claim to use the data for public health surveillance. However, since DMV did not maintain a codebook for their data set, it took several conversations with DMV staff to identify all the right variables to use for data cleaning and analysis.
Accounting for multiple cardholders per residence, and also for numerous cardholders at a given site per multiunit housing, about 1.7 million addresses were available for geocoding. Nearly 1.5 million addresses were geocoded; 90% of the addresses were successfully geocoded to a tax lot and 9% to a street. Based on the geocoding, we assigned each address to a census tract. Most of the remaining records were geocoded to city or postal code area centers and not used for this study.
Body mass index (BMI, in units of kg/m2) is the standard measure used to track the weight status of populations. General BMI classifications do not always accurately reflect an individual’s body composition (e.g., body builders have high BMI scores because of their muscle mass, but do not have excess body fat, like most people with a similar BMI). However, BMI remains the best measure for tracking population trends, as it correlates strongly with clinical assessments and health outcomes and is easily computed [
There are several different schemes for defining urban and rural areas throughout the United States [
Rural-urban commuting areas in Oregon.
Data on population characteristics came from the American Community Survey (ACS), an ongoing survey conducted by the Census Bureau since 2005. We used five-year ACS data (2006–2010) for the most stable and reliable Census tract level estimates. We analyzed ACS data on race, ethnicity, education, median home value, income, and commuting to work.
Estimates of the average intersection density (a simple measure of neighborhood walkability) were calculated from the 2010 Census Topologically Integrated Geographic Encoding and Referencing (TIGER) street network. Excluding intersections that led to cul-de-sacs, we counted the intersections within a one-square mile grid and used zonal statistics to create average intersection density per square mile for each census tract and block group. Intersection density is a widely used measure that is feasible to calculate for large areas [
Measures of proximity to various businesses were calculated using 2010 data from the Oregon Employment Department. Using North American Industry Classification System (NAICS) codes, we classified businesses as fast food restaurants, restaurants (including fast food), convenience stores, grocery stores or produce stands. We supplemented the produce stand list from a 2013 directory maintained by the Oregon Farmers’ Markets Association [
This study used an ecological, cross-sectional design. We computed summary statistics by RUCA type, weighting estimates by population counts where appropriate. We then fit linear regression models to the dependent variable of age-adjusted mean BMI at the tract level. Out of a total of 834 Census tracts, nine tracts with fewer than 100 people, according to the 2010 Census, were excluded from analysis, as were three additional tracts with missing home values in the ACS data. Variables were added to the model one at a time and retained if they improved model fit (as indicated by reduction in AIC scores or a statistically significant increase in adjusted R2) and beta coefficients were statistically significantly distinct from zero. We tested different classifications of the same constructs and kept the one that performed best. For example, models that described educational attainment using multiple categories performed better than models using a single “years of education” variable. The model was developed for the dependent variable of age-adjusted mean BMI for all adults, then for men and women separately. The tract-level models for men and women used sex-specific educational attainment rates.
Sixty percent of Oregon’s population resides in metropolitan core areas. These areas have the highest educational attainment rates and the largest non-white populations. Metropolitan core areas contain both the tracts with the highest median incomes and the lowest median incomes in the state. They also contain the tracts with the lowest and highest median age. Home values are lower and vary less across micropolitan and metropolitan commuting areas. Statewide, average age-adjusted BMI was 25.7 kg/m2 for women and 27.3 kg/m2 for men. Not controlling for other factors, BMI was highest in micropolitan and small town areas and lowest in metropolitan core areas (
Metropolitan core areas and rural areas had the lowest percent of workers commuting by car. In rural areas 13.7% of people worked at home, while in metropolitan core areas 8.9% commuted by bicycle or public transit. Walking to work was most common in small town and rural areas in Oregon (5.8%); micropolitan and metropolitan commuting areas had the smallest percent of workers commuting on foot (3.5%).
People in metropolitan areas in Oregon are less likely to live in a USDA-designated food desert. About 85% of households in metropolitan core areas are within a mile of a restaurant and most are within a mile of a supermarket. Only about one in five houses in metropolitan commuting and rural areas are within one mile of a supermarket, but about twice as many are within a mile of a restaurant. People in rural and micropolitan areas were more likely to live in a food desert, but even so, a quarter of the population in rural areas lives within one mile of a supermarket. Intersection density, a measure of walkability, was twice as high in metropolitan core areas than in micropolitan areas.
The statewide regression models explained most of the observed variation in BMI (R2 = 0.85 for women, R2 = 0.81 for men) (
Tract level statistics by RUCA category.
| Census Tract Characteristic | Metropolitan Area Core (N = 508) | Metropolitan Area Commuting (N = 99) | Micropolitan Area (N = 138) | Small Town (N = 32) | Rural Area (N = 45) |
|---|---|---|---|---|---|
| 2011 Population | 2,455,876 | 419,537 | 632,057 | 143,854 | 150,661 |
| Education, % of adults with | |||||
| No high school diploma | 10.2 | 10.3 | 13.7 | 14.5 | 14.2 |
| High school diploma | 22.0 | 29.7 | 30.7 | 32.4 | 33.2 |
| Some college | 33.4 | 38.7 | 36.5 | 35.1 | 35.3 |
| 4-year college degree | 21.6 | 14.0 | 12.5 | 11.8 | 11.5 |
| Graduate degree | 12.8 | 7.3 | 6.7 | 6.1 | 5.8 |
| Race/ethnicity | |||||
| % Asian | 5.6 | 1.0 | 1.2 | 0.8 | 0.8 |
| % Black | 2.3 | 0.3 | 0.5 | 0.5 | 0.6 |
| % Native American | 1.3 | 1.0 | 2.0 | 1.9 | 4.3 |
| % Hispanic | 12.0 | 6.2 | 13.4 | 11.1 | 9.7 |
| % White | 82.8 | 93.1 | 88.2 | 90.7 | 89.0 |
| Median age range | 19–61 | 27–57 | 22–69 | 32–54 | 20–56 |
| Median home value range ($ thousands) | 17–815 | 162–546 | 89–419 | 115–422 | 97–419 |
| Median household income range ($ thousands) | 10–147 | 32–96 | 20–70 | 23–59 | 25–64 |
| Age-adjusted mean BMI, all adults (Standard deviation) | 26.1 (1.1) | 26.7 (0.7) | 27.1 (0.6) | 27.1 (0.6) | 26.8 (0.9) |
| Females only (Standard deviation) | 25.2 (1.4) | 25.7 (0.8) | 26.4 (0.8) | 26.4 (0.8) | 25.9 (1.0) |
| Males only (Standard deviation) | 27.0 (0.9) | 27.6 (0.5) | 27.9 (0.6) | 27.8 (0.4) | 27.6 (0.8) |
| % population living in USDA defined food desert | 8.8 | 6.0 | 20.4 | 11.5 | 32.2 |
| Percent of households within 1 mile | |||||
| Convenience store | 65.3 | 24.1 | 41.2 | 39.6 | 12.8 |
| Fast food restaurant | 74.6 | 27.0 | 52.7 | 49.2 | 17.7 |
| Any restaurant | 82.0 | 33.7 | 62.7 | 62.8 | 34.3 |
| Supermarket | 56.0 | 22.4 | 43.3 | 46.9 | 25.7 |
| Produce stand | 23.9 | 12.7 | 21.0 | 27.8 | 11.5 |
| Workers commuting by car (%) | 80.3 | 87.7 | 88.3 | 86.5 | 77.6 |
| Intersections per square mile | 100.8 | 24.2 | 50.1 | 36.4 | 9.1 |
(
Linear regression results on age-adjusted mean BMI, statewide analysis by census tract (N = 822).
| Census Tract Characteristic | All Adults | Women | Men |
|---|---|---|---|
| Constant | |||
| Median Home Value (per $100,000) | |||
| Education | |||
| % with no high school diploma | |||
| % with high school diploma | |||
| % with some college | |||
| % with a bachelor or graduate degree (reference) | |||
| Race | |||
| % White (reference) | |||
| % African American | |||
| % Native American | |||
| % Asian or Pacific Islander | |||
| % Other | 0.006 | 0.000 | |
| % Hispanic ethnicity | 0.004 | 0.004 | 0.003 |
| % of Workers commuting by car | |||
| Rural urban commuting area | |||
| Metropolitan core (reference) | |||
| Metro commuting | 0.010 | ||
| Micropolitan | |||
| Small town | 0.000 | 0.038 | |
| Rural | 0.001 | ||
| Adj. R2 | 0.86 | 0.85 | 0.81 |
BMIs in micropolitan areas were statistically significantly higher than in metropolitan core areas, while rural areas had significantly lower mean BMIs than metropolitan areas for all adults and women. Metro commuting and small town tracts were not statistically significantly different from metropolitan cores. These results challenge previous research, which found higher obesity rates in rural areas, but are not necessarily inconsistent. Tract-level BMI data enabled us to use a more detailed urban/rural classification scheme than other studies, and as a result we may have revealed differences between metropolitan, micropolitan, small town and rural areas that were previously hidden. Other studies that identified disparities in health risks and outcomes between urban and rural areas used simpler classification schemes to compare urban and rural areas at the county level, thus the resolution of their analysis was not as precise as ours [
Race and ethnicity explained a small but statistically significant portion of the variance in BMI, with BMI consistently higher in areas with larger Native American populations and lower in areas with larger Asian populations. These finding are consistent with Oregon BRFSS data on obesity disparities by race. In contrast to the BRFSS, Hispanic ethnicity was not significantly associated with tract-level BMI estimates. It could be that the association was obscured after controlling for education and home values.
For every 10% of the population that commutes to work by car, mean BMI was higher by about 0.1 kg/m2. After controlling for other factors, intersection density and measures of proximity to businesses did not account for a meaningful amount of variation in BMI in our models. Other studies have found similar results [
Even though reporting bias is evident in the DMV data, the patterns observed in our statewide analysis were still consistent with other studies. We conclude DMV data are a powerful tool to guide public health obesity prevention work, and strongly encourage other public health practitioners in other states to collect and analyze DMV data. Not every state collects both height and weight information, but enough do to warrant building a nationwide dataset of tract-level BMI estimates from DMV data [
There are many factors related to a person’s weight status, but measures of affluence and education were the strongest predictors of average weight status at the population level. Almost all of the observed variation in average BMI for Census tracts in Oregon could be explained with a few variables relating to socioeconomic status, location of residence, and commuting behavior. Our study expanded on other research focused on metropolitan areas [
Population-level associations between socioeconomic conditions and BMI were stronger for women than for men. This result may reflect economic realities in America, where women and people of color earn lower wages for full-time work and are more likely than men to live in poverty [
We did not find meaningful associations between weight status and intersection density or average proximity to restaurants or groceries. This result is not a judgment on the efficacy of interventions that make it easier for people to walk or bike safely, or to access fresh fruits and vegetables It could very well be that the effects of these environmental factors are captured by the median home prices or RUCA variables. Also, associations may simply be obscured because our statewide measures did not describe neighborhood environments in enough detail. We did not do any ground truthing or other validity testing of the food environment measures. Examining the influence of community factors within a smaller study area, and measuring the environment in more detail, may produce different results. Based on our findings, additional research in micropolitan areas seems especially warranted.
Racial disparities are evident in our study but perhaps less pronounced than they would be elsewhere because Oregon’s population is predominantly white. Higher BMIs in areas with larger Native American populations highlight disparities on and around tribal lands. Replicating this study in states with more diverse populations may yield different results.
To our knowledge, this is the first statewide study on disparities in weight status to be conducted at a sub-county level. Our large sample of geocoded BMI data permitted a robust analysis across the state.
Data on height and weight were self-reported to the Oregon DMV and not verified by direct measurement. Though estimates do correlate well with those from the BRFSS, it is not possible to predict the error in any individual’s data [
One major limitation is that eight years of DMV records were aggregated for this study. Aggregating eight years of data provided stable estimates for less populated areas, but also prevented us from comparing time trends, which could provide stronger evidence of causal effects. Though many people likely moved residences during that time, we do not know how many updated their address information on file with the DMV. Information on both current and former addresses may explain more of the observed variation in BMI. If higher BMIs are associated with lower incomes, and lower income families are moving rapidly to areas with fewer amenities, the relationship between community factors and BMI will be obscured. With a single snapshot of the DMV database, we were not able to ascertain whether people were more likely to update their address or their weight. Future studies may resolve this issue by tracking people through multiple years of DMV data.
Measures of the food and physical environments are based on a snapshot at a single point in time, but we compared them to DMV records issued in multiple years. With access to more data on the changing environment, sequential years of DMV records could be used for time series analysis. As public health practitioners build those datasets with community assessments, opportunities for future research grow.
People without driver licenses or state-issued ID cards are not included in this study. This may influence results, especially in agricultural areas where many undocumented workers live. Legislative efforts to issue driver licenses to undocumented workers may enhance data completeness for future studies. However, since nearly every adult resident in Oregon has a state-issued ID card, we consider the data to be representative of the state resident population [
This study used a simple ecological, cross-sectional design and is therefore susceptible to the ecological fallacy. We therefore urge caution when attributing causal effects. Using BMI data from individual DMV records for a multi-level or longitudinal analysis may yield stronger results. Many factors known to be associated with obesity, such as physical activity and dietary patterns, were not included in this study because reliable population level estimates are unavailable for areas smaller than counties. As more local data become available, this study may be replicated elsewhere with different results.
State-issued identification cards are a promising new data source for BMI surveillance and may offer new insight into the association between environment and weight status. Public health agencies should explore options for developing their own obesity estimates from identification card data. DMV data can reveal previously-unseen variations in weight status between small geographic areas, challenging some assumptions about drivers of obesity. With access to DMV data, more robust studies can be done on the associations between weight status, population characteristics, and the environment.
This journal article was supported in part by Cooperative Agreement 2 U38 EH000950-04, funded by the Centers for Disease Control and Prevention. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the Centers for Disease Control and Prevention or the Department of Health and Human Services. For their thoughtful reviews of the manuscript, we thank David Hopkins, Sarah Bartelmann, Tia Henderson, Bruce Gutelius and Heather Gramp. We also thank the three anonymous reviewers who provided specific suggestions for improving the article.
Daniel Morris designed the study and was the primary author of the manuscript. Eric Main conducted geospatial data analysis. Jenine Harris designed the statistical analysis plan. Abraham Moland conducted research and helped write the manuscript. Curtis Cude helped design the study and write the manuscript. All authors contributed to the study and preparation of the manuscript.
The authors declare no conflict of interest.