Poverty, Educational Attainment and Health Among America’s Children: Current and Future Effects of Population Diversification and Associated Socioeconomic Change

A complex of interrelated factors including minority status, poverty, education, health status, and other factors determine the general welfare of children in America, particularly in heavily diverse states such as Texas. Although racial/ethnic status is clearly only a concomitant factor in that determination it is a factor for which future projections are available and for which the relationships with the other factors in the complex can be assessed. After examining the nature of the interrelationships between these factors we utilize direct standardization techniques to examine how the future diversification of the United States and Texas will affect the number of children in poverty, the educational status of the householders in households in which children in poverty live and the health status of children in 2040 assuming that the current relationships between minority status and these socioeconomic factors continue into the future. In the results of the analyses, data are compared with the total population of the United States and Texas in 2040 assumed in the first simulation scenario, to have the race/ethnicity characteristics of 2008 and in the second those projected for 2040 by the U.S. Census Bureau for the nation and by the Texas State Data Center for Texas in 2040. The results show that the diversification of the population could increase the number of children in poverty in the United States by nearly 1.8 million more than would occur with the lower levels of diversification evident in 2008. In addition, poverty would become increasingly concentrated among minority children with minority children accounting for 76.2 percent of all children in poverty by 2040 and with Hispanic children accounting for nearly half of the children in poverty by 2040. Results for educational attainment show an increasing concentration of minority children in households with householders with very low levels of education such that by 2040, 85.2 percent of the increase in the number of children in poverty would be in households with a householder with less than a high school level of education. Finally, the results related to several health status factors show that children in poverty will have a higher prevalence of nearly all health conditions. For example, the number of children with untreated dental conditions could increase to more than 4 million in the United States and to nearly 500,000 in Texas. The results clearly show that improving the welfare of children in America will require concerted efforts to change the poverty, educational, and health status characteristics associated with minority status and particularly Hispanic status. Failing to do so will lead to a future in which America’s children are increasingly impoverished, more poorly educated, and less healthy and which, as a result, is an America with a more tentative future.

Individual, Family, and Neighborhood Characteristics and Children's Food Insecurity

A growing body of work documents the influence of neighborhood environments on child health and well-being. Food insecurity is likely linked to neighborhood characteristics via mechanisms of social disadvantage, including access to and availability of healthy foods and the social cohesion of neighbors. In this paper, we utilize restricted, geo-coded data from the Early Childhood Longitudinal Study, which allows us to link individual children with their neighborhood's census characteristics, to assess how the neighborhoods of food secure and food insecure children differ at both the kindergarten level and in third grade. The average food insecure child lives in a neighborhood with a higher proportion of black and Hispanic residents, a higher proportion of residents living in poverty, and a higher proportion of foreign-born and linguistically isolated residents. After accounting for individual and household-level characteristics, children living in neighborhoods with a high proportion of Hispanic and foreign-born residents have a significantly increased risk of food insecurity compared to children living in neighborhoods which are predominantly white and have high socioeconomic status. We argue that interventions which take neighborhood context into account may be most efficacious for curbing child food insecurity.

An examination of the rate of high school completion experienced by handicapped and nonhandicapped students in the state of Texas

This research examines the graduation rate experienced by students receiving public education services in the state of Texas. Special attention is paid to that subgroup of Texas students who meet Texas Education Agency criteria for handicapped status. The study is guided by two research questions: What are the high school completion rates experienced by handicapped and nonhandicapped students attending Texas public schools? and What are the predictors of graduation for handicapped and nonhandicapped students?^ In addition, the following hypotheses are explored. Hypothesis 1: Handicapped students attending a Texas public school will experience a lower rate of high school completion than their nonhandicapped counterparts. Hypothesis 2: Handicapped and nonhandicapped students attending school in a Texas public school with a budget above the median budget for Texas public schools will experience a higher rate of high school completion than similar students in Texas public schools with a budget below the median budget. Hypothesis 3: Handicapped and nonhandicapped students attending school in large Texas urban areas will experience a lower rate of high school completion than similar students in Texas public schools in rural areas. Hypothesis 4: Handicapped and nonhandicapped students attending a Texas public school in a county which rates above the state median for food stamps and AFDC recipients will experience a lower rate of high school completion than students living in counties below the median.^ The study will employ extant data from the records of the Texas Education Agency for the 1988-1989 and the 1989-1990 school years, from the Texas Department of Health for the years of 1989 and 1990, and from the 1980 Census.^ The study reveals that nonhandicapped students are graduating with a two year average rate of.906, while handicapped students following an Individualized Educational Program (IEP) achieve a two year average rate of.532, and handicapped students following the regular academic program present a two year average graduation rate of only.371. The presence of other handicapped students, and the school district's average expense per student are found to contribute significantly to the completion rates of handicapped students. Size groupings are used to elucidate the various impacts of these variables on different school districts and different student groups.^ Conclusions and implications are offered regarding the need to reach national consensus on the definition and computation of high school completion for both handicapped and nonhandicapped students, and the need for improved statewide tracking of handicapped completion rates. ^

The location of majority and minority populations with respect to commercial treatment, storage, and disposal facilities in Texas

There is a growing interest in the location of Treatment, Storage, and Disposal (TSDF) sites in relation to minority communities. A number of studies have been completed, and the results of these studies have been varied. Some of the studies have shown a strong positive correlation between the location of TSDF sites and minority populations, while a few have shown no significance in that relationship. The major difference between these studies has been in the areal unit used.^ This study compared the minority populations of Texas census tracts and ZIP codes containing a TSDF using the associated county as the comparison population. The hypothesis of this study was that there was no difference between using census tracts and ZIP codes to analyze the relationship of minority populations and TSDF's. The census data used was from 1990, and the initial list of TSDF sites was supplied by the Texas Natural Resource Conservation Commission. The TSDF site locations were checked using graphical information systems (GIS) programs, in order to increase the accuracy of the identity of exposed ZIP codes and census tracts. The minority populations of the exposed areal units were compared using proportional differences, crosstables, maps, and logistic regression. The dependent variable used was the exposure status of the areal units under study, including counties, census tracts, and ZIP codes. The independent variables used included minority group proportion and grouping of the proportions, educational status, household income, and home value.^ In all cases, education was significant or near significant at the.05 level. Education rather than minority proportion was therefore the most significant predictor of the exposure status of a census tract or ZIP code. ^

An examination of the relationship between maternal and neighborhood-level characteristics and intrauterine growth retardation in Harris County, Texas, 1999--2001

The persistence of low birth weight and intrauterine growth retardation (IUGR) in the United States has puzzled researchers for decades. Much of the work that has been conducted on adverse birth outcomes has focused on low birth weight in general and not on IUGR. Studies that have examined IUGR specifically thus far have focused primarily on individual-level maternal risk factors. These risk factors have only been able to explain a small portion of the variance in IUGR. Therefore, recent work has begun to focus on community-level risk factors in addition to the individual-level maternal characteristics. This study uses Social Ecology to examine the relationship of individual and community-level risk factors and IUGR. Logistic regression was used to establish an individual-level model based on 155, 856 births recorded in Harris County, TX during 1999-2001. IUGR was characterized using a fetal growth ratio method with race/ethnic and sex specific mean birth weights calculated from national vital records. The spatial distributions of 114,460 birth records spatially located within the City of Houston were examined using choropleth, probability and density maps. Census tracts with higher than expected rates of IUGR and high levels of neighborhood disadvantage were highlighted. Neighborhood disadvantage was constructed using socioeconomic variables from the 2000 U.S. Census. Factor analysis was used to create a unified single measure. Lastly, a random coefficients model was used to examine the relationship between varying levels of community disadvantage, given the set of individual-level risk factors for 152,997 birth records spatially located within Harris County, TX. Neighborhood disadvantage was measured using three different indices adapted from previous work. The findings show that pregnancy-induced hypertension, previous preterm infant, tobacco use and insufficient weight gain have the highest association with IUGR. Neighborhood disadvantage only slightly further increases the risk of IUGR (OR 1.12 to 1.23). Although community level disadvantage only helped to explain a small proportion of the variance of IUGR, it did have a significant impact. This finding suggests that community level risk factors should be included in future work with IUGR and that more work needs to be conducted. ^

A survey of the implementation status of environmental management systems in U.S. colleges and universities, and, An environmental management system implementation model for U.S. colleges and universities

A census of 925 U.S. colleges and universities offering masters and doctorate degrees was conducted in order to study the number of elements of an environmental management system as defined by ISO 14001 possessed by small, medium and large institutions. A 30% response rate was received with 273 responses included in the final data analysis. Overall, the number of ISO 14001 elements implemented among the 273 institutions ranged from 0 to 16, with a median of 12. There was no significant association between the number of elements implemented among institutions and the size of the institution (p = 0.18; Kruskal-Wallis test) or among USEPA regions (p = 0.12; Kruskal-Wallis test). The proportion of U.S. colleges and universities that reported having implemented a structured, comprehensive environmental management system, defined by answering yes to all 16 elements, was 10% (95% C.I. 6.6%–14.1%); however 38% (95% C.I. 32.0%–43.8%) reported that they had implemented a structured, comprehensive environmental management system, while 30.0% (95% C.I. 24.7%–35.9%) are planning to implement a comprehensive environmental management system within the next five years. Stratified analyses were performed by institution size, Carnegie Classification and job title. ^ The Osnabruck model, and another under development by the South Carolina Sustainable Universities Initiative, are the only two environmental management system models that have been proposed specifically for colleges and universities, although several guides are now available. The Environmental Management System Implementation Model for U.S. Colleges and Universities developed is an adaptation of the ISO 14001 standard and USEPA recommendations and has been tailored to U.S. colleges and universities for use in streamlining the implementation process. In using this implementation model created for the U.S. research and academic setting, it is hoped that these highly specialized institutions will be provided with a clearer and more cost-effective path towards the implementation of an EMS and greater compliance with local, state and federal environmental legislation. ^

Bexar County's dental safety net for children: An estimate of our current capacity and need

Introduction. It has been well established that poor uninsured children lack access to dental care and have greater dental needs than their insured counterparts. ^ Objective. To assess the capacity of Bexar County's dental safety net to treat children. To assess the dental needs of Bexar County children ages 0-18 who are uninsured or are Medicaid or SCHIP recipients. ^ Methods. Information was requested from dental safety net clinics that treat children ages 0-18. Data from the census, NHANES and other sources was used to estimate the dental needs. ^ Results. The capacity of the current safety net to treat children is 33,537 patient encounters per year. The dental needs of the community are 227,124 patient encounters per year. ^ Conclusion. The results of the study indicate that Bexar County is not prepared to treat the dental needs of the underserved children in San Antonio.^

The local safety net, health insurance, preventable hospitalizations and costs: Harris County (Houston), Texas

The aim of this study was to examine the association between determinants of access to healthcare and preventable hospitalizations, based on Davidson et al.'s framework for evaluating the effects of individual and community determinants on access to healthcare. The study population consisted of the low income, non-elderly, hospitalized adults residing in Harris County, Texas in 2004. The objectives of this study were to examine the proportion of the variance in preventable hospitalizations at the ZIP-code level, to analyze the association between the proximity to the nearest safety net clinic and preventable hospitalizations, to examine how the safety net capacity relates to preventable hospitalizations, to compare the relative strength of the associations of health insurance and the proximity to the nearest safety net clinic with preventable hospitalizations, and to estimate and compare the costs of preventable hospitalizations in Harris County with the average cost in the literature. The data were collected from Texas Health Care Information Collection (2004), Census 2000, and Project Safety Net (2004). A total of 61,841 eligible individuals were included in the final data analysis. A random-intercept multi-level model was constructed with two different levels of data: the individual level and the ZIP-code level. The results of this study suggest that ZIP-code characteristics explain about two percent of the variance in preventable hospitalizations and safety net capacity was marginally significantly associated with preventable hospitalizations (p= 0.062). Proximity to the nearest safety net clinic was not related to preventable hospitalizations; however, health insurance was significantly associated with a decreased risk of preventable hospitalization. The average direct cost was $6,466 per preventable hospitalization, which is significantly different from reports in the literature. ^

A space-time analysis of reported shigellosis and salmonellosis cases in Texas from 2000 to 2004

Background. The purpose of this study was to describe the risk factors and demographics of persons with salmonellosis and shigellosis and to investigate both seasonal and spatial variations in the occurrence of these infections in Texas from 2000 to 2004, utilizing time series analyses and the geographic information system digital mapping methods. ^ Methods. Spatial Analysis: MapInfo software was used to map the distribution of age-adjusted rates of reported shigellosis and salmonellosis in Texas from 2000–2004 by zip codes. Census data on above or below poverty level, household income, highest level of educational attainment, race, ethnicity, and urban/rural community status was obtained from the 2000 Decennial Census for each zip code. The zip codes with the upper 10% and lower 10% were compared using t-tests and logistic regression to determine whether there were any potential risk factors. ^ Temporal analysis. Seasonal patterns in the prevalence of infections in Texas from 2000 to 2003 were determined by performing time-series analysis on the numbers of cases of salmonellosis and shigellosis. A linear regression was also performed to assess for trends in the incidence of each disease, along with auto-correlation and multi-component cosinor analysis. ^ Results. Spatial analysis: Analysis by general linear model showed a significant association between infection rates and age, with young children aged less than 5 and those aged 5–9 years having increased risk of infection for both disease conditions. The data demonstrated that those populations with high percentages of people who attained a higher than high school education were less likely to be represented in zip codes with high rates of shigellosis. However, for salmonellosis, logistic regression models indicated that when compared to populations with high percentages of non-high school graduates, having a high school diploma or equivalent increased the odds of having a high rate of infection. ^ Temporal analysis. For shigellosis, multi-component cosinor analyses were used to determine the approximated cosine curve which represented a statistically significant representation of the time series data for all age groups by sex. The shigellosis results show 2 peaks, with a major peak occurring in June and a secondary peak appearing around October. Salmonellosis results showed a single peak and trough in all age groups with the peak occurring in August and the trough occurring in February. ^ Conclusion. The results from this study can be used by public health agencies to determine the timing of public health awareness programs and interventions in order to prevent salmonellosis and shigellosis from occurring. Because young children depend on adults for their meals, it is important to increase the awareness of day-care workers and new parents about modes of transmission and hygienic methods of food preparation and storage. ^

The association of area socioeconomic status and cancer screening: A systematic review and multilevel study

Background. The gap between actual and ideal rates of routine cancer screening in the U.S., particularly for colorectal cancer screening (CRCS) (1;2), is responsible for an unnecessary burden of morbidity and mortality, particularly for disadvantaged groups. Knowledge about the effects of individual and area influences is being advanced by a growing body of research that has examined the association of area socioeconomic status (SES) and cancer screening after controlling for individual SES. The findings from this emerging and heterogeneous research in the cancer screening literature have been mixed. Moreover, multilevel studies in this area have not yet adequately explored the possibility of differential associations by population subgroup, despite some evidence suggesting gender-specific effects. ^ Objectives and methods. This dissertation reports on a systematic review of studies on the association of area SES and cancer screening and a multilevel study of the association between area SES and CRCS. The specific aims of the systematic review are to: (1) describe the study designs, constructs, methods, and measures; (2) describe the association of area SES and cancer screening; and (3) identify neglected areas of research. ^ The empiric study linked a pooled sample of respondents aged ≥50 years without a personal history of colorectal cancer from the 2003 and 2005 California Health Interview Surveys with a comprehensive set of census-tract level area SES measures from the 2000 U.S. Census. Two-level random intercept models were used to test 2 hypotheses: (1) area SES will be associated with adherence to two modalities of CRCS after controlling for individual SES; and (2) gender will moderate the relationship between area socioeconomic status and adherence to both modalities of CRCS. ^ Results. The systematic review identified 19 eligible studies that demonstrated variability in study designs, methods, constructs, and measures. The majority of tested associations were either not statistically significant or significant and in the positive direction, indicating that as area SES increased, the odds of CRCS increased. The multilevel study demonstrated that while multiple aspects of area SES were associated with CRCS after controlling for individual SES, associations differed by screening modality and in the case of endoscopy, they also differed by gender. ^ Conclusions. Conceptual and methodologic heterogeneity and weaknesses in the literature to date limit definitive conclusions about the underlying relationships between area SES and cancer screening. The multilevel study provided partial support for both hypotheses. Future research should continue to explore the role of gender as a moderating influence with the aim of identifying the mechanisms linking area SES and cancer prevention behaviors. ^

Hazardous air pollutants and childhood lymphohematopoietic cancer in Southeast Texas, 1995--2004

Southeast Texas, including Houston, has a large presence of industrial facilities and has been documented to have poorer air quality and significantly higher cancer rates than the remainder of Texas. Given citizens’ concerns in this 4th largest city in the U.S., Mayor Bill White recently partnered with the UT School of Public Health to determine methods to evaluate the health risks of hazardous air pollutants (HAPs). Sexton et al. (2007) published a report that strongly encouraged analytic studies linking these pollutants with health outcomes. In response, we set out to complete the following aims: 1. determine the optimal exposure assessment strategy to assess the association between childhood cancer rates and increased ambient levels of benzene and 1,3-butadiene (in an ecologic setting) and 2. evaluate whether census tracts with the highest levels of benzene or 1,3-butadiene have higher incidence of childhood lymphohematopoietic cancer compared with census tracts with the lowest levels of benzene or 1,3-butadiene, using Poisson regression. The first aim was achieved by evaluating the usefulness of four data sources: geographic information systems (GIS) to identify proximity to point sources of industrial air pollution, industrial emission data from the U.S. EPA’s Toxic Release Inventory (TRI), routine monitoring data from the U.S. EPA Air Quality System (AQS) from 1999-2000 and modeled ambient air levels from the U.S. EPA’s 1999 National Air Toxic Assessment Project (NATA) ASPEN model. Further, once these four data sources were evaluated, we narrowed them down to two: the routine monitoring data from the AQS for the years 1998-2000 and the 1999 U.S. EPA NATA ASPEN modeled data. We applied kriging (spatial interpolation) methodology to the monitoring data and compared the kriged values to the ASPEN modeled data. Our results indicated poor agreement between the two methods. Relative to the U.S. EPA ASPEN modeled estimates, relying on kriging to classify census tracts into exposure groups would have caused a great deal of misclassification. To address the second aim, we additionally obtained childhood lymphohematopoietic cancer data for 1995-2004 from the Texas Cancer Registry. The U.S. EPA ASPEN modeled data were used to estimate ambient levels of benzene and 1,3-butadiene in separate Poisson regression analyses. All data were analyzed at the census tract level. We found that census tracts with the highest benzene levels had elevated rates of all leukemia (rate ratio (RR) = 1.37; 95% confidence interval (CI), 1.05-1.78). Among census tracts with the highest 1,3-butadiene levels, we observed RRs of 1.40 (95% CI, 1.07-1.81) for all leukemia. We detected no associations between benzene or 1,3-butadiene levels and childhood lymphoma incidence. This study is the first to examine this association in Harris and surrounding counties in Texas and is among the first to correlate monitored levels of HAPs with childhood lymphohematopoietic cancer incidence, evaluating several analytic methods in an effort to determine the most appropriate approach to test this association. Despite recognized weakness of ecologic analyses, our analysis suggests an association between childhood leukemia and hazardous air pollution.^

Using small area estimation and geographic information systems technology to target health services for the uninsured

The three articles that comprise this dissertation describe how small area estimation and geographic information systems (GIS) technologies can be integrated to provide useful information about the number of uninsured and where they are located. Comprehensive data about the numbers and characteristics of the uninsured are typically only available from surveys. Utilization and administrative data are poor proxies from which to develop this information. Those who cannot access services are unlikely to be fully captured, either by health care provider utilization data or by state and local administrative data. In the absence of direct measures, a well-developed estimation of the local uninsured count or rate can prove valuable when assessing the unmet health service needs of this population. However, the fact that these are “estimates” increases the chances that results will be rejected or, at best, treated with suspicion. The visual impact and spatial analysis capabilities afforded by geographic information systems (GIS) technology can strengthen the likelihood of acceptance of area estimates by those most likely to benefit from the information, including health planners and policy makers. ^ The first article describes how uninsured estimates are currently being performed in the Houston metropolitan region. It details the synthetic model used to calculate numbers and percentages of uninsured, and how the resulting estimates are integrated into a GIS. The second article compares the estimation method of the first article with one currently used by the Texas State Data Center to estimate numbers of uninsured for all Texas counties. Estimates are developed for census tracts in Harris County, using both models with the same data sets. The results are statistically compared. The third article describes a new, revised synthetic method that is being tested to provide uninsured estimates at sub-county levels for eight counties in the Houston metropolitan area. It is being designed to replicate the same categorical results provided by a current U.S. Census Bureau estimation method. The estimates calculated by this revised model are compared to the most recent U.S. Census Bureau estimates, using the same areas and population categories. ^

Spatial analysis of water body distance and West Nile Virus human cases in Houston, TX

West Nile Virus (WNV) is an arboviral disease that has affected hundreds of residents in Harris County, Texas since its introduction in 2002. Persistent infection, lingering sequelae and other long-term symptoms of patients reaffirm the need for prevention of this important vector-borne disease. This study aimed to determine if living within 400m of a water body increases one’s odds of infection with WNV. Additionally, we wanted to determine if one’s proximity to a particular water type or water body source increased one’s odds of infection with WNV.^ 145 cases’ addresses were abstracted from the initial interview and consent records from a cohort of patients (Epidemiology of Arboviral Encephalitis in Houston study, HSC-SPH-03-039). After applying inclusion criteria, 140 cases were identified for analysis. 140 controls were selected for analysis using a population proportionate to size model and US Census Bureau data. MapMarker USA v14 was used to geocode the cases’ addresses. Both cases’ and controls’ coordinates were uploaded onto a Harris County water shapefile in MapInfo Professional v9.5.1. Distance in meters to the closest water source, closest water source type, and closest water source name were recorded.^ Analysis of Variance (p=0.329, R2 = 0.0034) indicated no association between water body distance and risk of WNV disease. Living near a creek (x2 = 11.79, p < 0.001), or the combined group of creek and gully (x 2 = 14.02, p < 0.001) were found to be strongly associated with infection of WNV. Living near Cypress Creek and its feeders (x2 = 15.2, p < 0.001) was found to be strongly associated with WNV infection. We found that creek and gully habitats, particularly Cypress Creek, were preferential for the local disease transmitting Culex quinquefasciatus and reservoir avian population.^

The contribution of change in ethnic composition to prevalence of obesity and overweight

Objective. The prevalence of overweight and obesity differs substantially among children of different ethnic origin in the United States. The objective of this project is to estimate to what extent changes in ethnic composition since 1980 have contributed to the current general “obesity epidemic” in the childhood population of the United States.^ Methods. Populations by single year of age, 0 to 19, male and female, for Hispanics, non-Hispanic whites, and non-Hispanic blacks, from the US Census’ July estimates for 1985, 1990, 1995, 2000 and 2005 were taken and compared to the population and percentage of those groups from 1980. Age, sex, and ethnicity specific prevalence rates for overweight in 1980 were then applied to the populations by age for the specified year and differences in expected and actual overweight populations were assessed.^ Result. The results from this investigation provide estimates of the contribution that different ethnic groups have made to the overall prevalence of overweight and obesity in the childhood population of the United States. Assuming that the 1976-1980 prevalence rates had remained unchanged, and then comparing the population had there been no change in ethnic composition with the population given the actual change in ethnicity, the percentage increase was 1.06% in 1985, 1.72% in 1990, 2.57% in 1995, 3.95% in 2000, and 4.39% in 2005.^ Conclusion. The changes in ethnic composition of the population, independent of changes in ethnicity-specific prevalence, have contributed substantially to the current overall prevalence of obesity in the United States childhood population. There are a number of factors that may be responsible for the apparent susceptibility of Mexican-Americans and non-Hispanic blacks to overweight and obesity. Further research is needed on specific characteristics of those populations.^

The relationship between supermarket nutrient movements and the mortality rates from fat-associated and non-fat-associated cardiovascular disease: Development of an ecologic method for determining community food patterns

Supermarket nutrient movement, a community food consumption measure, aggregated 1,023 high-fat foods, representing 100% of visible fats and approximately 44% of hidden fats in the food supply (FAO, 1980). Fatty acid and cholesterol content of foods shipped from the warehouse to 47 supermarkets located in the Houston area were calculated over a 6 month period. These stores were located in census tracts with over 50% of a given ethnicity: Hispanic, black non-Hispanic, or white non-Hispanic. Categorizing the supermarket census tracts by predominant ethnicity, significant differences were found by ANOVA in the proportion of specific fatty acids and cholesterol content of the foods examined. Using ecological regression, ethnicity, income, and median age predicted supermarket lipid movements while residential stability did not. No associations were found between lipid movements and cardiovascular disease mortality, making further validation necessary for epidemiological application of this method. However, it has been shown to be a non-reactive and cost-effective method appropriate for tracking target foods in populations of groups, and for assessing the impact of mass media nutrition education, legislation, and fortification on community food and nutrient purchase patterns. ^