Misspecification of confounder-exposure and confounder-outcome associations leads to bias in effect estimates.

BACKGROUND: Confounding is a common issue in epidemiological research. Commonly used confounder-adjustment methods include multivariable regression analysis and propensity score methods. Although it is common practice to assess the linearity assumption for the exposure-outcome effect, most researchers do not assess linearity of the relationship between the confounder and the exposure and between the confounder and the outcome before adjusting for the confounder in the analysis. Failing to take the true non-linear functional form of the confounder-exposure and confounder-outcome associations into account may result in an under- or overestimation of the true exposure effect. Therefore, this paper aims to demonstrate the importance of assessing the linearity assumption for confounder-exposure and confounder-outcome associations and the importance of correctly specifying these associations when the linearity assumption is violated. METHODS: A Monte Carlo simulation study was used to assess and compare the performance of confounder-adjustment methods when the functional form of the confounder-exposure and confounder-outcome associations were misspecified (i.e., linearity was wrongly assumed) and correctly specified (i.e., linearity was rightly assumed) under multiple sample sizes. An empirical data example was used to illustrate that the misspecification of confounder-exposure and confounder-outcome associations leads to bias. RESULTS: The simulation study illustrated that the exposure effect estimate will be biased when for propensity score (PS) methods the confounder-exposure association is misspecified. For methods in which the outcome is regressed on the confounder or the PS, the exposure effect estimate will be biased if the confounder-outcome association is misspecified. In the empirical data example, correct specification of the confounder-exposure and confounder-outcome associations resulted in smaller exposure effect estimates. CONCLUSION: When attempting to remove bias by adjusting for confounding, misspecification of the confounder-exposure and confounder-outcome associations might actually introduce bias. It is therefore important that researchers not only assess the linearity of the exposure-outcome effect, but also of the confounder-exposure or confounder-outcome associations depending on the confounder-adjustment method used.

Effect of Partial Sick Leave on Sick Leave Duration in Employees with Musculoskeletal Disorders.

Objective This study determined if partial sick leave was associated with a shorter duration of sick leave due to musculoskeletal disorders (MSD) based on routinely collected health data in Dutch sick-listed employees. Furthermore, the effect of timing of partial sick leave on sick leave duration was determined. Methods This cohort study consisted of 771 employees with partial sick leave and 198 employees with full-time sick leave who participated in an occupational health check, and had sick leave due to MSD for minimally 4 weeks and were diagnosed by an occupational physician. Multivariable linear regression models were performed to determine the effects of partial sick leave (unadjusted and adjusted for confounders and MSD diagnosis) and Kaplan-Meier curves were presented for visualization of return to work for different timings of starting partial sick leave. Furthermore, linear regression analysis were done in subsets of employees with different minimal durations of sick leave to estimate the effects of timing of partial sick leave. Results Initial results suggest that partial sick leave was associated with longer sick leave duration, also when adjusted for confounders and sick leave diagnosis. Secondary results which accounted for the timing of partial sick leave suggest that partial sick leave had no effect on the duration of sick leave. Conclusion Partial sick leave does not influence MSD sick leave duration in this study when accounting for the timing of partial sick leave.

Development of Prediction Models for Sick Leave Due to Musculoskeletal Disorders.

Purpose The aim of this study was to develop prediction models to determine the risk of sick leave due to musculoskeletal disorders (MSD) in non-sick listed employees and to compare models for short-term (i.e., 3 months) and long-term (i.e., 12 months) predictions. Methods Cohort study including 49,158 Dutch employees who participated in occupational health checks between 2009 and 2015 and sick leave data recorded during 12 months follow-up. Prediction models for MSD sick leave within 3 and 12 months after the health check were developed with logistic regression analysis using routinely assessed health check variables. The performance of the prediction models was evaluated with explained variance (Nagelkerke's R-square), calibration (Hosmer-Lemeshow test) and discrimination (area under the receiver operating characteristic curve, AUC) measures. Results A total of 376 (0.8%) and 1193 (2.4%) employees had MSD sick leave within 3 and 12 months after the health check. The prediction models included similar predictor variables (educational level, musculoskeletal complaints, distress, supervisor social support, work-home interference, intrinsic motivation, development opportunities, and work pace). The explained variances were 7.6% and 8.8% for the model with 3 and 12 months follow-up, respectively. Both prediction models showed adequate calibration and discriminated between employees with and without MSD sick leave 3 months (AUC = 0.761; Interquartile range [IQR] 0.759-0.763) and 12 months (AUC = 0.740; IQR 0.738-0.741) after the health check. Conclusion The prediction models could be used to determine the risk of MSD sick leave in non-sick listed employees and invite them to preventive consultations with occupational health providers.

Development of Prediction Model for the Prognosis of Sick Leave Due to Low Back Pain.

OBJECTIVE: The aim of this study was to develop a prediction model for the prognosis of sick leave due to low back pain (LBP). METHODS: This is a cohort study with 103 employees sick-listed due to non-specific LBP and spinal disc herniation. A prediction model was developed based on questionnaire data and registered sick leave data with follow up of 180 days. RESULTS: At follow up 31 (30.1%) employees were still sick-listed due to LBP. Forward selection procedure resulted in a model with: catastrophizing, musculoskeletal work load, and disability. The explained variance was 27.3%, calibration was adequate and discrimination was fair with area under the ROC-curve (AUC) = 0.761 (interquartile range [IQR]: 0.755-0.770). CONCLUSION: The prediction model of this study can adequately predict LBP sick leave after 180 days and could be used for employees sick listed due LBP.

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

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.