Dramatic increases in childhood obesity over the past decades suggest the predominance of influential behavioural factors over genetic factors. It is generally believed that today’s environment enforces an inactive lifestyle that is likely to contribute to a positive energy balance and childhood obesity(Reference Brettschneider and Naul1). Despite the important role of nutrition (energy input) and physical activity (PA; energy output) in the prevention of childhood obesity, there seems to be limited research addressing the relationship between multiple health behaviours or clusters based on this energy-balance approach(Reference Dodd2–Reference Prochaska, Spring and Nigg4). Previous literature in this area has mainly focused on describing bivariate associations between PA and sedentary behaviour(Reference Motl, McAuley and Birnbaum5), PA and food choices(Reference Sanchez, Norman and Sallis6, Reference Kremers, De Bruijn and Schaalma7) or sedentary behaviour and food choices(Reference Boynton-Jarrett, Thomas and Peterson8, Reference Vereecken, Todd and Roberts9). Although significant relationships may exist between two health behaviours, the reported correlations are predominantly weak and often too small to be of much clinical or practical significance. Nevertheless, the low correlations between health-related behaviours does not exclude the existence of clusters or groups of young people who show negative health profiles due to factors such as low PA levels, high televison (TV) watching and unhealthy food choices(Reference Driskell, Dyment and Mauriello10). In addition, previous research has shown that we cannot assume that all children who eat healthily are also sufficiently physically active and vice versa(Reference Sabbe, De Bourdeaudhuij and Legiest11). Consequently, a better understanding of clustering of multiple health behaviours in young people is needed to inform effective interventions(Reference Prochaska, Spring and Nigg4, Reference Pronk, Anderson and Crain12).
Cluster analysis is frequently used in health research to identify mutually exclusive, homogeneous groups of individuals sharing numerous characteristics. The rationale behind studying the clustering of health behaviours lies in the fact that it can identify high-risk groups and thus provide a first step towards effective primary prevention strategies(Reference Schuit, van Loon and Tijhuis13). Previous studies examining the clustering of multiple health-related behaviours mainly focused on patterns in sedentary and exercise behaviours(Reference Marshall, Biddle and Sallis14–Reference te Velde, De Bourdeaudhuij and Thorsdottir19). From an energy regulation perspective, however, it is clear that we need to focus on both energy input and energy output(Reference Dodd2). Only a few prior studies have addressed the clustering of PA and food habits in youth(Reference Sabbe, De Bourdeaudhuij and Legiest11, Reference Kremers, De Bruijn and Schaalma20), but they failed to incorporate sedentary behaviours. Therefore, the first aim of the present study was to explore the existence of clusters in multiple lifestyle behaviours, including PA, sedentary pursuits and food choices, in a sample of seventh graders aged 11–12 years. In addition, we investigated whether the identified clusters could be characterised by gender, weight status and cardiorespiratory fitness (CRF). CRF is a key component of overall physical fitness and higher levels of CRF are associated with a healthier cardiovascular profile in children and adolescents(Reference Anderssen, Cooper and Riddoch21).
Experimental methods
Sampling procedure
All middle schools in the region of Leuven, a medium-sized city of Flanders (Belgium), were eligible for participation in the study. Six schools were randomly selected out of the forty-three eligible schools. Two of the six school principals approached by telephone declined participation because of lack of time. In each of the four participating schools, all of the seventh graders were invited to participate. An informed consent was obtained from the school board and all parents. Owing to illness, ten (3·1 %) of the 327 students were absent on the day of data collection. Unfortunately, no further information on these non-participants regarding their health-related behaviours was available. The final sample consisted of 317 students (137 boys and 180 girls) with a mean age of 11·7 (sd 0·5) years. Data were collected in October and November 2007 during a physical education lesson at school lasting 50 min. The study protocol was approved by the ethical committee of the Katholieke Universiteit Leuven.
Measures
Leisure-time physical activity
Despite many PA opportunities during school hours (e.g. school recess), numerous studies have shown that a substantial proportion of young people’s daily PA occurs outside the school hours(Reference Gidlow, Cochrane and Davey22, Reference Mota, Santos and Guerra23). In the present study, the pupils’ PA behaviour during leisure time was assessed making use of the Flemish Physical Activity Questionnaire, whose reliability and validity has been described elsewhere(Reference Philippaerts, Matton and Wijndaele24). Leisure-time physical activity (LTPA) was assessed by asking for children’s main sports practised in leisure time (compulsory school physical education, school sports and active transportation did not count for PA in leisure time). For each activity, they reported the frequency and the usual time they spend on that activity. For describing our sample, the level of LTPA was dichotomised and a low-activity profile was defined as 0–3 h/week involvement in LTPA, which is consistent with previous research(Reference te Velde, De Bourdeaudhuij and Thorsdottir19). This is different from the internationally well-accepted recommendation that children should engage in moderate-to-vigorous PA for at least 1 h/d(Reference Strong, Malina and Blimkie25), because in the present study we were not able to assess engagement in physical activities during school hours or active transportation. For the cluster analyses, the total time spent in LTPA (expressed as h/week) was used.
Sedentary behaviour
There are specific domains of sedentary behaviour that have both content and face validity. During after school hours, screen-based sedentary behaviours (TV, videos and computers) and school homework have been identified as the most prevalent sedentary activities(Reference Atkin, Gorely and Biddle26). Therefore, this part of the questionnaire was developed to measure a range of sedentary activities common among adolescents reported in the literature(Reference Atkin, Gorely and Biddle26). With respect to the screen-based sedentary behaviours, the participants were asked to report the time (h/d) they usually spent in TV viewing, computer use (excluding computer-based homework) and playing electronic computer games on a typical school day and at the weekend. For describing our sample, high TV viewing was defined as >2 h/d and for computer use (including electronic computer games) the cut-off was set at >1 h/d, which is consistent with cut-offs used in other study(Reference Boynton-Jarrett, Thomas and Peterson8, Reference Gorely, Marshall and Biddle18). For the cluster analyses, a composite measure of screen-based media use (expressed as h/week) was calculated by summing the time spent watching TV, using the computer and playing computer games. In addition, children were asked to report the weekly time (h/week) they spent doing study or homework, as an index of time spent on educational activities.
Dietary behaviour
A Dutch version of a validated FFQ(Reference Mast, Körtzinger and Müller27) was used to assess children’s eating practices. In the present study, five healthy food items (wholemeal bread, fruits, vegetables, fish and milk/diary products) and five risk-related (energy-dense) food items (white bread, meat, soft drinks, fast food and sweets/chips) were chosen to characterise dietary pattern. Food items were classified in categories as follows: 1 = never, 2 = 1–2 times/week, 3 = 3–5 times/week and 4 = daily. Although overestimation of the mean consumption frequency of food items must be considered using the FFQ, a satisfactory test–retest reliability of the FFQ items (intra-class correlation coefficient = 0·68) has been demonstrated elsewhere(Reference Landsberg, Plachta-Danielzik and Much28). On the basis of children’s responses, two dietary indices were calculated. A ‘healthy food index’ and a ‘risk-related food index’ were calculated by summing the frequency of the consumption of the five ‘healthy’ and the five ‘risk-related’ food items, respectively. Children indicating ‘never’ to eat these food products scored 5, while score 20 was allocated to children consuming all considered food items ‘daily’.
Overweight
Self-reported weight and height were used to calculate the BMI (kg/m2). The International Obesity Task Force criteria were used to define overweight and obesity(Reference Cole, Bellizzi and Flegal29). These criteria identify gender- and age-adjusted BMI values associated with a predicted BMI of 25 kg/m2 (overweight) and 30 kg/m2 (obesity) at the age of 18 years.
Cardiorespiratory fitness
CRF was assessed indirectly using the 20 m progressive shuttle-run test with 1 min stages and expressed as the running speed in km/h at the final completed stage(Reference Léger, Mercier and Gadoury30). All measurements were carried out under standardised conditions. All students were familiar with the test because this is a commonly used fitness tests in the physical education curriculum in Belgium.
Statistical analyses
Data management and computations of descriptive statistics and prevalence rates were performed using the Statistical Package for the Social Sciences statistical software package version 15·0 for Windows (SPSS Inc., Chicago, IL, USA) and the significance level was set at P < 0·05. Means, standard deviations and/or proportions were calculated for the key variables. Spearman rank correlations were calculated as a measure of association between the different behaviours. Gender differences in proportions of the behaviours were tested with χ 2 tests.
Consistent with previous research(Reference Sabbe, De Bourdeaudhuij and Legiest11, Reference te Velde, De Bourdeaudhuij and Thorsdottir19), κ-means cluster analysis was used to identify specific patterns in five health behaviour indices. κ-means clustering places children in mutually exclusive groups characterised by their PA behaviour, sedentary pursuits and food habits. Before clustering, all variables were transformed into standardised Z-scores. This prevents variables measured in larger ranges from contributing to the distance to a larger extent than variables with smaller ranges do. A possible instability of the cluster solution could be one of the cluster analysis’ drawbacks. Therefore, the reliability of the created cluster solution was examined on a random sample of 50 % of the total sample. In addition, a kappa degree of agreement in cluster membership was calculated by comparing memberships of the subsample with these of the total sample. Following these procedures, we specified three-, four- and five-cluster solutions. The four-cluster solution was found to be the most reliable and meaningful regarding the different patterns found and is reported in the Results section. The different clusters are given labels that characterise the behavioural pattern. These labels may reflect a value judgement, which was not intended as such, but was used as an efficient way to further describe our results regarding the clusters.
Cluster differences on gender distribution and pupils’ weight status were examined with χ 2 tests. By means of the ANOVA and a Bonferroni post hoc test, cluster differences on the selected behavioural indices and CRF levels were investigated. Partial η 2 were calculated as a measure of effect size. For overweight and CRF, these tests were controlled for gender.
Results
Pupil’s characteristics
Sample characteristics are presented in Table 1. Gender differences in proportions of health-related behaviours were observed for LTPA, computer use and some risk-related food habits. In particular, a significantly higher proportion of boys was found for high LTPA and high computer use compared to girls. Some gender differences were also observed regarding type of sports, with boys being more involved in endurance sports and soccer but less involved in dance. With respect to their food habits, boys reported more daily consumption of meat and soft drinks than girls. Approximately 9·5 % of our sample was categorised as overweight (including obesity), but no gender difference in overweight prevalence was observed. With respect to educational activities, girls spent significantly more time on school homework than boys (girls: 14·7 (sd 6·9) h/week; boys: 11·3 (sd 7·3) h/week; P < 0·01).
Gender-specific prevalence rates (%) of the self-reported health-related behaviours, sports type preference and weight status in 317 Flemish adolescents
Frequencies were significantly different between boys and girls: *P<0·01, **P < 0·05.
Description of the clusters
Spearman’s correlation coefficients between the five health behaviour indices were very low to moderate (r-values ranged from −0·433 to 0·319; see Table 2). On the basis of five health behaviour indices, the applied κ-means cluster analysis resulted in four reliable and meaningful lifestyle clusters. The stability of this cluster solution was guaranteed by repeating the procedure in a randomly selected subsample of the total sample and good agreement was found (Cohen’s κ = 0·79, P < 0·001). Distinguishing characteristics of each cluster are identified by high or low Z-scores. As shown in Fig. 1, cluster 1 (‘Sporty media-oriented mixed eaters’, n 50) shows high Z-scores (above 0) on all behaviours, except for time spent on educational activities (Z-score = −0·60). Cluster 2 (‘Academic healthy eaters’, n 88) shows an average level of LTPA (Z-score = 0·10), but has high values on homework and healthy food index combined with low Z-scores (below 0) for screen-based media use and consumption of risk-related food items. Cluster 3 (‘Inactive healthy eaters’, n 95) has low Z-scores on all behaviours, except for the consumption of healthy food items (Z-score = 0·50). Cluster 4 (‘Inactive media-oriented unhealthy eaters’, n 84) is characterised by high Z-scores on screen-based media use and risk-related food consumption but low Z-scores on LTPA and healthy food index. Mean absolute values of the behavioural factors by cluster solution are presented in Table 3.
Inter-correlations of indices indicating health-related lifestyle behaviors in 317 Flemish adolescents
Spearman’s ρ correlation coefficients were significantly different from 0: *P < 0·01.
Patterns of the final cluster centres (expressed as standardised Z-values) of the selected lifestyle behaviour factors, identified using κ-means cluster analysis; 
, Sporty media-oriented mixed eaters; 
, Academic healthy eaters; 
, Inactive healthy eaters; 
, Inactive media-oriented unhealthy eaters
Mean absolute scores of the specific health-related behaviours by lifestyle cluster in 317 Flemish adolescents
Mean values were significantly different between lifestyle cluster groups: **P < 0·01.
†κ-means cluster analysis.
‡Significant difference compared with ‘Sporty media-oriented mixed eaters’.
§Significant difference compared with ‘Academic healthy eaters’.
||Significant difference compared with ‘Inactive healthy eaters’.
¶Significant difference compared with ‘Inactive media-oriented unhealthy eaters’.
Association between lifestyle clusters and gender, weight status and cardiorespiratory fitness
Table 4 presents the association between cluster solution and gender distribution, overweight prevalence (including obesity) and CRF. For overweight and CRF, the results were controlled for gender. With respect to gender distribution, boys formed the majority of the first cluster, whereas the girls formed the majority of the other clusters. After controlling for gender distribution, the different clusters could not be characterised by pupils’ weight status, although the lowest prevalence rate for overweight (including obesity) was observed for the boys and girls in the ‘Sporty media-oriented mixed eaters’ group. With respect to CRF, a significant relationship between the shuttle-run test performance and lifestyle cluster membership was observed for boys only (P < 0·05; partial η 2 = 0·10). Post hoc tests revealed that the final speed of boys belonging to the ‘Sporty media-oriented mixed eaters’ group was approximately 1 km/h or two stages higher than their counterparts in clusters 3 and 4.
Association between lifestyle clusters and gender, overweight prevalence and cardiorespiratory fitness (CRF) in 317 Flemish adolescents
Frequencies or mean values were significantly different between lifestyle cluster groups: *P < 0·05.
†κ-means cluster analysis.
‡Significant difference compared with ‘Sporty media-oriented mixed eaters’.
§Significant difference compared with ‘Inactive healthy eaters’.
||Significant difference compared with ‘Inactive media-oriented unhealthy eaters’.
Discussion
The first objective of the present study was to explore the existence of clusters in multiple health behaviours including LTPA, sedentary behaviour and food habits in a sample of Flemish pupils aged 11–12 years. Despite the lack of strong correlations between the single behavioural indices, the applied κ-means cluster analysis showed that co-occurrence of specific PA, sedentary and dietary components may exist within certain groups of adolescents. For instance, the ‘Sporty media-oriented mixed eaters’ (comprising 15·8 % of our sample) represent a lifestyle that could be called a sporting lifestyle, although sports participation was not their only interest. Even though their LTPA index was very high, adolescents in this group also reported high amounts of screen-based media use. Comparable lifestyle profiles were also reported in other studies and were labelled as ‘techno-actives’(Reference Marshall, Biddle and Sallis14), ‘general activists’(Reference Telama, Nupponen and Pieron15) or ‘skaters and gamers’(Reference Nelson, Gordon-Larsen and Adair16, Reference Nelson and Gordon-Larsen17). However, the relationship between PA and screen-based sedentary behaviour is the subject of debate(Reference Biddle31). Whereas some have suggested that those who engage less in PA are more likely to watch a lot of TV(Reference Eisenmann, Bartee and Wang32), more recent studies have suggested a weak relationship between screen-based sedentary activities (i.e. TV viewing) and PA in youth(Reference Gorely, Marshall and Biddle18, Reference te Velde, De Bourdeaudhuij and Thorsdottir19, Reference Atkin, Gorely and Biddle26, Reference Marshall, Gorely and Biddle33). This suggests that there is time for both sedentary and active behaviours through the day and provides further support for the need to study more than one sedentary behaviour(Reference Biddle31). Apart from the screen-based sedentary behaviours, our adolescents in the ‘Academic healthy eaters’ group (comprising 27·8 % of our sample) were highly occupied with educational activities but also demonstrated an above average level of LTPA behaviour (Z-score above 0). This confirms previous research indicating that engagement in PA does not appear to displace the time spent for doing homework(Reference Atkin, Gorely and Biddle26).
From a public health perspective, the question appears to be whether or not a large volume of sedentary behaviour is compensated by a large enough volume of PA. It may therefore be more appropriate to speak of a sedentary behaviour compensation effect. Evidence for the compensation hypothesis has recently been found in a prospective cohort study among Flemish adolescents and showed that the relative risk of being overweight 2 years after the baseline measurement was considerably higher in those respondents who did not balance each 7 h of TV viewed per week by at least two and a half hours of PA per week(Reference Van den Bulck and Hofman34). Although high PA levels and high sedentary-related behaviours may coexist in some groups of adolescents at a particular stage in their life, longitudinal research has reported that PA declines throughout adolescence and is often replaced with less-active pursuits(Reference Nelson, Gordon-Larsen and Adair16). Increasing PA and decreasing sedentary behaviours are therefore key factors in promoting health in youth(Reference Nelson, Gordon-Larsen and Adair16). On the basis of aforementioned considerations, the adolescents in our ‘Inactive media-oriented unhealthy eaters’ group (comprising 26·5 % of our sample) may be at particular risk for developing chronic diseases later in life, not only because of their unbalanced screen time–LTPA ratio but also due to their low consumption of healthy food items in combination with a high consumption of energy-dense food items such as soft drinks and sweets/chips. Nonetheless, prospective cohort studies are needed to examine the stability of the identified cluster solution and to investigate the impact of this behavioural clustering on health outcomes later in life.
Our results pertaining to the second goal of the present study – the analysis of the association between lifestyle cluster membership and health-related physical fitness indices – revealed that the lifestyle clusters could not be characterised by the pupils’ weight status. Some other research has found an association between lifestyle clusters and overweight in European 9- to 14-year-old boys, with the most unhealthy behaviour pattern having the highest risk of being overweight(Reference te Velde, De Bourdeaudhuij and Thorsdottir19). Nevertheless, the detected association was less evident in girls and their cluster analysis did not include food habits. Our finding is however consistent with another study reporting no significant association between lifestyle clusters including PA and food habits and overweight prevalence rates among 10-year-old Flemish children(Reference Sabbe, De Bourdeaudhuij and Legiest11). During adolescence, the pupil’s body weight and body height is rapidly changing due to the adolescent growth spurt(Reference Rogol, Roemmich and Clark35). In fact, physiological demands of growth and maturation dictate that young people must necessarily be in a slight positive energy balance(Reference Dodd2). Consequently, this natural process of physical growth and developmental of (young) adolescents may have masked the relationship between lifestyle clusters and weight status in our sample. Nevertheless, it could be assumed that a continuous pattern of low PA, high media use and unhealthy eating might result in a higher overweight prevalence in cluster 4 (‘Inactive media-oriented unhealthy eaters’) in the long term. Longitudinal follow-up is needed to solve this question.
Our analysis regarding the association between cluster membership and CRF, on the other hand, revealed significant results, particularly among the boys. Boys in cluster 1 performed significantly better than boys in clusters 3 and 4. These latter groups should therefore be encouraged to increase their PA levels. Although physical fitness is in part genetically determined, it can also be greatly influenced by environmental factors. PA seems to be a key element in physical fitness enhancement(Reference Ortega, Ruiz and Castillo36). The present study clearly showed the well-known unequal PA pattern between boys and girls(Reference Sallis, Prochaska and Taylor37, Reference Van der Horst, Paw and Twisk38) and gender differences in cluster allocation mainly reflected this higher LTPA engagement in boys. This may be one of the reasons for the lack of association between lifestyle cluster membership and CRF among girls. Moreover, boys as a group performed significantly better on the 20 m shuttle run test than the girls in present study (boys: 12·0 (sd 1·2) km/h; girls: 10·9 (sd 0·9) km/h, P < 0·01). Our results also showed that boys practised different kinds of sports than girls (Table 1). On the basis of the compendium of energy expenditures for youth (expressed in metabolic equivalents or MET)(Reference Ridley, Ainsworth and Olds39), it was found that our boys participated to a greater extent in sports involving higher levels of energy expenditure such as soccer (on average 8·8 MET), bicycle touring (on average 6 MET) and running/athletics (on average 8·5 MET). This is an important finding, because previous research has observed that increased levels of vigorous PA (>6 MET), rather than light/moderate physical activity, are associated with higher CRF levels in children and adolescents(Reference Ruiz, Rizzo and Hurtig-Wennlof40). In order to enhance physical fitness levels, intervention programmes might need to focus on physical activities with appropriate intensity(Reference Ortega, Ruiz and Castillo36).
Limitations
The present study is subject to some limitations. First, the cross-sectional nature of the present study limits our ability to determine any causality in the results. Second, a relatively small sample was used. Although direct comparison between studies is hindered due to difference in year of data collection and/or recruited age groups, the frequencies of health-related behaviours and proportions of overweight in the current study were more or less in line with observations in other studies involving nationally representative samples(Reference te Velde, De Bourdeaudhuij and Thorsdottir19, Reference Janssen, Katzmarzyk and Boyce41). Despite the speculated resemblance of the current sample to nationwide representative cohorts, generalisation of the results might not be appropriate.
A third limitation of the study was related to the fact that only a few factors were available for characterising the clusters and future studies should investigate other factors in order to better characterise these subgroups. Factors such as socio-economic status (SES) or environmental characteristics were not assessed but may also be relevant correlates of adolescents’ lifestyles. Nonetheless, other research found no clear distinction between lifestyle clusters with respect to SES or environmental factors(Reference Sabbe, De Bourdeaudhuij and Legiest11, Reference Gorely, Marshall and Biddle18). This lack of association may have been caused by the different impact of SES and environmental factors on single health behaviours in youth(Reference Van Lenthe, Boreham and Twisk42). The association between lifestyle cluster allocation and SES or environmental factors may therefore disappear when multiple health behaviours are combined. Furthermore, it is important to keep in mind that school-related PA and active commuting to school were not assessed in the present study. Owing to the multifaceted character of PA, definitive statements regarding the association between our clusters and CRF cannot be made. With respect to active travel to school, for instance, previous research concluded that cycling to school was related to higher CRF levels in adolescents(Reference Cooper, Wedderkopp and Jago43, Reference Andersen, Lawlor and Cooper44). Cooper et al.(Reference Cooper, Wedderkopp and Wang45) also reported a positive association between active school travel and PA levels for adolescent girls. This may reflect the importance of PA accumulated while travelling to or from school and should therefore be taken into account in further research.
A final limitation of the present study is that prevalence data for overweight and obesity were derived from self-reported anthropometrics. Using self-reported height and weight instead of measured values is popular in many studies for practical or financial reasons. In the present study, all data (shuttle run test and completion of questionnaire) were collected during one physical education lesson that lasted 50 min. Owing to this limited time frame, it was not feasible to measure all children’s body dimensions. Nevertheless, using BMI values from self-reports may lead to an underestimation of the actual prevalence of overweight(Reference Elgar, Roberts and Tudor-Smith46) and should be taking into consideration while interpreting the results of the present study. However, correlation coefficients between self-reported and directly measured height and weight are generally high and self-reported anthropometrics are believed to be valuable if these are the only available source of data(Reference Sherry, Jefferds and Grummer-Strawn47). Previous research also showed that 96 % of a large sample was correctly classified as obese on the basis of self-reported measures(Reference Tsigilis48).
Conclusions
A strength of the present study was the application of a cluster analysis, identifying patterns of behaviours relating to both energy input and energy output. From the cluster analysis, we can conclude that risk-related behaviours such as high levels of screen-based media use or high consumption of energy dense food items can be components of more than one lifestyle. For instance, in one lifestyle (‘Sporty media-oriented mixed eaters’) the association between screen-based media use and LTPA was positive, but was negative in another lifestyle cluster (‘Inactive media-oriented unhealthy eaters’). Consequently, health-enhancing and risk-related behaviours may coexist in some groups of adolescents, which imply that health-related behaviours do not always discriminate in the same direction. Moreover, high levels of screen-based media use or high consumption of energy dense food items are not necessarily related with negative health outcomes, on the condition that these risk-related behaviours are compensated with healthy behaviours such as PA. This assumption can be supported by the significant higher CRF levels of the boys in the ‘Sporty media-oriented mixed eaters’ group, observed in the present study. Even though the practical implications are rather preliminary, the postulated compensation hypothesis may offer interesting pathways for prevention and intervention. However, this implies that children need selective advice on PA and eating habits in combination with clear reinforcements of healthy behaviours. As our clusters could not be identified by pupil’s weight status, screening of health behaviour is needed. If young people could be assessed in terms of the cluster they belong to, effective ways of directing them to relevant advice could be developed.
Acknowledgements
The present study received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. The authors have no conflict of interest to declare. J.S. developed and co-ordinated the overall study, carried out the statistical analysis and drafted the manuscript. C.R. supported in the statistical analysis and helped to draft the manuscript. The authors are grateful for all pupils, Physical Education teachers and school principals who collaborated in the study.
