The inter-relationship between physical activity, sedentary behaviour and sleep (collectively defined as physical behaviours) is of interest to researchers from different fields. Each of these physical behaviours has been investigated in epidemiological studies, yet their codependency and interactions need to be further explored and accounted for in data analysis. Modern accelerometers capture continuous movement through the day, which presents the challenge of how to best use the richness of these data. In recent years, analytical approaches first applied in other scientific fields have been applied to physical behaviour epidemiology (eg, isotemporal substitution models, compositional data analysis, multivariate pattern analysis, functional data analysis and machine learning). A comprehensive description, discussion, and consensus on the strengths and limitations of these analytical approaches will help researchers decide which approach to use in different situations. In this context, a scientific workshop and meeting were held in Granada to discuss: (1) analytical approaches currently used in the scientific literature on physical behaviour, highlighting strengths and limitations, providing practical recommendations on their use and including a decision tree for assisting researchers’ decision-making; and (2) current gaps and future research directions around the analysis and use of accelerometer data. Advances in analytical approaches to accelerometer-determined physical behaviours in epidemiological studies are expected to influence the interpretation of current and future evidence, and ultimately impact on future physical behaviour guidelines.
- physical activity
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Physical activity (PA), sedentary behaviour (SB) and sleep, collectively described as physical behaviours,1 are of interest to many researchers from previously separate fields. Accounting for the inter-relations of these behaviours is important because: (1) they share the 24 hours of the day (ie, closure), so change in one behaviour results in change in others; and (2) the relation of a specific behaviour with health depends on other behaviours (eg, SB and mortality relation depends on moderate-to-vigorous PA (MVPA)2 3). These inter-relations should be considered in modelling their association with health, with attention to collinearity issues potentially leading to spurious findings.
Accelerometers are increasingly being used to estimate different constructs/dimensions of physical behaviours (eg, types (walking, cycling, dancing), intensities (light, moderate, vigorous), and postures (reclining, sitting, standing)). Other constructs focus more on the description of the acceleration signal (eg, time spent within acceleration bands with no energy expenditure interpretation (eg, min/day between 0 and 100 mg), or the acceleration above which the most active 30 min of the day occur). The data-analytical approach usually includes: (1) reduction of the acceleration signal into meaningful behaviours/descriptors; (2) mathematical treatment of the descriptors if needed and (3) selection of the statistical model. Multiple choices are available for each step, and decisions should be adapted to the research question and account for potential collinearity issues arising from behaviours’ inter-relationships. However, there are currently no consensus/recommendations to help to choose the most appropriate approach. Online supplemental appendix 1 presents the different choices for descriptors, mathematical treatments and statistical models discussed in the Analytical approaches section.
The ‘International Workshop: A focus on statistical methods to analyse accelerometer-measured PA’ was held in Granada on 21–22 October 2019. This event brought together a panel of researchers to discuss, reach consensus, and provide recommendations about the most frequently used analytical approaches in the field and about future research directions in physical behaviour epidemiology. The focus was on modelling physical behaviour constructs (mainly related to PA and SB, although we also included sleep to cover the 24-hour continuum) as exposure variables and health indicators as outcomes. We covered time-use descriptors as those quantified in time over the day, and acceleration-based as those quantified as acceleration magnitude.
Data collection decisions are outside our focus, although decisions on body attachment site,4–8 number of days recorded,9 treatment of weekdays and weekend days,10 11 seasonality,12 among others, affect the ability of accelerometer data to identify specific constructs/aspects of physical behaviours. For example, attaching the accelerometer to the hip, wrist or thigh may be considered depending on the constructs of interest (eg, PA intensity,13–15 postures16 or sleep patterns,17 18 among others). A recent consensus report discussed best practices on these decisions.19
Analytical approaches: discussion and practical considerations
This section discusses different analytical approaches’ applicability in various situations (or research questions). Analytical approaches include the combination of accelerometer descriptors (table 1, online supplemental appendix 1 (Section 1)) and statistical models (table 2, online supplemental appendix 1 (Section 3)) with and without mathematical (compositional) transformation (online supplemental appendix 1 (Section 2)). We provide practical considerations on (1) informativeness of each analytical approach for public health messaging and (2) appropriateness of the analytical approaches for certain research questions. Additionally, table 3 shows the performance of these approaches regarding closure or collinearity, relationship assumptions and interpretation for PA guidelines.
Total PA and linear regression
Average acceleration (or steps per day) provides the simplest estimate of the overall movement and proxy for total daily PA-related energy expenditure. Statistical interpretation of findings using linear regression is straightforward since there is a single variable representing the overall activity volume. Thus, codependence with other explanatory variables is not usually a concern and linear regression models are an option for the analysis. The opinion of the consensus group is that the average acceleration is useful for reducing the confounding effect of PA in a given association analysis (eg, is the association of sugar consumption with body mass index dependent on overall PA?), or as the main exposure in cases where it explains a large proportion of the PA-related energy expenditure in a certain cohort (eg, is PA-related energy expenditure associated with protein intake?). Beyond this, the average acceleration alone is not very informative relative to associations of specific physical behaviours with health outcomes, limiting its applicability for public health messaging. A recent study proposed the minimum clinically informative difference for average acceleration from wrist data,20 but further studies are needed. Although these descriptors cannot be interpreted in terms of meeting or not meeting the PA guidelines, they may be the best descriptor to test the ‘move more’ message reported in several guidelines.
Time-use behaviours or intensity spectrum and linear regression
Among time-use constructs, time spent in PA intensities is the most frequently used, while PA types and postures have gained momentum recently. These descriptors are often introduced in linear regression models to test the association of time spent in a certain intensity/behaviour with health outcomes. As it is widely used, it is useful for comparing estimates with other cohorts. The intensity spectrum is an extension of PA intensities with higher resolution energy bands. When using such time-use behaviours, requirements for bouts in these behaviours should be considered. We observe a lack of consensus in the literature on how a bout should be calculated, including the definition of both acceptable allowance drop period without terminating the bout and minimum and maximum duration. Bouts of 30 min for SB and 10 min for MVPA, often allowing short time intervals outside the behaviour of interest, are frequently used.21 It is unclear how much these choices are driven by a desire for harmonisation, by public health guidelines, or by evidence. Although observational data based on 1-minute or longer epoch lengths suggest that any bout duration can produce health benefits, randomised controlled trials investigating the effects of differing bout durations on health outcomes are lacking.21 Based on the observational studies, the recommendation about accumulating PA in certain bout durations has been excluded from recent guidelines.22
Time-use behaviours (or intensity spectrum) include a set of codependent variables, and linear regression does not handle closure and collinearity among explanatory variables. When using these descriptors, linear regression adjusted for all physical behaviour components may be affected by multicollinearity.23 Variance inflation factors are unable to explain inconsistencies between linear regression models sequentially excluding a behaviour from the explanatory variables24 and might not be an acceptable diagnostic indicator for the interdependency between time-use descriptors.24 Additionally, the assumption of linearity of the association between these descriptors and health outcome might not be sufficiently met for analyses to yield valid results. This consensus group recommends moving towards other analytical approaches more suitable for studying the codependencies among time-use behaviours. In this regard, transforming time-use behaviours using the compositional data transformation (isometric log ratio (ILR), online supplemental appendix 1 (Section 2)) represents an option. Using the ILR transformation, each variable indicates the time spent in a given behaviour relative to the time spent in the rest of behaviours of the composition (eg, SB, light PA, MVPA and sleep). In other words, it quantifies the effect of increasing the time in a behaviour while proportionally reducing the time in the rest. Pair-wise reallocations of time can be interpreted from linear regression predictions on specific time compositions arising from hypothetical reallocations of time rather than from regression coefficients (as in isotemporal substitution models, online supplemental appendix 1 (Section 3.2)). By transforming the variables, the codependency among the time-use descriptors relative to their time closure is solved (ie, it accounts for the codependency of time among variables). However, transformed variables can still be collinear, and collinearity should be investigated because linear regression cannot handle collinearity, regardless of its source. This is especially problematic when analysing the intensity spectrum since it provides a wide range of variables (usually more than 10) that are highly correlated, even if using ILR-transformed variables.25 As such, we recommend testing the correlations and risk of collinearity among the explanatory variables (even when compositionally transformed). In absence of high correlations and collinearity, linear regression can be appropriate.
The opinion of the consensus group is that physical behaviour epidemiology should move to studying the combined effects and interactions of physical behaviours on health, and a feasible option is using ILR-transformed time-use descriptors and linear regression.24 26 This approach is informative for public health messaging as it provides information on combinations of behaviours (considering every behaviour that occurs in the day) which are beneficial for health. Clustering groups of people based on their behaviours is also an alternative to investigating the interactions between behaviours, although compositional analyses allow the variables to be studied on a continuous scale. With the intensity spectrum, the use of linear regression models is not possible because of collinearity issues in the variables (either transformed or not).25 The collinearity problem, however, can be solved by using partial least square (PLS) regression. Regression models can be used in different study designs, including longitudinal studies, either with absolute27–29 or compositional data.30 Linear regression with compositional data may need appropriate graphical representation of the results to interpret the magnitude of the association.24
Time-use behaviours or intensity spectrum and isotemporal substitution models
Isotemporal substitution modelling carries forward the main limitations of linear regression, that is, multicollinearity and assumption of linearity (as the magnitude of the association is derived from regression coefficients). These important limitations preclude us from recommending the use of isotemporal models with time-use descriptors. However, it is notable that this approach provides broadly similar findings to compositional ILR transformation of time-use descriptors and linear regression.31 Public health messaging can be complemented with information on the effect of reallocating daily time across behaviours (either with isotemporal substitution models or with linear regression with compositional data, the Time-use behaviours or intensity spectrum and linear regression section). The intensity spectrum has not been analysed with isotemporal substitution models thus far. We do not recommend such an analysis since the large number of variables in the intensity spectrum would complicate the interpretation.
Time-use behaviours or intensity spectrum and multivariate pattern analysis
Multivariate pattern analysis fully handles the collinearity among explanatory variables using latent variable modelling. Collinearity is approached as a dimensionality reduction problem in which the variance of the explanatory variables shared with the outcome is retained. Multivariate pattern analysis describes the pattern of associations for the descriptors with the outcome, accounting for the correlated structure of the data. Associations with health are interpreted for each descriptor (each PA intensity or band in the intensity spectrum) considering its codependency with the rest, but without quantifying time exchange between descriptors. A limitation of this analytical approach is that PLS regression models cannot be adjusted as usually done in linear regression. If covariates are included in the PLS model, they will contribute their shared variance with PA and the outcome. Aadland et al proposed obtaining residuals for the outcome from a linear regression model including confounders as explanatory variables, prior to entering the outcome variable in the PLS model.25 32 33 This challenge remains for the analysis of categorical or time-to-event outcomes (eg, mortality).
Likewise, time-use behaviours or the intensity spectrum could be transformed as compositional data before performing multivariate pattern analysis. Since multivariate pattern analysis can handle singular data, the use of ILR coordinates is not necessary. Aadland et al recently compared the use of centred log ratio (CLR)-transformed time-use and intensity spectrum descriptors with respect to associations with metabolic health using multivariate pattern analysis.25 While associations appeared to differ, the interpretation of associations, considering the absolute and relative interpretation, were partly equivalent. The interpretation of CLR-transformed variables may not be very informative for public health messaging as they represent the effect of time exchange from the geometric mean of the time-use descriptors distribution to a specific time-use descriptor (eg, MVPA or any intensity spectrum band).
Other similar alternatives to reduce dimensionality of the data while retaining relevant information by increasing covariance among descriptors (rather than with the outcome) include factor analysis, principal component analysis, or joint and individual variation explained.34 This consensus group recommends considering these approaches to analyse many explanatory variables (eg, intensity spectrum) in relation to health.35 There is no clear recommendation on the number of bands (or number of explanatory variables) to generate for this analytical approach, though previous studies have used from 16 (uniaxial data)32 to 102 (triaxial data)33 intensity bands. Resolution (number of bands) may influence the relationship with the outcome and depend on the sample characteristics; thus, further research is needed.
Intensity gradient and linear regression
The intensity gradient describes the straight line negative slope of the natural logs of time and acceleration intensity.36 The intensity gradient was developed to: (1) capture the entire intensity distribution, (2) avoid reliance on population and protocol-specific calibration protocols, and (3) provide information that complements average acceleration. The intensity gradient can be used alongside average acceleration to more fully describe the 24-hour movement profile by capturing both volume and intensity of PA. Using the intensity gradient and average acceleration together in linear regression models allows investigation of independent, additive and interactive associations of volume and intensity of PA with health. More work is needed to interpret the intensity gradient relative to the adherence to PA guidelines.
Intensity gradient or intensity distribution and functional data analysis
The acceleration distribution over time of the day, the acceleration density or the intensity gradient function can be used in functional data analysis. Using scalar-on-function data analysis,37 these acceleration functions can be used as an explanatory variable in regression models including linear,38 logistic or Cox regression models. For example, in the case of the acceleration density function as explanatory variable, the association with the event of interest is described along the acceleration range.38 This shows acceleration sections that are associated with the outcome by accounting for the full acceleration distribution, allowing identification of a cut-point such that proportion of time spent above this acceleration cut-point is associated with the outcome. Once these cut-points are identified, it is possible to estimate differences in the outcome by allocating time below to time above this cut-point.39
Functional data analysis has several advantages: (1) it is not affected by multicollinearity since it handles the data continuity; (2) it can test the effects of time reallocation and thus consider closure; and (3) it detects sections of the accelerometer data that are important for a certain health outcome, thus relaxing assumptions of linearity in particular behaviours made by other statistical models. Among its main limitations, acceleration functions usually carry much information that may be irrelevant to the outcome, but is considered in the analysis. Its main drawback is difficulty translating the findings into useful, straightforward public health messages. Investigation of how to make the conclusions of functional data analysis relevant to public health guidelines is highly encouraged by this consensus group (see39 for an example).
MX metrics and linear regression
MX metrics represent the acceleration above which a person’s most active non-consecutive X minutes over the day are spent. An advantage of using MX metrics is that analysis is not affected by cut-point assumptions on energy expenditure, while cut-points may be post-hoc applied to enable public health messaging. For example, if the M60 of a child is 230 mg, this can be compared with an MVPA cut-point, for example, 200 mg, 14 showing that the child meets the 60 min daily MVPA recommendation. However, if compared with a more stringent 250 mg MVPA cut-point, the child does not reach the recommendation. The post-hoc application of cut-points can, therefore, be skipped and keep the interpretation to the descriptive MX values instead. Another advantage of this approach is that the intensity of PA for the specified duration is captured regardless of how inactive a person is. Regarding statistical modelling, the MX metrics usually include a wide range of variables (table 1, online supplemental appendix 1 (Section 1)). These MX metrics are likely to be codependent as they are time-use descriptors, which may increase the multicollinearity risk. Likewise, each MX metric would carry partial and relative information on the pattern, and compositional transformation would also be interesting, although this approach has not been tested yet.
The usefulness of MX metrics with multivariate pattern analysis has not been investigated yet. However, since one of the limitations of MX metrics with linear regression is collinearity among the explanatory variables, multivariate pattern analysis could provide new insights by overcoming collinearity.
Multiple descriptors and machine learning
Machine learning (ML) describes a broad range of techniques to automate finding patterns in data with a focus on predictive ability (supervised ML) or data reduction (unsupervised ML). ML methods have been widely applied to derive accelerometer descriptors,40 yet rarely applied to study health associations.41 42 Different ML approaches have different strengths and limitations. In general, strengths include their usefulness for data-driven hypothesis generation, their capacity to handle multidimensional data, their ability to find non-linear patterns, and the possibility of training a model in one dataset and updating it in another. However, it can be difficult to interpret how results are obtained and their significance for public health guidelines. ML methods can also be data hungry and computationally intensive. Overfitting and sensitivity to (potentially unknown) biases in the training data are risks.
In some ways, multivariate pattern analyses and other dimension reduction methods can be considered ML methods. The Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis (TRIPOD) Initiative developed a set of recommendations for reporting of studies developing, validating, or updating ML-based prediction models for diagnostic or prognostic purposes.43 The TRIPOD statement should be considered when developing or applying ML-based prediction models in physical behaviour epidemiology.
Future research directions
The workshop in Granada, as well as the later meeting and the work developed in the following months by this author group, has initiated a discussion on analytical approaches and their usefulness for public health guidelines. Currently, 150 min/week of MVPA in adults and older adults, and 60 min/day of MVPA in children are recommended by different agencies.22 44 Recent guidelines removed the 10-minute bout requirement for MVPA in adults, and included the importance of replacing SB for PA.22 The Canadian PA guidelines are the first attempt to promote the combined effects of behaviours on health,45 although the evidence used was not based on the 24-hour paradigm (and so, appropriateness of statistical approach can be discussed).
We propose future research directions based on the research gaps identified, that is, the uncertainty regarding the accelerometer data descriptors to use and what analytical models are the most appropriate given the research question being addressed. The authors of this consensus article agree that investigations determining associations between physical behaviours and health should be extended to understand the interplay of physical behaviours (PA, SB and sleep) in their relationship with health. Measurement and processing capacity is increasing and offers an opportunity to provide further information on how different intensities and types of PA interact to improve health. At the same time, the focus on translating findings to meaningful information for interpretation in practice cannot be lost when using advanced analytical models. It is notable that most of the information presented comes from the PA and SB fields; thus, the relevance for the sleep research filed can be further discussed. The main implications for the analysis of accelerometer data proposed and agreed by the authors of this consensus manuscript are presented below.
Clear communication on the rationale for the use and limitations of each analytical approach in studies is important for a meaningful interpretation of the findings. Practical recommendations for this are provided in the Analytical approaches section of this document and a decision tree was developed (figure 1) to assist researchers’ decision-making.
Investigation of the associations of physical behaviours with health using different analytical approaches is encouraged. Ideally, physical behaviour epidemiology would draw consistent conclusions independently of analytical approach. To do so, clear reporting on the interpretation of findings derived from each analytical approach is crucial to understand ‘a priori’ inconsistencies across methods. Triangulation of results from different analytical approaches is currently the best solution to quantify associations of physical behaviours with health. Additionally, using the best-suited analytical approaches for a given research question is crucial (see figure 1).
Although little explored so far, ML-based approaches for diagnostic/prognostic purposes are worth implementing. We encourage transparent reporting of the resulting tools (TRIPOD initiative checklist).
Translating findings to meaningful information for guidelines should be a priority. Accurate reporting of study findings, interpretation and practical implications is highly encouraged.
How to conveniently adjust for confounders in multivariate pattern analysis should be investigated and its application extended to time-dependent outcomes (eg, survival analysis with mortality).
Further efforts are needed to translate functional data analysis and other advanced analytical approaches’ outputs into meaningful information for public health guidelines.
To evaluate whether the information gathered from the analytical approaches discussed herein can result in complementary information for public health guidelines. Such complementary information may result in more specific recommendations for certain health outcomes or populations, or even in their implementation at population level through movement sensors using evidence-based goals on PA intensity, duration, timing or type, among others.
This group agreed on several consensus points and research needs for physical behaviour epidemiology (see box 1 and figure 1). This consensus article will increase researchers’ understanding of different analytical approaches used in recent epidemiological studies of physical behaviours. This article and the decision tree provided aim to assist researchers in selecting analytical approaches based on their research questions and data. This will ultimately have an impact on the scientific evidence and, therefore, on future public health guidelines on physical behaviours.
Consensus points from the GRANADA report on analytical approaches to assess associations with accelerometer-determined physical behaviours (physical activity (PA), sedentary behaviour (SB) and sleep) in epidemiological studies
The study of the association between physical behaviours (ie, PA, SB and sleep) and health should move to a more thorough investigation of the interactions and codependencies between different behaviours (or PA intensities) and health. Several analytical approaches are provided in this consensus document, although none of them is free from limitations.
We recommend investigating more detailed PA intensities than the typically studied (ie, SB and moderate-to-vigorous PA). Examples include light PA of different intensities or the more fine-grained intensity bands as described in this document.
Public health guidelines on physical behaviours should acknowledge that behaviours are codependent and this may affect the guidelines as traditionally understood.
Further investigation in functional data analysis and machine learning is needed concerning the associations of physical behaviours with health.
There is not a gold standard able to test which analytical approach is the best for a given research question. Thus, we cannot make a strong recommendation on a single analytical approach. Instead, we provide some practical recommendations to select analytical approaches well suited for a given research question. Triangulation across findings from different analytical approaches is currently the best solution.
Patient consent for publication
Because of logistic and funding reasons, the researchers participating in the workshop and authors of this manuscript are mainly based in Europe and do not represent governing body or public health organisation. Thus, the authors of this consensus manuscript are a group of European researchers selected from the speaker and participant list of the GRANADA workshop: 'a focus on statistical methods to analyse accelerometer-measured physical activity'. Additionally, the manuscript was sent for review to a set of external experts in the field from different countries, including Europe, USA and Australia. Ten reviewers provided comments (more than 150 comments received) out of the 15 initially contacted. Therefore, it has an international dimension beyond Europe given the important input that the comments from the international review group had on the manuscript. The authors jointly accumulate reasonable expertise around methodology and application of methodology in physical behaviour research. This manuscript contains the authors’ opinions and it is not written under the name of any health association or governing body. This is part of a PhD Thesis conducted in the Biomedicine Doctoral Studies at the University of Granada, Spain.
Twitter @Jairohm8, @epiAgeing
Correction notice This article has been corrected since it published Online First. The co-corresponding author has been added.
Collaborators Alexander Burchartz (Institute for Sports and Sports Science, Karlsruhe Institute of Technology, Karlsruhe, Germany), Cain Clark (Centre for Sport, Exercise and Life Sciences, Coventry University, Coventry, UK), Paddy Dempsey (MRC Epidemiology Unit, Institute of Metabolic Science, University of Cambridge, Cambridge, UK; Physical Activity & Behavioural Epidemiology Laboratories, Baker Heart and Diabetes Institute, Melbourne, Australia), Aiden Doherty (Nuffield Department of Population Health, Big Data Institute, University of Oxford, Oxford, UK; NIHR Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, John Radcliffe Hospital, Oxford, UK), Ulf Ekelund (Department of Sports Medicine, Norwegian School of Sport Sciences, Oslo, Norway), Timothy Olds (Alliance for Research in Exercise, Nutrition and Activity, University of South Australia, Adelaide, Australia), Eric J Shiroma (Laboratory of Epidemiology and Population Science, National Institute on Aging, Baltimore, Maryland), Emmanuel Stamatakis (Charles Perkins Centre, Prevention Research Collaboration, Sydney School of Public Health, Faculty of Medicine and Health, University of Sydney, Sydney, NSW, Australia), Richard P Troiano (Division of Cancer Control and Population Sciences, National Cancer Institute, NIH, HHS, Rockville, Maryland, USA), Stewart Trost (Institute of Health and Biomedical Innovation at Queensland Centre for Children's Health Research, Queensland University of Technology, South Brisbane, Australia; Faculty of Health, School of Exercise and Nutrition Sciences, Queensland University of Technology, Brisbane, Australia) and Vadim Zipunnikov (Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University).
Contributors JHM and FBO conceived the idea, organised and led the workshop, the meeting, and the GRANADA consensus design and writing. All authors contributed to the idea, reviewed the manuscript and provided detailed suggestions for revisions.
Funding This study was conducted under the umbrella of the ActiveBrains and the SmarterMove projects supported by the MINECO/FEDER (DEP2013‐47540, DEP2016‐79512‐R, RYC‐2011‐09011) and the CoCA project supported by the European Union’s 2020 research and innovation programme (667302). JHM is supported by a grant from the Spanish Ministry of Education, Culture and Sport (FPU15/02645). AR is supported by the NIHR Leicester Biomedical Research Centre, and the Collaboration for leadership in Applied Health Research and Care (CLAHRC) East Midlands. SS is supported by the French National Research Agency (ANR-19-CE36-0004-01). RW is supported by a Medical Research Council Industrial Strategy Studentship (MR/S502509/1). Additional funding was obtained from the Andalusian Operational Programme supported with European Regional Development Funds (ERDF in English, FEDER in Spanish, project ref: B-CTS-355-UGR18), the University of Granada, Plan Propio de Investigación 2016, Excellence actions: Units of Excellence; Scientific Excellence Unit on Exercise and Health (UCEES), Junta de Andalucía, Consejería de Conocimiento, Investigación y Universidades and European Regional Development Funds (ref. SOMM17/6107/UGR). In addition, funding was provided by the SAMID III network, RETICS, funded by the PN I+D+I 2017‐2021 (Spain), ISCIII‐Sub‐Directorate General for Research Assessment and Promotion, the European Regional Development Fund (ERDF) (Ref. RD16/0022), the EXERNET Research Network on Exercise and Health in Special Populations (DEP2005‐00046/ACTI).
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.