Article Text

Download PDFPDF

A worksite vitality intervention to improve older workers' lifestyle and vitality-related outcomes: results of a randomised controlled trial
  1. Jorien E Strijk1,
  2. Karin I Proper1,2,
  3. Allard J van der Beek1,2,
  4. Willem van Mechelen1,2
  1. 1Department of Public and Occupational Health, EMGO+ Institute for Health and Care Research, VU University Medical Center, Amsterdam, The Netherlands
  2. 2Body@Work, Research Center Physical Activity, Work and Health, TNO-VUmc, Amsterdam, The Netherlands
  1. Correspondence to Professor Allard J van der Beek, Occupational Epidemiology, VU University Medical Center, Department of Public and Occupational Health, EMGO+ Institute, Amsterdam 1081 BT, The Netherlands; a.vanderbeek{at}vumc.nl

Abstract

Objective To evaluate the effectiveness of a worksite vitality intervention on vigorous physical activity (VPA), fruit intake, aerobic capacity, mental health and need for recovery after work among older hospital workers (ie, 45 years and older).

Methods The 6-month intervention was evaluated using a randomised controlled trial design. Workers who were randomised to the intervention group (n=367; control: n=363) received the Vital@Work intervention containing (1) a Vitality Exercise Program (VEP) combined with (2) three visits to Personal Vitality Coach. The VEP consisted of a weekly yoga session, a weekly workout session and weekly unsupervised aerobic exercising. Free fruit was provided at the VEP. Data on the outcome measures were collected (ie, year 2009–2010) at baseline (n=730) and 6 months of follow-up after baseline (n=575) using questionnaires, accelerometers and 2 km walk tests. Effects were analysed according to the intention-to-treat principle with complete cases (n=575) and imputed data (n=730) using linear regression analyses. Additional analyses were performed for high yoga and workout compliance (ie, >mean number of sessions).

Results Effects were found for sports activities (β=40.4 min/week, 95% CI 13.0 to 67.7) and fruit intake (β=2.7 pieces/week, 95% CI 0.07 to 4.7) and were stronger for workers with high compliance to yoga (sport: β=49.6 min/week, 95% CI 13.9 to 85.2; fruit: β=3.8 pieces/week, 95% CI 1.1 to 6.4) and workout sessions (sport: β=72.9 min/week, 95% CI 36.1 to 109.8; fruit: β=4.0 pieces/week, 95% CI 1.1 to 6.4). The intervention group lowered their need for recovery, when compared to controls (β=−3.5, 95% CI −6.4 to −0.54), with stronger effects for high workout compliance (β=−5.3, 95% CI −9.3 to −1.3). No effects were found on VPA, aerobic capacity or mental health.

Conclusions Implementation of worksite yoga and workout facilities and minimal fruit interventions should be considered by employers to promote transitions into healthier lifestyles and thereby health.

  • Lifestyle
  • older workers
  • randomised controlled trial
  • workplace yoga
  • workplace aerobic exercising
  • ageing
  • workplace
  • health promotion
  • randomised trials

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/3.0/ and http://creativecommons.org/licenses/by-nc/3.0/legalcode

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Introduction

The baby boom after the Second World War, longer life expectancies and lower birth rates are leading to an ageing society and subsequently a shrinkage of the workforce. Hence, ageing workers are required in the near future. To do so, most Western countries, at which the official retirement age is around the age of 65 years,1 ,2 are considering to raise or have already raised the official retirement age. For instance, in the Netherlands, the government recently decided to raise the retirement age from 65 to 67 years in the upcoming decades.3 But also other solutions are needed to maintain ageing workers, who are defined as workers aged 45 years and older.4 This definition is based on the period at which major changes occur in functional capacities (eg, decline in aerobic capacity, higher prevalence of chronic diseases)5–7 and relevant work-related outcomes (eg, higher need for recovery (NFR) from work, lower work ability and more sick leave).8 ,9 Also, previous studies have shown that an age-related decline in health is a major contributor to early exit from work. Thus, in order to prolong the working life of older workers and increase their employability, it is important to promote and maintain good health.

A concept assumed to be closely related to and therefore may influence health is vitality. Vitality is related to both mental and physical factors of health.10–15 Regarding the mental factors, vitality reflects well-being, lower levels of fatigue, higher levels of emotional energy, mental resilience and perseverance.10–14 With respect to the physical factors, vitality is characterised by high energy levels and feeling ‘strong and fit’.14 In the field of occupational health, vitality has been described as one of the three dimensions of work engagement and is characterised by ‘feeling full of energy, strong and fit and being able to keep on working indefatigable’.14 ,15 In the Vital@Work study, a worksite lifestyle intervention was developed aiming at improving both mental (ie, by yoga sessions) and physical (ie, by aerobic exercising) factors of vitality.16

Healthy lifestyle choices, such as sufficient physical activity (PA) and healthy dietary habits (eg, sufficient fruit intake), contribute to better health outcomes, for example, mental health and lower risk for chronic diseases (ie, cardiovascular diseases, diabetes, cancer).17–21 Although the beneficial effects of yoga are not widely reported yet, studies among patient populations showed favourable effects on mental health and well-being.22–25 Therefore, a worksite health promotion (WHP) programme containing physical exercising (ie, aimed at improving aerobic capacity) and yoga is considered a potentially effective tool to keep older workers vital, promote their health and thereby prolong their labour participation. The beneficial effects of WHP programmes on work-related outcomes, such as sick leave and productivity, have indeed been reported.26–28 Also, positive effects of WHP programmes on health29 and aforementioned lifestyle behaviours have been reported.30 ,31

The objective of the present study was to evaluate the effectiveness of the Vital@Work intervention on (1) lifestyle behaviours, that is, sports, vigorous physical activities (VPA), fruit intake and (2) vitality-related outcomes, that is, aerobic capacity, mental health and the NFR after a day of work.

Methods

Study population and design

All workers aged 45 years and older from two academic hospitals in the Netherlands were invited to participate between April and October 2009. A worker was considered eligible when working at least 16 h a week, giving written informed consent and having no risk for developing adverse health effects when becoming physically active as assessed using the Physical Activity Readiness Questionnaire.32 The Medical Ethics Committee of VU University Medical Center approved the study protocol. Details on the randomised controlled trial design have been described extensively elsewhere.16

The workers who consented to participate were, after baseline measurements, individually randomised to the intervention or control group using Random Allocation Software (V.1.0, May 2004; Isfahan University of Medical Sciences, Isfahan, Iran). The research assistant notified each worker to which group he or she had been allocated and did not reveal the group allocation to the investigator responsible for data analyses. Blinding of participants or intervention providers was impossible. The sample size calculation is described extensively elsewhere16 but showed that 189 participants per group were needed at follow-up.

After randomisation, workers of both the intervention and control group received written information about a healthy lifestyle in general (ie, diet, PA and relaxation). Additionally, the intervention group received a 6-month lasting intervention consisting of (1) a Vitality Exercise Program (VEP) with (2) provision of free fruit and combined with (3) three visits to a Personal Vitality Coach (PVC). The VEP consisted of a weekly 45 min: (1) yoga session, (2) workout session and (3) unsupervised aerobic exercise session. Yoga was guided by a qualified yoga instructor and included relaxation and preparation postures for the joints, series of standing and forward bending postures and twists, light back bending postures and total relaxation and meditation. Workout sessions were guided by certified fitness instructors and consisted of a warming-up followed by aerobic exercises, resistance training and cooling-down. The intensity of the workout had to be 65%–90% of the age-predicted maximum heart rate (HRmax).33 ,34 The resistance training was progressive in nature and provided stimulus to all major muscle groups. Besides the yoga and workout sessions, workers were prescribed to perform weekly 45 min unsupervised vigorous PA (eg, fitness, running, spinning) with a similar intensity as the guided workout sessions. At the VEP group sessions, there was free provision of fruit. As to the PVC visits, the first PVC visit was scheduled at the start of the intervention and was followed by two consecutive visits at 4–6 weeks and 10–12 weeks after the first visit. During the 30-min PVC visits, five items were discussed, namely goal setting, confidence in achieving formulated goals, feedback on formulated goals, discussing barriers for formulated goals and problem solving. At the first visit, the items goal setting and confidence in achieving formulated goals were discussed. At the second and third visit, the following items were discussed, namely feedback on formulated goals, discussing barriers for formulated goals and problem solving.

Outcome measures

Data on outcome measures were collected (ie, year 2009–2010) at baseline (n=730) and 6 months of follow-up after baseline (n=575) using questionnaires, accelerometers and 2 km walk tests.

Lifestyle behaviours

The level of PA was measured both subjectively, using the Short QUestionnaire to ASses Health-enhancing physical activity (SQUASH),35 and objectively among a random subsample of 196 workers using Computer Science Application (CSA) accelerometers (Type GTM1 and ActiTrainer; ActiGraph, Pensacola, Florida, USA).

The SQUASH questionnaire35 measures habitual PA levels referring the past week of four PA domains, that is, commuting, occupational, household and leisure time. For each domain, workers were asked to indicate the frequency (times per week), self-reported intensity (light, moderate or vigorous) and average duration of the activity per day.35 As part of leisure time PA, workers could report up to four sports, with accompanying frequency, intensity and duration. For each domain, activities were subdivided into three age-dependent intensity categories (ie, light/moderate/vigorous), corresponding to the metabolic equivalents (METs) derived from Ainsworth's compendium of PAs.36 ,37 Since the Vital@Work intervention was aimed at improving sports and vigorous physical activities (VPA), the outcome measures used for this study were total minutes per week of: (1) sports activities, (2) VPA (aged<55 years: ≥6.5 METs; aged ≥55 years: ≥5.0 METs) and (3) total moderate-to-vigorous physical activities (MVPA). Total minutes per week of VPA and MVPA (aged <55 years: ≥4.0 METs; aged ≥55 years: ≥3.0 METs) were calculated by summing the time spent on at least moderate intensity activities across all domains.38

Among 214 older workers (nintervention=102; ncontrol=112), that is, a random sample out of the 730 workers included at baseline, PA was measured using accelerometers, which registered actual PA by detecting the magnitude of vertical accelerations and decelerations.39 Workers were asked to wear the accelerometers on the right hip during waking hours for 7 days. Activity counts were measured over a time interval of 60 s (epoch). Data were cleaned and scored using Meterplus40 according to the following criteria: (1) a wearing period of three or more valid days, (2) the minimal wearing time for a valid day was 10 h and (3) Freedson's cut-off points41 were used to quantify PA intensities. Since the wearing period varied between workers (range 3–7 days, the total minutes per intensity category were divided by the valid wearing days (resulting in minutes per day) multiplied by 7 (in minutes per week). The outcome measures obtained from the accelerometers were the total minutes per week of VPA: ≥5725 counts/min and MVPA: ≥1953 counts/min.

Weekly fruit intake was assessed using the validated Short Fruit and Vegetable Questionnaire,42 of which only the questions about fruit consumption were included, by verifying how many days per week they usually eat fruit and the number of pieces of fruit per day.

Vitality-related outcomes

Aerobic capacity (VO2max) was estimated using the UKK 2 km walk test. This test has shown to be a feasible and accurate method for predicting VO2max among healthy populations43 ,44 and was performed in a public park near the workplace. Workers were asked to fill out a form with their name, age and body height and weight and put on a heart rate monitor (type S610I; Polar Electro, Lake Success, New York, USA). Workers had to walk 2 km individually at a pace as brisk as possible, but without running. At the finish, the research assistant noted heart rate and performance time. VO2max was estimated using gender-specific equations including age, body mass index, performance time (in minutes) and heart rate at finish (HR).45

Mental health was measured using the mental health scale of the RAND-36 general health questionnaire,46 which consists of five questions that refer to the past 4 weeks: “Did you feel …?” (1) nervous, (2) down in dumps, (3) peaceful, (4) sad and (5) happy. The RAND-36 mental health scale, which has shown to be sufficiently reliable, ranges from 0 to 100 points (higher scores indicating better subjective mental health), with a score of >76.8 considered as good.46

NFR from work was assessed with a scale from the Dutch Questionnaire on the Experience and Evaluation of Work (VBBA)47 consisting of 11 statements (yes/no) concerning the recovery period after a day work. Examples of negative statements are: “My job causes me to feel rather exhausted at the end of a working day” and “I find it hard to relax at the end of a working day”. An illustration of a positive statement is: “After the evening meal, I generally feel good in shape”. The NFR scale ranges from 0 to 100 (higher scores being more unfavourable), with a score ≤54 being indicated as a sufficient NFR.47

Potential confounders and effect modifiers

At baseline, data on potential confounders and effect modifiers were assessed by questionnaire including age, gender (male/female), education (low=elementary school or less, medium=secondary education and high=college/university), chronic disease status (yes/no), smoking (yes/no), intervention location (Amsterdam/Leiden), type of work (blue/white collar) and marital status (married/cohabitating/single/divorced/widowed).

Statistical analysis

Differences in baseline characteristics (presented in table 1) between the intervention and control group and differences in outcome measures between completers and non-completers were tested using independent t test for continuous variables (eg, age, working hours per week) and Pearson's χ2 tests for categorical (ie, education) and dichotomous variables (eg, chronic diseases, smoking). To analyse the intervention effects (table 2), the differences in change over time between the intervention and control group (ie, β=the regression coefficient) were analysed using linear regression analyses.48 In this analyses, the outcome measures over 6-month follow-up (ie, PA, fruit intake, aerobic capacity, mental health and NFR) were regressed onto the baseline values of these outcomes. All analyses were performed according to the intention-to-treat principle. As adding potential confounders to crude models did not change intervention effects, >10% and no effect modifiers were found; only crude effect estimates are presented in this paper.

Table 1

Mean and SD for continuous measures and numbers and percentages (%) for dichotomous and categorical measures for the intervention (n=367) versus control group (n=363) and for completers (n=575) versus non-completers (n=155)

Table 2

Mean and SD for complete cases and imputed data for missing values on PA (ie, sport, VPA and MVPA), fruit intake, aerobic capacity (estimated VO2max), mental health and need for recovery for the intervention and control group at baseline and after 6 months of follow-up after baseline

As not all the individuals included in the intervention group had participated in yoga sessions and workout sessions in the same way,49 additional analyses were performed taking this into account. Additionally, data analyses were performed to determine significant relationships between the compliance of workers to the guided yoga and workout group sessions and the study outcomes (table 3). The compliance to the guided group sessions was defined based on the mean number of followed yoga and workout group sessions, which were distinguished from the process evaluation of the Vital@Work intervention and were 10.4 and 11.1 sessions per 24 weeks, respectively.49 Compliance categories defined were (1) control group workers (ie, reference), (2) intervention group workers, who did not follow a guided session, (2) low compliance: ≤mean number of sessions and (3) high compliance: >mean number of sessions. To test differences between these compliance groups, linear regression analyses were used with dummy variables for each compliance category with the control group as reference category.

Table 3

Dose–response relation for sports, fruit intake and need for recovery among yoga and workout subgroups for both complete cases and imputed data for missing values analyses

As the possible effects of the missing participants should be considered,50 ,51 it is recommended to perform, in addition to complete case analyses, sensitivity analyses with imputed data.52 For the sensitivity analyses (tables 2 and 3), missing data were imputed using multiple imputations based on Multivariate Imputation by Chained Equations.53 ,54 The multiple imputations procedure was performed in PASW (V.18.0) by generating 40 different data sets. By using Rubin's rules, PASW enabled to pool effects from these 40 data sets.55 All statistical analyses were performed using PASW.

Results

In figure 1, a flow diagram of the study population is presented. Enrolment of the study population took place between April and September 2009. In total, 730 older workers were included as they completed the baseline questionnaire and were subsequently randomised to the intervention (n=367) or control group (n=363). Of those workers randomised to the intervention group, 329 (89.6%) followed PVC visits, 259 (70.6%) and 234 (63.8%) guided yoga and workout sessions, respectively. After 6 months of follow-up, 575 workers (intervention n=293, control n=282) completed the questionnaire 6 months after baseline and were therefore used for complete cases analyses. Also, sensitivity analyses with imputed data among the total population (n=730) were performed. No adverse events of the intervention were reported. In table 1, baseline characteristics of the study population are presented with no significant differences between the groups in any of these variables. Employees with regular working hours (n=96) were more likely to complete participation (p=0.042) when compared to workers with irregular working hours (n=634).

Figure 1

Flow diagram Vital@Work study. PVC, Personal Vitality Coach.

Effectiveness on lifestyle and vitality-related outcomes

Complete cases analyses, as shown in table 2, revealed effectiveness on sports activities (β: 40.4 min/week, 95% CI 13.0 to 67.7). The control group workers increased their sports activities with 35.1 min/week, but when compared to the intervention group, this increase was statistically higher (75.3 min/week). As for the subjectively measured VPA, in total 134 workers (intervention n=63, control n=73) completed these measures. It appeared that both the intervention and control group increased their VPA from baseline to 6 months later (+159.5 vs +110.3 min/week, respectively), with no significant differences between groups (β=48.5 min/week, 95% CI −81.0 to 178.1) (table 2). Also based on the accelerometer data, there were no significant differences between groups (β=8.5 min/week, 95% CI −0.34 to 17.3). No effects were found on total weekly MVPA (SQUASH: β=−1.4 min/week, 95% CI −126.0 to 123.2; CSA: β=13.8 min/week, 95% CI −25.9 to 53.5). Regarding fruit intake, the intervention group workers improved their fruit intake significantly more when compared to the control group (+5.7 vs +2.7 pieces/week), resulting in an intervention effect on increasing fruit intake (β=2.7 pieces/week, 95% CI 0.63 to 4.7) (table 2). As for the vitality-related outcomes, no significant effects were found on aerobic capacity or mental health (table 2). As for NFR, the intervention group significantly decreased their NFR more when compared to the control group (−3.2 vs 0.6 points). Hence, the intervention was effective in decreasing workers' NFR (β=−3.5 points, 95% CI −6.4 to −0.54) (table 2).

Additional analyses

In table 3, only the relationships between the outcomes measures and guided group compliance that appeared to be significant are presented. A significant relationship was found between sports and high compliance to the guided yoga (β=49.6 min, 95% CI 13.9 to 85.2) and workout sessions (β=72.9 min/week, 95% CI 36.1 to 109.8) when compared to the control group (table 3). Also for fruit intake, effects were stronger in the high compliance group of both the yoga (β=3.8 pieces, 95% CI 1.1 to 6.4) and the workout sessions (β=4.0 pieces/week, 95% CI 1.1 to 6.4).

Sensitivity analyses

Sensitivity analyses with imputed data for missing values showed similar findings when compared to the complete cases analyses (tables 2 and 3). However, the effect sizes, derived from the analyses with imputed data, were consistently smaller when compared to the complete cases.

Discussion

This study showed that intervention group workers significantly increased their weekly sports activities and fruit intake when compared to control group workers. Also, the intervention favourably affected the NFR after a day of work. No effects were observed for VPA, aerobic capacity and mental health.

Limitations and strengths

Some limitations of this study can be indicated. First, we studied a relatively healthy older workers population, mainly consisting of female workers, making it more difficult to generalise the study results. Another limitation is that we failed to ensure vigorous intensity physical activity compliance during the guided workout session, which was required to improve aerobic capacity. Also, sensitivity analyses showed similar but smaller estimates of effects, when compared to the complete cases analyses, indicating that the risk of bias is minimal. This is a commonly seen consequent of imputation52 but could indicate a potentially biased estimation obtained from complete cases. Nevertheless, both complete cases and imputed data analyses showed comparable intervention effects (table 2). Presenting sensitivity analyses with imputed data, as done in our study, can also be indicated as strength, as complete cases analyses are mostly used in trials.52 Another strength was that the intervention was specifically tailored to the needs of this specific group of older workers.16 Furthermore, to our knowledge, this study is the first one investigating a worksite intervention consisting of both yoga and aerobic exercising aiming to promote the mental and physical components of health. Another strength was that, both subjective (ie, questionnaires) and objective measurements (ie, accelerometers and 2 km walk tests for measuring of PA and aerobic capacity, respectively) were used to obtain data, although the latter due to practical reasons only among a random subsample. As structural labour shortages are expected in the near future,2 it is essential to extend the labour participation of older workers. The results of this study are therefore innovative and provided valuable information.

Comparison with other studies

Reviews have shown that workplace PA interventions have positive effects on PA behaviour,29 ,56 ,57 but effects on VPA are less clear. Our VPA findings were supported by a previous study evaluating a workplace PA intervention consisting of resistance training and physical exercise, which reported no differences on VPA.58 Interestingly, our study showed a positive effect on weekly sports activities but not on VPA. A likely explanation for this could be that the SQUASH is not able to detect small changes over time (ie, poor responsiveness) regarding PA intensities59 but is more useful for detecting changes in clearly defined activities (ie, sports). As the responsiveness of the SQUASH has not been evaluated before, it is recommended to do so. Furthermore, objectively measured could be indicated as more suitable to study changes in PA over time, sample sizes in the present study were too small for detecting significant differences between the intervention and control group. To illustrate, a sample size calculation based on data from our study (α=0.05; power=0.90; two-sided testing, mean baseline VPA: 13.3±34.5 min/week, loss to follow-up of 25%) showed that 199 older workers per group, and 498 in total, were needed. In contrast, our study managed to collect data regarding VPA among 134 workers. To detect changes over time and consequently effectiveness in future studies using accelerometers, it is therefore advised to use sample size calculations based on the specific PA outcome measure, for instance VPA.

The positive finding regarding fruit intake were supported by two other studies.60 ,61 Both offered free fruit at the workplace aimed to promote fruit intake by increasing availability and accessibility. In most Europe countries, daily fruit recommendations are not met by the majority of the adult population,62 and a worldwide trend for decreased fruit consumption has been seen.63 Thus, to promote population-wide fruit intake, effective strategies are needed. As the majority of the adults spend most their time at work and the positive effects shown in this study, we recommend implementation of minimal worksite fruit interventions.

The lack of impact on aerobic fitness and mental health may be caused by the relatively healthy, fit and physically active group of highly educated middle-aged workers we studied. This could indicate a healthy worker effect.64 ,65 Healthier workers are more likely to stay in the workforce than those who are sick or physically unfit. This may be especially true for older hospital workers, a population that has to deal with higher physical workloads than an average Dutch worker. Specific for aerobic fitness, a review of Proper et al (2003)56 showed inconclusive evidence for worksite PA programmes on aerobic capacity. This is probably due to the fact that quite intensive PA, required for enhancement of aerobic capacity, are often not reached.56 As the mean number of attended workout sessions was 11.1 (SD=7.2) during a 24-week intervention period,49 this could indeed be the case in our study. To ensure exercising with certain intensity in future research, it is recommended to objectively monitor exercise intensity compliance using heart rate monitors or accelerometers.

Effects of PA or yoga programmes on NFR have not been investigated before. The improvement found in NFR can be indicated as a relevant effect. It is known that the NFR worsens during ageing, and this seems to be especially true for highly educated women, who are largely represented in our study.66 Besides, high NFR predicts sickness absence duration,67 which is an important predictor for early retirement.68 ,69

Conclusion and implications

The Vital@Work intervention was successful in increasing sports participation and fruit intake, and positively affected NFR after a day work, but was not effective in improving VPA, aerobic capacity or mental health. To increase involvement in sports among older workers, implementation of worksite yoga and workout facilities might be considered by employers. Moreover, minimal worksite fruit interventions are recommended to promote transitions into better lifestyles and, on the long run, health.

What is already known on this subject

  • A structural labour shortage is expected in the near future, making it important to extend the labour participation of older workers. By doing so, it is essential to promote and maintain their good health.

  • Previous research has shown that worksite physical activity programmes has positive effects on health- and work-related outcomes.

What this study adds

  • This study is the first one investigating a worksite intervention containing both yoga and aerobic exercising aiming to promote both the mental and the physical components of health.

  • As the labour participation of specifically older workers need to be extended, the results of this study are innovative and provide valuable information for both (occupational) health epidemiology and employers.

  • Implementation of worksite yoga and workout facilities might be considered by employers.

  • Minimal worksite fruit interventions are recommended to promote healthier transitions into better lifestyles and, on the long run, health.

References

Footnotes

  • Funding The Vital@Work study is financially supported by the ‘Foundation Institute GAK’. The trial is registered at the Dutch Trial Register (NTR) under trial registration number: NTR1240 (http://www.trialregister.nl/trialreg/admin/rctview.asp?TC=1240).

  • Competing interests None.

  • Patient consent Obtained.

  • Ethics approval The study was approved by Medical Ethics Committee of VU University Medical Center.

  • Provenance and peer review Not commissioned; externally peer reviewed.