Next Article in Journal
A Bibliometric Analysis of Research Publications of the Bucharest University of Economic Studies in Time of Pandemics: Implications for Teachers’ Professional Publishing Activity
Next Article in Special Issue
Factors Affecting Mothers’ Adherence to Home Exercise Programs Designed for Their Children with Cerebral Palsy
Previous Article in Journal
Leveraging Citizen Science and Low-Cost Sensors to Characterize Air Pollution Exposure of Disadvantaged Communities in Southern California
Previous Article in Special Issue
The Impact of Stand-Biased Desks on Afterschool Physical Activity Behaviors of Elementary School Children
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 19
Issue 14
10.3390/ijerph19148778
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Physical Behaviour Intensity Spectrum and Body Mass Index in School-Aged Youth: A Compositional Analysis of Pooled Individual Participant Data

by
Stuart J. Fairclough
1,*,
Liezel Hurter
2,
Dorothea Dumuid
3,4,
Ales Gába
5,
Alex V. Rowlands
6,7,
Borja del Pozo Cruz
8,
Ashley Cox
1,
Matteo Crotti
2,
Lawrence Foweather
2,
Lee E. F. Graves
2,
Owen Jones
2,
Deborah A. McCann
2,
Robert J. Noonan
9,
Michael B. Owen
10,
James R. Rudd
11,
Sarah L. Taylor
2,
Richard Tyler
1 and
Lynne M. Boddy
2
1
Movement Behaviours, Nutrition, Health & Wellbeing Research Group, Department Sport & Physical Activity, Edge Hill University, St Helens Road, Ormskirk L39 4QP, UK
2
Research Institute of Sport and Exercise Science, Liverpool John Moores University, Liverpool L3 3AF, UK
3
Alliance for Research in Exercise, Nutrition and Activity, Allied Health and Human Performance, University of South Australia, Adelaide, SA 5001, Australia
4
Centre for Adolescent Health, Murdoch Children’s Research Institute, Melbourne, VIC 3052, Australia
5
Faculty of Physical Culture, Palacký University Olomouc, CZ 779 00 Olomouc, Czech Republic
6
Diabetes Research Centre, Leicester General Hospital, University of Leicester, Gwendolen Road, Leicester LE5 4PW, UK
7
National Institute for Health Research (NIHR) Leicester Biomedical Research Centre (BRC), University Hospitals of Leicester NHS Trust and the University of Leicester, Leicester LE5 4PW, UK
8
Department of Sports Science and Clinical Biomechanics, University of Southern Denmark, 5230 Odense, Denmark
9
Department of Psychology, University of Liverpool, Liverpool L69 7ZA, UK
10
Department of Applied Health and Social Care and Social Work, Faculty of Health, Social Care and Medicine, Edge Hill University, Ormskirk L39 4QP, UK
11
Department of Teacher Education and Outdoor Studies, Norwegian School of Sport Sciences, 0863 Oslo, Norway
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(14), 8778; https://doi.org/10.3390/ijerph19148778
Submission received: 14 June 2022 / Revised: 4 July 2022 / Accepted: 15 July 2022 / Published: 19 July 2022
(This article belongs to the Special Issue Physical Activity and Health in Children)
Browse Figure
Review Reports Versions Notes

Abstract

:
We examined the compositional associations between the intensity spectrum derived from incremental acceleration intensity bands and the body mass index (BMI) z-score in youth, and investigated the estimated differences in BMI z-score following time reallocations between intensity bands. School-aged youth from 63 schools wore wrist accelerometers, and data of 1453 participants (57.5% girls) were analysed. Nine acceleration intensity bands (range: 0–50 mg to ≥700 mg) were used to generate time-use compositions. Multivariate regression assessed the associations between intensity band compositions and BMI z-scores. Compositional isotemporal substitution estimated the differences in BMI z-score following time reallocations between intensity bands. The ≥700 mg intensity bandwas strongly and inversely associated with BMI z-score (p < 0.001). The estimated differences in BMI z-score when 5 min were reallocated to and from the ≥700 mg band and reallocated equally among the remaining bands were −0.28 and 0.44, respectively (boys), and −0.39 and 1.06, respectively (girls). The time in the ≥700 mg intensity band was significantly associated with BMI z-score, irrespective of sex. When even modest durations of time in this band were reallocated, the asymmetrical estimated differences in BMI z-score were clinically meaningful. The findings highlight the utility of the full physical activity intensity spectrum over a priori-determined absolute intensity cut-point approaches.

1. Introduction

Overweight and obesity in children and adolescents (hereafter referred to as youth) continue to increase in prevalence [1] and are significant risk factors for health and wellbeing [2,3,4]. Consequently, there has long been a research focus on the relationships between youth physical behaviours and adiposity [5]. This research demonstrates a favourable relationship between physical activity and adiposity, with the most consistent associations observed for moderate-to-vigorous physical activity (MVPA) and vigorous physical activity (VPA) [6,7]. Conversely, unfavourable associations between adiposity and sedentary screen time (television viewing in particular) are commonly reported [8].
Accelerometers are used in health research to estimate the duration and mode of different physical behaviours [9]. The focus of epidemiological and intervention studies has typically been on the total volume or bouted time spent in specific physical behaviours and intensities. Traditional accelerometer data reduction methods rely heavily on cut-points whereby the same absolute intensity thresholds are applied to all participants’ data to generate estimates of time spent above these thresholds. These resultant free-living physical activity estimates are prone to intensity misclassification and bias because the cut-points used are specific to the original calibration protocols and sample populations [10]. In addition, the condensing of continuous data into only a small number of pre-specified intensity categories (e.g., sedentary time, light physical activity (LPA), MVPA) causes an important loss of information from the captured accelerometer data because a very limited portion of the data is used [11]. In recognition of these factors, emerging analytical approaches for assessing the associations between accelerometer-determined physical behaviours and health outcomes were recently advocated in the GRANADA consensus [9].
One such approach is to use a wide range of incremental acceleration intensity bands to describe daily physical behaviour patterns across the full intensity spectrum [9]. This gives high-resolution descriptions of the full physical activity intensity pattern and allows for an examination of the relationships between health outcomes and a wider range of intensities [11]. The small number of intensity spectrum studies conducted to date in school-aged youth have observed higher intensity activity to be most beneficially associated with health outcomes [11,12,13,14]. The inverse relationship between physical activity intensity and duration dictates that levels of higher intensity physical activity are typically low (e.g., 0.6% to 3.9% of waking hours [15,16]). Arguably, this is one reason why higher intensity physical activity has been less emphasised in physical activity promotion messaging, which focuses more on physical activity of at least a moderate intensity [17], which is considered more accessible and attainable for population health. However, even small amounts of higher intensity physical activity are beneficial for cardiorespiratory fitness in youth [18], and it has been shown in adults that similarly low doses may be favourably associated with other outcomes related to cardiometabolic health [19], such as adiposity [7]. Moreover, a range of structured and unstructured opportunities exist for youth to accumulate time in higher intensity physical activity, including sports participation, physical education classes, and time-efficient exercise modalities such as high intensity interval training (HIIT) [18].
Compositional data analysis with isotemporal substitution accommodates the complexity of analysing the intensity spectrum and allows the theoretical effects of displacing fixed durations of time between mutually exclusive incremental intensity bands to be investigated [9]. Almost all of the previous accelerometer-derived compositional analysis studies in school-aged youth have used three- or four-part compositions as the physical behaviour exposure and have relied on the traditional cut-points approach [20,21,22,23,24,25,26,27]. These studies described waking-hours physical behaviours as sedentary time, LPA, and MVPA and applied various accelerometer metrics and data reduction methods to produce time-use estimates of activity. An ActiGraph count cut-point for an MVPA of 2296 counts·min−1 was applied to describe the compositional associations between physical behaviours and adiposity across 24 h [25] and the school day [20]. Other studies used a much lower MVPA count cut-point of 1499 counts·min−1 [24,27], while raw acceleration cut-points (e.g., 200 mg; [21]) have also been applied to describe the compositional associations between physical behaviours and adiposity. Each of these studies reported significant associations between activity compositions and the adiposity outcomes, which included BMI z-score, percent and absolute body fat, and waist circumference. Moreover, a common finding in each study was that the greatest predicted differences in adiposity were when the time was reallocated to and from MVPA relative to the other intensities.
In contrast, only one dataset, using a sample of 10-year-old children with ActiGraph counts data, has examined the full intensity spectrum in school-aged youth [11,14]. Focusing on 5-to-16-year-olds, our novel study extends this work by employing compositional data analysis across the physical activity intensity spectrum, using raw acceleration data obtained from different accelerometer brands. The aims of this study were to (i) examine the compositional associations between the intensity spectrum derived from multiple incremental raw acceleration intensity bands and body mass index (BMI) z-score in youth and (ii) investigate the estimated differences in BMI z-score when durations of time were reallocated between incremental intensity spectrum bands.

2. Materials and Methods

2.1. Data Acquisition and Study Eligibility

Ethically approved wrist accelerometry studies led or supervised by the first or last authors were identified for inclusion within this pooled individual participant data analysis. Eligible studies involved school-aged youth who participated in observational or interventional physical activity research studies during school term time. Inclusion in the analysis studies required non-intervention assessments of wrist accelerometer-derived physical behaviours; thus, baseline data were used for contributing intervention studies. In addition, studies needed to provide anthropometric and demographic data including age, sex, and area-level socioeconomic status. Where published, details of these studies can be found elsewhere [22,28,29,30,31,32,33]. Investigators with a major involvement in the eligible studies (e.g., past PhD students, co-supervisors) were approached by email and invited to contribute individual participant data to allow for data harmonisation and subsequent pooled analysis. On receipt of the signed data transfer agreements, all the contributing investigators transferred their de-identified data via a secure file sharing system. Ethical approval for this pooled analysis study was granted by Edge Hill University’s Science Research Ethics Committee (#ETH2021-0034). The data were available from ten studies conducted in 63 schools between 2015 and 2019 in the Merseyside, Lancashire, and Greater Manchester regions of northwest England.

2.2. Outcomes

2.2.1. Outcome Variable and Covariates

BMI was calculated from stature and body mass measured using standard procedures, with participants wearing light clothing and no shoes [34]. BMI z-scores were assigned using British 1990 growth reference data [35], and International Obesity Task Force BMI cut-points were applied to classify participants by weight status [36]. Socioeconomic status (SES) was measured at the neighbourhood level using the English Indices of Multiple Deprivation (EIMD) [37,38] based on home postcodes. EIMD deciles were generated, where smaller values indicated a lower SES. Group (school) mean centering was applied to the participants’ ages to aid the model interpretation and to reduce the risk of multicollinearity.

2.2.2. Physical Behaviour Acceleration Exposure Variables

In the contributing studies, ActiGraph GT9X (ActiGraph, Pensacola, FL, USA; seven studies), GENEActiv Original (Activinsights, Cambs, UK; two studies), and Axivity AX3 (Axivity Ltd., Newcastle-Upon-Tyne, UK; one study) triaxial accelerometers were used. The devices have a dynamic range of ±8 g and were requested to be worn for up to seven consecutive days on the non-dominant wrist using either 24 h (seven studies) or waking hours wear protocols (two studies), with the sampling frequency set at 100 Hz (seven studies) or 30 Hz (two studies). The devices were initialised and data were downloaded using the latest releases of the respective ActiLife (versions 6.13.1 to 6.13.4), GENEActiv (versions 2.2 to 3.1), and OMGUI software (version 43) available at the time of data collection. Physical behaviour metrics were generated from the raw accelerometer data files (ActiGraph: gt3x followed by conversion to .csv format; GENEActiv: .bin format; Axivity: .wav format) and were processed in R using the package GGIR version 1.11-0 [39].
The data were harmonised by re-processing the raw accelerometer files using a standardised waking hours ‘day’ of 07:00 to 23:00 h (16 h or 960 min) to enable the inclusion of those studies that did not use a 24 h wear protocol. Signal processing included autocalibration, using local gravity as a reference [40], the detection of implausible values, and the detection of non-wear. Non-wear was imputed by default in GGIR, whereby invalid data were imputed by the average at similar time points on other days of the week [41]. The participants’ accelerometer data were excluded from the analyses if the post-calibration error was >10 milligravitational units (mg) and/or if <3 days of valid wear (i.e., ≥600 min·day−1) were available. We calculated the average magnitude of dynamic acceleration (i.e., average acceleration) as the Euclidean norm of the three accelerometer axes, with 1 g subtracted and negative values truncated to zero (ENMO) [42] averaged over 1-s epochs and expressed in mg. The average acceleration from the three devices worn on the non-dominant wrist has demonstrated equivalence in adults [43]. We generated nine acceleration bands using 50 mg increments for the first seven bands and wider bands for the remaining two (0–50 mg, 50–100 mg, 100–150 mg, 150–200 mg, 200–250 mg, 250–300 mg, 300–350 mg, 350–700 mg, and ≥700 mg). Based on the empirical evidence, the acceleration bands that are 350–700 mg and ≥700 mg reflect the upper range of MVPA [44] and higher intensity activities such as jogging and running [45], respectively. They were combined to reflect the anticipated very short duration (or absence) of accumulated accelerations within 50 mg bands above 350 mg [10]. The average minutes per day spent in each intensity band were calculated to create a nine-part composition per participant. The time in the nine bands was summed to 960 min.

2.3. Data Analysis

In these analyses, the outcome variable was the BMI z-score, and the exposure variables were the proportions of time spent in each of the nine acceleration intensity bands. Data were available for 1803 participants. Of these, 67 cases with missing BMI z-scores were removed because under a missing at a random assumption, there is no advantage in multiple imputation for missing data on the outcome variable [46]. Of the remaining 1736 participants, 233 did not achieve the accelerometer minimum wear criteria. These cases were also removed because the imputation of these summary activity estimates would rely on too many unknown assumptions about the pattern of missingness and would thus introduce random variation. No significant differences were observed between the included and excluded participants for centred age (p = 0.89), BMI z-score (p = 0.89), school type (p = 0.28), and sex (p = 0.15). A higher proportion of participants who did not achieve the accelerometer wear criteria were low SES (i.e., EIMD deciles 1–3, 66.5% vs. 55.1% who met the criteria), while fewer were high SES (EIMD deciles 7–10, 9.3% vs. 23.3%). Home postcodes were missing for 50 of the remaining 1503 participants (3.3% of the data), which prevented the calculation of EIMD deciles; these cases were also removed, leaving an analytical sample of N = 1453.
Compositional analyses were conducted using the R package compositions (v. 1.40-5, Mathsoft, Cambridge, MA, USA) [47]. Nine-part time-use compositions were expressed as nine specific sets of eight isometric log-ratios (ILRs) [48], which were used in multivariate linear regression models. BMI z-score was the dependent variable, and the intensity band composition ILRs were the explanatory variables. Each set of ILRs contained one ILR1 (i.e., the first pivot coordinate), which captured the time in one specific intensity band relative to all the remaining bands (i.e., the geometric mean of the remaining intensity bands), ensuring that each of the nine intensity bands were considered against all the remaining bands. The models were adjusted for sex, centred age, SES, accelerometer model, and accelerometer recording frequency. The influence of school on the models was trivial (ICC = 0.03), and the model parameter VIF values ranged from 1.01 to 1.96, indicating no multicollinearity. Where sex was significantly associated with BMI z-score, follow-up sex-stratified analyses were performed to assess the moderating influence of sex. The sex-specific models included the same covariates as the full sample models, with sex removed.
If the ANOVA table of the model fit showed that the set of intensity band ILRs was significantly associated with BMI z-score (p < 0.05), follow-up analyses were performed. In these, nine subsequent models were computed, examining the associations between BMI z-score and each of the ILR sets, including adjustment for the covariates. The co-efficient of the ILR1 was extracted from each of the nine models (and therefore each of the intensity bands) to examine the relationship of each intensity band (relative to all other intensity bands) with the BMI z-score [49]. Regression analyses were performed with the lmtest (version 0.9-40) [50] and car (version 3.0-12) [51] R packages, and model diagnostics were undertaken using the performance package (version 0.9.0) [52].
Where an intensity band (relative to the remaining bands) was significantly associated with BMI z-score (p < 0.05), this intensity band was then the focus of ‘one-for-remaining’ compositional isotemporal substitution analyses. These used predicted BMI z-scores for the initial ‘baseline’ average intensity band composition and compared them to the new BMI z-score predicted for subsequent compositions created using hypothetical time reallocations (1 to 20 min) between that intensity band and all the others. The estimated differences in BMI z-score for the reallocations of time to and from the selected average composition intensity band, and equally between the remaining intensity bands, were calculated by finding the difference between the two predicted BMI z-scores [53]. Adjustment for covariates was included in all the time reallocation analyses. Ninety-five percent confidence intervals (CIs) for the estimated differences in predicted BMI z-score were generated using the deltacomp R package (version 0.2.2, Mathsoft, Cambridge, MA, USA) [54]. The estimated difference in BMI z-score was considered significant when the 95% CI did not cover zero.

3. Results

3.1. Descriptive Results

Of the 1453 participants in the analytical sample, 57.5% were girls, 26.6% were overweight or obese, and two-thirds attended primary school (Table 1). Among boys, 71.3% lived in lower SES neighbourhoods (EIMD deciles 1–5) compared to 64.5% of girls. The participants were highly compliant in wearing the accelerometers, averaging 5.3 days of valid wear for 15.3 ± 1.0 h·day−1.
The geometric means of the intensity bands (linearly adjusted to collectively sum to 960 min) are presented in Table 2. Around 75% of waking hours were spent in the ‘inactive’ acceleration band of 0–50 mg. A linearly decreasing pattern of accumulated time in the intensity bands was observed, with fewer minutes spent in the higher intensity bands. The exceptions to this were the 350–700 mg and ≥700 mg bands, which were the sum of the combined 50 mg bands. The variation matrices of the time-use compositions representing the variability of the compositional dataset are included in Supplementary File S1 in Tables S1–S3.

3.2. Compositional Regression Analyses

BMI z-score was significantly associated with the intensity spectrum composition ILR1 coordinates (F(8,1430) = 11.9, p < 0.001) and sex (F(1,1430) = 10.6, p < 0.001) after the adjustment for covariates. Subsequent sex-stratified adjusted analyses revealed that BMI z-score was significantly associated with the intensity spectrum compositions for boys (F(8,1430) = 5.0, p < 0.001) and girls (F(8,1430) = 9.6, p < 0.001) (Supplementary File S2, Tables S4–S6). The unstandardised beta coefficients for each intensity band ILR1 determined which bands were the most dominant in the relationship with BMI z-score relative to time in all the remaining intensity bands. The ≥700 mg intensity band ILR1 was most strongly and inversely associated with BMI z-score (boys: βILR1 = −0.77, p < 0.001; girls: βILR1 = −0.71, p < 0.001) relative to the other ILRs in the composition (Table 3). Further, for girls, there were significant associations between BMI z-score and the 50–100 mg (βILR1 = 1.39, p = 0.002) and 100–150 mg (βILR1 = −2.55, p = 0.002) intensity band ILR1s relative to the remaining ones.

3.3. Compositional Isotemporal Substitution Analyses: One-to-Remaining Reallocations

The BMI z-scores for the baseline intensity band compositions (Table 2) were 0.72 units (boys) and 0.56 units (girls). Figure 1a shows the estimated differences in the boys’ BMI z-scores when 1 to 20 min were added to, and when 1 to 10 min (10 min was the maximum duration that could be reallocated from the ≥700 mg intensity band ILR1 because its geometric mean in the average composition was 11.1 min) were subtracted from—the baseline composition duration of the ≥700 mg intensity band (i.e., the most influential intensity band in Table 3) and redistributed equally among the remaining intensity bands (Supplementary File S3, Tables S7–S10). Reallocating time from the ≥700 mg intensity band was reflected by unfavourable estimated differences in BMI z-score that were greater than when the time was reallocated to this intensity band. For example, the estimated differences in BMI z-score when 10 min were reallocated to and from the ≥700 mg intensity band were −0.48 units (95% CI = −0.67, −0.29) and 1.71 units (95% CI = 1.03, 2.39), respectively.
Figure 1b–d present the estimated differences in the girls’ BMI z-scores following time reallocations involving the most influential intensity bands (50–100 mg, 100–150 mg, and ≥700 mg, respectively) and the remaining intensity bands. The reallocated time durations were 1 to 20 min, except for the ≥700 mg intensity band, where 1 to 5 min (5 min was the maximum duration that could be reallocated from the ≥700 mg intensity band ILR1 because its geometric mean in the average composition was 6.3 min) were subtracted from the baseline composition duration (Supplementary File S3). Adding 20 min to the 50–100 mg band (Figure 1b) reflected an estimated difference of 0.26 BMI z-score units (95% CI = 0.10, 0.32) compared to –0.30 (95% CI = −0.49, −0.12) when 20 min were substituted from this band. The estimated differences in BMI z-score when 20 min were reallocated to and from the 100–150 mg intensity band (Figure 1c) were –0.81 (95% CI = −1.34, −0.28) and 1.18 (95% CI = 0.41, 1.96), respectively. Reallocating 5 min to the ≥700 mg intensity band yielded a predicted difference in BMI z-score of −0.39 units (95% CI = −0.51, −0.28) compared to 1.06 (95% CI = 0.76, 1.36) when 5 min were substituted from the ≥700 mg band and reallocated equally among the remaining bands (Figure 1d).

4. Discussion

In this study we examined the associations between the physical behaviour intensity spectrum composition and BMI z-score across a wide age range of English youth. This is the first pooled individual participant data analysis in which harmonised raw acceleration data derived from different accelerometer models have been analysed in this way. For boys and girls, the highest intensity band (≥700 mg) was the most dominant in the relationship with BMI z-score (relative to the time spent in the remaining intensity bands). Respectively, time spent by girls in the 50–100 mg and 100–150 mg intensity bands was also positively and inversely associated with BMI z-score. Reallocating time between these dominant intensity bands and equally among the remaining bands resulted in predicted increases in BMI z-score when minutes were taken from the 100–150 mg and ≥700 mg intensity bands and decreases in BMI z-score when minutes were substituted from the 50–100 mg band. The predicted BMI z-score increases were larger than the predicted decreases when corresponding time reallocations were made to the 100–150 mg and ≥700 mg intensity bands.
Our main finding that the strongest associations with BMI z-score were from the highest intensity band concurs to some extent with the previous compositional analysis of the intensity spectrum in school-aged children, although it should be noted that this study did not focus on BMI z-score as the health outcome. Aadland et al.’s analysis in Norwegian 10-year-olds used hip-worn ActiGraph counts data and reported that the time in the 7500–7999 counts·min−1 intensity band was significantly and negatively associated with cardiometabolic risk [11]. The lack of direct comparability between proprietary ActiGraph hip counts data and ActiGraph, GENEActiv, and Axivity raw acceleration wrist data makes it difficult to translate Aadland et al.’s findings to ours. Although not specific to BMI z-score or other adiposity outcomes, these authors have also examined the accelerometer data intensity spectrum using multivariate pattern analysis [9]. They found that time spent in the 5000–7000 counts·min−1 intensity bands (ActiGraph hip counts data) was most strongly associated with children’s metabolic health [55,56]. Similar conclusions about the higher intensity bands (4000–5000 counts·min−1) were made in a more recent pooled analysis of data from over 11,000 children [12]. Although we also applied an intensity spectrum approach, it is difficult to make direct comparisons with the findings from these multivariate pattern analysis studies because of the different health outcomes studied and the aforementioned methodological and analytical differences between them. However, notwithstanding these differences, the common conclusion is that, when using the full range of available acceleration data to examine the relationships with indicators of adiposity or metabolic risk in youth (of which adiposity is a contributory factor), the magnitude of associations is stronger for time spent in the highest, rather than lower, intensity ranges.
In contrast to our analysis of the physical behaviour intensity spectrum, almost all previous compositional analyses of youth physical behaviours have used three or four component compositions defined by published cut-points. These have consistently reported MVPA to be most strongly and negatively associated with adiposity indicators relative to other physical behaviours [20,21,24,25,27,53]. However, when using MVPA as the highest intensity component in physical behaviour compositions, important parts of the accelerometer data are lost between the lower MPA threshold and the upper levels of VPA. This might lead to a loss of information that may be a hallmark of childhood obesity, making it unclear which specific intensities are most strongly related to adiposity. Some light has been shed on this question by compositional analyses of waking hours sedentary time, LPA, MPA, and VPA among Czech [57,58] and American youth [59], which observed VPA to be negatively and significantly associated with adiposity indicators relative to the other behaviours. These results support ours and demonstrate how the commonly reported finding that MVPA has the strongest influence on adiposity-related indicators relative to other physical behaviours actually masks the potentially more important contribution of higher intensity activities, while simultaneously amplifying the influence of lower intensity activities that fall within the MVPA range, as determined by the chosen absolute intensity cut-points. The use of absolute intensity cut-points derived from specific calibration study samples and protocols increases the likelihood of misclassifying time estimates in each physical behaviour [60], which can result in the over- and under-estimation of intensity-specific activity patterns [61]. By using multiple smaller incremental intensity bands rather than three-to-four broad cut-point categories, we allowed for a more nuanced description and analysis of participants’ activity patterns. Moreover, our analysis demonstrates the utility of compositional analysis for raw accelerometer data presented as an intensity spectrum, which further supports the integration of multiple physical behaviours and intensities to promote health and wellbeing among youth [59].
Among girls, BMI z-score was positively associated with time in the 50–100 mg intensity band, which may reflect a combination of low energy expenditure seated postures with arm movements and/or standing stationary postures [62]. In contrast, time spent in the hypothesised ‘sedentary’ 0–50 mg intensity band [62], relative to the other bands, was not significantly associated with BMI z-score. We used a standardised waking hours ‘day’ (i.e., 07.00 to 23.00); thus, it was not possible to differentiate between sedentary time and sleep within the 0–50 mg intensity band. As a consequence, some participants may have been asleep after 07:00 or before 23:00, and, as sleep is favourably associated with obesity risk [63], this may have confounded the strength of any positive associations between BMI z-score and sedentary time captured in the 0–50 mg band. It may also partly explain why accumulated very low intensity time in the 50–100 mg band had a significant and stronger positive association with girls’ BMI z-scores. Significant negative associations were observed between girls’ BMI z-scores and time in the 100–150 mg intensity band. Based on the available published raw acceleration data from wrist-worn devices, this lower intensity band may reflect slow-medium paced walking, akin to LPA [10]. Compositional analysis studies using count cut-points have reported both positive associations [24,27] and negative associations between LPA and BMI z-scores [57]. Moreover, an earlier non-compositional study observed negative relationships between LPA and children’s fat mass [64]. In most free-living studies, LPA represents the longest accumulated duration of waking hours physical activity and typically encompasses ActiGraph counts in the 100 to 2000+ counts·min−1 range. This likely includes some misclassification of sedentary time and MPA, which, in addition to differences in analytical approaches (i.e., non-compositional vs. compositional), provides insight into why the findings from previous cut-points studies lack agreement. Applying compositional analysis, which accounts for the mutually exclusive relationships between activity intensities, with higher-resolution intensity bands rather than a limited number of cut-points can advance this area of research by allowing for a more detailed examination of activity intensity associations with health outcomes.
The significant associations between lower intensity bands and BMI z-scores were only apparent in girls, who are typically less active and less fit than boys [21]. For our pooled sample, the 100–150 mg intensity band, in which girls accrued more time than boys, may have represented a higher relative intensity for some girls that was favourably associated with BMI z-score. This is consistent with previous research showing volume of physical activity to be more strongly associated with health outcomes in lower active and less fit groups, while among higher active and more fit groups, associations were strongest with increasing physical activity intensity [65]. Furthermore, the associations may have been influenced by the environmental contexts in which the physical behaviours occurred. For example, significant negative associations were recently reported between out-of-school LPA and adiposity in girls but not boys, while associations with school time LPA were positive and non-significant irrespective of sex [66].
We found that when hypothetical time reallocations were modelled between one intensity band and equally among the remaining bands, the most favourable predicted differences in boys’ and girls’ BMI z-scores were when time was added to the ≥700 mg intensity bands. Consistent with previous compositional analyses, the predicted changes in BMI z-score were asymmetrical [20,21,24,25,27,53]. This demonstrates that the potential detrimental health effects of taking time away from activity above the 700 mg intensity band and redistributing it equally to the remaining bands were greater than the beneficial effects of adding time at this higher intensity. The relative contributions of each intensity band to the waking hours day provide some insight into these asymmetrical relationships. Taking time away from the ≥700 mg intensity band, which contributed to 1.2% (boys) and 0.7% (girls) of the 16 h waking day, is a substantially larger relative change than taking time from the 50–100 mg intensity band, for example, which contributed to 10.7% (boys) and 11.2% (girls) of the day [21]. The observed predicted differences in BMI z-score were also not linearly related to the durations of reallocated time in the ≥700 mg and 100–150 mg bands (girls only). This has been consistently observed in previous compositional analysis studies [21,25,53] and reflects findings from experimental research whereby diminishing health benefits are predicted by marginal increases in physical activity [67]. Moreover, the wider 95% confidence intervals corresponding to the girls’ 100–150 mg time reallocations indicate greater variability and relatively less precision in the predicted BMI z-score differences compared to those from the ≥700 mg band reallocations. This complements the consistent relationship between the higher intensity activity and health outcomes reported in compositional analysis studies [20,21,24,25,27,53,57,59] and underscores the importance of promoting, providing, and not withholding developmentally appropriate and enjoyable opportunities for youth to be active at these intensities, even for short accumulated durations (e.g., through active play, sports, physical education, HIIT, etc.). In our sample, as little as 3 min (girls) and 5 min (boys) of additional time spent in the ≥700 mg intensity band relative to the remaining bands were associated with predicted decreases in BMI z-score of 0.26 and 0.28 units, respectively. These values are greater than the BMI z-score mean differences between the intervention and control groups for adolescents reported in a recent Cochrane review of obesity prevention interventions [5] and those observed in an effective school-based healthy weight intervention for children [68]. Hence, accumulating time spent being active at these higher intensities is advocated as being beneficial for health in youth.
This study employed a large pooled individual participant dataset comprising primary and secondary school youth. Assessments of the physical behaviours from different devices were used to generate raw acceleration data that enabled the analysis of the physical behaviour intensity spectrum rather than a small number of activity intensities derived from absolute intensity cut-points. Compositional data analysis with ‘one-for-remaining’ isotemporal substitutions enabled the associations with BMI z-score and predicted changes to be presented. There were, however, limitations that should be considered when interpreting the results. The contributing studies were cross-sectional, so causality between the intensity spectrum compositions and BMI z-score associations cannot be inferred, and the possibility of bi-directional associations is acknowledged. Further, the use of a cross-sectional analysis meant that the estimated differences in BMI z-score reflected more of a sample shift in intensity spectrum time allocations than actual differences for individual participants [69]. As the resolution of intensity spectrum bands can result in different interpretations of the relationship between activity intensity and health [14], our choice of nine intensity spectrum bands may have influenced the associations with BMI z-score. In future, using a higher resolution intensity spectrum with emerging analytical techniques such as functional data analysis may help improve our understanding of the relationships between specific acceleration ranges and health outcomes [9]. Although the analyses were adjusted for sociodemographic and methodological variables, it is possible that residual confounding from other non-measured variables such as sleep may have influenced the results. Lastly, the regional sample limits the generalisation of the findings to the wider UK and beyond.

5. Conclusions

This innovative compositional data analysis of the physical activity intensity spectrum in 5-to-16-year-olds is the first to use pooled raw acceleration data from different devices. Time in the highest intensity band (≥700 mg), relative to the remaining intensity bands, was significantly and negatively associated with BMI z-score, irrespective of sex. ‘One-for-remaining’ time reallocations involving the ≥700 mg intensity band indicated that the asymmetrical estimated differences in BMI z-score were meaningful at even modest volumes of reallocated time. The consistency of the results with previous findings suggests that pooling raw acceleration data from different devices was appropriate, which highlights the utility of this approach. Our novel results highlight the utility of the full physical activity intensity spectrum over a priori-determined absolute intensity cut-point approaches and further emphasise the benefits of promoting higher intensity physical activity for health in youth. Furthermore, they can provide researchers, public health professionals, and physical activity deliverers with important insights to inform obesity prevention intervention design and physical activity programming.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19148778/s1, Supplementary File S1: Table S1: Compositional variation matrix of the time spent by all participants in each intensity band; Table S2: Compositional variation matrix of the time spent by boys in each intensity band; Table S3: Compositional variation matrix of the time spent by girls in each intensity band; Supplementary File S2: Table S4: ANOVA results of the adjusted association analysis between the overall activity spectrum composition and BMI z-score in all participants; Table S5: ANOVA results of the adjusted association analysis between the overall activity spectrum composition and BMI z-score in boys; Table S6: ANOVA results of the adjusted association analysis between the overall activity spectrum composition and BMI z-score in girls; Supplementary File S3: Table S7: One-for-remaining compositional isotemporal substitutions to predict changes in girls’ BMI z-scores resulting from hypothetical time reallocations involving the 50–100 mg intensity band; Table S8: One-for-remaining compositional isotemporal substitutions to predict changes in girls’ BMI z-scores resulting from hypothetical time reallocations involving the 100-150 mg intensity band; Table S9: One-for-remaining compositional isotemporal substitutions to predict changes in girls’ BMI z-scores resulting from hypothetical time reallocations involving the ≥700 mg intensity band; Table S10: One-for-remaining compositional isotemporal substitutions to predict changes in boys’ BMI z-scores resulting from hypothetical time reallocations involving the ≥700 mg intensity band.

Author Contributions

Conceptualisation: S.J.F. and L.M.B.; Methodology: S.J.F., L.M.B., L.H., A.C., M.C., L.F., O.J., D.A.M., R.J.N., M.B.O., J.R.R., S.L.T., L.E.F.G. and R.T. contributed to the data collection in the original studies; Data analysis: S.J.F. with support from D.D. and L.H.; Writing—original draft preparation: S.J.F., D.D., A.G., A.V.R., B.d.P.C. and L.M.B.; Writing—review and editing: all authors commented on previous versions of the manuscript and read and approved the final manuscript; Funding acquisition: where applicable, funding for the original studies was acquired by S.J.F. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for the selected contributing studies was provided by the Waterloo Foundation (#1669/3509), West Lancashire Sport Partnership, West Lancashire Leisure Trust, Edge Hill University, and Wigan Council. Alex Rowlands is supported by the Lifestyle Theme of the Leicester NHR Leicester Biomedical Research Centre and NIHR Applied Research Collaborations East Midlands (ARC-EM). Dorothea Dumuid is supported by the Australian National Health and Medical Research Council (NHMRC) Early Career Fellowship APP1162166 and by the Centre of Research Excellence in Driving Global Investment in Adolescent Health funded by NHMRC APP1171981.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helinski and approved by the Research Ethics Committee of Edge Hill University (#ETH2021-0034; 2 December 2020).

Informed Consent Statement

Informed consent was obtained from all participants involved in the original studies.

Data Availability Statement

The dataset analysed during the current study is available from https://osf.io/kae5r/. The R code used for the data analysis is available from https://osf.io/ze4nk/.

Acknowledgments

The authors thank the participating children and schools for their participation in the contributing studies and West Lancashire Sport Partnership and Liverpool School Sport Partnership for their assistance. Full or partial funding support for the studies from the Waterloo Foundation is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design, execution, interpretation, or writing of the study.

References

  1. Office for National Statistics. National Child Measurement Programme, England 2020–21 School Year. Available online: https://digital.nhs.uk/data-and-information/publications/statistical/national-child-measurement-programme/2020-21-school-year# (accessed on 13 May 2022).
  2. Jago, R.; Drews, K.L.; McMurray, R.G.; Baranowski, T.; Galassetti, P.; Foster, G.D.; Moe, E.; Buse, J.B. BMI change, fitness change and cardiometabolic risk factors among 8th grade youth. Pediatr. Exerc. Sci. 2013, 2, 52–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kolsgaard, M.L.P.; Joner, G.; Brunborg, C.; Anderssen, S.A.; Tonstad, S.; Andersen, L.F. Reduction in BMI z-score and improvement in cardiometabolic risk factors in obese children and adolescents. The Oslo Adiposity Intervention Study—A hospital/public health nurse combined treatment. BMC Pediatr. 2011, 11, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Halfon, N.; Larson, K.; Slusser, W. Associations between obesity and comorbid mental health, developmental, and physical health conditions in a nationally representative sample of US children aged 10 to 17. Acad. Pediatr. 2013, 13, 6–13. [Google Scholar] [CrossRef] [PubMed]
  5. Brown, T.; Moore, T.H.; Hooper, L.; Gao, Y.; Zayegh, A.; Ijaz, S.; Elwenspoek, M.; Foxen, S.C.; Magee, L.; O’Malley, C.; et al. Interventions for preventing obesity in children. Cochrane Database Syst. Rev. 2019, 23, CD001871. [Google Scholar] [CrossRef] [Green Version]
  6. Poitras, V.J.; Gray, C.E.; Borghese, M.M.; Carson, V.; Chaput, J.-P.; Janssen, I.; Katzmarzyk, P.T.; Pate, R.R.; Connor Gorber, S.; Kho, M.E.; et al. Systematic review of the relationships between objectively measured physical activity and health indicators in school-aged children and youth. Appl. Physiol. Nutr. Metab. 2016, 41 (Suppl. 3), S197–S239. [Google Scholar] [CrossRef]
  7. Hamer, M.; Stamatakis, E. Relative proportion of vigorous physical activity, total volume of moderate to vigorous activity, and body mass index in youth: The Millennium Cohort Study. Int. J. Obes. 2018, 42, 1239–1242. [Google Scholar] [CrossRef] [Green Version]
  8. Carson, V.; Hunter, S.; Kuzik, N.; Gray, C.E.; Poitras, V.J.; Chaput, J.-P.; Saunders, T.J.; Katzmarzyk, P.T.; Okely, A.D.; Connor Gorber, S.; et al. Systematic review of sedentary behaviour and health indicators in school-aged children and youth: An update. Appl. Physiol. Nut.r Metab. 2016, 41 (Suppl. 3), S240–S265. [Google Scholar] [CrossRef]
  9. Migueles, J.H.; Aadland, E.; Andersen, L.B.; Brønd, J.C.; Chastin, S.F.; Hansen, B.H.; Konstabel, K.; Kvalheim, O.M.; McGregor, D.E.; Rowlands, A.V.; et al. GRANADA consensus on analytical approaches to assess associations with accelerometer-determined physical behaviours (physical activity, sedentary behaviour and sleep) in epidemiological studies. Br. J. Sports Med. 2022, 56, 376–384. [Google Scholar] [CrossRef]
  10. Rowlands, A.V. Moving forward with accelerometer-assessed physical activity: Two strategies to ensure meaningful, interpretable, and comparable measures. Pediatr. Exerc. Sci. 2018, 30, 450–456. [Google Scholar] [CrossRef] [Green Version]
  11. Aadland, E.; Kvalheim, O.M.; Anderssen, S.A.; Resaland, G.K.; Andersen, L.B. Multicollinear physical activity accelerometry data and associations to cardiometabolic health: Challenges, pitfalls, and potential solutions. Int. J. Behav. Nutr. Phys. Act. 2019, 16, 74. [Google Scholar] [CrossRef]
  12. Aadland, E.; Kvalheim, O.M.; Hansen, B.H.; Kriemler, S.; Ried-Larsen, M.; Wedderkopp, N.; Sardinha, L.B.; Møller, N.C.; Hallal, P.C.; Anderssen, S.A. The multivariate physical activity signature associated with metabolic health in children and youth: An International Children’s Accelerometry Database (ICAD) analysis. Prev. Med. 2020, 141, 106266. [Google Scholar] [CrossRef] [PubMed]
  13. Aadland, E.; Nilsen, A.K.O.; Andersen, L.B.; Rowlands, A.V.; Kvalheim, O.M. A comparison of analytical approaches to investigate associations for accelerometry-derived physical activity spectra with health and developmental outcomes in children. J. Sports Sci. 2021, 39, 430–438. [Google Scholar] [CrossRef] [PubMed]
  14. Aadland, E.; Andersen, L.B.; Migueles, J.H.; Ortega, F.B.; Kvalheim, O.M. Interpretation of associations between the accelerometry physical activity spectrum and cardiometabolic health and locomotor skills in two cohorts of children using raw, normalized, log-transformed, or compositional data. J. Sports Sci. 2020, 38, 2708–2719. [Google Scholar] [CrossRef] [PubMed]
  15. Verloigne, M.; Van Lippevelde, W.; Maes, L.; Yildirim, M.; Chinapaw, M.; Manios, Y. Levels of physical activity and sedentary time among 10- to 12-year-old boys and girls across 5 European countries using accelerometers: An observational study within the ENERGY-project. Int J Behav Nutr. Phys. Act. 2012, 9, 34. [Google Scholar] [CrossRef] [Green Version]
  16. Beltran-Valls, M.R.; Janssen, X.; Farooq, A.; Adamson, A.J.; Pearce, M.S.; Reilly, J.K.; Basterfield, L.; Reilly, J.J. Longitudinal changes in vigorous intensity physical activity from childhood to adolescence: Gateshead Millennium Study. J. Sci. Med. Sport. 2019, 22, 450–455. [Google Scholar] [CrossRef] [Green Version]
  17. DHSC. UK Chief Medical Officers’ Physical Activity Guidelines; DHSC: London, UK, 2019.
  18. Costigan, S.A.; Eather, N.; Plotnikoff, R.C.; Taaffe, D.R.; Lubans, D.R. High-intensity interval training for improving health-related fitness in adolescents: A systematic review and meta-analysis. Br. J. Sports Med. 2015, 49, 1253–1261. [Google Scholar] [CrossRef]
  19. Janssen, I.; Ross, R. Vigorous intensity physical activity is related to the metabolic syndrome independent of the physical activity dose. Int. J. Epidemiol. 2012, 41, 1132–1140. [Google Scholar] [CrossRef]
  20. Fairclough, S.J.; Dumuid, D.; Mackintosh, K.A.; Stone, G.; Dagger, R.; Stratton, G.; Davies, I.; Boddy, L.M. Adiposity, fitness, health-related quality of life and the reallocation of time between children’s school day activity behaviours: A compositional data analysis. Prev. Med. Rep. 2018, 11, 254–261. [Google Scholar] [CrossRef]
  21. Fairclough, S.J.; Dumuid, D.; Taylor, S.; Curry, W.; McGrane, B.; Stratton, G.; Maher, C.; Olds, T. Fitness, fatness and the reallocation of time between children’s daily movement behaviours: An analysis of compositional data. Int. J. Behav. Nutr. Phys. Act. 2017, 14, 64. [Google Scholar] [CrossRef]
  22. Fairclough, S.J.; Tyler, R.; Dainty, J.R.; Dumuid, D.; Richardson, C.; Shepstone, L.; Atkin, A.J. Cross-sectional associations between 24-hour activity behaviours and mental health indicators in children and adolescents: A compositional data analysis. J. Sports Sci. 2021, 39, 1602–1614. [Google Scholar] [CrossRef]
  23. Burns, R.D.; Kim, Y.; Byun, W.; Brusseau, T.A. Associations of school day sedentary behavior and physical activity with gross motor skills: Use of compositional data analysis. J. Phys. Act. Health 2019, 16, 811–817. [Google Scholar] [CrossRef] [PubMed]
  24. Talarico, R.; Janssen, I. Compositional associations of time spent in sleep, sedentary behavior and physical activity with obesity measures in children. Int. J. Obes. 2018, 42, 1508–1514. [Google Scholar] [CrossRef] [PubMed]
  25. Dumuid, D.; Stanford, T.E.; Pedišić, Ž.; Maher, C.; Lewis, L.K.; Martín-Fernández, J.-A.; Katzmarzyk, P.T.; Chaput, J.-P.; Fogelholm, M.; Standage, M.; et al. Adiposity and the isotemporal substitution of physical activity, sedentary time and sleep among school-aged children: A compositional data analysis approach. BMC Public Health 2018, 18, 311. [Google Scholar] [CrossRef]
  26. Dumuid, D.; Maher, C.; Lewis, L.K.; Stanford, T.E.; Martin Fernandez, J.A.; Ratcliffe, J.; Katzmarzyk, P.T.; Barreira, T.V.; Chaput, J.P.; Fogelholm, M.; et al. Human development index, children’s health-related quality of life and movement behaviors: A compositional data analysis. Qual. Life Res. 2018, 27, 473–482. [Google Scholar] [CrossRef] [PubMed]
  27. Carson, V.; Tremblay, M.S.; Chaput, J.-P.; Chastin, S.F.M. Associations between sleep duration, sedentary time, physical activity, and health indicators among Canadian children and youth using compositional analyses. Appl. Physiol. Nutr. Metab. 2016, 41 (Suppl. 3), S294–S302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Taylor, S.L.; Noonan, R.J.; Knowles, Z.R.; Owen, M.B.; McGrane, B.; Curry, W.B.; Fairclough, S.J. Predictors of segmented school day physical activity and sedentary time in children from a northwest England low-income community. Int. J. Environ. Res. Public Health 2017, 14, 534. [Google Scholar] [CrossRef]
  29. Taylor, S.L.; Noonan, R.J.; Knowles, Z.R.; Owen, M.B.; McGrane, B.; Curry, W.B.; Fairclough, S.J. Evaluation of a pilot school-based physical activity clustered randomised controlled trial-Active Schools: Skelmersdale. Int. J. Environ. Res. Public Health 2018, 15, 1011. [Google Scholar] [CrossRef] [Green Version]
  30. Hurter, L.; Rowlands, A.V.; Fairclough, S.J.; Gibbon, K.C.; Knowles, Z.R.; Porcellato, L.A.; Cooper-Ryan, A.M.; Boddy, L.M. Validating the Sedentary Sphere method in children: Does wrist or accelerometer brand matter? J. Sports Sci. 2019, 37, 1910–1918. [Google Scholar] [CrossRef]
  31. Foweather, L.; Crotti, M.; Foulkes, J.D.; O’Dwyer, M.V.; Utesch, T.; Knowles, Z.R.; Fairclough, S.J.; Ridgers, N.D.; Stratton, G. Foundational movement skills and play behaviors during recess among preschool children: A compositional analysis. Children 2021, 8, 543. [Google Scholar] [CrossRef]
  32. Noonan, R.J.; Boddy, L.M.; Kim, Y.; Knowles, Z.R.; Fairclough, S.J. Comparison of children’s free-living physical activity derived from wrist and hip raw accelerations during the segmented week. J. Sports Sci. 2017, 35, 2067–2072. [Google Scholar] [CrossRef]
  33. Owen, M.; Kerner, C.; Taylor, S.; Noonan, R.; Newson, L.; Kosteli, M.-C.; Curry, W.; Fairclough, S.J. The feasibility of a novel school peer-led mentoring model to improve the physical activity levels and sedentary time of adolescent girls: The Girls Peer Activity (G-PACT) project. Children 2018, 5, 67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lohman, T.G.; Roche, A.F.M.; Martorell, R. Anthropometric Standardization Reference Manual; Human Kinetics Books: Champaign, IL, USA, 1991. [Google Scholar]
  35. Cole, T.J.; Freeman, J.V.; Preece, M.A. Body mass index reference curves for the UK, 1990. Arch. Dis. Child 1995, 73, 25–29. [Google Scholar] [CrossRef] [Green Version]
  36. Cole, T.J.; Bellizzi, M.; Flegal, K.; Dietz, W. Establishing a standard definition for child overweight and obesity worldwide: International survey. Br. Med. J. 2000, 320, 1240–1243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Department for Communities and Local Government. The English Indices of Deprivation. 2015. Available online: https://www.gov.uk/government/statistics/english-indices-of-deprivation-2015 (accessed on 12 December 2017).
  38. Ministry of Housing, Communities, and Local Government. English Indices of Deprivation 2019. Available online: https://www.gov.uk/government/statistics/english-indices-of-deprivation-2019 (accessed on 24 April 2022).
  39. Migueles, J.H.; Rowlands, A.V.; Huber, F.; Sabia, S.; Van Hees, V.T. GGIR: A Research community–driven open source R package for generating physical activity and sleep outcomes from multi-day raw accelerometer data. J. Meas. Phys. Behav. 2019, 2, 188–196. [Google Scholar] [CrossRef] [Green Version]
  40. van Hees, V.T.; Fang, Z.; Langford, J.; Assah, F.; Mohammad, A.; da Silva, I.C.; Trenell, M.I.; White, T.; Wareham, N.J.; Brage, S. Autocalibration of accelerometer data for free-living physical activity assessment using local gravity and temperature: An evaluation on four continents. J. Appl. Physiol. 2014, 11, 738–744. [Google Scholar] [CrossRef] [Green Version]
  41. van Hees, V.T.; Renstrom, F.; Wright, A.; Gradmark, A.; Catt, M.; Chen, K.Y.; Lof, M.; Bluck, L.; Pomeroy, J.; Wareham, N.J. Estimation of daily energy expenditure in pregnant and non-pregnant women using a wrist-worn tri-axial accelerometer. PLoS ONE 2011, 6, e22922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. van Hees, V.T.; Gorzelniak, L.; Dean León, E.C.; Eder, M.; Pias, M.; Taherian, S.; Ekelund, U.; Renström, F.; Franks, P.W.; Horsch, A. Separating movement and gravity components in an acceleration signal and implications for the assessment of human daily physical activity. PLoS ONE 2013, 8, e61691. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Rowlands, A.V.; Plekhanova, T.; Yates, T.; Mirkes, E.M.; Davies, M.; Khunti, K.; Edwardson, C.L. Providing a basis for harmonization of accelerometer-assessed physical activity outcomes across epidemiological datasets. J. Meas. Phys. Behav. 2019, 2, 131–142. [Google Scholar] [CrossRef]
  44. Phillips, L.R.; Parfitt, G.; Rowlands, A.V. Calibration of the GENEA accelerometer for assessment of physical activity intensity in children. J. Sci. Med. Sport 2012, 16, 124–128. [Google Scholar] [CrossRef]
  45. Hildebrand, M.; Van Hees, V.T.; Hansen, B.H.; Ekelund, U. Age-group comparibility of raw accelerometer output from wrist- and hip-worn monitors. Med. Sci. Sports Exerc. 2014, 46, 1816–1824. [Google Scholar] [CrossRef]
  46. Sterne, J.A.C.; White, I.R.; Carlin, J.B.; Spratt, M.; Royston, P.; Kenward, M.G.; Wood, A.M.; Carpenter, J.R. Multiple imputation for missing data in epidemiological and clinical research: Potential and pitfalls. Br. Med. J. 2009, 338, b2393. [Google Scholar] [CrossRef] [PubMed]
  47. van den Boogaart, K.G.; Tolosana-Delgado, R. ‘Compositions’: A unified R package to analyze compositional data. Comput. Geosci. 2008, 34, 320–338. [Google Scholar] [CrossRef]
  48. Dumuid, D.; Stanford, T.E.; Martin-Fernandez, J.A.; Pedisic, Z.; Maher, C.A.; Lewis, L.K.; Hron, K.; Katzmarzyk, P.T.; Chaput, J.P.; Fogelholm, M. Compositional data analysis for physical activity, sedentary time and sleep research. Stat. Methods Med. Res. 2018, 27, 3736–3738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Dumuid, D.; Pedišić, Ž.; Palarea-Albaladejo, J.; Martín-Fernández, J.A.; Hron, K.; Olds, T. Compositional data analysis in time-use epidemiology: What, why, how. Int. J. Environ. Res. Public Health 2020, 17, 2220. [Google Scholar] [CrossRef] [Green Version]
  50. Zeileis, A.; Hothorn, T. Diagnostic checking in regression relationships. R News 2002, 2, 7–10. [Google Scholar]
  51. Fox, J.; Weisberg, S. An R Companion to Applied Regression; Sage: Thousand Oaks, CA, USA, 2019. [Google Scholar]
  52. Ludecke, D.; Ben-Shachar, M.S.; Patil, I.; Waggoner, P.; Makowski, D. performance: An R package for assessment, comparison and testing of statistical models. J. Open Sci. Softw. 2021, 60, 3139. [Google Scholar] [CrossRef]
  53. Dumuid, D.; Wake, M.; Clifford, S.; Burgner, D.; Carlin, J.B.; Mensah, F.K.; Fraysse, F.; Lycett, K.; Baur, L.; Olds, T. The association of the body composition of children with 24-hour activity composition. J. Pediatr. 2019, 208, 43–49. [Google Scholar] [CrossRef]
  54. Stanford, T.E. Tystan/Deltacomp: Deltacomp Version 0.2.2. 2022. Available online: https://zenodo.org/record/6340017 (accessed on 5 May 2022).
  55. Aadland, E.; Andersen, L.B.; Anderssen, S.A.; Resaland, G.K.; Kvalheim, O.M. Associations of volumes and patterns of physical activity with metabolic health in children: A multivariate pattern analysis approach. Prev. Med. 2018, 115, 12–18. [Google Scholar] [CrossRef] [Green Version]
  56. Aadland, E.; Kvalheim, O.M.; Anderssen, S.A.; Resaland, G.K.; Andersen, L.B. The multivariate physical activity signature associated with metabolic health in children. Int. J. Behav. Nutr. Phys. Act. 2018, 15, 77. [Google Scholar] [CrossRef]
  57. Štefelová, N.; Dygrýn, J.; Hron, K.; Gába, A.; Rubín, L.; Palarea-Albaladejo, J. Robust compositional analysis of physical activity and sedentary behaviour data. Int. J. Environ. Res. Public Health 2018, 15, 2248. [Google Scholar] [CrossRef] [Green Version]
  58. Rubín, L.; Gába, A.; Pelclová, J.; Štefelová, N.; Jakubec, L.; Dygrýn, J.; Hron, K. Changes in sedentary behavior patterns during the transition from childhood to adolescence and their association with adiposity: A prospective study based on compositional data analysis. Arch. Public Health 2022, 80, 1–9. [Google Scholar] [CrossRef]
  59. Carson, V.; Tremblay, M.S.; Chaput, J.-P.; McGregor, D.; Chastin, S. Compositional analyses of the associations between sedentary time, different intensities of physical activity, and cardiometabolic biomarkers among children and youth from the United States. PLoS ONE 2019, 14, e0220009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Pedišić, Ž.; Bauman, A. Accelerometer-based measures in physical activity surveillance: Current practices and issues. Br. J. Sports Med. 2015, 49, 219–223. [Google Scholar] [CrossRef] [PubMed]
  61. Trost, S.G. Population-level physical activity surveillance in young people: Are accelerometer-based measures ready for prime time? Int. J. Behav. Nutr. Phys. Act. 2020, 17, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Hurter, L.; Fairclough, S.J.; Knowles, Z.; Porcellato, L.; Cooper-Ryan, A.; Boddy, L. Establishing raw acceleration thresholds to classify sedentary and stationary behaviour in children. Children 2018, 5, 172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Chaput, J.-P.; Gray, C.E.; Poitras, V.J.; Carson, V.; Gruber, R.; Olds, T.; Weiss, S.K.; Connor Gorber, S.; Kho, M.E. Systematic review of the relationships between sleep duration and health indicators in school-aged children and youth. Appl. Physiol. Nutr. Metab. 2016, 41 (Suppl. 3), S266–S282. [Google Scholar] [CrossRef]
  64. Kwon, S.; Janz, K.F.; Burns, T.L.; Levy, S.M. Association between light-intensity physical activity and adiposity in childhood. Pediatr Exerc Sci. 2011, 23, 218–229. [Google Scholar] [CrossRef] [Green Version]
  65. Rowlands, A.V.; Fairclough, S.J.; Yates, T.O.M.; Edwardson, C.; Davies, M.; Munir, F.; Khunti, K.; Stiles, V. Activity intensity, volume, and norms: Utility and interpretation of accelerometer metrics. Med. Sci. Sports Exerc. 2019, 51, 2410–2422. [Google Scholar] [CrossRef]
  66. Gába, A.; Dygrýn, J.; Štefelová, N.; Rubín, L.; Hron, K.; Jakubec, L. Replacing school and out-of-school sedentary behaviors with physical activity and its associations with adiposity in children and adolescents: A compositional isotemporal substitution analysis. Environ. Health Prev. Med. 2021, 26, 16. [Google Scholar] [CrossRef]
  67. Sattelmair, J.; Pertman, J.; Ding, E.L.; Kohl, H.W.; Haskell, W.; Lee, I.M. Dose response between physical activity and risk of coronary heart disease: A meta-analysis. Circulation 2011, 124, 789–795. [Google Scholar] [CrossRef] [Green Version]
  68. Fairclough, S.J.; Hackett, A.; Davies, I.; Gobbi, R.; Mackintosh, K.; Warburton, G.; Stratton, G.; van Sluijs, E.; Boddy, L. Promoting healthy weight in primary school children through physical activity and nutrition education: A pragmatic evaluation of the CHANGE! randomised intervention study. BMC Public Health. 2013, 13, 626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Chastin, S.F.M.; Palarea-Albaladejo, J.; Dontje, M.L.; Skelton, D.A. Combined effects of time spent in physical activity, sedentary behaviors and sleep on obesity and cardio-metabolic health markers: A novel compositional data analysis approach. PLoS ONE 2015, 10, e0139984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijerph 19 08778 g001 550
Figure 1. (a–d) Predicted differences in BMI z-score for the time reallocations between the most dominant intensity bands and the remaining intensity bands. Note: (a) ≥700 mg intensity band (boys); (b) 50–100 mg (girls); (c) 100–150 mg (girls); (d) ≥700 mg (girls). Note: coloured lines represent the lower and upper boundaries of the 95% confidence intervals.
Figure 1. (a–d) Predicted differences in BMI z-score for the time reallocations between the most dominant intensity bands and the remaining intensity bands. Note: (a) ≥700 mg intensity band (boys); (b) 50–100 mg (girls); (c) 100–150 mg (girls); (d) ≥700 mg (girls). Note: coloured lines represent the lower and upper boundaries of the 95% confidence intervals.
Ijerph 19 08778 g001
Table
Table 1. Participants’ descriptive characteristics (M (SD) or %).
Table 1. Participants’ descriptive characteristics (M (SD) or %).
All (n= 1453)Boys (n = 624)Girls (n = 829)
Age (years)10.5 (2.6)10.0 (2.6)10.8 (2.5)
Height (cm)142.1 (16.2)139.8 (16.8)143.7 (15.6)
Weight (kg)39.6 (14.8)37.3 (14.4)41.3 (15.0)
BMI (kg·m−2)19.0 (3.9)18.4 (3.6)19.4 (4.1)
BMI z-score0.51 (1.24)0.53 (1.29)0.49 (1.21)
Weight status
  Normal weight (%)73.474.272.9
  Overweight/obese (%)26.625.827.1
EIMD decile
  Deciles 1–5 (%)67.471.364.5
  Decile 6–10 (%)32.628.735.5
School type
  Primary (%)66.372.461.8
  Secondary (%)33.727.638.2
Note: BMI = body mass index; EIMD = English Indices of Multiple Deprivation.
Table
Table 2. Geometric means of the time spent in the intensity spectrum bands.
Table 2. Geometric means of the time spent in the intensity spectrum bands.
All (n = 1453)Boys (n = 641)Girls (n = 862)
Intensity Bandmin·day−1%min·day−1%min·day−1%
0–50 mg 722.975.3719.974.972475.5
50–100 mg 105.711102.610.710811.2
100–150 mg 51.15.349.65.2525.4
150–200 mg 27.82.927.82.9282.9
200–250 mg 15.41.616.11.7151.5
250–300 mg 9.10.99.8190.9
300–350 mg 5.80.66.40.750.6
350–700 mg 14.31.516.71.7131.3
≥700 mg 8.00.811.11.160.6
Note: mg = milligravitational units.
Table
Table 3. Regression model results for boys and girls in assessing the compositional association between each intensity spectrum band ILR1 and BMI z-score, relative to the remaining intensity bands, with adjustment for SES, centred age, accelerometer model, and accelerometer sampling frequency.
Table 3. Regression model results for boys and girls in assessing the compositional association between each intensity spectrum band ILR1 and BMI z-score, relative to the remaining intensity bands, with adjustment for SES, centred age, accelerometer model, and accelerometer sampling frequency.
BoysGirls
Intensity Band ILR1 (mg)βILR195% CIpβILR195% CIp
0–50 mg vs. remaining−0.20−0.61, 0.210.340.03−0.28, 0.340.87
50–100 mg vs. remaining−0.74−1.92, 0.440.221.39 0.53, 2.250.002
100–150 mg vs. remaining−0.34−2.65, 1.970.76−2.55−4.18, −0.920.002
150–200 mg vs. remaining0.88−1.96, 3.720.551.87−0.05, 3.790.06
200–250 mg vs. remaining−0.50−3.24, 2.240.72−1.64−3.56, 0.280.11
250–300 mg vs. remaining−1.75−4.14, 0.640.151.04−0.51, 2.590.19
300–350 mg vs. remaining1.39−0.61, 3.390.17−0.09−1.38, 1.200.89
350–700 vs. remaining0.57−0.41, 1.550.260.39−0.39, 0.990.26
≥700 mg vs. remaining−0.77−1.08, −0.46<0.001−0.71−0.91, −0.51<0.001
Note: BMI = body mass index; CI = confidence interval; ILR = isometric log-ratio; mg = milligravitational units.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fairclough, S.J.; Hurter, L.; Dumuid, D.; Gába, A.; Rowlands, A.V.; Cruz, B.d.P.; Cox, A.; Crotti, M.; Foweather, L.; Graves, L.E.F.; et al. The Physical Behaviour Intensity Spectrum and Body Mass Index in School-Aged Youth: A Compositional Analysis of Pooled Individual Participant Data. Int. J. Environ. Res. Public Health 2022, 19, 8778. https://doi.org/10.3390/ijerph19148778

AMA Style

Fairclough SJ, Hurter L, Dumuid D, Gába A, Rowlands AV, Cruz BdP, Cox A, Crotti M, Foweather L, Graves LEF, et al. The Physical Behaviour Intensity Spectrum and Body Mass Index in School-Aged Youth: A Compositional Analysis of Pooled Individual Participant Data. International Journal of Environmental Research and Public Health. 2022; 19(14):8778. https://doi.org/10.3390/ijerph19148778

Chicago/Turabian Style

Fairclough, Stuart J., Liezel Hurter, Dorothea Dumuid, Ales Gába, Alex V. Rowlands, Borja del Pozo Cruz, Ashley Cox, Matteo Crotti, Lawrence Foweather, Lee E. F. Graves, and et al. 2022. "The Physical Behaviour Intensity Spectrum and Body Mass Index in School-Aged Youth: A Compositional Analysis of Pooled Individual Participant Data" International Journal of Environmental Research and Public Health 19, no. 14: 8778. https://doi.org/10.3390/ijerph19148778

APA Style

Fairclough, S. J., Hurter, L., Dumuid, D., Gába, A., Rowlands, A. V., Cruz, B. d. P., Cox, A., Crotti, M., Foweather, L., Graves, L. E. F., Jones, O., McCann, D. A., Noonan, R. J., Owen, M. B., Rudd, J. R., Taylor, S. L., Tyler, R., & Boddy, L. M. (2022). The Physical Behaviour Intensity Spectrum and Body Mass Index in School-Aged Youth: A Compositional Analysis of Pooled Individual Participant Data. International Journal of Environmental Research and Public Health, 19(14), 8778. https://doi.org/10.3390/ijerph19148778

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop