Original article

Scand J Work Environ Health 2018;44(2):156-162    pdf

https://doi.org/10.5271/sjweh.3703 | Issue date:

Prediction models to identify workers at risk of sick leave due to low-back pain in the Dutch construction industry

by Bosman LC, Dijkstra L, Joling CI, Heymans MW, Twisk JWR, Roelen CAM

Objective The aim of this study was to develop a prediction model based on variables measured in occupational health checks to identify non-sick listed workers at risk of sick leave due to non-specific low-back pain (LBP).

Methods This cohort study comprised manual (N=22 648) and non-manual (N=9735) construction workers who participated in occupational health checks between 2010 and 2013. Occupational health check variables were used as potential predictors and LBP sick leave was recorded during 1-year follow-up. The prediction model was developed with logistic regression analysis among the manual construction workers and validated in non-manual construction workers. The performance of the prediction model was evaluated with explained variances (Nagelkerke’s R-square), calibration (Hosmer-Lemeshow test), and discrimination (area under the receiver operating curve, AUC) measures.

Results During follow-up, 178 (0.79%) manual and 17 (0.17%) non-manual construction workers reported LBP sick leave. Backward selection resulted in a model with pain/stiffness in the back, physician-diagnosed musculoskeletal disorders/injuries, postural physical demands, feeling healthy, vitality, and organization of work as predictor variables. The Nagelkerke’s R-square was 3.6%; calibration was adequate, but discrimination was poor (AUC=0.692; 95% CI 0.568–0.815).

Conclusions A prediction model based on occupational health check variables does not identify non-sick listed workers at increased risk of LBP sick leave correctly. The model could be used to exclude the workers at the lowest risk on LBP sick leave from costly preventive interventions.

This article refers to the following texts of the Journal: 2009;35(4):245-260  2004;30(6):459-467
The following article refers to this text: 2019;45(2):101-102

The prevalence of low-back pain (LBP) is high, especially in Western countries where it ranks the highest cause of disability (1). The costs of productivity loss and both sick leave and disability pension benefits are higher than the medical costs of LBP diagnosis and treatment (2). LBP sick leave is therefore a considerable societal burden and an important topic of research.

LBP sick leave is a multidimensional problem considering that it is associated with various individual, psychosocial, physical and work-related factors (35), which complicates the identification of workers at risk of LBP sick leave. Combining these factors in prediction models would improve the ability to determine the risk of LBP sick leave of workers. Several studies have developed prediction models for the duration or outcome of LBP sick leave among sick-listed workers (611). The ability of these models to discriminate between workers with and without long-term LBP sick leave was fair to good. After identifying these high risk workers, healthcare providers can refer them to interventions aimed at facilitating return to work shortly after reporting sick.

Prediction models that identify workers with an increased risk of LBP sick leave before they report sick would enable healthcare providers to refer workers at high risk of future LBP sick leave to preventive interventions. However, only few studies have examined models for predicting the risk of LBP sick leave in non-sick-listed workers (12, 13). Melloh et al (12) developed a prediction model that predicted the risk of LBP sick leave during 6-week follow-up of workers experiencing LBP, with job control, depression and functional limitations as predictor variables. The authors reported that the prediction model failed to predict LBP sick leave during longer follow-up periods. Salathé et al (13) succeeded in developing a prediction model for LBP sick leave during 1-year follow-up of workers experiencing LBP. The strongest predictors were a history of sick leave, work-related fear-avoidance behavior, and life satisfaction. There is still a need for more research on predicting the risk of LBP sick leave among healthy workers as the aforementioned studies only included workers who reported LBP at baseline, and it is well conceivable that these models do not apply to healthy workers without LBP. Furthermore, only the overall predictive performance of these prediction models was presented in terms of explained variances. Melloh et al (12) reported 47% explained variance for the 6-week model and Salathé et al (13) found 25–28% explained variance for the 1-year model.

The aim of the present study was to develop a prediction model specifically for the risk of LBP sick leave in non-sick-listed workers based on variables commonly measured in occupational health checks.

Methods

Study design and population

In the Netherlands, employers are obliged to enable their workers to participate in occupational health checks addressing health and working conditions every 2–4 years. Occupational health checks consist of a questionnaire followed by a consultation with an occupational physician (OP). Participation in occupational health checks is voluntary for workers. The present study was designed as a cohort study with the occupational health checks as baseline. The occupational health check variables were used as predictors of LBP sick leave during 1-year follow-up. We retrieved the occupational health checks performed among 32 383 construction workers between 2010 and 2013 because LBP is a major cause of sick leave among construction workers (14). The study population consisted of 22 648 manual workers (eg, stone masons, carpenters and plasterers) and 9735 non-manual (eg, technical workplace and office) workers in the construction industry as shown in Figure 1. Workers gave informed consent during the occupational health check for the use of the questionnaire data and sick leave data for research purposes. We wanted to investigate the performance of the LBP sick leave prediction model in different occupational settings. Therefore, the prediction model was developed among the manual workers and validated in the non-manual workers. Sick leave data recorded in the year following the occupational health check was obtained from an occupational health register.

Figure 1

Flow chart of the study population.

SJWEH-44-156-g001.tif

Outcome variable

The outcome variable of this study was non-specific LBP sick leave, which was diagnosed by an OP. In the Netherlands, sick-listed workers consult with an OP within six weeks of reporting sick. OP certify sick leave with a diagnostic code derived from the 10th version of the International Classification of Diseases (ICD-10). Sick leaves OP-certified within the ICD-10 categories M54.3 (sciatica), M54.4 (lumbago with sciatica) and M54.5 (low-back pain) were defined as non-specific LBP sick leave. Since we used OP-certified sick leaves, the study only included long-term LBP sick leave as outcome variable; LBP sick leave had a median duration of 52 days (range of 3–730 days).

Predictor variables

Variables included in the occupational health check questionnaire were selected as potential predictors based on the literature on LBP sick leave (1520) and based on the distributions of the variables. Variables with a narrow distribution (eg, gender 99.8% men) were discarded as predictors (21).

Health factors

The health check questionnaire measured age, pain/stiffness in the back (no=0, yes=1), pain/stiffness in neck, upper, and lower extremities (range 0–3 areas), number of physician-diagnosed musculoskeletal disorders/injuries (range 0–3 diagnoses), health complaints caused or worsened by work (no=0, yes=1), feeling healthy (no=0, yes=1), and sport activities (ranging from <1 time/month=0 to ≥3 times/week=4).

The hereafter mentioned variables were scored in accordance with a previous study on work stress among Dutch construction workers (22). The health check questionnaire items on “feeling irritated, agitated and/or depressed”, “difficulties with relaxing after work”, and “not feeling able to cope with the situation” were answered with “no” (=0) or “yes” (=1) and summed to the variable stress (Cronbach’s α=0.79).

Fatigue was measured with three items “feeling tired”, “feeling fit after dinner”, and “feeling exhausted at the end of a workday” (α=0.73). Items were answered with “no” (=0) or “yes” (=1) and were summed. Vitality included four items on “being active”, “being healthy”, “enjoying daily activities”, and “having confidence in the future” (α=0.74). Items were answered on a frequency scale ranging from “never” (=0) to “always” (=4), and were subsequently summed.

Work factors

Work factors included work satisfaction, which consisted of seven items on job insecurity, appreciation, rewards, and enjoying work (α=0.69). The organization of work included four items on work planning, communication at work, information on work tasks, and supervisor support (α=0.77). Psychological demands consisted of four items on mentally demanding work, working under time pressure, large amounts of work, and the number of working hours/week (α=0.61). Autonomy was measured with two items: “Can you decide how you do your work?” and “Can you influence your work pace?”. The six occupational health check questionnaire items on sitting, standing, and working in one position for a long time were clustered together with working in bended, kneeled, squatted, and awkward positions into the variable postural physical work demands (α=0.73). The health check items on lifting/pushing/pulling/carrying loads, repetitive movements, exerting large forces, and physically demanding work were clustered into the variable dynamic physical work demands (α=0.79). The items of the aforementioned work factors were answered with “no” (=0) or “yes” (=1) and were summed for each separate work factor.

The occupational health check questionnaire also included the items of the work ability index (WAI), which were summed to a workability score (range 7–49) according to the WAI guidelines (23). Furthermore, the number of years employed in the construction industry was assessed.

Statistical analysis

Statistical analyses were performed using SPSS for Windows, version 22 (released 2013. IBM Corp, Armonk, NY, USA) and R (statistical computing) version 3.2.5.

Prediction model development

The prediction model was developed in the sample of manual workers (N=22 648). Based on the criterion of ten events per variable (24) and the fact that 178 of the manual workers had LBP sick leave during follow-up, we decided to include a maximum of 18 variables in a multivariable logistic regression model with LBP sick leave during 1-year follow-up (no=0, yes=1) as dependent outcome variable. The number of predictors in the 18-predictor model was then reduced by statistical backward selection to develop an easy-to-use prediction model for healthcare practice. The Akaike’s information criterion (AIC) was used as stopping rule for the exclusion of predictors, which meant that all variables with P<0.157 remained in the model (25).

The Nagelkerke’s R-square was regarded as a measure for overall predictive performance (26). Calibration (ie, the agreement between predicted and observed LBP sick leave risks) was assessed by the Hosmer-Lemeshow goodness-of-fit test; P≥0.05 indicates that predicted risks do not deviate significantly from observed risks, meaning that calibration is adequate. The ability of the prediction model to discriminate between initially non-sick listed workers with and without LBP sick leave during follow-up was investigated with receiver operating characteristic (ROC) analysis. The area under the ROC curve (AUC) reflects discrimination.

The performance of prediction models is generally better for subjects used to develop the model than for new subjects, a phenomenon called over-optimism. This caused by the adaptation of the prediction model to the data in the development sample. The prediction model was adjusted for over-optimism by drawing 250 bootstrap samples from the manual workers and comparing the performance of the prediction model in the bootstrap samples with the performance in the original sample (21).

Performance of the prediction model across occupational settings

The over-optimism adjusted final model was then applied to the sample of non-manual workers (N=9735) to compare the performance from both groups.

Based on the 0.8% incidence of LBP sick leave among manual workers, we investigated specificities and sensitivities among non-manual workers for predicted probabilities of 0.2–3.6% (0.25–4× incidence).

Results

The characteristics of the study populations are shown in table 1. The manual workers reported more physician-diagnosed musculoskeletal disorders/injuries, experienced higher (postural and dynamic) physical work demands, and had more LBP sick leaves than the non-manual workers.

Table 1

Study population characteristics. [SD=standard deviation.]

Manual construction workers (N=22 648) Non-manual construction workers (N=9735)


N % Mean SD N % Mean SD
Gender
 Men 22 603 99.8 7942 81.6
 Women 45 0.2 1793 18.4
Age (years) 44.5 12.2 44.4 10.7
Employed (years)
 Construction industry 24.3 13.0 19.7 12.7
 Current company 14.2 11.0 12.3 10.6
 Current position 19.1 13.1 10.3 8.8
Work (hours/week) 39.8 6.5 41.1 10.0
Pain/stiffness (back) 9057 40.0 3004 30.9
Pain/stiffness (other) 13 143 58.1 4788 49.2
Physician-diagnosed musculoskeletal disorders/injuries 9888 43.9 3426 35.5
Health complaints caused/worsened by work 6537 28.9 1191 12.2
Feeling healthy 2765 12.2 991 10.2
Sport activities
 <1 per month 7022 31.0 1964 20.2
 Once a month 1901 8.4 790 8.1
 2–3 times/month 2855 12.6 1390 14.3
 1–2 times/week 6643 29.3 3892 40.0
 ≥3 times/week 4178 18.4 1681 17.3
 Missing 49 18
Stress 0.77 1.46 0.88 1.51
Fatigue 0.94 1.10 0.92 1.09
Vitality 3.38 0.99 3.48 0.94
Work satisfaction 1.22 1.50 0.82 1.24
Work organization 0.5 1.00 0.4 0.9
Psychological work demands 1.31 1.16 2.41 1.32
Autonomy 0.45 0.98 0.39 0.89
Workability 39.54 6.48 41.77 5.24
Postural physical work demands 1.59 1.66 0.52 0.99
Dynamic physical work demands 2.08 1.37 0.29 0.77
Sick leave due to low-back pain 178 0.79 17 0.17

Development of the LBP sick leave prediction model

When all 18 predictors were included in the logistic regression model (Nagelkerke’s R-square 3.8%), pain/stiffness in the back was the strongest predictor of LBP sick leave, followed by physician-diagnosed musculoskeletal disorders/injuries (table 2). After backward statistical reduction, six predictor variables (pain/stiffness in the back, number of physician-diagnosed musculoskeletal disorders/injuries, postural physical demands, feeling healthy, organization of work, and vitality) remained in the final model (Nagelkerke’s R-square 3.5%; over-optimism adjusted Nagelkerke’s R-square 2.9%). Calibration of the 6-predictor model was adequate, ie, the predicted LBP sick leave risks did not deviate significantly (model Chi2=11.168; df=8; P=0.192) from the observed LBP sick leave risks. Discrimination between manual workers with and without LBP sick leave during follow-up was poor (AUC=0.659; 95% CI 0.615–0.702); the over-optimism adjusted AUC was 0.648.

Table 2

Regression coefficients of the 18-predictor and the reduced 6-preditor models in the manual construction workers sample. [SE=standard error.]

18-predictor model 6-predictor model


B SE Wald B SE Wald
Age 0.003 0.012 0.046
Years employed in the construction industry -0.005 0.011 0.234
Pain/stiffness (back) 0.619 0.180 11.823 0.632 0.173 13.329
Pain/stiffness (other) -0.061 0.091 0.451
Musculoskeletal disorders/injuries 0.257 0.091 7.958 0.232 0.075 9.538
Health complaints caused/worsened by work 0.195 0.194 1.012
Feeling healthy 0.308 0.227 1.851 0.341 0.215 2.530
Sports activities 0.030 0.051 0.345
Stress 0.032 0.061 0.271
Fatigue 0.073 0.080 0.813
Vitality 0.201 0.106 3.552 0.237 0.086 7.556
Work satisfaction -0.075 0.061 1.518
Work organization 0.161 0.076 4.420 0.126 0.067 3.470
Psychological work demands -0.118 0.072 2.718
Autonomy 0.086 0.111 0.604
Workability 0.014 0.019 0.558
Postural physical work demands 0.120 0.058 4.326 0.140 0.048 8.598
Dynamic physical work demands 0.073 0.071 1.045

Performance of the LBP sick leave prediction model

When applied to the non-manual workers, the explained variance of the 6-predictor model was 3.6%. The model showed adequate calibration (Chi2=8.455 df=7; P=0.294), but poor discrimination (AUC=0.692; 95% CI 0.568–0.815) (Figure 2) between non-manual workers with and without LBP sick leave during 1-year follow-up.

Figure 2

Receiver operating characteristic (ROC) curve of the 6-predictor model among manual construction workers; the diagonal indicated no discrimination above chance.

SJWEH-44-156-g002.tif

Specificities were higher and sensitivities lower in non-manual workers as compared with manual workers for all cut-off probabilities (table 3). Specificities ranged between 1.0–99.9%, and sensitivities between 5.9–100% in non-manual workers.

Table 3

Cut-off probabilities for the 6-predictor model.

Manual construction workers (N=22 648) Non-manual construction workers (N=9735)


Cut off probability Specificity / true negative False positive Sensitivity / true positive False negative Specificity / true negative False positive Sensitivity / true positive False negative









% N % N % N % N % N % N % N % N %
0.2 145 0.6 22 325 99.4 178 100 0 0 94 1.0 9624 99.0 17 100 0 0
0.4 5253 23.4 17 217 76.6 158 88.8 20 11.2 4070 41.9 5648 58.2 14 82.4 3 17.6
0.8 14 740 65.6 7730 34.4 106 59.6 72 40.4 7794 80.2 1924 19.8 7 41.2 10 58.8
1.6 20 660 91.9 1810 8.1 46 25.8 132 74.2 9465 97.4 253 2.6 2 11.8 15 88.2
2.4 22 146 98.6 324 1.4 13 7.3 165 92.7 9680 99.6 38 0.4 1 5.9 16 94.1
3.2 22 390 99.6 80 0.4 4 2.2 174 97.8 9709 99.9 9 0.1 1 5.9 16 94.1

Discussion

In this study, we developed a prediction model to identify non-sick listed workers at risk of LBP sick leave with variables measured in occupational health checks. A 6-predictor model including pain/stiffness in the back, physician-diagnosed musculoskeletal disorders/injuries, postural physical demands, feeling healthy, organization of work, and vitality as predictors showed adequate calibration, but low explained variance and poor discrimination between construction workers with and without LBP sick leave during 1-year follow up.

Predictors of LBP sick leave

Pain/stiffness in the back was the strongest predictor of LBP sick leave, in agreement with previous studies (1518). Furthermore, musculoskeletal complaints have previously been reported to predict LBP sick leave (15, 19), although these were not physician-diagnosed. In agreement with previous studies, postural physical work demands was a strong predictor in the 6-predictor model (17, 19). However, dynamic physical work demands were excluded from the 6-predictor model, while frequent bending and lifting/carrying of loads are risk factors of LBP according to the guidelines for OP (27). Dynamic physical demands were also significant predictors of LBP sick leave in a study among nurses’ aides (15) and among blue- and white-collar workers (20). However other studies confirmed our present finding that dynamic physical work demands did not predict LBP sick leave (17, 18).

LBP sick leave prediction model

The performance of the prediction model was poor in terms of explained variance and discrimination. Previously, Salathé et al (13) reported that prediction models including gender, age, sickness absence at baseline, fear avoidance behavior, life satisfaction, job satisfaction, and doing sports explained 25–28% of the variance of LBP sick leave risk during 1-year follow-up. Melloh et al (12) reported that a prediction model with job control, depression and functional limitations explained 47% of the risk of LBP sick leave after 6 weeks, but none of the predictor variables predicted LBP sick leave during longer follow-up. Both studies included workers who experienced LBP at baseline. It is conceivable that predicting LBP sick leave is more difficult among workers without LBP complaints. Furthermore, the fact that the prediction model presented by Melloh et al (12) failed to predict LBP sick leave occurring >6 weeks after baseline underlines the difficulty of predicting sick leave during longer periods of follow-up. It can also be debated whether LBP is too heterogeneous and therefore more difficult to predict. However no clear consensus is reached on subgrouping non-specific LBP patients to decrease the heterogeneity (28).

Strengths and limitations of the study

The prospective design with different data sources (health check questionnaire for baseline predictor variables and occupational health register for LBP sick leave) are strengths of the present study. LBP sick leave was OP-diagnosed and registered instead of worker-reported, which minimized self-report bias.

The study population consisted of older workers who have been employed in the construction industry for a long time. In other words, they were healthy enough to stay in the construction industry (healthy worker bias). Furthermore, healthy workers are more likely to participate in occupational health checks than workers with health complaints (healthy volunteer bias) (29). Both biases may have under-estimated the associations between risk factors and LBP sick leave. The low incidence of LBP sick leave in the study population could also be explained by these types of bias. This restricted the number of variables to be included in the full prediction model and could well explain the low predictive performance in terms of explained variance (21).

The low explained variance also indicates that important predictors could have been missed. This study only used data collected by the occupational health check questionnaire. Consequently, not all known risk factors of LBP were included in the prediction model. Furthermore, the answer options of the items of the factors were often limited to “yes” and “no”, which provides less information and may have reduced the predictive performance of the factors. Additionally, it can be debated whether physical demands can be measured with self-report questionnaires as its validity is often found uncertain (30). However, physical demands measured by a questionnaire (15) and video recordings (20) were both predictive of LBP sick leave.

Practical implications and future research

The prediction model presented in this study is based on occupational health check data and could therefore be readily implemented on a large scale as a screening tool for LBP sick leave if the performance would be sufficient. Although the manual and non-manual workers differed on baseline characteristics, the performance of the model was similar in both populations. The AUC is slightly higher among the non-manual than manual workers, but considering the overlapping confidence intervals, this difference was not statistically significant. The similar performance would indicate that the prediction model performs equally in different populations.

The prediction model based on occupational health check variables poorly discriminated between workers with and without LBP sick leave during follow-up. However, the model’s risk predictions were adequate in the sense that the risks predicted by the model did not deviate significantly from observed LBP sick leave frequencies. The model could be used to estimate a worker’s risk of LBP sick leave. For example, if a construction company considers an ergonomic training aimed at preventing LBP, but is unable or unwilling to train all staff, then the workers at lowest risk of LBP sick leave could be excluded from the training. In the present study, there were no false negatives when the cut-off risk was set at 0.2% (ie, a quarter of the population incidence). In other words, none of the construction workers with a predicted risk of LBP sick leave lower than 0.2% had LBP sick leave during 1-year follow-up.

The occupational health check variables did not provide sufficient information to identify workers at risk of LBP sick leave. It should be noted that occupational health checks are conducted to address health and working conditions, not specifically for predicting (LBP) sick leave. Additional predictors of LBP sick leave are needed to improve the discriminative ability of the prediction model. The present study investigated a prediction model for risk of LBP sick leave among healthy workers. Since it is difficult to predict LBP sick leave among healthy workers, it might be better to focus on prediction models for sick leave durations and outcomes among workers who have reported sick with LBP sick leave because previous studies have shown higher performance measures for such models (611).

Acknowledgements

The authors wish to thank Jan Golsteijn, director of Volandis knowledge centre for the construction industry, for supplying the data on workplace settings.

The authors declare no conflicts of interests.

References

1 

Hoy, D, March, L, Brooks, P, Blyth, F, Woolf, A, Bain, C, et al. (2014). The global burden of low back pain:estimates from the Global Burden of Disease 2010 study. Ann Rheum Dis, 73(7), 1309-15, https://doi.org/10.1136/annrheumdis-2013-204431.

2 

Dagenais, S, Caro, J, & Haldeman, S. (2008). A systematic review of low back pain cost of illness studies in the United States and internationally. Spine J, 8(1), 8-20, https://doi.org/10.1016/j.spinee.2007.10.005.

3 

Campbell, P, Wynne-Jones, G, Muller, S, & Dunn, KM. (2013). The influence of employment social support for risk and prognosis in nonspecific back pain:a systematic review and critical synthesis. Int Arch Occup Environ Health, 86(2), 119-37, https://doi.org/10.1007/s00420-012-0804-2.

4 

Ramond, A, Bouton, C, Richard, I, Roquelaure, Y, Baufreton, C, Legrand, E, et al. (2010). Psychosocial risk factors for chronic low back pain in primary care—a systematic review. Fam Pract, 28(1), 12-21, https://doi.org/10.1093/fampra/cmq072.

5 

Steenstra, I, Verbeek, J, Heymans, M, & Bongers, P. (2005). Prognostic factors for duration of sick leave in patients sick listed with acute low back pain:a systematic review of the literature. Occup Environ Med, 62(12), 851-60, https://doi.org/10.1136/oem.2004.015842.

6 

Lindell, O, Johansson, SE, & Strender, LE. (2010). Predictors of stable return-to-work in non-acute, non-specific spinal pain:low total prior sick-listing, high self prediction and young age. A two-year prospective cohort study. BMC Fam Pract, 11(1), 53, https://doi.org/10.1186/1471-2296-11-53.

7 

McIntosh, G, Frank, J, Hogg-Johnson, S, Bombardier, C, & Hall, H. (2000). Prognostic factors for time receiving workers'compensation benefits in a cohort of patients with low back pain. Spine, 25(2), 147-57, https://doi.org/10.1097/00007632-200001150-00003.

8 

Heymans, MW, De Vet, HC, Knol, DL, Bongers, PM, Koes, BW, & Van Mechelen, W. (2006). Workers' beliefs and expectations affect return to work over 12 months. J Occup Rehabil, 16(4), 685-95, https://doi.org/10.1007/s10926-006-9058-8.

9 

Jensen, OK, Stengaard-Pedersen, K, Jensen, C, & Nielsen, CV. (2013). Prediction model for unsuccessful return to work after hospital-based intervention in low back pain patients. BMC Musculoskelet Disord, 14(1), 140, https://doi.org/10.1186/1471-2474-14-140.

10 

Storheim, K, Brox, JI, Holm, I, & Bo, K. (2005). Predictors of return to work in patients sick listed for sub-acute low back pain:A 12-month follow-up study. J Rehabil Med, 37(6), 365-71, https://doi.org/10.1080/16501970510040344.

11 

Truchon, M, Cote, D, Schmouth, ME, Leblond, J, Fillion, L, & Dionne, C. (2010). Validation of an Adaptation of the Stress Process Model for Predicting Low Back Pain Related Long-term Disability Outcomes A Cohort Study. Spine, 35(13), 1307-15, https://doi.org/10.1097/BRS.0b013e3181c03d06.

12 

Melloh, M, Elfering, A, Salathe, CR, Kaser, A, Barz, T, Roder, C, et al. (2012). Predictors of Sickness Absence in Patients with a New Episode of Low Back Pain in Primary Care. Ind Health, 50(4), 288-98, https://doi.org/10.2486/indhealth.MS1335.

13 

Salathé, CR, Melloh, M, Mannion, A, Tamcan, Ö, Müller, U, Boos, N, et al. (2012). Resources for preventing sickness absence due to low back pain. Occup Med (Lond), 62(4), 273-80, https://doi.org/10.1093/occmed/kqs024.

14 

Afrian, K. (2010). Ziekteverzuim in de Bouw [Absenteeism in construction.]. Amsterdam: Stichting Economisch Instituut voor de Bouw.

15 

Eriksen, W, Bruusgaard, D, & Knardahl, S. (2004). Work factors as predictors of intense or disabling low back pain;a prospective study of nurses'aides. Occup Environ Med, 61(5), 398-404, https://doi.org/10.1136/oem.2003.008482.

16 

Morken, T, Riise, T, Moen, B, Hauge, SH, Holien, S, Langedrag, A, et al. (2003). Low back pain and widespread pain predict sickness absence among industrial workers. BMC Musculoskelet Disord, 4(1), 21, https://doi.org/10.1186/1471-2474-4-21.

17 

Tubach, F, Leclerc, A, Landre, M-F, & Pietri-Taleb, F. (2002). Risk factors for sick leave due to low back pain:a prospective study. J Occup Environ Med, 44(5), 451-8, https://doi.org/10.1097/00043764-200205000-00013.

18 

van den, Heuvel SG, Ariëns, GA, Boshuizen, HC, Hoogendoorn, WE, & Bongers, PM. (2004). Prognostic factors related to recurrent low-back pain and sickness absence. Scand J Work Environ Health, 30(6), 459-67, https://doi.org/10.5271/sjweh.835.

19 

Alexopoulos, EC, Konstantinou, EC, Bakoyannis, G, Tanagra, D, & Burdorf, A. (2008). Risk factors for sickness absence due to low back pain and prognostic factors for return to work in a cohort of shipyard workers. Eur Spine J, 17(9), 1185-92, https://doi.org/10.1007/s00586-008-0711-0.

20 

Hoogendoorn, W, Bongers, P, De, Vet H, Ariens, G, Van, Mechelen W, & Bouter, L. (2002). High physical work load and low job satisfaction increase the risk of sickness absence due to low back pain:results of a prospective cohort study. Occup Environ Med, 59(5), 323-8, https://doi.org/10.1136/oem.59.5.323.

21 

Steyerberg, E. (2009). Clinical prediction models:a practical approach to development, validation, and updating, New York, Springer Science & Business Media, https://doi.org/10.1007/978-0-387-77244-8.

22 

Dijkstra, L, Houtkamp, R, & Hamming, C. (2016). Diagnose Werkdruk/Stress:Ontwikkeling van een regel voor bedrijfsartsen in de bouw [Diagnosis work pressure / Stress:Development of a rule for occupational physicians in the construction industry]. Netherlands: Harderwijk.

23 

Tuomi, K, & Oja, G. (1998). Work ability index. Helsinki: Finnish Institute of Occupational Health.

24 

Peduzzi, P, Concato, J, Kemper, E, Holford, TR, & Feinstein, AR. (1996). A simulation study of the number of events per variable in logistic regression analysis. J Clin Epidemiol, 49(12), 1373-9, https://doi.org/10.1016/S0895-4356(96)00236-3.

25 

Ambler, G, Brady, AR, & Royston, P. (2002). Simplifying a prognostic model:a simulation study based on clinical data. Stat Med, 21(24), 3803-22, https://doi.org/10.1002/sim.1422.

26 

Heinzl, H, Waldhör, T, & Mittlböck, M. (2005). Careful use of pseudo R2measures in epidemiological studies. Stat Med, 24(18), 2867-72, https://doi.org/10.1002/sim.2168.

27 

Gezondheidsraad. (2005). Verzekeringsgeneeskundige protocollen:Aspecifieke lage rugpijn, Hartinfarct. [Insurance medicine protocols:A specific low-back pain, myocardial infarction.] (Vol. 15). The Hague: Gezondheidsraad. Publication.

28 

Wand, BM, & O'Connell, NE. (2008). Chronic non-specific low back pain–sub-groups or a single mechanism? BMC Musculoskelet Disord, 9(1), 11, https://doi.org/10.1186/1471-2474-9-11.

29 

Martikainen, P, Laaksonen, M, Piha, K, & Lallukka, T. (2007). Does survey non-response bias the association between occupational social class and health? Scand J Public Health, 35(2), 212-5, https://doi.org/10.1080/14034940600996563.

30 

Barrero, LH, Katz, JN, & Dennerlein, JT. (2009). Validity of self-reported mechanical demands for occupational epidemiologic research of musculoskeletal disorders. Scand J Work Environ Health, 35(4), 245, https://doi.org/10.5271/sjweh.1335.