Modelling the relationship between the WOMAC osteoarthritis index and EQ-5D
- Allan Wailoo^{1}Email author,
- Monica Hernandez Alava^{1} and
- Antonio Escobar Martinez^{2}
DOI: 10.1186/1477-7525-12-37
© Wailoo et al.; licensee BioMed Central Ltd. 2014
Received: 1 July 2013
Accepted: 5 March 2014
Published: 12 March 2014
Abstract
Objective
Economic evaluation typically is conducted using health state utilities to estimate treatment benefits. However, such outcomes are often missing from studies of clinical effectiveness. This study aims to bridge that gap by providing appropriate methods to link values from the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) to the EQ-5D utility instrument.
Method
Patients from a large registry of Spanish patients (n = 7072 observations) with knee or hip osteoarthritis who completed both WOMAC and EQ-5D was used. A mixture model approach was used based on distributions bespoke to the EQ-5D UK value set to estimate EQ-5D as a function of WOMAC pain, stiffness and function subscores.
Results
A five class mixture model provides very close fit to the observed data at all levels of disease severity. The overall mean (0.542 vs 0.542), median (0.620 vs 0.636) and the percentage of observations at full health (15 vs 14.8) were very similar between the observed data and the estimated model respectively. Stiffness has limited relationship to EQ-5D, whereas functional disability and pain are strong predictors.
Conclusion
EQ-5D can be reliably estimated from WOMAC subscale scores without any systematic bias using these results.
Keywords
EQ-5D WOMAC Mapping Osteoarthritis Statistical methodsBackground
Osteoarthritis is one of the leading causes of pain and limitation in function, typically affecting the hip, knee and small hand joints. It is the most common form of arthritis and is characterised by extreme variability in clinical presentation, rate of disease progression and patient outcomes.
One of the most widely reported outcome measures in Knee and Hip osteoarthritis is the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) [1]. WOMAC comprises 24 questions covering the dimensions of pain, stiffness and physical functioning. With substantial evidence of psychometric validity it is in widespread use in clinical studies designed to assess the efficacy of technologies in this disease area.
However, WOMAC is not a tool suitable for use in economic evaluation. In particular, whilst analysts typically seek to evaluate interventions in terms of costs and Quality Adjusted Life Years (QALYs), the scores from the WOMAC instrument are not appropriate since they do not provide either a cardinal or preference based scale. Generic instruments that would facilitate the calculation of benefits in terms of QALYs, such as the EQ-5D, Health Utilities Index Mark 3 and SF-6D, have rarely been included in clinical studies. There is therefore often a gap in the evidence.
One method of bridging this evidence gap is by using statistical methods to estimate the relationship between clinical outcomes, such as WOMAC, and measures suitable for use in economic evaluation, using separate datasets of patients that have completed both instruments. This allows the benefits of treatment to be converted to an appropriate scale for economic evaluation. This practice has been variously labelled “mapping”, “cross-walking” and “transfer to utility”.
However, utility data exhibit peculiar characteristics that raise statistical challenges: they are limited at full health (1) and tend to have a mass of observations at that upper bound, they are limited below at the worst health state described. EQ-5D tends to be bi or tri-modal with a gap between full health and the next feasible health state. Standard regression methods are not suited to data with such distributions and the resulting bias has been demonstrated in many disease areas [2, 3] and is evidenced in the published evidence on WOMAC and osteoarthritis [4, 5]. This bias leads to the systematic undervaluing of the benefits of health care technologies in economic evaluations, and ultimately to inappropriate decisions for patients and their treating clinicians.
New methods have been developed specifically for modelling EQ-5D and been demonstrated to overcome the limitations of standard methods in other disease areas [8]. This paper reports on the application of such methods to modelling EQ-5D from WOMAC and other explanatory variables using data from patients with osteoarthritis of the hip or knee.
Methods
Dataset
A total of 1768 patients were recruited to this registry study conducted in 15 hospitals from three regions of Spain: 3 in Andalusia, 3 in the Canary Islands, and 9 in the Basque Country. The Institutional Review Boards of each hospital approved the study.
Consecutive patients scheduled to undergo primary joint replacement surgery due to knee or hip osteoarthritis in one of the participating hospitals between October 2005 and October 2006 and who received postoperative management in the hospitals were eligible for the study. Patients with cancer or severe organic or psychiatric diseases were excluded because these conditions could prevent them from completing all the questionnaires included in the study.
All patients were sent a letter informing them about the study and asking for their voluntary participation. Questionnaires were mailed to each patient at baseline prior to surgery, 3, 6 and 12 months after surgery. Reminder letters were sent 15 days after each mailing to patients who had not replied promptly. The baseline questionnaire included items about expectations for the operation, along with the SF-12, WOMAC, and EQ-5D-3L questionnaires, plus questions requesting sociodemographic information like gender, age, hospital, weight and height (for the calculation of body mass index [BMI]).
A total of 7072 observations were obtained.
We applied the EQ-5D UK tariff for use in this analysis. The Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) is one of the most widely reported outcome measured in clinical studies of osteoarthritis. It comprises 24 questions with responses to each given on a Likert scale: none (0), mild (1), moderate (2), severe (3) or extreme (4) problems. It is possible to construct a single index score between 0–96 by summing the scores to each question or one may consider subscale scores based on pain (5 questions yielding a range of 0–20), stiffness (2 questions range 0–8) and functional impairment (17 questions range 0–68).
Statistical methods
Health state utility values in general, and the EQ-5D instrument in particular, are characterised by several features that render standard statistical models redundant. Feasible values are limited both above and below at full health and the worst health state described by the instrument. There tend to be a large number of individuals that indicate they are at full health. There is a relatively large gap from full health to the next feasible value (0.883). The typical distribution of EQ-5D tends to be tri-modal with distinctly non-normal elements within each section. In this situation it is obvious that standard multiple regression models, which assume conditional normality, are inappropriate. There now exists a battery of studies that consistently demonstrate this to be the case in a range of disease areas. Linear regression leads to estimates that are consistently biased downwards when considering the value of change in disease level on the EQ-5D scale.
We have previously developed a method specifically designed to produce appropriate estimates of EQ-5D from clinical and other predictors. It has been demonstrated to overcome the limitations of other commonly used direct and indirect methods in several different settings and full details of the model are provided elsewhere [8, 9]. We apply this method here in order to estimate EQ-5D from WOMAC responses, incorporating other socio-demographic and clinical covariate information.
The method applied here uses multiple regression with two novel features. First, we use mixture models which identify a number of different classes within the data which when combined form a new distribution. Whilst the individual classes may estimate EQ-5D as a function of WOMAC and other variables based on standard distributional forms, when these are combined they can form very complex, non standard distributions such as the distribution of EQ-5D. Whilst the relationship between EQ-5D and each of the explanatory variables will differ for each class it is not the identification of different classes that leads us to use this framework, though this may be a valuable insight. It is the flexibility of the modelling framework that makes it particularly useful for modelling EQ-5D values ([10]). The overall estimate of EQ-5D, predicted from any set of WOMAC pain, stiffness and disability scores, is a weighted function of these individual components.
The second novel feature of the analysis is that, in this case, instead of basing each separate class of the mixture on a standard normal distribution, we base it instead on a distribution specific to the characteristics of EQ-5D, namely, limited above at full health (1), below at −0.594 and adjusted to reflect the gap in feasible values between 1 and 0.883.
Explanatory variables may enter the model in two ways: either as predictors of the relationship with EQ-5D within each of the individual classes, as in standard regression, or as predictors of component membership. We compared several different variants of using the explanatory variables in these two ways. Within each class we adopted a consistent approach to the inclusion of variables, opting not to exclude non-significant terms on an individual basis.
Results
Characteristics of the patient sample
Variable | Obs | Mean | SD | Min | Max |
---|---|---|---|---|---|
Age (yrs) | 1593 | 69.14 | 9.81 | 25.64 | 100.04 |
Body Mass Index | 1455 | 29.65 | 4.67 | 14.17 | 50.22 |
EQ-5D | 1520 | 0.29 | 0.35 | −0.594 | 1 |
WOMAC pain scale | 1660 | 11.37 | 3.78 | 0 | 20 |
WOMAC stiffness scale | 1657 | 4.70 | 1.95 | 0 | 8 |
WOMAC function scale | 1664 | 43.85 | 11.83 | 0 | 68 |
Proportion - Osteoporosis of the hip | 1768 | 0.61 | 0.49 | ||
Proportion - Male | 1764 | 0.39 | 0.49 |
WOMAC scale scores indicate that patients are of moderate severity on average at baseline, as expected from a hospital recruited cohort. Similarly mean EQ-5D of 0.29 (sd 0.35) indicates a significant reduction in quality of life with patients spanning the entire feasible range for the scale (−0.594 to 1). This is reflected in Figure 1 which provides a histogram of EQ-5D scores collected over the 4 time points. This illustrates a distribution with the classic characteristics of EQ-5D data: there is a mass of observations at full health (1), a gap between these observations and the rest, and at least two additional modes. In total there are 598 observations, 8.5% of the sample, at this upper limit.
Parameter estimates from the 5 class mixture model
Class 1 | Class 2 | Class 3 | Class 4 | Class 5 | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Parameter | p-value | p-value | p-value | p-value | p-value | |||||||||||
Within class | Pain/10 | 0.0166 | 0.260 | −0.2357 | 0.000 | *** | 0.1475 | 0.000 | *** | 0.001 | 0.285 | 0.0116 | 0.600 | |||
Stiffness/10 | 0.0118 | 0.583 | −0.0287 | 0.760 | −0.0182 | 0.662 | −0.0029 | 0.057 | * | 0.0037 | 0.914 | |||||
Function/10 | −0.0356 | 0.000 | *** | −0.0274 | 0.103 | −0.07 | 0.000 | *** | −0.0003 | 0.087 | * | −0.0403 | 0.000 | *** | ||
Variance | 0.0024 | 0.000 | *** | 0.0406 | 0.000 | *** | 0.0047 | 0.001 | *** | 0 | 0.000 | *** | 0.0058 | 0.002 | *** | |
Random effects terms | Intercept | 0.7636 | 0.000 | *** | 0.3762 | 0.000 | *** | 0.9662 | 0.000 | *** | 0.555 | 0.000 | *** | 0.2118 | 0.000 | *** |
Age^{1} | 0.005 | 0.000 | *** | |||||||||||||
Male | 0.0083 | 0.000 | *** | |||||||||||||
Variance | 0.0016 | 0.000 | *** | |||||||||||||
Class probabilities | Intercept | 8.8821 | 0.000 | *** | 7.0881 | 0.000 | *** | 11.4076 | 0.000 | *** | 5.1497 | 0.0004 | *** | |||
Pain/10 | −5.1772 | 0.006 | *** | −7.9363 | 0.000 | *** | −13.4885 | 0.000 | *** | −6.4654 | 0.0006 | *** | ||||
Stiffness/10 | −1.0227 | 0.625 | −6.1566 | 0.004 | *** | −0.7275 | 0.751 | −1.0918 | 0.6023 | |||||||
Function/10 | −1.3706 | 0.026 | ** | −0.9434 | 0.162 | −2.0612 | 0.002 | *** | 0.37 | 0.5729 | ||||||
(Pain/10)^{2} | 1.2266 | 0.114 | 2.3464 | 0.001 | *** | 4.5149 | 0.000 | *** | 1.4951 | 0.0568 | * | |||||
(Stiffness/10)^{2} | −1.158 | 0.604 | 3.7049 | 0.156 | −1.483 | 0.570 | −1.3521 | 0.5465 | ||||||||
(Function/10)^{2} | 0.0797 | 0.286 | 0.1946 | 0.008 | *** | 0.2545 | 0.001 | *** | −0.0626 | 0.4179 | ||||||
EQ-5D estimates by class | Mean | 0.678 | 0.124 | 0.881 | 0.557 | 0.106 | ||||||||||
Median | 0.679 | 0.129 | 0.865 | 0.557 | 0.107 | |||||||||||
Min | 0.307 | −0.594 | 0.335 | 0.382 | −0.377 | |||||||||||
Max | 1 | 1 | 1 | 0.748 | 0.598 | |||||||||||
Proportion of ones | 0.004 | 0.001 | 0.431 | 0.000 | 0.000 | |||||||||||
Prob of class membership | 0.333 | 0.109 | 0.240 | 0.146 | 0.172 |
Pain, stiffness, function and their quadratic terms were entered as predictors of class membership. Estimates are also shown in Table 2. The probability of class 5 is 1 minus the sum of the probability of being in classes 1 to 4. Pain and function are important predictors of class membership whilst stiffness is statistically significant only as a determinant of class 2 membership.
Differences both in the determinants of class membership and the impact of explanatory variables within each class are partly due to the fact that each class is centred around different elements of the disease severity spectrum. The lower element of Table 2 provides statistics describing the location and spread of each class. Classes 2 and 5 are centred around the most severe patients with mean EQ-5D of approximately 0.1. This peak in the overall distribution of EQ-5D scores is clearly visible in Figure 1. Class 2 has a much greater variance however, spanning the entire feasible range of EQ-5D (−0.594 to 1). Note that because the model is designed bespoke to the features of EQ-5D it does not generate values outside the feasible range. Classes 1, 3 and 4 jointly reflect the peaks in the distribution seen around 0.7 and 1. Class 1 is the largest of the groups and has a mean EQ-5D of 0.68. Class 3 has a mean of 0.88 (the upper limit of EQ-5D values below full health) but is the only class that contributes a substantial proportion of full health observations.
Summary measures of performance
Data | Mixture model | Linear model | |
---|---|---|---|
Mean | 0.5422 | 0.5421 | 0.5402 |
Median | 0.620 | 0.636 | 0.552 |
SD | 0.363 | 0.361 | 0.361 |
Minimum | −0.594 | −0.594 | −0.946 |
Maximum | 1 | 1 | 1.942 |
Proportion of ones | 15.02 | 14.81 | 1.121^{*} |
The figures also illustrate the expected values for the individual classes. There is clear variation in the slopes for each class. Class 4 is the least responsive to changes in disease severity whilst class 2 is the most responsive. Note that this is not simply a reflection of the within class coefficients since patients with greater levels of disease severity on one measure also tend to have greater levels of impairment on other measures. None of the slopes are as steep as for the combined relationship which is formed as the weighted average of the five classes with the weights themselves a function of disease severity. The overall slope at low levels of impairment is formed predominantly from classes 1 and 3, whereas at the most severe levels of impairment it is formed from classes 2 and 5.
Discussion
Whilst several previous studies have been published that look at the relationship between some measure of health utility and the WOMAC instrument, they are all based on linear regression models (Xie et al. [6] include a comparison with the Censored Least Absolute Deviations (CLAD) model). Theoretically these are unlikely to be appropriate, have been shown to be biased in many other disease areas, and also appear to be biased in the osteoarthritis setting [4, 5]. Studies have also been based on relatively small samples of patients.
In this paper we have applied a method based on mixture models that suffers from none of these problems, using data of over 7,000 patient observations that cover the full spectrum of disease. It therefore provides an appropriate method to express treatment benefits from clinical studies which use the WOMAC as an outcome measure in terms of QALYs. This link between clinical and economic outcomes is increasingly important for reimbursement authorities. This is therefore not just of academic interest. Patients and their treating clinicians are directly affected by the methods used to estimate QALYs. To assist those wishing to use these results we provide an online calculator in Excel that provides predicted EQ-5D scores for any user-defined set of WOMAC inputs, age and sex (see Additional file 1).
We find that the stiffness component of WOMAC is relatively unimportant in determining EQ-5D values, consistent with Barton et al. [5] and Grootendorst et al. [7]. Pain and functional disability are much more strongly associated with changes in health utilities but these relationships are different between five classes of patients, identifiable from their disease characteristics. Whether disease was of the hip or the knee did not impact the relationships we identified.
Mitchell et al. [11] recently report a “mapping” study linking to the recently developed ICECAP-0 capabilities measure which, whilst anchored around a 0–1 scale, is not a health utility based instrument. The study only considers linear models on the basis that there are only a small number of observations at the maximum value. Whether the spread of values in the dataset truly represents the spread of severity requires consideration but in any case, the performance of the mixture model approach advocated here in the context of health state utility values should also be considered in the context of other outcome measures such as ICECAP.
Our analysis is based on single sample of Spanish patients receiving a variety of treatments within the broad spectrum of standard practice. In using the results, consideration must be given as to whether this sample includes these for the decision problem at hand and whether the specific treatments under consideration might in some way change the relationship between pain, stiffness, function and health utility.
Conclusions
EQ-5D can be reliably estimated from WOMAC subscale scores without any systematic bias by using the results presented here.
Declarations
Authors’ Affiliations
References
- Bellamy N, Buchanan W, Goldsmith CH, Campbell J, Stitt L: Validation study of WOMAC: a health status instrument for measuring clinically important patient relevant outcomes to antirheumatic drug therapy in patients with osteoarthritis in the knee or the hip. J Rheumatol 1988, 30: 167–178.Google Scholar
- Brazier J, Yang Y, Tsuchiya A: A review of studies mapping (or cross walking) non-preference based measures of health to generic preference-based measures. Eur J Health Econ 2010, 11: 21525.View ArticleGoogle Scholar
- Hernandez Alava M, Wailoo A, Wolfe F, Michaud K: The relationship between EQ-5D, HAQ and pain in patients with rheumatoid arthritis. Rheumatology 2013. doi:10.1093/rheumatology/kes400Google Scholar
- Marshall D, Pericak D, Grootendorst P, Gooch K, Faris P, Frank C, Bellamy N, Torrance G, Feeny D: Validation of a prediction model to estimate health utilities index mark 3 utility scores fromWOMAC index scores in patients with osteoarthritis of the Hip. Value Health 2008, 11(3):470–477. 10.1111/j.1524-4733.2007.00258.xView ArticlePubMedGoogle Scholar
- Barton G, Sach T, Jenkinson C, Avery A, Doherty M, Muir K: Do estimates of cost-utility based on the EQ-5D differ from those based on the mapping of utility scores? Health Qual Life Outcomes 2008, 6(1):51. 10.1186/1477-7525-6-51PubMed CentralView ArticlePubMedGoogle Scholar
- Xie F, Pullenayegum EM, Li S, Hopkins R, Thumboo J, Lo N: Use of a disease-specific instrument in economic evaluations: mapping WOMAC onto the EQ-5D utility index. Value Health 2010, 13(8):873–878. 10.1111/j.1524-4733.2010.00770.xView ArticlePubMedGoogle Scholar
- Grootendorst P, Marshall D, Pericak D, Bellamy N, Feeny D, Torrance GW: A model to estimate health utilities index mark 3Utility scores from WOMAC index scores in patients with osteoarthritis of the knee. J Rheumatol 2007, 34: 534–542.PubMedGoogle Scholar
- Hernandez Alava M, Wailoo AJ, Ara R: Tails from the peak district: adjusted limited dependent variable mixture models of EQ-5D health state utility values. Value Health 2012. doi:10.1016/j.jval.2011.12.014Google Scholar
- Hernandez Alava M, Wailoo A, Wolfe F, Michaud K: “A comparison of direct and indirect methods for the estimation of health utilities from clinical outcomes”. Med Decis Mak in pressGoogle Scholar
- McLachlan GJ, Peel D: Finite Mixture Models. New York: Wiley; 2000.View ArticleGoogle Scholar
- Mitchell PM, Roberts TE, Barton PM, Pollard BS, Coast J: Predicting the ICECAP-O capability index from the WOMAC osteoarthritis index: is mapping onto capability from condition-specific health status questionnaires feasible? Med Decis Mak 2013. doi:10.1177/0272989X12475092Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.