Skip to main content

Day-to-day pattern of work and leisure time physical behaviours: are low socioeconomic status adults couch potatoes or work warriors?

Abstract

Background

Most studies on day-to-day patterns of physical behaviours (i.e. physical activities and sedentary behaviour) are based on adults with high socioeconomic status (SES) and without differentiating between work and leisure time. Thus, we aimed to characterise the day-to-day leisure time physical behaviours patterns among low SES adults and investigate the influence of work physical behaviours.

Methods

This cross-sectional study included 963 adults from low SES occupations (e.g. manufacturing, cleaning and transportation). The participants wore accelerometers for 1–7 days to measure physical behaviours during work and leisure time, expressed as time-use compositions consisting of time spent sedentary, standing or being active (walking, running, stair climbing, or cycling). Compositional multivariate multilevel models were used to regress daily leisure time-use composition against work time-use compositions. Interaction between weekday and (1) type of day, (i.e., work/non-work) and (2) the work time-use composition were tested. Compositional isotemporal substitution was used to interpret the estimates from the models.

Results

Each weekday, workers consistently spent most leisure time being sedentary and most work time standing. Leisure time physical behaviours were associated with type of day (p < 0.005, more sedentary on workdays vs. non-workdays), weekday (p < 0.005, more sedentary on Friday, Saturday and Sunday), standing work (p < 0.005, more sedentary and less standing and active leisure time on Sunday), and active work (p < 0.005, less sedentary and more standing and active leisure time on Sunday). Sedentary leisure time increased by 18 min, while standing and active leisure time decreased by 11 and 7 min, respectively, when 30 min were reallocated to standing at work on Sunday. Conversely, sedentary leisure time decreased by 25 min, and standing and active leisure time increased by 15 and 10 min, respectively, when 30 min were reallocated to active time at work on Sunday.

Conclusions

While low SES adults’ leisure time was mostly sedentary, their work time was predominantly standing. Work physical behaviours differently influenced day-to-day leisure time behaviours. Thus, public health initiatives aiming to change leisure time behaviours among low SES adults should consider the influence of work physical behaviours.

Peer Review reports

Background

Leisure time physical behaviours (i.e. physical activities (PA) and sedentary behaviour) are well-known determinants of non-communicable diseases and mortality [1,2,3,4]. Worldwide efforts are made to increase PA [5, 6]. Nevertheless, current PA promoting interventions and policies seem to fail in reaching those in most need of health-enhancing leisure time PA: individuals of lower socioeconomic status (SES). Current data on global leisure time PA levels show a persistent, steadily widening SES gap, favouring high SES groups [7,8,9]. Clearly, there is an urgent need to strengthen PA promoting strategies to narrow this gap and thus, decrease inequalities in health.

One strategy to increase PA levels among low SES adults could be to recommend health enhancing physical activities on non-workdays days (during the week or weekends) with more energy and free time available than on workdays. This pattern of higher PA, e.g. during non-working weekends, than on workdays has been found among high SES adults (termed “weekend warriors”) [10,11,12]. An opposite pattern, with more PA on weekdays than during the weekend has been found among low SES adults [13, 14]. However, none of the previous studies differentiated between work and leisure PA. As low SES adults are often in blue collar positions with high levels of PA as part of their job [15, 16], occupational PA is likely the driving factor for their high accumulation of PA on workdays.

The few studies assessing determinants of physical behaviour patterns over the week have primarily focused on individual factors, such as age [17], income [11] and education [13]. However, at best, individual-level factors explain 20–40% of the variance in PA levels [18]. Thus, research and policies on physical behaviours has increasingly adopted a broader approach which also considers environmental determinants, such as work factors [19, 20]. Understanding how work factors could affect daily leisure time physical behaviours is considered particularly important among low SES groups, as such factors are modifiable and amenable to change with interventions [21, 22]. Nevertheless, to our knowledge, no study has investigated how work factors, such as work physical behaviours, influence day-to-day leisure time physical behaviours. Such insight into domain-specific physical behaviours over the week and potentially modifiable determinants is needed to inform approaches aiming to increasing health-enhancing leisure time PA. As highlighted by the World Health Organization (WHO), this knowledge on domain-specific physical behaviour patterns is lacking, yet, urgently needed among less privileged adults [23].

When investigating the relationship between work and leisure time physical behaviours, it is essential to consider the potential co-dependency between these behaviours. That is, increasing time spent in one behaviour may take away time from another behaviour. Consequently, work and leisure time physical behaviours are likely to be co-dependent and collinear. These type of data are known as compositional data and they convey relative rather than absolute information [24]. Accordingly, assessment of patterns and associations between domain-specific physical behaviours can be made in terms of relative information [25, 26]. This is done by using compositional data analysis (CoDA), which is based on a log-ratio methodology. To our knowledge, no study has assessed the day-to-day pattern of domain specific physical behaviours while taking the co-dependency between time-use behaviours into account.

The aim of this study was to characterize the day-to-day pattern of work and leisure time physical behaviours among low SES adults. Moreover, we investigated the influence of day of the week and its interaction with type of day (i.e. workday vs. non-workday) and work physical behaviours on day-to-day leisure time physical behaviours using a CoDA approach.

Methods

This study used baseline data from two Danish studies with identical data procedure and collection: the Danish PHysical ACTivity cohort with Objective measurements (DPhacto) [27] and the New Method for Objective Measurements of Physical Activity in Daily Living (NOMAD) study [28]. The study population consisted of low SES workers recruited from Danish workplaces within cleaning, transportation, manufacturing, construction, road maintenance, garbage disposal, assembly, mobile plant operator, and health care. Workers were eligible if they were employed in one of the mentioned sectors for at least 20 h/week; between 18 and 65 years old; and had given voluntary consent to participate. Workers were excluded if they had band-aid allergy, fever on the date of data collection or were pregnant.

Data collection

Details of data collection and related instruments used for the DPhacto and NOMAD studies have been described previously [27, 28]. In short, data collection included questionnaires, health checks, and accelerometer-based measurements. Workers eligible for participation were invited to complete a questionnaire and to participate in a health check, which consisted of anthropometric measurements and a physical health examination. Moreover, participants were asked to wear accelerometers for a minimum of two consecutive workdays and to complete a daily diary reporting time at work; time going to bed at night and getting up in the morning; and non-wear time.

Accelerometer measurements of physical Behaviours

Physical behaviours at work and leisure time was assessed using data from one tri-axial ActiGraph GT3X+ accelerometer (Actigraph, Pensacola, FL, USA). The accelerometer was placed on the right thigh using double-sided adhesive tape (3 M, Hair-Set, St. Paul, MN, USA) and Fixomull (Fixomull BSN medical GmbH, Hamburg, Germany). Accelerometer data were downloaded using the Actilife Software version 5.5 (Actigraph, Pensacola, FL, USA) [29] and analysed using the custom-made MATLAB program Acti4 (The National Research Centre for the Working Environment, Copenhagen, Denmark) [30]. In brief, the Acti4 program derive instantaneous average parameters (i.e. mean and standard deviation of acceleration, inclination) from thigh-worn accelerometer signal using overlapping 2-s intervals. Based on these parameters and a rule-based decision tree, physical behaviours (i.e., cycling, stair climbing, running, walking, standing, sitting and lying) were classified [19]. The Acti4 program has been shown to separate physical behaviours with high sensitivity and specificity under both semi-standardized [30] and non-standardized conditions [31].

Day of the week, daily work hours, leisure time and time-in-bed were defined from the participants’ daily diary. Specifically, leisure time was defined as waking time not at work. Non-workdays were defined as days where the participants had not reported a working period. Only workers having at least 1 day of valid accelerometer measurements of both work and leisure time periods were included in this study. A valid day was defined as having at least 4 h of accelerometer-derived work and leisure time or what corresponded to at least 75% of the individual’s average work and leisure time. To decrease the risk of reverse causality between work periods and leisure time, we only considered leisure time following work. Moreover, time-in-bed was not considered in this study.

Figure 1 shows the flowchart of the study population. A total of 1207 eligible workers answered the questionnaire and/or participated in the physical health check. Of these workers, 40 were excluded as they were department leaders or students, on holiday, pregnant or because they did not want to participate. A total of 204 workers were excluded from the study as they did not valid have leisure time accelerometer measurements on at least one weekday. Therefore, a total of 963 workers were included in the study.

Fig. 1
figure 1

Flow chart of the study population

Covariates

Sex and age were determined from each worker’s unique Danish civil registration number. BMI (Body Mass Index) was calculated as weight (kg) divided by height (m) squared (kg/m2). Information on smoking-status was obtained by the question: “Do you smoke?” with four response categories: daily smoking, occasionally smoking, formerly smoked, and never smoked. The variable was dichotomized into smokers and non-smokers (including former smokers). Information on shift work was obtained using the question: “At what time(s) of the day do you usually work in your main occupation?” with three response categories: fixed day work, night/varying work hours with night, and other. The variable was dichotomized into workers with and without fixed day work. Work duration was calculated as the log of total accelerometer-derived work time [32].

Statistical analysis

Compositional descriptive

Time use of daily work and leisure time behaviours was treated as two compositions of activities performed within a 24-h day. Work and leisure time were defined as a 3-part composition, both consisting of time spent on sedentary (i.e. sitting or lying), standing and active (i.e. walking, running, stair climbing or cycling).

Compositional means were used to describe the day-to-day pattern of work and leisure time physical behaviours [24, 33]. These were obtained by calculating the geometric mean of each physical behaviour of the respective compositions and then normalising the geometric means to the workers’ average accelerometer-derived daily work and leisure time (i.e. 450 min and 450 min, respectively). On non-workdays, the leisure time composition consisted of daily waking time, normalised to the workers’ average accelerometer-derived daily time spent awake (i.e. 960 min).

Calculation of pivot isometric log-ratios (ilrs) and model development

Daily work and leisure time-use compositions were expressed using pivot isometric log-ratio (ilr) coordinates [34]. The first pivot-coordinate was calculated as the normalised log-ratio of the first compositional part (i.e. behaviour), relative to the geometric mean of the remaining parts within the work and leisure time composition, respectively. The work and leisure time behaviours were sequentially rearranged to place each behaviour in the first position, where after the corresponding ilr-coordinate sets were computed. This way, the relative importance of each behaviour was sequentially represented in the first ilr-coordinate (ilr1) and used in the regression analysis. A detailed description of how the pivot-coordinates were calculated and model development is provided in Additional file 1.

The analysis was performed using two multivariate multilevel models. In both models, the outcome variables were the ilr-coordinates expressing the leisure time-use composition.

In Model 1, we investigated if leisure time physical behaviours differed between each day of the week (e.g. Monday, Tuesday, Wednesday, etc.) and between workday and non-workdays. Thus, day of the week and an interaction between day of the week and type of day (reference = non-workday) were entered as predictors in Model 1. Model 1 was fitted three times. This was done to isolate the association with one of the leisure time behaviours in relation to the others in the first ilr-coordinate (denoted by ilr1).

In Model 2, only workdays were considered as we investigated if work behaviours influenced day-to-day leisure time behaviours. The following predictors were entered in Model 2: the work time-use composition (i.e. work time sent standing, active and sedentary, expressed as ilr-coordinates) and an interaction term between day of the week and the work time-use composition. Model 2 was fitted six times to investigate the association between each part of the leisure time and work compositions, respectively. Of note, only results of the associations between the relative work time spent active and standing (as a proxies of physical work demands) and leisure time physical behaviours are shown. Results on the association between relative work time spent sedentary and day-to-day leisure time physical behaviours are shown in Additional file 2.

Both Model 1 and 2 were adjusted for the following covariates (reference in parenthesis for categorical variables): sex (men), smoking-status (smoker), BMI, and age. Model 2 was further adjusted for work duration. These covariates were chosen as potential confounders based on theoretical assumptions concerning their possible influence on day-to-day pattern of leisure time and work behaviours and work duration [13, 17]. In all models, Monday was selected as the reference category when entering the “type of the day” variable, as this is considered as the first day of the week in Denmark, where the International Standard ISO 8601 is followed.

Compositional isotemporal substitution analysis

Compositional isotemporal substitution analysis was used to provide meaningful interpretation of the expected change (in min/day) in leisure time-use compositions when time was reallocated between behaviours during work on workdays. This was done using the multivariate regression Model 2, stratified on each workday of the week. First, a “reference” leisure time-use composition (average daily leisure time spent sedentary, standing and active) was estimated for the workers’ mean work time-use composition (average min of work time spent sedentary, standing and active for that particular day). Second, new work time-use compositions were calculated where time (15, 30 and 45 min) had been reallocated between behaviours. This enabled us to express effect sizes as expected changes in leisure time behaviours in min/day. Note that results are only shown for workdays where the work behaviours were significantly associated with the leisure time behaviours. A detailed description of this method based on ilr linear regression with non-compositional and compositional outcomes can be found in Dumuid et al. [35] and Lund Rasmussen et al. [36].

All analyses were performed in R version 1.1.3 [37], using the compositions [38] and MCMCglmm [39] packages. We used the MCMCglmm package to conduct the multivariate multilevel analysis, following the guide provided by Baldwin et al. [40], by which a Bayesian approach with uninformative priors were used. The assumptions of normality and homoscedasticity of the residuals were assessed for all models by visual inspection of residuals versus predicted values and quantile-quantile plots.

Sensitivity analysis

A sensitivity analysis was conducted in which only workers with at least 2 days of measurements were included (N = 831). Results are shown in Additional file 3.

Assessment of potential selection bias

To identify potential selection bias, we compared the characteristics of the blue-collar workers included and excluded from the study. Differences between groups were investigated by calculating means and standard deviations or frequencies and percentages. Group differences were tested using t-test and Chi-squared statistics and a 5% significance threshold. Results are shown in Additional file 4.

Results

Study population characteristics

Table 1 shows the basic characteristics of the study population. Mean age was 44.9 (SD = 10.0) years, mean BMI was 27.2 kg/m2 (SD = 4.9), 45% were women and the majority were working within manufacturing (59%). On average, the participants had 4.1 (SD = 1.3) days with valid accelerometer measurements. Across all weekdays, the mean accelerometer-derived valid leisure and work hours were 7.4 (SD = 2.6) and 7.6 (SD = 2.1).

Table 1 Characteristics of the study population (n = 963)

Compositional descriptive of day-to-day patterns in work and leisure time physical behaviours

Table 2 shows the compositional means of the day-to-day work and leisure time-use compositions on workdays and non-workdays. On workdays, the workers were predominantly sedentary during leisure time throughout the week. In contrast, most work time was spent standing on all working weekdays. On non-workdays, most time was spent sedentary throughout the week.

Table 2 Compositional means of leisure time and work physical behaviours on workdays and non-workdays

Results of multilevel models

Leisure time physical behaviours on workdays vs. non-workdays (model 1)

Table 3 shows the results of model 1, investigating if the workers’ day-to-day leisure time physical behaviours differed between workdays and non-workdays. On Friday, Saturday and Sunday, the workers spent significantly more leisure time being sedentary on workdays compared with non-workdays (Table 3, ilr1(Sedentary_leisure) βinteraction = 0.26, 95% CI = (0.07; 0.44), ilr1(Sedentary_leisure) βinteraction = 0.34, 95% CI = (0.13; 0.55), and ilr1(Sedentary_leisure) βinteraction = 0.27, 95% CI = (0.04; 0.46), respectively). Moreover, the workers spent significantly less leisure time standing on Friday and Saturday on workdays compared with non-workdays (Table 3, ilr1(Standing_leisure) βinteraction = − 0.16, 95% CI = (− 0.32; − 0.02), and ilr1(Standing_leisure) βinteraction = − 0.21, 95% CI = (− 0.34; − 0.07), respectively).

Table 3 Association between day-to-day leisure time physical behaviours and weekday and workday

These results suggested that the workers day-to-day leisure time behaviours differed between workdays and non-workdays and thus, further supported the investigation of whether day-to-day work physical behaviours caused this difference.

Standing work time and leisure time behaviours (model 2)

On all working weekdays, we observed a trend of relative standing work time to be positively associated with relative sedentary leisure time and negatively associated with relative standing and active leisure time. However, the associations were only statistically significant on Tuesday (Table 4; ilr1(Standing_leisure) βinteraction = − 0.17, 95% CI = (− 0.33; − 0.02)) and Sunday (ilr1(Sedentary_leisure) βinteraction = 0.93, 95% CI = (0.23; 1.60), ilr1(Standing_leisure) βinteraction = − 0.50, 95% CI = (− 0.93; − 0.12) and ilr1(Active_leisure) βinteraction = − 0.45, 95% CI = (− 0.91; − 0.01)). Specifically, the positive ilr1 beta-coefficient of 0.93 for Sedentary_leisure indicated that more standing work was associated with more sedentary leisure time on Sunday. In contrast, the negative ilr1 beta-coefficient of − 0.50 and − 0.45 for Standing_leisure and Active_leisure, respectively, showed that, on Sunday, more standing work was associated with less standing and active leisure time.

Table 4 Association between day-to-day leisure time physical behaviours, weekday, and relative standing work time

Active work time and leisure time behaviours (model 2)

For relative active work time, we observed a negative association with relative sedentary leisure time and positive associations with relative standing and active leisure time throughout the week. These associations were only statistically significant on Tuesday (Table 5; ilr1(Standing_leisure) βinteraction = 0.20, 95% CI = (0.05; 0.42)), Friday (ilr1(Standing_leisure) βinteraction = 0.16, 95% CI = (0.01; 0.34)), and Sunday (ilr1(Sedentary_leisure) βinteraction = − 1.02, 95% CI = (− 1.67; − 0.33), ilr1(Standing_leisure) βinteraction = 0.56, 95% CI = (0.14; 1.04), and ilr1(Active_leisure) βinteraction = 0.47, 95% CI = (0.01; 0.97)).

Table 5 Association between day-to-day leisure time physical behaviours, weekday and relative active work time

Results of compositional isotemporal substitutions

Tables 6 and 7 show the results of the compositional isotemporal substitution. Note that only results for reallocations on Tuesday and Sunday are shown, as we observed the strongest relationship between work and leisure time behaviours on these working weekdays.

Table 6 Expected difference in leisure time behaviours following reallocation between work physical behaviours on Tuesday
Table 7 Expected difference in leisure time behaviours following reallocation between work physical behaviours on Sunday

On Tuesday (Table 6), reallocating 30 min work time to standing from the remaining work behaviours was associated with 1 min less standing leisure time. Reallocating 30 min to active work time was associated with 3 min less sedentary leisure time and 1 min more standing leisure time.

On Sunday (Table 7), reallocating 30 min work time to standing from the remaining work behaviours was associated with 18 min more sedentary time, 11 min less standing time, and 7 min less active time during leisure. Reallocating 30 min to active work time from the remaining work behaviours was associated with 25 min less sedentary time, a 15 min increase in standing time and a 10 min increase active time during leisure.

Sensitivity analysis

Results of sensitivity analysis among workers with at least 2 days of valid accelerometer data corresponded to those from the primary analyses (results shown in Additional file 3).

Assessment of potential selection bias

Comparison of the characteristics of blue-collar workers excluded and included in the current study is shown in Additional file 4. Blue-collar workers excluded from the study sample (n = 244) had, on average, a lower seniority (mean = 11.6 years, SD = 10.2) and aerobic capacity (mean = 30.0, SD = 7.0) and a higher proportion was women and cleaners (24%) or transporters (9%) compared to those included (n = 963).

Discussion

In this study, we investigated the day-to-day pattern of leisure time physical behaviours among blue-collar workers on workdays and non-workdays. Furthermore, we assessed the association between day-to-day work physical behaviours and leisure time physical behaviours. The workers were primarily sedentary during leisure throughout the week on both workdays and non-workdays. Moreover, the workers were more sedentary at leisure during the weekend on workdays compared to non-workdays. Regarding the association between work and leisure time physical behaviours, standing work time was positively associated with sedentary leisure time and negatively associated with standing and active leisure time. The opposite direction was found for the association between active work time and leisure time physical behaviours.

The overall finding of our study is that over the course of a week, low SES adults were predominantly sedentary during leisure, while most work time was spent standing or active. To our knowledge, this is the first study to investigate the domain-specific pattern of day-to-day physical behaviours among low SES adults, while taking the inter-dependency between daily time-use into account using CoDA. Nevertheless, our results correspond to those reported in reviews finding leisure time PA to be less prevalent and work PA to be more prevalent among low SES adults compared with high SES adults when considering weekly averages [15, 41]. Moreover, accelerometer-based studies have consistently found low SES adults to be more active on workdays compared to non-workdays. One study found Swiss workers with manual jobs (e.g. craftsmen, machine operators and labourers) to be more active on workdays than non-workdays [42]. Another study reported adults with low educational level as less likely to be active during the weekend (i.e. non-workdays) and more likely to be active during the weekdays (i.e. workdays) compared with those with high educational level [13]. Finnish low-level occupational groups (e.g. cleaners, plumbers and construction workers) and Australian blue-collar workers have been found to be more sedentary and take less steps on weekends (i.e. non-workdays) compared with weekdays (i.e. workdays) [14, 43]. Taken together, these findings suggest that low SES adults derive the greatest proportion of day-to-day PA from work activities. Consequently, when only considering leisure time physical behaviours, low SES adults could erroneously be perceived to be “couch potatoes”. However, taking work physical behaviours into consideration, this population group might also be “work warriors”.

We further investigated whether the workers’ work physical behaviours was associated with their high amount of sedentary leisure time throughout the week. For example, on Sunday, reallocating 30 min to standing work time was associated with an 18 min increase in sedentary leisure time, an 11 min decrease in standing leisure time, and a 7 min decrease in active leisure time. This finding could be explained by the fact that standing for prolonged, uninterrupted periods increase blood pooling in the legs which can cause swelling, pain and muscle fatigue in the lower extremities [44, 45]. Consequently, workers with much standing work time could perceive to have an increased need to compensate with more sedentary time and less physical activities when coming home from work [46]. We have previously assessed the association between percentage work hours spent standing and percentage leisure time spent sedentary over consecutive workdays within the same population [47]. However, contrary to the present study, no relationship between occupational standing and sedentary leisure time was observed. One explanation for the discrepancy between the two studies could be that we did not consider the co-dependency between work and leisure time physical behaviours as done in the current study.

Counter to standing work time, active work was associated with less sedentary and more standing and active leisure time. For example, reallocating 30 min to active work time was associated with a 25 min decrease in sedentary leisure time, a 15 min increase in standing leisure time, and a 10 min increase in active leisure time on Sunday. This result is in line with one study finding a positive relationship between accelerometer-derived time spent on moderate-to-vigorous physical activity at work and leisure among 233 adults [48]. In contrast, another study based on accelerometer measurement from 112 adults found no difference in leisure physical activities between those with low and high work activity levels [49]. Thus, the relationship between work and leisure time physical behaviours among low SES adults remains unclear. We encourage future research on this topic based on technical measurements and considering the inter-dependency between time spent in daily behaviours.

Practical implications

Most physical activity interventions focus on individual factors such as motivation and self-efficacy [50]. This could be a result of cognitive social theories dominating the behavioural research on physical activity for decades. Accordingly, this research field has been shaped by the assumption that being active or inactive are deliberate choices. However, motivating individuals to be active without considering barriers is likely to be ineffective, particularly for low SEP adults. In fact, interventions that solely rely on individual agency to increase physical activity levels risk widening social inequalities due to higher reach and effectiveness among high SEP adults [51]. Instead, interventions should both address modifiable barriers, such as work factors, that limit opportunities for being active, while protecting and enhancing factors that enable and encourage this behaviour [52].

Based on the findings of this study, standing at work could be a barrier for low SES adults to be physically active during leisure time. Thus, we encourage future intervention studies to investigate the effectiveness of reducing this physical work demand on leisure time behaviours. For workers for whom work time spent standing cannot be altered, it will be important to investigate the health effect and feasibility of different combinations of leisure time spent on sedentary (to promote post-work recovery) and health-enhancing activities over the course of a week.

Although the estimated changes in leisure time physical behaviours might appear small (e.g. a 7 min decrease in active leisure time when increasing standing work time), the found changes in leisure time should be considered in relation to the overall leisure time physical behaviours. As this group of adults was predominantly sedentary during leisure time, any additional trade-off between sedentary and active leisure time could have important long-term health implications. Moreover, it should be noted that this study was conducted in a high-income country. Considering that work has been found to be a greater contributor to moderate-to-vigorous physical activity levels in lower compared with higher income countries [53], we would expect our findings to be more prominent among lower income countries.

Strengths and limitations

This study was based on technical measurements, which enabled accurate information on daily time spent in physical behaviours at work and during leisure time. Particularly, the use of the Acti4 program to measure physical behaviours was a strength of the study as it enabled distinction between specific physical behaviours, which fall within the same intensity, such as standing and walking, with high sensitivity and specificity [30, 31]. The use of multivariate multilevel models was another strength by allowing analyses which included the repeated measurements as well as multiple outcomes for each participant.

Given the close link between SES and work and leisure time physical behaviours [15, 41, 54], there is always a risk of socioeconomic confounding if SES is not appropriately accounted for. While there is no single best indicator of SES, occupational class is well-acknowledged to reflect social standing among working adults [55, 56]. Accordingly, the use of a study population consisting of adults who were homogenous in occupational class was a strength as it limited the possibility of socioeconomic confounding. Nevertheless, the lack of information on other measures of socioeconomic status such as education and income was a limitation of our study.

This study was based on blue-collar workers recruited through low profit workplaces from a range of different sectors and locations in Denmark to ensure a study population representative of typical blue-collar workers [57, 58]. Accordingly, the characteristics of the study population in the DPhacto study has been considered as representative for the target population [57]. Nevertheless, we observed that the blue-collar workers excluded in the study had a lower aerobic capacity compared to those included. A high aerobic capacity can enhance the workers’ ability to sustain their physical work, and thereby decrease levels of post-work fatigued [59]. Accordingly, the found differences could have underestimated the relationship between work and leisure time physical behaviours.

The cross-sectional design of this study hinders causal inference. The workers included in the NOMAD and DPhacto studies were asked to wear the accelerometers for at least two workdays. This was to enhance the feasibility of the workers to participate in the study but might have compromised the representativeness of the measurements. Most of the workers had the accelerometers mounted on Monday or Tuesday. Thus, it is possible that the workers behaved differently these days because they wore accelerometers. Moreover, participants with heavy physical work or a very active leisure time could have experienced issues with wearing the accelerometer, causing them to take off the devices at an earlier stage than less active participants. However, the sensitivity analysis excluding those workers with only 1 day of valid accelerometer data revealed results similar to that of the primary analysis.

We lacked information about the context in which day-to-day physical behaviours were performed. Thus, we are unable to state if active leisure time spent was spent on active transport; health-promoting planned activities such as sports; or on domestic work. The adults included in this study spent less than 1% of daily work and leisure time spent in activities which could be considered of higher intensities (i.e. running, stair climbing, and biking). Thus, we decided to combine these activities with walking, which could have attenuated some of the tested associations. Finally, we did not include time in bed at night in the analysis, which we encourage future studies on day-to-day patterns of daily physical behaviours to consider.

Conclusions

We found low SES adults to be primarily sedentary during leisure time, while work time was spent mostly standing or active. Thus, this group of adults could be characterised as ‘work warriors’. This highlights the need of differentiating between work and leisure time when assessing daily behavioural patterns among low SES adults. Standing work time was associated with a more sedentary leisure time, whereas the opposite was found for the association between active work time and leisure time physical behaviours. Accordingly, we recommend strategies for increasing health-enhancing leisure time PA among low SES adults to consider the potential influence of work physical behaviours.

Availability of data and materials

The datasets supporting the conclusions of this article are available at the Danish National Archives, https://www.sa.dk/en/k/about-us.

Abbreviations

DPhacto:

Danish PHysical ACTivity cohort with Objective measurements

NOMAD:

New Method for Objective Measurements of Physical Activity in Daily Living

PA:

physical activity

SES:

socioeconomic status

WHO:

World Health Organization

References

  1. Sofi F, Capalbo A, Marcucci R, Gori AM, Fedi S, Macchi C, et al. Leisure time but not occupational physical activity significantly affects cardiovascular risk factors in an adult population. Eur J Clin Investig. 2007;37(12):947–53. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-2362.2007.01884.x.

    Article  CAS  Google Scholar 

  2. Richard A, Martin B, Wanner M, Eichholzer M, Rohrmann S. Effects of leisure-time and occupational physical activity on total mortality risk in NHANES III according to sex, ethnicity, central obesity, and age. J Phys Act Health. 2015;12(2):184–92. https://0-doi-org.brum.beds.ac.uk/10.1123/jpah.2013-0198.

    Article  PubMed  Google Scholar 

  3. Smith AD, Crippa A, Woodcock J, Brage S. Physical activity and incident type 2 diabetes mellitus: a systematic review and dose-response meta-analysis of prospective cohort studies. Diabetologia. 2016;59(12):2527–45. https://0-doi-org.brum.beds.ac.uk/10.1007/s00125-016-4079-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Stamatakis E, Rogers K, Ding D, Berrigan D, Chau J, Hamer M, et al. All-cause mortality effects of replacing sedentary time with physical activity and sleeping using an isotemporal substitution model: a prospective study of 201,129 mid-aged and older adults. Int J Behav Nutr Phys Act. 2015;12(1):121. https://0-doi-org.brum.beds.ac.uk/10.1186/s12966-015-0280-7.

    Article  PubMed  PubMed Central  Google Scholar 

  5. American Center for Disease Control. 2008 Physical activity guidelines for Americans. 2008.

    Google Scholar 

  6. World Health Organization. More active people for a healthier world: global action plan on physical activity 2018-2030. 2018.

    Google Scholar 

  7. Scholes S, Bann D. Education-related disparities in reported physical activity during leisure-time, active transportation, and work among US adults: repeated cross-sectional analysis from the National Health and nutrition examination surveys, 2007 to 2016. BMC Public Health. 2018;18(1):926. https://0-doi-org.brum.beds.ac.uk/10.1186/s12889-018-5857-z.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Szilcz M, Mosquera PA, Sebastián MS, Gustafsson PE. Time trends in absolute and relative socioeconomic inequalities in leisure time physical inactivity in northern Sweden. Scand J Public Health. 2018;46(1):112–23. https://0-doi-org.brum.beds.ac.uk/10.1177/1403494817713123.

    Article  PubMed  Google Scholar 

  9. Van Cauwenberg J, Loyen A, Lakerveld J, Cardon G, De Craemer M, Gheysen F, et al. Differential influences of population densification and economic growth on Europeans’ physical activity and sitting time. Cities. 2018;82:141–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cities.2018.07.006.

    Article  Google Scholar 

  10. Evenson KR, Wen F, Metzger JS, Herring AH. Physical activity and sedentary behavior patterns using accelerometry from a national sample of United States adults. Int J Behav Nutr Phys Act. 2015;12(1):20. https://0-doi-org.brum.beds.ac.uk/10.1186/s12966-015-0183-7.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Shuval K, Li Q, Gabriel KP, Tchernis R. Income, physical activity, sedentary behavior, and the “weekend warrior” among U.S. adults. Prev Med. 2017;103:91–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ypmed.2017.07.033.

    Article  PubMed  Google Scholar 

  12. Metzger JS, Catellier DJ, Evenson KR, Treuth MS, Rosamond WD, Siega-Riz AM. Patterns of objectively measured physical activity in the United States. Med Sci Sports Exerc. 2008;40(4):630–8. https://0-doi-org.brum.beds.ac.uk/10.1249/MSS.0b013e3181620ebc.

    Article  PubMed  Google Scholar 

  13. Gubelmann C, Vollenweider P, Marques-Vidal P. Of weekend warriors and couch potatoes: socio-economic determinants of physical activity in Swiss middle-aged adults. Prev Med. 2017;105:350–5. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ypmed.2017.10.016.

    Article  PubMed  Google Scholar 

  14. Pulakka A, Stenholm S, Bosma H, Schaper NC, Savelberg HHCM, Stehouwer CDA, et al. Association between employment status and objectively measured physical activity and sedentary behavior-the Maastricht study. J Occup Environ Med. 2018;60(4):309–15. https://0-doi-org.brum.beds.ac.uk/10.1097/JOM.0000000000001254.

    Article  PubMed  Google Scholar 

  15. Beenackers MA, Kamphuis CB, Giskes K, Brug J, Kunst AE, Burdorf A, et al. Socioeconomic inequalities in occupational, leisure-time, and transport related physical activity among European adults: a systematic review. Int J Behav Nutr Phys Act. 2012;9(1):116. https://0-doi-org.brum.beds.ac.uk/10.1186/1479-5868-9-116.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Cusatis R, Garbarski D. Different domains of physical activity: the role of leisure, housework/care work, and paid work in socioeconomic differences in reported physical activity. SSM - Population Health. 2019;7:100387. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ssmph.2019.100387.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Jones SA, Wen F, Herring AH, Evenson KR. Correlates of US adult physical activity and sedentary behavior patterns. J Sci Med Sport. 2016;19(12):1020–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jsams.2016.03.009.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Spence JC, Lee RE. Toward a comprehensive model of physical activity. Psychol Sport Exerc. 2003;4(1):7–24. https://0-doi-org.brum.beds.ac.uk/10.1016/S1469-0292(02)00014-6.

    Article  Google Scholar 

  19. Rhodes RE, Mark RS, Temmel CP. Adult sedentary behavior: a systematic review. Am J Prev Med. 2012;42(3):e3–28. https://0-doi-org.brum.beds.ac.uk/10.1016/j.amepre.2011.10.020.

    Article  PubMed  Google Scholar 

  20. Kohl HW, Craig CL, Lambert EV, Inoue S, Alkandari JR, Leetongin G, et al. The pandemic of physical inactivity: global action for public health. Lancet. 2012;380(9838):294–305. https://0-doi-org.brum.beds.ac.uk/10.1016/S0140-6736(12)60898-8.

    Article  PubMed  Google Scholar 

  21. Owen N, Leslie E, Salmon J, Fotheringham MJ. Environmental determinants of physical activity and sedentary behavior. Exerc Sport Sci Rev. 2000;28(4):153–8.

    CAS  PubMed  Google Scholar 

  22. Kirk MA, Rhodes RE. Occupation correlates of adults’ participation in leisure-time physical activity: a systematic review. Am J Prev Med. 2011;40(4):476–85. https://0-doi-org.brum.beds.ac.uk/10.1016/j.amepre.2010.12.015.

    Article  PubMed  Google Scholar 

  23. Dempsey PC, Biddle SJH, Buman MP, Chastin S, Ekelund U, Friedenreich CM, et al. New global guidelines on sedentary behaviour and health for adults: broadening the behavioural targets. Int J Behav Nutr Phys Act. 2020;17:1–12.

    Article  Google Scholar 

  24. Aitchison J. The statistical analysis of compositional data. London: Chapman and Hall; 1986. https://0-doi-org.brum.beds.ac.uk/10.1007/978-94-009-4109-0.

    Book  Google Scholar 

  25. Gupta N, Rasmussen CL, Holtermann A, Mathiassen SE. Time-based data in occupational studies: the whys, the Hows, and some remaining challenges in compositional data analysis (CoDA). Ann Work Expo Health. 2020;64(8):778–85. https://0-doi-org.brum.beds.ac.uk/10.1093/annweh/wxaa056.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Dumuid D, Pedišić Ž, Palarea-Albaladejo J, Martín-Fernández JA, Hron K, Olds T. Compositional data analysis in time-use epidemiology: what, why, how. Int J Environ Res Public Health. 2020;17(7). https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17072220.

  27. Jørgensen MB, Korshøj M, Lagersted-Olsen J, Villumsen M, Mortensen OS, Skotte J, et al. Physical activities at work and risk of musculoskeletal pain and its consequences: protocol for a study with objective field measures among blue-collar workers. BMC Musculoskelet Disord. 2013;14:1.

    Article  Google Scholar 

  28. Gupta N, Jensen BS, Søgaard K, Carneiro IG, Christiansen CS, Hanisch C, et al. Face validity of the single work ability item: comparison with objectively measured heart rate reserve over several days. Int J Environ Res Public Health. 2014;11(5):5333–48. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110505333.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Software | ActiGraph [Internet]. Available from: http://actigraphcorp.com/support/software/. Cited 2016 Dec 20.

  30. Skotte J, Korshøj M, Kristiansen J, Hanisch C, Holtermann A. Detection of physical activity types using triaxial accelerometers. J Phys Act Health. 2014;11(1):76–84. https://0-doi-org.brum.beds.ac.uk/10.1123/jpah.2011-0347.

    Article  PubMed  Google Scholar 

  31. Stemland I, Ingebrigtsen J, Christiansen CS, Jensen BR, Hanisch C, Skotte J, et al. Validity of the Acti4 method for detection of physical activity types in free-living settings: comparison with video analysis. Ergonomics. 2015;58(6):953–65. https://0-doi-org.brum.beds.ac.uk/10.1080/00140139.2014.998724.

    Article  PubMed  Google Scholar 

  32. Pawlowsky-Glahn V, Egozcue JJ, Lovell D. Tools for compositional data with a total. Stat Model. 2015;15(2):175–90. https://0-doi-org.brum.beds.ac.uk/10.1177/1471082X14535526.

    Article  Google Scholar 

  33. Pawlowsky-Glahn V, Egozcue JJ, Tolosana-Delgado R. Modeling and analysis of compositional data. Chichester: Wiley; 2015.

  34. Hron K, Filzmoser P, de Caritat P, Fišerová E, Gardlo A. Weighted pivot coordinates for compositional data and their application to geochemical mapping. Math Geosci. 2017;49(6):797–814. https://0-doi-org.brum.beds.ac.uk/10.1007/s11004-017-9684-z.

    Article  CAS  Google Scholar 

  35. Dumuid D, Pedišić Ž, Stanford TE, Martín-Fernández J-A, Hron K, Maher CA, et al. The compositional isotemporal substitution model: a method for estimating changes in a health outcome for reallocation of time between sleep, physical activity and sedentary behaviour. Stat Methods Med Res. 2019;28(3):846–57.

  36. Lund Rasmussen C, Palarea-Albaladejo J, Bauman A, Gupta N, Nabe-Nielsen K, Jørgensen MB, et al. Does physically demanding work hinder a physically active lifestyle in low socioeconomic workers? A Compositional Data Analysis Based on Accelerometer Data. Int J Environ Res Public Health. 2018;15:1306.

    Article  Google Scholar 

  37. R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017. Available from: https://www.R-project.org/

    Google Scholar 

  38. Boogaart K, Tolosana R, Bren M. compositions: Compositional Data Analysis. R (version ≥ 220); 2014.

    Google Scholar 

  39. Hadfield JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J Stat Softw. 2010;33:1–22.

    Article  Google Scholar 

  40. Baldwin SA, Imel ZE, Braithwaite SR, Atkins DC. Analyzing multiple outcomes in clinical research using multivariate multilevel models. J Consult Clin Psychol. 2014;82(5):920–30. https://0-doi-org.brum.beds.ac.uk/10.1037/a0035628.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Stalsberg R, Pedersen A. Are differences in physical activity across socioeconomic groups associated with choice of physical activity variables to report? Int J Behav Nutr Phys Act. 2018;15:922.

    Google Scholar 

  42. Brighenti-Zogg S, Mundwiler J, Schüpbach U, Dieterle T, Wolfer DP, Leuppi JD, et al. Physical workload and work capacity across occupational groups. PLoS One. 2016;11(5):e0154073. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0154073.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Miller R, Brown W. Steps and sitting in a working population. Int J Behav Med. 2004;11(4):219–24. https://0-doi-org.brum.beds.ac.uk/10.1207/s15327558ijbm1104_5.

    Article  PubMed  Google Scholar 

  44. Balasubramanian V, Adalarasu K, Regulapati R. Comparing dynamic and stationary standing postures in an assembly task. Int J Ind Ergon. 2009;39(5):649–54. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ergon.2008.10.017.

    Article  Google Scholar 

  45. Bahk JW, Kim H, Jung-Choi K, Jung M-C, Lee I. Relationship between prolonged standing and symptoms of varicose veins and nocturnal leg cramps among women and men. Ergonomics. 2012;55(2):133–9. https://0-doi-org.brum.beds.ac.uk/10.1080/00140139.2011.582957.

    Article  PubMed  Google Scholar 

  46. Bláfoss R, Micheletti JK, Sundstrup E, Jakobsen MD, Bay H, Andersen LL. Is fatigue after work a barrier for leisure-time physical activity? Cross-sectional study among 10,000 adults from the general working population. Scand J Public Health. 2019;47(3):383–91.

  47. Lund Rasmussen C, Nabe-Nielsen K, Jørgensen MB, Holtermann A. The association between occupational standing and sedentary leisure time over consecutive workdays among blue-collar workers in manual jobs. Int Arch Occup Environ Health. 2018;92:481–90.

  48. JaKa MM, Haapala JL, Wolfson J, French SA. Describing the relationship between occupational and non-occupational physical activity using objective measurement. Prev Med Rep. 2015;2:213–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.pmedr.2015.03.003.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Tigbe WW, Lean MEJ, Granat MH. A physically active occupation does not result in compensatory inactivity during out-of-work hours. Prev Med. 2011;53(1-2):48–52. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ypmed.2011.04.018.

    Article  PubMed  Google Scholar 

  50. Rhodes RE, Janssen I, Bredin SSD, Warburton DER, Bauman A. Physical activity: health impact, prevalence, correlates and interventions. Psychol Health. 2017;32(8):942–75. https://0-doi-org.brum.beds.ac.uk/10.1080/08870446.2017.1325486.

    Article  PubMed  Google Scholar 

  51. Adams J, Mytton O, White M, Monsivais P. Why are some population interventions for diet and obesity more equitable and effective than others? The role of individual agency. PLoS Med. 2016;13(4):e1001990. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pmed.1001990.

  52. Craike M, Wiesner G, Hilland TA, Bengoechea EG. Interventions to improve physical activity among socioeconomically disadvantaged groups: an umbrella review. Int J Behav Nutr Phys Act. 2018;15(1):43. https://0-doi-org.brum.beds.ac.uk/10.1186/s12966-018-0676-2.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Strain T, Wijndaele K, Garcia L, Cowan M, Guthold R, Brage S, et al. Levels of domain-specific physical activity at work, in the household, for travel and for leisure among 327 789 adults from 104 countries. Br J Sports Med. 2020;54(24):1488–97. https://0-doi-org.brum.beds.ac.uk/10.1136/bjsports-2020-102601.

    Article  PubMed  Google Scholar 

  54. Gidlow C, Johnston LH, Crone D, Ellis N, James D. A systematic review of the relationship between socio-economic position and physical activity. Health Educ J. 2006;65(4):338–67. https://0-doi-org.brum.beds.ac.uk/10.1177/0017896906069378.

    Article  Google Scholar 

  55. Galobardes B, Shaw M, Lawlor DA, Lynch JW, Davey Smith G. Indicators of socioeconomic position (part 1). J Epidemiol Community Health. 2006;60(1):7–12. https://0-doi-org.brum.beds.ac.uk/10.1136/jech.2004.023531.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Galobardes B, Shaw M, Lawlor DA, Lynch JW, Davey Smith G. Indicators of socioeconomic position (part 2). J Epidemiol Community Health. 2006;60(2):95–101. https://0-doi-org.brum.beds.ac.uk/10.1136/jech.2004.028092.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Jørgensen MB, Gupta N, Korshøj M, Lagersted-Olsen J, Villumsen M, Mortensen OS, et al. The DPhacto cohort: an overview of technically measured physical activity at work and leisure in blue-collar sectors for practitioners and researchers. Appl Ergon. 2019;77:29–39. https://0-doi-org.brum.beds.ac.uk/10.1016/j.apergo.2019.01.003.

    Article  PubMed  Google Scholar 

  58. Gupta N, Christiansen CS, Hallman DM, Korshøj M, Carneiro IG, Holtermann A. Is objectively measured sitting time associated with low back pain? A cross-sectional investigation in the NOMAD study. PLoS One. 2015;10:e0121159.

    Article  Google Scholar 

  59. Kruger J, Yore MM, Ainsworth BE, Macera CA. Is participation in occupational physical activity associated with lifestyle physical activity levels? J Occup Environ Med. 2006;48(11):1143–8. https://0-doi-org.brum.beds.ac.uk/10.1097/01.jom.0000245919.37147.79.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the DPhacto and NOMAD research groups and personnel who contributed to the data collection.

Funding

The work was supported by a grant from the Danish government (Satspulje) and by The Danish Working Environment Research Fund (grant number 20175100213). Dorothea Dumuid is supported by the National Health and Medical Research Council (NHMRC) [Early Career Fellowship 1162166] and the National Heart Foundation of Australia [102084].

Author information

Authors and Affiliations

Authors

Contributions

CLR conceived the research idea of the study, conducted the analyses and was responsible for interpretation of the results and drafting the first draft of the manuscript. DD and KH advised on planning the study’s analytic strategy. AH and MBJ were principal investigators and responsible for the dataset used in this paper. All authors (i.e. CLR, DD, KH, NG, MBJ, KNN and AH) contributed to prepare the introduction and discussion sections, interpretation of the results, and have read and commented on the draft version and approved the final version of the entire manuscript.

Corresponding author

Correspondence to Charlotte Lund Rasmussen.

Ethics declarations

Ethics approval and consent to participate

The DPhacto and NOMAD studies were approved by the Danish Capital Region Committee on Health Research Ethics (Den Videnskabsetiske Komite for Region Hovedstaden) (reference number H-2-2012-011 [27] and reference number H-2-2011-047 [28], respectively). Both studies were conducted according to the Helsinki declaration and all data were anonymized in relation to individuals and workplaces. Written consent to participate was obtained from the study participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Calculation of pivot-coordinates and model development. Detailed description of how pivot-coordinates were calculated and models developed.

Additional file 2.

Association between day-to-day leisure time physical behaviours, weekday and relative sedentary work time. Results on the association between day-to-day relative work time spent sedentary and leisure time physical behaviours.

Additional file 3.

Sensitivity analyses. Results of sensitivity analyses in which only workers with at least two days of measurements were included.

Additional file 4.

Assessment of potential selection bias. Comparison of baseline characteristics between blue-collar workers excluded from and included in the study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rasmussen, C.L., Dumuid, D., Hron, K. et al. Day-to-day pattern of work and leisure time physical behaviours: are low socioeconomic status adults couch potatoes or work warriors?. BMC Public Health 21, 1342 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12889-021-11409-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12889-021-11409-0

Keywords