Objectives: Sickness absence is the outcome in many epidemiologic studies and is often based on summary measures such as the number of sickness absences per year. In this study the use of modern statistical methods was examined by making better use of the available information. Since sickness absence data deal with events occurring over time, the use of statistical models for survival data has been reviewed, and the use of frailty models has been proposed for the analysis of such data.
Methods: Three methods for analyzing data on sickness absences were compared using a simulation study involving the following: (i) Poisson regression using a single outcome variable (number of sickness absences), (ii) analysis of time to first event using the Cox proportional hazards model, and (iii) frailty models, which are random effects proportional hazards models. Data from a study of the relation between the psychosocial work environment and sickness absence were used to illustrate the results.
Results: Standard methods were found to underestimate true effect sizes by approximately one-tenth [method i] and one-third [method ii] and to have lower statistical power than frailty models.
Conclusions: An uncritical use of standard methods may underestimate the effect of work environment exposures or leave predictors of sickness absence undiscovered.