|Home | About | Journals | Submit | Contact Us | Français|
School environments may contribute to adolescent behavior, reproductive physiology and body composition (BC). The Nutritional and Health Survey in Taiwan (2010) for 1458 junior high school students was geo-mapped for 30 school environs. Facilities for physical activity (fitness centers, gymnasia and sports stadiums, activity centers and parks), sedentary activities (reading material rental shops (RMRS), internet cafes) and food and beverage outlets (FBOs) were calculated as weighted numbers within 1000m of schools. Multiple linear regressions were used to predict BC variable z-scores. For boys, higher fitness center densities and, for girls, gymnasia and sports stadiums were associated with less abdominal fatness. For girls, body mass index, waist circumference (WC) and triceps skinfold thickness (TSF) were greater when RMRS density was higher as was TSF with internet café density. Where there were no FBOs, boys’ WC and TSF were less with more parks, but girls were shorter and WC more adverse. With greater RMRS density and no FBOs, girls still had increased WC/Hip ratio, and less mid-arm muscle circumference. Boys’ findings were more evident after considering puberty. Physical activity facilities (differently by gender), food and beverage outlets absence for boys and low reading material rental shop density for girls increase the likelihood of healthier body composition.
The pubertal transition to adolescence and with continuing growth is one of considerable body compositional change and vulnerability. The changes may have life-long health consequences, albeit on a genetic and epigenetic background, whose trajectory may not be inevitable but, rather, modifiable during this time[2, 3]. Personal behaviors and the home, school and recreational environments are factors to be considered[4–11].
While the increase of childhood obesity is global, occurring even in the face of nutritional deprivation through poor dietary quality[12, 13], the overweight and obesity prevalence in some Asia Pacific regional countries are particularly high, in ascending order (for boys and girls): Japan (23 and 17%), China (24 and 16%), South Korea (25 and 20%), Chile (26 and 27%), Mexico (28 and 29%), the US (30 and 30%) and New Zealand (34 and 34%)[12–14]. In China, overweight (including obesity) prevalence increased over the past 20 years from 14% and 7.4% in 1985 to 34.2% and 30.3% in 2005, for 7-18y males and females, respectively. From 1991 to 2003, the prevalence of adolescent overweight and obesity in Taiwan went from 16.4% to 25.2% in boys and 8.9% to 15.2% in girls. A national cross sectional study in Taiwan between 2006 and 2007 found that the prevalence of overweight and obese adolescents was 19.2% and 10.3% for boys and 12.7% and 3.9% for girls, respectively. The Organization for Economic Co-operation and Development (OECD) report and other evidence from Australia indicates that food insecurity and social disparities with the global financial crisis contributed to increases in obesity prevalence within a short time. However, in France, where there has been active whole-of-community intervention, childhood obesity prevalence has plateaued. Thus, we need to look beyond the simple energy equation of input and output to understand the basis of body compositional disorders affecting fat and lean mass and how to minimize their occurrence. It is known, for example, that emotional health is associated with childhood obesity. For children, especially adolescents, various phenomena may be reflected in personal insecurity, family and personal relationships, educational and sporting achievements, recreational activities and diverse personal behaviors[20, 21]. In an appraisal of support to US food insecure families, girls were more responsive than boys to obesity prevention; gender differences in environmental susceptibility and prevention of body compositional disorders require further elucidation. The gender susceptibility may be more evident at puberty given that being overweight may lead to earlier puberty[23–25]. In Taiwan, the onset of puberty has ranged from 9–14 years of age for boys and 8–13 for girls.
It has become clear that, although urbanization does not necessarily pre-date obesity, the built environment may assume characteristics which make it more likely. The school environs are one such setting, negotiated by children in various states of food and other need at least twice a day, alone, with school-mates or with family.
Most attention has been focused on walkability to school and the food outlet exposures which are associated with body composition[7, 28–30]. We have been interested in how other facilities might interact or synergize with food outlets to affect body composition, with implications for later health patterns. In Taiwan, these may include fitness centers and gymnasia, internet and reading material rental shops (for books and magazines directed at adolescents) and public parks. It has been possible to map these across the country and locate them in relation to schools. We hypothesized that these might act co-operatively or detrimentally with food outlets for adolescent body composition. We further hypothesized that the advent of puberty would alter any such associations. To these ends, we have studied Junior High School students and their school environs across Taiwan.
Participants were adolescents aged 11–16 years (grade 7–9) who participated in the Nutrition and Health Survey in Taiwan (NAHSIT) 2010–2011. All of Taiwan’s 358 townships/districts were classified into 5 strata (northern 1, northern 2, central, southern, and eastern area) by geographical location and the population density. The survey used the PPS (probabilities proportional to sizes) sampling method to select 1620 students from 30 junior high schools (6 schools from each stratum) randomly. Their demographic and anthropometric information was obtained through interviews and health examinations. The present study excluded participants without body composition information, leaving 1458 junior high school students eligible for analysis. The study was approved by the Institutional Review Board of the National Health Research Institutes, Taiwan.
The NAHSIT 2010–2011 was funded by the Department of Health to provide assessment of health and nutrition of junior high school students in Taiwan. It was a cross-sectional and nationally representative survey. The questionnaire-based interviews and anthropometry measurements were conducted by trained personnel. The survey was conducted face-to-face at school with students and their teachers and through household interviews to obtain information regarding family members, socio-demographics, nutritional attitudes and behaviors, physical activity, and diet. Physical examinations were also conducted at school to provide body compositional information and allow venesection for fasting blood samples. All personal information was obtained through the survey with informed consent of parents or guardians.
The point of interest (POI) dataset for specific point locations was purchased from RITI Technology Inc., and corresponded to the NAHSIT 2010–2011 time frame. We gathered geographic points for food, physical or recreational and sedentary locations in the navigation system and added this layer onto the school environment data sets (Fig 1).
Junior high school addresses were obtained from the NAHSIT data and transferred to a geocoded database with Geo Gadget designed by the Center for Geographic Information Science (GIS), Academia Sinica, Taiwan. We applied Geo Gadget to transfer 30 schools’ geocoded nominal data into spatial coding with XY coordination for analysis. Obesogenic environments around the school were classified into 3 groups, (1) those to do with food (convenience, fast food, and beverage shops), (2) physical or recreational activity (fitness centers, gymnasia and sports stadiums, activity centers, and parks), and (3) sedentary activity (book and magazine rental shops and internet cafes). All locations of food store and activity facilities information were part of the POI datasets. Convenience stores included 7–11, Family Mart, OK, Hi-Life and Niko Mart. Fast food stores included McDonald’s, KFC and Burger King. Beverage stores were identified as providers of sweet beverages. Fitness centers were identified as a place of business with equipment and facilities for exercising and improving physical fitness. Gymnasia and sports stadiums were identified as places designed and equipped for indoor or outdoor sports, exercise, or physical education. Activity Centers were defined as a place that attracts people for dancing, working, studying, recreation or socializing. Parks can be used for pleasure and exercise. People can rent magazines, novels (e.g. classic, historical, mystery, romance and science-fiction), and comic books and so on in book and magazine rental shops. The internet cafe is a place which provides internet access to the public, usually for a fee, and these businesses usually provide snacks and drinks.
A circular buffer of 1000 meters was calculated around each school by ArcGIS (ArcInfo, version 10.2; ESRI Inc., Redlands, CA, USA), corresponding to an approximate 12–15 minutes’ walk distance at 4–5 km/hour[32–34]. Food store and activity facility availabilities were computed within this buffer area. For a measure of accessibility of food store and physical activity facilities, the weighted number of facilities was calculated and summed according to their distances from each school, with a weight of 1 for those with distances smaller than 200 m, 0.75 for those between 200–500 m, and 0.5 for those between 500 m and 1000 m.
Covariates considered in the present study included age, father’s ethnicity (Fukienese, Hakka, Mainlander, and Indigenous), mother’s education (up to primary school, secondary education, and above), household income (0–30,000, 30,000–50,000, 50,000–80,000, >80,000 NTD/month), ever smoking (no, yes), drinking alcohol (no, yes), dietary quality (the Youth Healthy Eating Index-Taiwan; YHEI-TW), physical activity (moderate or heavy physical activity ≥ 30 min/d), read on weekdays (0–1, 1–3, ≥3 hr/d), and development of puberty (boys: beard growth (not yet, at first, in progress, completed); girls: menarche (yes, no)).
Personal behaviors. Habits like reading, watching TV, playing computer games, were recorded.
Anthropometrics including height, weight, body mass index (BMI), waist circumference (WC) and triceps skinfold thickness (TSF) were collected during the NAHSIT physical examinations conducted at the schools. Muscle mass was determined by anthropometric measurement of the mid upper arm where circumference and TSF allowed its computation as Mid Arm Muscle Circumference (MAMC = MAC–(π*TSF)). We also calculated the ratio of WC/Height and WC/Hip. Adolescents were in the pubertal stage of physical development. Adolescents’ growth rates differ greatly in this age group, and therefore we calculated gender-age-specific z-scores for height, weight, BMI, WC, TSF, and MAMC. For purposes of classification of weight disorders, we used the Taiwanese Growth Charts (2002 year) where, for age and gender, underweight and overweight were BMI centiles <15% or ≥ 85%, respectively, and obesity was a BMI centile ≥ 95%. By these criteria, we defined “body compositional disorders (BCD)” such as underweight, overweight and overfatness, the latter taking into account BMI, WC, and TSF with lean mass indicator by MAMC. Evidence that a putative environmental risk factor is associated with BMI, WC and TSF, as the dependent variables, on multivariable regression analysis, was provided by the significance of the β coefficient (Fig 2). Blood pressure (systolic blood pressure (SBP) and diastolic blood pressure (DBP)) was taken.
Fasting blood was obtained by venesection for plasma metabolic analytes (glucose, total cholesterol, triglycerides, HDL, LDL, and uric acid).
Statistical analyses were performed using SAS software version 9.1 for Windows, and SUDAAN (Survey Data Analysis). SUDAAN was used to adjust for the design effects of cluster sampling and to obtain unbiased estimates of standard errors. The distributions of basic characteristics among adolescents according to gender were assessed by chi-square and t tests. Multiple linear regression analysis was used to compare schoolchildren’s body composition and the number of food outlets among regions. Models were adjusted for father's ethnicity, mother’s education, household income, smoking, drinking alcohol, dietary quality (YHEI-TW), physical activity (moderate or heavy physical activity ≥ 30 min/d), reading during weekdays and development of puberty (boys: beard growth; girls: menarche). Statistical significance was defined as P<0.05.
The characteristics of junior high school students by gender are shown in Tables Tables11 and and2.2. Their mean age was 13.5 years and ranged from 11 to 16 years. Socio-economic information about parents and students is provided. In particular, the prevalence of post-secondary education for both mothers and fathers was high at around or above 60%, which may represent a selective advantage for their children in high school attendance, since these education profiles exceed the national averages (40%) .
For potentially adverse personal behaviors, smoking was admitted by 10% of students (13% for boys and 6% for girls) and drinking alcohol by 12% (11% for boys and 14% for girls).
Time spent per day on physical and sedentary activities (reading, watching TV, and playing computer games) are also shown. Dietary patterns are described by way of the YHEI-TW score and in regard to consumption frequency (times/week) of convenience foods, sugary beverages and dairy products. YHEI-TW scores of 0 to 90 are regarded as acceptable, with higher scores indicating better dietary quality . It should be noted that the average use of fast foods in this age group across Taiwan is relatively low at less than once in 2 weeks. By contrast, the use of sugary beverages exceeds 5 times per week and that of liquid and flavored milk combined more than 4 times per week. This amounts to some 2 energy-replete beverages per school-day, if they were consumed on week-days.
In regard to pubertal development, we regarded the presence of beard hair as the onset in boys (32% with first indication, 33% in progress and 1% complete); and for girls, the menarche (94% had commenced menstruation).
Indices of body composition derived by anthropometry included BMI (kg/m2) (a mean of 21.1 for boys and 20.5 for girls); waist circumference (WC in cm) (a mean of 75.0 for boys and 73.4 for girls); mid arm muscle circumference (MAMC in cm, 21.4 for boys and 18.4 for girls); triceps skin fold thickness (TSF in cm) (14.2 for boys and 19.3 for girls);and WC/hip ratio (0.83 for boys and 0.81 for girls).
From BMI it can be seen that 20.1% of boys and 33.4% of girls were underweight and 30.2% of boys and 21.4% of girls overweight before the onset of puberty; after its onset, the respective prevalences, for underweight, were 13.9% for boys and 7.78% for girls and, for overweight, 32.4% for boys and 20.8% for girls.
Cardiovascular risk factors are shown by way of systolic and diastolic blood pressure (BP) (mean of 105/60) and metabolic variables (fasting plasma glucose (mg/dL) (mean of 95.4) and plasma lipids (mg/dL) (means for total cholesterol, triglycerides and HDL cholesterol were 158, 72 and 55 respectively).
Maps to show the representative junior high school environs facility densities studied are depicted by dot size in Fig 1. Gymnasia and sports stadiums, parks, book and magazine rental shops and fast food stores are shown in separate panels. It can be seen that the spectrum of densities differs markedly between major cities, provincial towns and coastal, rural or mountainous (Taiwan has a central mountain range) areas and also between the north-west (where Taipei and Taoyuan are situated) and the rest of the country.
There were no significant associations between food environments around junior high schools and student’s body compositions (Table 3).
The models for physical activity and sedentary activity facilities are shown in Tables Tables44 and and5.5. More favorable body composition is seen for boys with greater fitness center and also gymnasia and sports stadium density, but, for girls, only with gymnasia and sport stadium density. However, this applies to fat distribution as judged by waist circumference (WC) with a significantly negative β coefficient, and not necessarily fat mass with a non-significant β coefficient for BMI. Nevertheless, for girls, while WC/hip ratio is less with the greater availability of gymnasia, so also is muscle mass, judged by the MAMC z-score.
With more parks in the neighborhoods of the high schools, girls were more likely to have higher BMIs, even when extensive adjustments were made to the models. However, the significance of this association was lost with adjustment for the advent of puberty. This indicates that if park availability is a risk for overfatness in girls, that this might depend on the stage of physiological development.
The greater availability of book and magazine rental shops is associated with greater total, central and limb fat in girls, but not in boys. This applies in models which adjust for basic demography, personal behaviors, and reading itself, but not with adjustment for puberty. This suggests that puberty is a contributing factor to the risk for overfatness in this setting.
In view of a likely role of puberty in the associations between park and reading material rental shop availabilities, we modeled these in their own right, which is to say, there were 15 schools where no food outlet stores (convenience or fast food) or beverage shops were present in the school neighborhood, as shown in Table 6.
Here we find that, in this setting, boys are advantaged in their body composition by parks, but girls are not. But, girls’ height (indicative of growth velocity) is less. Likewise, boys’ body fatness and height are now susceptible to the reading material rental shops, dependent on puberty: and the susceptibility of girls is unchanged, irrespective of puberty.
In this study of the environs of junior high school students, we have found those facilities which may alter energy balance to act in different directions and that the net effect depends on whether they are present and their density. More particularly, ours and others’ previous findings of the risk that fast food outlets pose to body compositional disorders (BCD), may be countered by a sufficient density of physical activity facilities in the present age group. This was revealed by the modeling of BC with and without FBOs for parks whose association with less body fatness depended on the absence of FBOs, at least for boys.
Gender differences in the associations of facilities with BC are evident. Whereas boys appear to benefit in their BC by the availability of fitness centers as well as gymnasia and sports stadiums, girls only benefit from gymnasia and sports stadiums. For neither boys nor girls, do other activity centers, which comprise ping-pong, basketball, dancing, music, debating, or holding meetings, have any association with BC. The association with park density was unexpectedly adverse for girls’ BC and persistent through highly adjusted models, even when FBOs were excluded, and irrespective of puberty, although the fat distributions changed somewhat with these different models.
The adverse relationships for girls BC with both park density and RMRS may be related. It is known that much of the reading material rented by girls is in the form of fiction as novels and comic books, and of magazines, with an emphasis on story-telling and relationships and on personal care. Some of this is in the genre of romance comics[38, 39] as known in Western youth culture, but which, in Taiwan may take the form of “Yaoi” where adolescent girls create and read male homosexual comic books. Partly for these reasons, Taiwanese girls are the greater users of these shops than boys as well[39, 40]. It is possible and plausible that this reading material plays a role in forming body image and even BC. The counterpart is seen in aspiring sports people and in eating disorders with varying degrees of mis-match between body image and composition. Parks and the RMRS may serve a similar function if girls use them to read or share their experiences with others in these settings. In this event, unlike boys, they may not increase, but rather decrease their energy expenditure through park availability. This may be compounded if girls actually avoid parks on account of personal safety. There is some evidence that personal security is a factor which contributes to BC.
On the other hand, we find that the beneficial association of park availability to boys is evident once FBOs are absent (see above), although the fat distribution is puberty dependent in this setting.
Following rapid linear growth in the first 2 years of life, the onset of the growth spurt coincides with the onset of puberty; the earlier the onset the less the magnitude of the growth spurt, although not correlated with final height. But since almost 20% of final height is contributed by the pubertal growth spurt, this is a period of relative vulnerability to factors affecting stature. We found that, for girls, the greater the availability of reading rental shops, the shorter they were, when adjusted for puberty, which suggests that non-pubertal factors (possibly what the shops provided) were affecting growth at this age. However, this was not evident when considered in those school environs with no FBOs. Perhaps the FBOs were permissive in the manifestation of the reading rental shop and height association. Similarly, for the girls in school neighborhoods without FBOs, park accessibility increased the likelihood of shortness and of abdominal fatness, irrespective of puberty. We have previously found that elementary school girls exposed to more fast food outlets are more likely to be taller and discussed how this may be relevant to the increased risk of breast cancer in pre-menopausal Taiwanese women on account of rapid growth. Thus, it does not follow that the taller the better even though stunting can occur as a result of food deprivation and recurrent infection and even be adaptive to these circumstances[44–47]. Factors which are themselves not health adverse, but limit linear growth, as might be the case with reading and sedentary use of parks for social activity would not be reason to view shortness as unhealthy.
In the case of boys, where there were no FBOs, reading rental shop availability was associated with greater TSF thickness and also shortness irrespective of puberty; but while more park accessibility was not associated with height, it was with less fatness, with some dependence on puberty for its abdominal distribution.
In summary, reading rental shops encourage shortness in girls and boys, dependent on whether there are FBOs, and, for girls, in relation to puberty. Access to parks, probably through encouragement of physical play, reduce the risk of obesity in boys, but not in girls. Parks also have predisposition towards shortness in girls without regard to puberty, but not in boys.
The evaluation of the relevance of school environs facilities for energy balance and BC requires their constellation to be integrated in some way and related to personal behaviors. In the present study, we have used mapping, sub-analysis and adjustment to address this need.
As far as body compositional health is concerned, only those facilities that support physical activity (fitness centers, gymnasia and sports stadiums), depending on gender, and parks (for boys where there are no FBOs), play a beneficial role. The more sedentary activities (accruing from reading material rental shops-RMRS and internet cafes) are associated with less favorable BC, for girls in all settings and, for boys undergoing puberty where there are no FBOs. Therefore, it would be expected that, once sedentary opportunities increase, the availability of physical activity facilities becomes more critical for BC health. While we did not observe an explicit association between FBOs and BC, as previously reported for elementary school environs in Taiwan, this was evident by the findings in their absence, which showed parks to be of BC benefit to boys if they were not present. Further sensitivity analysis did not allow a distinction to be made between the various FBOs (data not shown).
Also relevant to the potential role of FBOs in student BC is their dietary information. On average, across Taiwan, the students bought fast foods less than once in 2 weeks, whereas they consumed sugary beverages 5 times per week, milk almost 4 times per week and total dairy almost 6 times per week; adjustments for these or a score of dietary quality (YHEI-TW) did not explain any of the models developed to predict BC by setting. This is not to say that dietary factors are unimportant in these students BC. Given the 10-fold greater use of sugary beverages than fast foods, we would expect this to be a potential risk factor for body compositional disorders. This would be in accord with the current concern about excessive consumption of liquid energy, without the constraint of food structure on appetite, in the growing epidemic of child obesity.
The recognition of factors affecting BC during puberty is complicated by the need to distinguish physiological growth and development from its pathology and also to take account of secular trends in the onset of puberty towards younger ages[1, 24–26]. There are probably bidirectional effects of energy imbalance and earlier reproductive maturity, further compounded by associated social and behavioral immaturity. Nevertheless, this is a vulnerable period for the development of BC disorders and one which merits attention to these children in their school environs, as elsewhere. We have, therefore, taken relatively unequivocal markers for the onset of puberty, the first of 4 stages in boys based on facial hair change, and one in girls, the menarche. We have adjusted our models accordingly.
We have found that puberty alters the associations between the school environs and BC and it, in turn, may be altered as part of an overall ecological disorder. Girls but not boys were less often underweight after than before puberty. This may be corrective and attributable to a relative decrease in growth velocity, to a more appropriate body image at this stage of development, to change in personal behaviors to do with eating and physical activity or to other ecological factors. Inevitably, any ecological disorder will include ones to do with relationships and emotional development, themselves intimately linked to growth and development[19, 50, 51]. In contemporary medical practice, the underlying basis of evolving health patterns is increasingly understood to be ecological. This is to say that people are an integral part of their environment, especially through sensory, microbiological and endocrinological pathways. Each of these conduits is now known to be pathogenetically linked to growth, development and body compositional disorders. This has extensive implications for present and future health[25–27].
Facilities indicative of opportunities for physical and sedentary activities have been mapped for their associations with BC, but these will not be a comprehensive set of environmental determinants of BC. However, they do show that there are synergies which may result in more or less healthful BC. This underscores the need for town and community planning to be attentive to a wide range of environmental factors in the health of schoolchildren, especially at or about puberty.
Facility availability does not necessarily represent usage- we have self-reported physical activity for which we have made adjustments in the models and which does not account for BC patterns in the different settings. Nevertheless, the reports may have had to do with activity undertaken outside the study environs.
Again, the focus on BC should not be taken as an isolated measure of student health. Thus, even though a facility may not be conducive to a healthy BC, it may provide other health benefits. A pertinent example from our work is the apparently disadvantageous associations in girls between settings with more available parks and RMRS. If these are a source of relationship-building, this may contribute to overall health.
Mapping of facilities in school environs which have to do with energy balance and body composition is increasingly feasible. The most plausibly relevant are food and beverage outlets and ones which provide for physical activity, purposefully built like fitness centers and gymnasia or recreational spaces like parks and gardens. Some facilities will provide for potentially healthful activities, not necessarily reflected in body composition, like meeting places, internet cafes and libraries or other providers of reading material. Students may be affected by this environment during puberty in ways different to earlier life and later adolescence. We have shown for a representative sample of junior high school students and their school environs that the associations with sedentary activities and recreational spaces food and beverage outlets with body composition are dependent on the presence or absence of food and beverage outlets, but with beverage outlets having the greater potential risk. The accessibility of reading material providers is associated with a strikingly adverse body composition in girls with greater total and abdominal fat, irrespective of puberty and even the presence of food and beverage outlets. On the other hand, boys’ body composition is more dependent on whether there are food and beverage shops and on puberty than is the case for girls.
The findings indicate that healthy school environments require attention to a composite of facilities to which students may have access. Girls are at risk of body compositional disorders from sedentary activities in ways which may be socio-behaviorally driven and that could be mediated by alterations in energy throughput, but might also be other energy-regulatory or metabolic pathways; there are, for example, recognized cortico-hypothalamic- neuro-endocrine pathways linking behavior to intra-abdominal omental fat metabolism[53–55]. Where the age group is pubertal, the vulnerability of boys to the aggregate environment is of added concern. To speak of the aggregate environment may be unusual, but it is very much the nub of what we have found. Namely, that several heretofore unconnected facilities in the school environment may be co-operative in their net effect on child health. The documentation of these probable collective environmental hazards to the body compositional health of young adolescent students in conjunction with regular health checks should allow risk mitigation.
Physical activity facilities in school environments favor healthier adolescent body composition, but differently for boys and girls (Fig 2). They act conjointly with facilities that provide sedentary activities. Food and beverage outlet absence increases the likelihood of healthier body composition for boys. The net association of facilities for physical and sedentary activity around junior high schools on body composition is dependent on stage of puberty.
Prof Wen-Harn Pan is the coordinator of the Nutrition and Health Surveys in Taiwan (NAHSIT) on which the project was based.
The present study was supported by the internal grant of National Health Research Institutes in Taiwan (04A1-PHPP21-014 and 103-2410-H-400-001-MY2).
Our study used the NAHSIT (Nutrition and Health Survey in Taiwan) dataset, which is third party data. The NAHSIT data set is only available on application to the Taiwanese office of National Nutrition Survey for researchers who meet the criteria for access to confidential data. Data requests may be sent to the Nutrition Survey office (wt.ude.acinis.smbi@oamgnauk). The data set is available for applicant use only. The present authors were granted access to these data. The data used in this paper were collected for the "Nutrition and Health Survey in Taiwan on Junior High School Students (NAHSIT 2010-2011)" sponsored by the Food and Drug Administration, Department of Health, Executive Yuan (99TFDA-FS-408).