- Research article
- Open Access
Large country differences in work outcomes in patients with RA – an analysis in the multinational study COMORA
Arthritis Research & Therapyvolume 19, Article number: 216 (2017)
We aimed to explore whether country of residence or specific country characteristics are associated with work outcomes in rheumatoid arthritis (RA).
Data from the 17 countries participating in the Comorbidities in RA (COMORA) study were used. Work outcomes were measured by the Work Productivity and Activity Impairment Questionnaire, addressing employment (yes/no), absenteeism (percentage of time; 3 categories) and presenteeism (percentage of at-work productivity restrictions; 4 categories). Contribution of country of residence, gross domestic product (GDP), Human Development Index (HDI), unemployment rate, social protection expenditures (SPE) or world region to work outcomes was investigated in adjusted (ordered) logistic regressions.
The patients (n = 2395) were younger than 60 years; mean age 48 (SD 9.2) years, 1972 (84%) female and 1065 (45%) employed. Large country differences were found. Taking the country with the best work outcome as reference, Moroccan patients had the lowest odds of being employed (OR 0.2 (95% CI 0.1; 0.3) vs. Germany) and highest odds of absenteeism (OR 13.2 (3.6; 48.3) vs. Japan). Patients in Taiwan had the highest odds of presenteeism (OR 13.0 (5.5; 30.9) vs. Venezuela). All country indices except SPE were associated with work outcomes. For example, patients in low-GDP countries had lower odds of employment (OR 0.6 (0.5; 0.8)), higher odds of absenteeism (OR 2.8 (2.0; 4.1)), but lower odds of presenteeism (OR 0.5 (0.4; 0.7)) compared to higher-GDP countries.
Substantial differences in work outcomes among patients with RA were observed between countries. Lower economic wealth and human development of countries were associated with worse employment and higher absenteeism, but lower presenteeism.
The worldwide economic crisis and the aging of the population pose challenges to employment perspectives. On the one hand, there is pressure on individuals to work longer and more efficiently, and on the other hand the number of available jobs is increasingly limited. Persons suffering from a chronic disease such as rheumatoid arthritis (RA) are already at increased risk for adverse work outcomes including restrictions in productivity while at work (presenteeism), sick leave (absenteeism) and eventually work disability or unemployment  and therefore are particularly vulnerable to a negative economic climate.
There is increasing evidence that environmental contextual factors such as nature of work or support from colleagues influence work outcome independent of biomedical and personal contextual factors [1, 2]. On this line, limited data suggest that system-level environmental factors, such as the prosperity or development of a country, play a decisive role. Dadoniene et al. observed lower employment rates among patients with RA in Lithuania, a country with an economy in transition, when compared to Norway, one of the wealthiest countries in Europe . Mau et al. reported lower employment in patients with RA residing in the former Eastern compared to Western German states, and attributed this to higher prevailing economic unemployment in former Eastern German states . Chung et al. observed a higher hazard of withdrawal from work due to RA-related work disability in Finland compared to the USA, for similar clinical and sociodemographic characteristics , and suggested the more supportive Finnish social security system with favorable disability allowance in case of disease-related work disability partly explained such differences. Finally, the QUEST RA study among 32 countries showed that patients in low-income countries continued to work with high levels of disease activity .
Better insight into system-related differences in work outcomes might be important when designing and implementing intervention programs aiming to improve short-term and long-term work outcomes of patients with RA, as the outcome may be strongly dependent on culture or regulations of the country from which participants originate. Moreover, such knowledge is crucial in multinational studies, as stratification by country or (groups of countries) based on a specific characteristic may be essential to avoid bias. As the focus of research on work outcome in rheumatology shifts towards presenteeism (i.e. productivity loss while at work) and sick leave, and as there is increasing evidence that the various work outcomes are differentially influenced by biomedical and contextual factors , it is important to address each of these outcomes in addition to employment status.
The Comorbidities in Rheumatoid Arthritis (COMORA) international study collected recent patient-level data from 17 countries across 5 continents including data on work outcomes. These data offer a unique opportunity to explore the impact of country of residence and specific country characteristics on employment, absenteeism and presenteeism. We hypothesized that higher economic wealth and greater human development of a country is associated with higher employment and lower absenteeism and presenteeism.
COMORA, an international, cross-sectional multicenter study in 17 countries, collected data in patients diagnosed with RA according to the 1987 American College of Rheumatology classification criteria for RA, who were older than 18 years, and had sufficient command of the questionnaire language. The study protocol was approved by one central and all local institutional ethics committees. Written informed consent was obtained from all subjects. Investigators and patients completed a questionnaire on demographics and self-reported health, including questions on work outcome .
Work outcome was measured by the Work Productivity Impairment Questionnaire (WPAI) [9,10,11]. Outcomes comprised current employment status (employed vs. not employed), and in those employed, the percentage of working hours absent (absenteeism) and percentage of at-work productivity restrictions (presenteeism), measured on a numeric rating scale (10 indicating health status completely prevented work) due to health problems during the past 7 days. Persons with long-term disability were therefore considered as not employed, and persons with short-term work disability (sick leave) were considered as employed and included in analyses with absenteeism or presenteeism as the outcome.
Sociodemographic and clinical characteristics
Sociodemographic and lifestyle information comprised age, gender, highest level of education achieved (primary school, secondary school, university), marital status (single, married, widowed, divorced), smoking status (past, current, never smoked), and weight (kg) and height (cm) from which the body mass index (BMI) was calculated (underweight, BMI < 18.5; normal weight, 18.5 ≥ BMI < 25; overweight, 25 ≥ BMI ≤ 30; and obese, BMI > 30) [12, 13].
Clinical characteristics comprised the modified 8-item Health Assessment Questionnaire (mHAQ) assessing the level of physical function and ranging from 0 to 3 (where 3 = worst function) [14, 15], the 28-joint Disease Activity Scale DAS28 ranging from 0 to (virtually) 10 , high rheumatoid factor (RF)/high positive anti-citrullinated protein antibody (ACPA). Physician-confirmed ischemic cardiovascular disease (myocardial infarction, stroke), cancer (colon, skin, lung, breast and uterus in women, prostate in men, and lymphoma), gastrointestinal diseases (diverticulitis, ulcers), infections (hepatitis), lung disease (chronic obstructive pulmonary disease and asthma) and psychiatric disorders (depression) were used to compute the Rheumatic Disease Comorbidity Index (RDCI) [17, 18]. As history of fractures was not collected, these were not included in the RDCI and therefore the score ranged from 0 to 8.
Data on socio-economic country indices were extracted from the International Monetary Fund (for gross domestic product (GDP) per capita), the International Labour Organization (for unemployment rates and social protection expenditures (SPE)) and the United Nations Development program (for Human Development Index (HDI)) [19,20,21,22]. The GDP per capita is an indicator of economic wealth and is expressed in international dollars (Intl$), to adjust local currency for purchasing power parity (PPP) allowing cross-country comparison . Based on the visual inspection of data, GDP was dichotomized as ≤ Int$20,000 or > Intl$$20,000. Unemployment rate, reflecting economic growth and efficiency, was dichotomized based on the median as ≤ 7.6% or > 7.6% . The HDI reflects a country’s level of development and includes life expectancy, education assessed by years of school attendance and the standard of living (GDP/capita) . HDI was dichotomized based on the United Nations Development Program classification resulting in two categories (moderate versus high and very high). The SPE reflects social expenditure for sickness, maternity, employment injury and disability, and was dichotomized using (a) the percentage SPE of the country’s GDP, dichotomized at the median (<2.15% vs. ≥ 2.15%) and (b) the absolute SPE in international dollars dichotomized at the median (SPE < Intl$801.50 vs. ≥ Intl$801.50) . For all indices we used the most recent year for which the data were available for all countries. For GDP per capita this was 2011, for HDI 2013 and for SPE 2014. Finally, countries were categorized by world region (Africa, Asia, North America, Europe, and Latin America) reflecting socio-cultural coherence. Additional file 1: Table S1 provides an overview of the countries according to the categories of the different indices and the overlap between categories.
Descriptive statistics were used to characterize the total sample. To avoid a strong influence of the age on retirement on the differences between countries in work outcomes, all further analyses were performed on the sample at age ≤ 60 years. Logistic regression was used to understand factors contributing to employment. When absenteeism or presenteeism was the outcome, the sample was limited to those employed and only countries with at least 30 employed persons were considered. A considerable proportion of participants indicated no absence or presenteeism. Based on the model fit (deviance), ordered logistic regression was preferred above zero-inflated regression approaches. Based on inspection and testing for the proportional odds assumption for ordered logistic regression, the categories chosen for absenteeism were 0%, > 0% to <100% and 100%. The categories for presenteeism were 0%, > 0% to 30%, > 30% to 50% and > 50% to 100%. To select the individual covariates associated with each outcome in multivariable models, manual forward selection (cutoff p value <0.05) was used. To understand the additional influence of country of residence, “country” was added to the multivariable model as a categorical variable, using the country with the highest employment, or the lowest absenteeism or presenteeism in the raw data as reference. Finally, “country” was replaced by the different country indices. For each model, the fit of models with the different country-related variables was compared to the model without any country variable using log-likelihood ratio tests (with a cutoff of p < 0.05 for significance). Statistical models were restricted to complete cases. Interactions between sociodemographic and clinical characteristics with GDP were explored for each work outcome. Analyses were performed using Stata 12 .
Overall, 3920 individuals (with a mean age of 56 years and of whom 81% were women) from 17 countries in 5 world regions comprising Europe (i.e. Austria (AT), France (FR), Germany (DE), Hungary (HU), Italy (IT), Netherlands (NL), Spain (ES), United Kingdom (UK)), Latin America (i.e. Uruguay (UR), Venezuela (VE), Argentina (AR)), Asia (i.e. Taiwan (TW), Japan (JP), Korea (KR)), Africa (i.e. Morocco (MA), Egypt (EG)) and North America (i.e. the United States of America (USA)) were included in COMORA. Clinical and disease characteristics are described in Table 1.
Overall 1,065/2359 (45.2%) of persons ≤ 60 years were employed with the highest percentage of employment in Germany (n = 71/113; 62.8%) and the lowest in Morocco (n = 33/181; 18.2%). As only 10 (55.6%) and 12 (70.6%) persons were employed in the UK and Uruguay respectively, these countries were excluded from analyses on absenteeism and presenteeism. Of the 1065 employed persons, 993 completed the question on absenteeism. The percentage of persons in each of the categories for absenteeism was 0% (n = 764 (76.9%)), > 0% and < 100% (n = 164 (16.5%)) and 100% (n = 65 (6.6%)). For the 914 completing the question on presenteeism the distribution was 0% (n = 283 (31.0%)), > 0% and ≤ 30% (n = 420 (46.0%)), > 30% and ≤ 50% (n = 203 (22.2%)) and > 50% and ≤ 100% (n = 8 (0.9%)). The mean absenteeism was 17.4% (SD ± 30.5) and was highest in Morocco 43.8% (SD 43.9) and lowest in Japan 2.6% (SD 9.2%). Among those that could experience presenteeism (employed and not 100% of the time absent), mean presenteeism was 25.3% (SD 25.6) and was highest in Morocco 34.0% (SD 22.3) and lowest in Venezuela 13.1% (SD 22.0). Additional file 2: Table S2 shows characteristics of the employed sample without age restriction.
Multivariable models with work outcomes
Missing values in covariates were present in 147/2359 (6.2%) cases when employment was the outcome, in 79/993 (8.0%) when absenteeism was the outcome and in 74/914 (8.1%) of those who had completed the presenteeism question (of note: persons with 100% absenteeism could not indicate presenteeism). An overview of missing values per variable in the total sample and differences in characteristics of the sample with and without missing values are presented in the Additional file 3: Table S3 and Additional file 4: Table S4, respectively.
Older age (OR 0.98 (95% CI 0.97; 0.99) per year), female gender (OR 0.36 (95% CI 0.28; 0.47)), lower educational level (low vs. high, OR 0.30 (95% CI 0.24; 0.47); medium vs. high, OR 0.63 (95% CI 0.50; 0.78)) were all significantly associated with lower odds of employment. In addition, each point increase in the HAQ score or increase in the DAS28 decreased the odds of employment by 0.64 (95% CI 0.53; 0.76) and 0.83 (95% CI 0.77; 0.89), respectively. However, comorbidities (RDCI) had no significant influence (OR 0.93 (95% CI 0.85; 1.02)).
In the next step, country of residence was significantly associated with employment. Germany, the country with the highest unadjusted employment was used as the reference category. In Argentina (OR 0.40 (95% CI 0.21; 0.75)), Morocco (OR 0.15 (95% CI 0.08; 0.30)), Egypt (OR 0.30 (95% CI 0.17; 0.55)), Venezuela (OR 0.34 (95% CI 0.18; 0.63)), Korea (OR 0.26 (95% CI 0.15; 0.46)), France (OR 0.56 (95% CI 0.32; 0.99)) and Taiwan (OR 0.40 (95% CI 0.22; 0.71)) there was significantly lower odds of being employed compared to Germany (Fig. 1a).
Patients from low-GDP countries had significantly lower odds of employment compared to high-GDP countries (OR 0.62 (95% CI 0.51; 0.77)), when country was replaced by GDP. As significant interactions were found between GDP and age, gender, level of education and the RDCI, the analyses were stratified by GDP. In high-GDP countries, older age, female gender, higher mHAQ score, higher DAS28 and presence of comorbidities were associated with the chance of not being employed. In lower-GDP countries lower education, higher DAS 28 and female sex reinforced the odds of not being employed (Table 2).
Lower HDI and higher SPE, but not higher unemployment, were significantly associated with lower odds of employment and each of the country indices significantly improved the model fit compared to the model without any country indices. Countries in Africa, Latin America and Asia had lower odds of employment compared to Europe. However, there was no significant difference between North America and Europe (Additional file 5: Table S5).
The odds of being in a higher category of absenteeism increased with each point increase in the HAQ score (OR 3.08 (95% CI 2.24; 4.23)), the DAS28 (OR 1.33 (95%CI 1.17; 1.51)) and the RDCI (OR 1.46 (95%CI 1.26; 1.69)). When adding country of residence, country was significantly associated with absenteeism. Since Japan had the lowest absenteeism we used it as the reference category.
The odds of being in the next higher category of absenteeism were significantly higher for Argentina (OR 4.17 (95% CI 1.18; 14.67)), Hungary (OR 7.60 (95% CI 2.14; 26.92)) and Morocco (OR 13.15 (95% CI 3.58; 48.29)) compared to Japan (Fig. 1b).
When replacing country of residence by GDP, the odds of being in the higher absenteeism group were higher in low-GDP countries (OR 2.85 (95% CI 1.98; 4.09)) than in high-GDP countries (Table 3). Interactions between GDP and the other covariates were not significant.
Similarly, lower HDI and higher unemployment rate, but not SPE, were associated with significantly increased odds of being in the higher absenteeism group, and each of the indices significantly improved the model fit compared to the model without any country indices. Countries in Europe, Africa and Latin America had significantly higher odds of absenteeism than countries in Asia (Additional file 5: Table S5).
The odds of being in a worse category of presenteeism increased for each point increase in the HAQ score (OR 4.45 (95% CI 3.22; 6.14)) and DAS28 (OR 1.56 (95% CI 1.40; 1.74)). Country of residence was significantly associated with presenteeism. When adding country, Venezuela had the lowest presenteeism and was used as the reference category. All countries except Spain had significantly higher odds of being in a higher presenteeism group compared to Venezuela. The highest odds were for Taiwan (OR 12.99 (95% CI 5.47; 30.85), Korea (OR 8.70 (95% CI 3.68; 20.58)) and Italy (OR 9.23 (95% CI 3.49; 24.31)) (Fig. 1c).
When replacing country of residence by GDP, low-GDP countries had significantly lower odds of being in a higher presenteeism group (0.49 ([0.35; 0.69)) compared to high-GDP countries (Table 3). Interactions between GDP and the other covariates were not significant.
Similarly, lower HDI and higher unemployment rate were associated with lower odds of being in a higher presenteeism group while lower SPE was associated with increased odds of being in the higher presenteeism group. Countries in Asia had significantly higher odds of being in a higher presenteeism group compared to Latin America (Additional file 5: Table S5).
Contribution of country indices to explanation of the work outcomes
When employment and presenteeism were the outcome, the largest improvement in log likelihood compared to the model without country index was found when the variable HDI was added (log likelihood chi-square test statistic 24.80 and 30.96, respectively; p < 0.001). When absenteeism was the outcome, the largest improvement in log likelihood compared to the model without country index was when the variable GDP was added (log likelihood chi-square test statistic 31.34; p < 0.001) (Additional file 5: Table S5). Naturally, when adding country of residence or continent the difference was largest for all outcomes due to the number of categories included in this variable, but also indicating the complexity of the construct “country of residence” comprising various cultural and system-related facets.
This study is the first to explore the impact of country of residence on sick leave (absenteeism) and presenteeism (at-work productivity loss) in addition to employment in patients with RA. We showed that country of residence was associated with relevant differences in employment status, absenteeism and presenteeism, independent of individual sociodemographic and clinical characteristics.
Specific-country socioeconomic indices could partially shed light on these differences. As hypothesized, higher economic wealth (GDP) and better human development (HDI) were associated with better employment rates and lower absenteeism among patients with RA. Paradoxically, the direction of association changed when at-work productivity, or presenteeism was the outcome: higher economic wealth and human development and lower unemployment were associated with more presenteeism. When grouping countries by world region, North America had the highest employment and lowest odds of absenteeism. The highest presenteeism was for Asia.
Country differences in work status in persons older than 50 years have been previously investigated in a European comparative household study. Analyses were adjusted for physician-confirmed chronic diseases such as musculoskeletal diseases. The chance of not participating in the labor force among persons with “poor self-perceived health” differed from 1.07 (Sweden) to 3.99 (Switzerland) . As most countries in this study guaranteed access to health care, differences were attributed to availability in regulations for employment protection across countries.
Also in RA, international differences in work status among patients with RA have been reported. Dadoniene et al. assumed better employment with higher economic prosperity, and Mau et al. assumed better employment with lower national rates of economic unemployment [3, 4]. Our findings provided evidence that all mechanisms suggested in the literature might indeed play a role. The HDI, accounting for educational level and life expectancy in addition to GDP, was at least as good and perhaps even a slightly better indicator of employment perspectives. Of course, we can only speculate that higher HDI or GPD is also associated with more efficient support systems for keeping the chronically ill in the labor force. While Sokka suggested that higher work disability among Finish compared to US patients with RA might be attributable to higher income substitution in the case of work disability , we failed to show an influence of expenditure on social protection. We need to realize that SPE comprises expenditures for all kind of social hazards. It would be interesting to have a national indicator of the level of income substitution for sick leave or work disability to understand how this influences worker outcomes.
The QUEST-RA study has already revealed differences in employment across a large number of countries,  but has not reported the independent magnitude of the effect of GDP. In addition to 1.5-fold higher odds of patients from countries with lower (compared to higher) GDP having no employment, we also showed that women and patients with low education were less likely to be employed. Of interest, in high-GDP countries higher age and more comorbidity had a stronger effect on not being employed. As in QUEST-RA, we saw a clear trend that in low-income countries, patients continue to work with poorer health compared to those in high-income countries, although the interaction was not significant. Overall, these findings suggest age and health restrictions (HAQ, comorbidities) are better acknowledged by the social security system in wealthier countries. Alternatively, increasing job demands and pressure for productivity in wealthier countries make it harder to sustain worker participation at the same level of disease severity in the current economic climate. International combined qualitative and quantitative research might provide more insight into the true explanations.
In the literature we found no studies reporting country differences in sick leave or presenteeism in persons with musculoskeletal disease, and also in RA we are the first to explore the role of country of residence on sick leave and presenteeism. While overall absenteeism followed the same pattern as employment, paradoxically there were higher rates of presenteeism in wealthier countries: likely, patients in high-income countries experience more pressure towards efficiency and productivity. This seems to confirm the suggestion that people in high-income countries experience high work pressure.
Several limitations have to be addressed. First, the national samples might not be representative of patients with RA in participating countries, especially when sample sizes were small. If participation bias was different in low-income compared to high-income countries, likely lower-educated and unemployed patients would have no access to health care in low-income countries and our data would even underestimate the country effect. Second, our study had a cross-sectional design, thus hindering causal inferences. With regard to exposure to country characteristics, it is of note that this factor is present before the outcome. Third, for the classification of countries following different socioeconomic indices the cutoff was chosen based on the distribution of the variable. For some variables (GDP and HDI) there are international classifications but these could not be applied to our sample as the distribution of countries would have been too uneven. For that reason, categories for the indices in the current analyses were based on data inspection that allowed identification of an arbitrary threshold that best discriminated between the two groups. In future research with a larger number of countries, it would be interesting to understand whether there is a specific threshold below which the effect of country wealth or human development adversely influences work outcomes or whether the relationship is more linear. Also, the country indices are not always assessed annually and at a similar point in time in all countries. However, the indices used in this study spanned an acceptable period.
Fourth, in the absence of data on levels of income substitution in the case of disability in sickness absence, we considered SPE as a remote proxy in the case of work restrictions due to illness. It must be noted that the personal insurance for work disability may be a potential confounder for which we were not able to adjust. Fifth, work status in the COMORA study only considered employment (yes/no), and therefore no further analyses could be performed to understand reasons for not participating (anymore) in the labor force. Also, information on presenteeism and absenteeism was collected through self-report from patients and comprised a short time span (7 days). However, it should be noted that the WPAI questionnaire and the recall periods are well-validated Specifically for absenteeism, the short time-frame is valid when the sample is sufficiently large, as sick leave likely occurs randomly between different persons over time [10, 25,26,27]. Seventh, our models were not adjusted for disease duration, as information in the COMORA study was restricted to identifying those with early disease (≤6 weeks). It is worth noting that in univariable analyses, having early disease was not associated with any of the work outcomes. Moreover in RA, collinearity between age and disease duration is common and often forces to chose on of both factors as confounder. Finally, in order to avoid a dominating impact of differences in the formal age of retirement on the results, the main analyses were restricted to persons age ≤60 years. Notwithstanding, it might also be of interest to compare work outcomes across countries while acknowledging an additional influence of the official retirement age (a well-known country-specific regulation). Additional file 6: Table S6, Additional file 7: Table S7 and Additional file 8: Table S8 and Additional file 9: Figures S1A-C present the results for the total sample (i.e. not excluding those age >60 years). Changes were seen in the raw and adjusted employment rates, but the overall conclusions of the regression models remained unaltered. An expected difference was a stronger impact of age on work status.
This study contributes to the understanding of influence of country of residence on three relevant work outcomes among patients with RA. This is important for international research, as outcomes of interventions aiming to improve health or work ability might not be simply transferrable between countries. Moreover, healthcare providers aiming to improve work outcomes in patients with RA should account for national socioeconomic system characteristics.
In conclusion, adverse work outcome in patients with rheumatoid arthritis is not only associated with clinical status and individual sociodemographic characteristics but also with the country of residence and country-specific features such as economic wealth, human development or geographical region. Comparison of work outcome across countries should account for system-level contextual factors.
Body mass index
Comorbidities in Rheumatoid Arthritis
28-Joint Disease Activity Scale
Gross domestic product
Human Development Index
Modified 8-item Health Assessment Questionnaire
Rheumatic Disease Comorbidity Index
Social protection expenditures
United States of America
Work Productivity Impairment Questionnaire
de Croon EM, Sluiter JK, Nijssen TF, Dijkmans BA, Lankhorst GJ, Frings-Dresen MH. Predictive factors of work disability in rheumatoid arthritis: a systematic literature review. Ann Rheum Dis. 2004;63(11):1362–7.
Lenssinck MLB, Burdorf A, Boonen A, Gignac MA, Hazes JMW, Luime JJ. Consequences of inflammatory arthritis for workplace productivity loss and sick leave: a systematic review. Ann Rheum Dis. 2013;72(4):493–505.
Dadoniene J, Uhlig T, Stropuviene S, Venalis A, Boonen A, Kvien TK. Disease activity and health status in rheumatoid arthritis: a case-control comparison between Norway and Lithuania. Ann Rheum Dis. 2003;62(3):231–5.
Mau W, Listing J, Huscher D, Zeidler H, Zink A. Employment across chronic inflammatory rheumatic diseases and comparison with the general population. J Rheumatol. 2005;32(4):721–8.
Chung CP, Sokka T, Arbogast PG, Pincus T. Work disability in early rheumatoid arthritis: higher rates but better clinical status in Finland compared with the US. Ann Rheum Dis. 2006;65(12):1653–7.
Sokka T, Kautiainen H, Pincus T, Verstappen SMM, Aggarwal A, Alten R, et al. Work disability remains a major problem in rheumatoid arthritis in the 2000s: data from 32 countries in the QUEST-RA Study. Arthritis Res Ther. 2010;12(2):R42.
Stolwijk C, Castillo-Ortiz JD, Gignac M, Luime J, Boonen A. Importance of contextual factors when measuring work outcome in ankylosing spondylitis: a systematic review by the OMERACT Worker Productivity Group. Arthritis Care Res. 2015;67(9):1316–27.
Dougados M, Soubrier M, Antunez A, Balint P, Balsa A, Buch MH, et al. Prevalence of comorbidities in rheumatoid arthritis and evaluation of their monitoring: results of an international, cross-sectional study (COMORA). Ann Rheum Dis. 2014;73(1):62–8.
Lambert J, Hansen BB, Arnould B, Grataloup G, Guillemin I, Hojbjerre L, et al. Linguistic validation into 20 languages and content validity of the rheumatoid arthritis-specific Work Productivity and Activity Impairment questionnaire. Patient. 2014;7(2):171–6.
Zhang W, Bansback N, Boonen A, Young A, Singh A, Anis AH. Validity of the work productivity and activity impairment questionnaire–general health version in patients with rheumatoid arthritis. Arthritis Res Ther. 2010;12(5):R177.
Reilly MC, Zbrozek AS, Dukes EM. The validity and reproducibility of a work productivity and activity impairment instrument. Pharmaco Econ. 1993;4(5):353–65.
WHO. Obesity: preventing and managing the global epidemic. Report of a WHO Consultation. WHO Technical Report Series 894. Geneva: World Health Organization; 2000.
WHO. Physical status: the use and interpretation of anthropometry. Report of a WHO Expert Committee. WHO Technical Report Series 854. Geneva: World Health Organization; 1995.
Uhlig T, Haavardsholm EA, Kvien TK. Comparison of the Health Assessment Questionnaire (HAQ) and the modified HAQ (MHAQ) in patients with rheumatoid arthritis. Rheumatology. 2006;45(4):454–8.
Bruce B, Fries JF. The Stanford Health Assessment Questionnaire: a review of its history, issues, progress, and documentation. J Rheumatol. 2003;30(1):167–78.
Prevoo MLL, Vanthof MA, Kuper HH, Vanleeuwen MA, Vandeputte LBA, Vanriel PLCM. Modified disease-activity scores that include 28-joint counts - development and validation in a prospective longitudinal-study of patients with rheumatoid-arthritis. Arthritis Rheum. 1995;38(1):44–8.
England BR, Sayles H, Mikuls TR, Johnson DS, Michaud K. Validation of the Rheumatic Disease Comorbidity Index. Arthritis Care Res (Hoboken). 2015;67(6):865-72. doi:10.1002/acr.22456.
Wolfe F, Michaud K, Li T, Katz RS. Chronic conditions and health problems in rheumatic diseases: comparisons with rheumatoid arthritis, noninflammatory rheumatic disorders, systemic lupus erythematosus, and fibromyalgia. J Rheumatol. 2010;37(2):305–15.
IMF. World Economic Outlook Database (http://www.imf.org/external/pubs/ft/weo/2011/02/weodata/index.aspx). Accessed Feb 2011.
Organization IL. Trends Econometric Models, http://www.ilo.org/global/research/global-reports/global-employment-trends/2014/WCMS_233936/lang%2D-en/index.html. Accessed Oct 2013.
PROGRAMME UND. Human Development Index and its components. 2013.
Organization IL. World Social Protection Report 2014/15-Building economic recovery, inclusive development and social justice. 2014.
StataCorp. Stata Statistical Software: Release 12. College Station: StataCorp LP. 2011.
Alavinia SM, Burdorf A. Unemployment and retirement and ill-health: a cross-sectional analysis across European countries. Int Arch Occup Environ Health. 2008;82(1):39–45.
Leggett S, van der Zee-Neuen A, Boonen A, Beaton D, Bojinca M, Bosworth A, et al. Content validity of global measures for at-work productivity in patients with rheumatic diseases: an international qualitative study. Rheumatology. 2016;55(8):1364–73.
Severens JL, Mulder J, Laheij RJ, Verbeek AL. Precision and accuracy in measuring absence from work as a basis for calculating productivity costs in The Netherlands. Soc Sci Med. 2000;51(2):243–9.
van Roijen L, Essink-Bot ML, Koopmanschap MA, Bonsel G, Rutten FF. Labor and health status in economic evaluation of health care. The Health and Labor Questionnaire. Int J Technol Assess Health Care. 1996;12(3):405–15.
The authors would like to thank the national principal investigators of the COMORA study: Gustavo Casado (Argentina), Josef Smolen (Austria), Bassel Kamal El-Zorkany (Egypt), Martin Soubrier (France), Gerd Burmester (Germany), Peter Balint (Hungary), Carlo Maurizio Montecucco (Italy), Masayoshi Harigai (Japan), Yeong-Wook Song (Korea), Najia Hajjaj-Hassouni (Morocco), Mart van de Laar (the Netherlands), Emilio Martin-Mola (Spain, Shue-Fen Luo (Taiwan), Jonathan Kay & Kevin Winthrop (USA), Gabriel Maciel (Uruguay), Anna Antunez (Venezuela); the scientific committee of COMORA study: Peter Balint, Gerd Burmester, Paul Emery, Jonathan Kay, Emilio Martin-Mola, Iain McInnes, Carlo Maurizio Montecucco and Josef Smolen; and the patients who agreed to participate in this study.
Data collection was supported by an unrestricted grant from Roche Ltd. Secondary analyses presented in the manuscript received no financial support.
Availability of data and materials
The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request. Permission was required to use data from the COMORA study. All data presented in the current manuscript are contained within the manuscript and Additional files 1, 2, 3, 4, 5, 6, 7, 8 and 9.
Ethics approval and consent to participate
Ethical approval was not required as this paper is based on secondary analyses of the previously published COMORA study (Ann Rheum Dis. 2014 Jan;73(1):62–8. doi:10.1136/annrheumdis-2013-204223. Epub 2013 Oct).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Overview of the countries in the categories of each country index and overlap between categories. (DOCX 1591 kb)
Sociodemographic and clinical characteristics for employed population by country and high and low gross domestic product (GDP). (DOCX 16 kb)
Missing values, n (%) in variables by country; n (%). (DOCX 15 kb)
Sociodemographic and lifestyle characteristics in subjects with and without missing data in the outcome variables. (DOCX 13 kb)
Overview of odds ratios for each country index (in fully adjusted models) and results of likelihood ratio chi-square (LR chi2) tests comparing models with and without country index variables for the sample age ≤60 years. (DOCX 17 kb)
Overview of odds ratios for each country index (in fully adjusted models) and results of likelihood ratio chi-square (LR chi2) tests comparing models with and without country index variables (for total sample). (DOCX 16 kb)
Association of individual sociodemographic and clinical characteristics with being employed stratified by gross domestic product (GDP) (total sample). (DOCX 14 kb)
Association of individual sociodemographic and clinical characteristics and gross domestic product (GDP) with absenteeism and presenteeism (total sample). (DOCX 13 kb)
Odds of being employed compared to reference country (total sample). Figure S1B. Odds of being in the higher absenteeism group per country compared to reference country Japan (total sample). Figure S1C. Odds of being in the higher presenteeism group per country compared to reference country Venezuela (total sample). (ZIP 77 kb)