 Research
 Open Access
 Published:
Improved RIDIT statistic approach provides more intuitive and informative interpretation of EQ5D data
Health and Quality of Life Outcomes volume 18, Article number: 63 (2020)
Abstract
Background
EQ5D is generic measure of healthrelated quality of life. Studies using EQ5D generate ordinal data that are interpreted as categories ordered by severity. New analytic approaches taking into account the ordinal nature of the health dimension severity and leading to a better interpretation of EQ5D data are needed to better elucidate differences in healthrelated quality of life. We propose utilizing the Improved RIDIT statistical method to analyze EQ5D outcomes.
Methods
556 Moroccan participants aged over 18 years representing four chronic diseases: back pain (n = 158), renal insufficiency (n = 56), diabetes (n = 82) or hypertension (n = 80) and healthy subjects (n = 180). All participants received the two EQ5D versions. Two other published data sets were included. The first was extracted from a diabetic Spain study and the second was extracted from a clinical trial study. The Improved RIDIT analyses were carried out using an R statistic program we developed.
Results
Applying the Improved RIDIT on the EQ5D data allowed estimating for the first time the ordinal odds, the Absolute Risk Reduction (ARR) or the Absolute Risk Increase (ARI) and the Number Needed to Treat. The ARI values estimated for Moroccan patients showed that (i) hypertension increased anxiety/depression by 66% and reduced mobility by 65%; (ii) back pain increased pain/discomfort by 69%; (iii) renal insufficiency impacts mobility (ARI = 57%, odds_{ordinal} = 9.95) and usual activities (ARI = 44%, odds_{ordinal} = 6.41) and (iv) diabetes acts only on anxiety/depression (ARI = 50%, odds_{ordinal} = 4.8). Also, we demonstrated that the approach works well in clinical trials.
Conclusions
Improved RIDIT provides more intuitive and informative interpretation of the EQ5D data by (1) taking into account the level severity; estimating (2) the odds ordinal, (3) the ARR/ARI and the NNT; (4) analyzing the five dimensions of the EQ5D separately, which gives clinical teams more precision in understanding the treatment/pathology impacts on the health status and completes the EQ5D data analysis based on score utilities.
Background
The EQ5D instrument is a generic measure of patient reported outcomes. It was developed for estimating utility scores based on population preferences for the different health states. It is used in observational studies to evaluate healthrelated quality of life in patients with different pathologies [1,2,3,4,5], to estimate normative values of health states or so called score utilities in different populations [6,7,8,9,10,11,12] and in clinical trials to describe treatment efficacy [13, 14].
The EQ5D exists in two versions called EQ5D3 L and EQ5D5 L and recently the childfriendly EQ5D version (EQ5DY) was introduced [15, 16]. The EQ5D consists of five dimensions (Mobility, Selfcare, Usual Activities, Pain/Discomfort and Anxiety/Depression). Each one has three (EQ5D3 L) or five (EQ5D5 L) levels of problems representing the degree of the health state severity. This descriptive system is a set of a number of health states, each coded by a number of 5 digits. The three levels version (EQ5D3 L) contains a set of 3^{5} = 243 health states, whereas the 5 L version contains 5^{5} = 3125 health states. To summarize the EQ5D responses the individual response is converted into a utility score. Also, the EQ5D consists of a Visual Analog Scale (VAS) to assess patients’ current health status (scale 0–100, where 0 = worse imaginable health state and 100 = best imaginable heath state). The VAS has been established as valid and reliable in a range of clinical and research applications [17, 18].
The EQ5D instrument has been expressly developed to yield two forms of data: score utilities and VAS values. Appropriate analysis of EQ5D data and consequent interpretation of results is of vital importance to the field. In studies whose objective is to compute qualityadjusted life years or in medical studies where the aim is to compare the quality of life between two groups of subjects (for example, healthy and ill; before and after receiving a treatment), the score utilities and the VAS values are most often used. When utility scores are missing, it would be acceptable to apply another country’s value set to estimate utilities. Utility scores and VAS are univariate continuous variables and the statistical inference can be done through the common Student or MannWhitney statistical tests. These statistics fail to reveal the difference between two samples of subjects at each dimension separately since the utilities and VAS measure the global health state by grouping the 5 dimensions in a single value. Measurement on the 5 distinct dimensions of the EQ5D can also be informative in order to determine which dimensions of the EQ5D are affected by disease or improved by treatment. Some studies [1, 2, 4, 14] have compared the impact of a treatment or disease on the quality of life by comparing the frequencies of “level 1” obtained from the two samples for each of the five dimensions. A chisquare test is often used to compare these frequencies. For example, in a clinical trial, the positive impact of treatment on the health status of patients is significant if the frequency of the “level 1” corresponding to “no problems” is observed to be higher after treatment [14]. Healthrelated quality of life studies using the EQ5D generate outcome ordinal data that are interpreted as categories ordered by severity. Nevertheless, the Chisquare test fails to incorporate the ordered nature of the random variable even when used at each dimension separately. By dichotomizing the responses into two categories: “no problem: level 1” and “other problems: the remaining levels together”, the Chisquare masks the importance of the severity of the levels, which is the fundamental property of EQ5D.
New analytic approaches taking into account the ordinal nature of the EQ5D data are needed to better elucidate differences in health related quality of life. Bross developed a nonparametric test RIDIT (Relative to an Identified Distribution) to analyze ordered data [19]. The mean RIDIT comparison of two samples from two conditions estimates the probability that a randomly selected subject from one sample will have lower health state than a randomly selected subject from the other one. RIDIT analysis works well when the reference population is known, as this analysis can lead to conflicting conclusions depending on the group used as reference. In some cases notably when the sample sizes are very different, interchanging the groups can affect the standard deviation and the statistical result. To overcome this limitation, Flora proposed an “Improved RIDIT” approach [20] based on the MannWhitney test [21] and theoretical bases of Conover [22]. Contrary to Bross’ RIDIT analysis, the Improved RIDIT does not require a reference group. The main difference between RIDIT and Improved RIDIT relates to the probabilities under consideration. The Bross method considers only the probability that a randomly selected subject from a test group is in more serious condition than a randomly selected subject from the reference population group. The Improved RIDIT procedure considers three probabilities: more serious, as serious as and less serious condition.
To our knowledge, the Improved RIDIT has not yet been used to analyze EQ5D data. The aim of our study is to apply the Improved RIDIT approach in the context of ordered data of healthrelated quality of life outcomes from EQ5D. We demonstrate the efficacy and advantages of the Improved RIDIT analytical tool on 3 sets of EQ5D data.
The improved RIDIT statistic method comes with an added value for the EQ5D questionnaire by revealing information on health status changes for each EQ5D dimension separately. Utility scores yield information on the change in overall health status. Together, utilities and improved RIDIT can provide total information on the change of healthrelated quality of life.
Methods
Data collection
We applied the Improved RIDIT method on three sets of EQ5D data. The first data set (Set 1) was extracted from the Spanish study on diabetes by Mateo et al. [23]. The second data set (Set 2) was extracted from clinical trial study of Devlin et al. [14]. The third set (Set 3) were data collected in our study representing five Moroccan citizens groups: a sample of healthy subjects and 4 samples of chronic diseases (back pain, hypertension, renal insufficiency and diabetes).
We chose Set 1 because the authors used the EQ5D5 L version to evaluate the healthrelated quality of life of individuals with diabetes, a chronic disease that we incorporated in our data set (Set 3). We chose Set 2 for two reasons: (i) in their study, Devlin et al. [14] used the EQ5D3 L version, which is the case also in our study (Set 3); (ii) Devlin et al. evaluated in their clinical trial the positive impact of medical treatment on the quality of life of patients. Devlin et al. data set (Set 2) enabled us to test the effectiveness of the Improved RIDIT method to evaluate the health quality of life in clinical trials.
The 556 subjects associated with the Moroccan data (Set 3) were recruited from the unit of hypertensive treatment (Settat, Morocco), hemodialysis center (Temara, Morocco), Association de Soutien des Diabetiques de Settat (Settat, Morocco) and University Hospital of Rabat (Rabat, Morocco). All subjects had one of the four following diagnoses: back pain (n = 158), renal insufficiency (n = 56), diabetes (n = 82) or hypertension (n = 80). The control sample consisted of healthy subjects (n = 180) working in the clinical care setting. Only individuals who have none of the 4 pathologies studied are included in the healthy sample. All study participants were over 18 years of age and signed an informed consent form after receiving information on study objectives and the instrument used to collect data. For the back pain, hypertension and diabetes samples patients with hypertension/renal insufficiency/diabetes, back pain/renal insufficiency/diabetes and back pain/renal insufficiency/hypertension were excluded, respectively. For the renal insufficiency sample patients with hypertension or back pain were excluded. For the renal insufficiency group, 78 patients gave their consent from a total of 90 patients followed in the hemodialysis center. 56/78 met the exclusion criteria. For the diabetes and the hypertension samples, 82/135 and 80/127 patients gave their consent and met the exclusion criteria, respectively. For the control and back pain samples, 3 interviewers collected the data for two months (one day per week) with a frequency of 10 participants per interviewer per day. All participants were selected randomly with no repeated measurement. In total, we had the consent of N = 3 * 8 * 10 = 240 participants. For the back pain sample 158 patients met the exclusion criteria. For the control group 180 had no chronic pathology. The ethics committee of the Hassan 1st University (Settat, Morocco) provided approval of the study protocol.
Participants’ responses to the EQ5D were collected according to the Scalone protocol [3]. Subjects from the clinical diagnostic groups completed the survey in the clinic setting; control subjects completed the survey at their place of employment. All subjects reported their age and gender and received the Arabic EQ5D5 L version followed by the Arabic 3 L version and finally the visual analogue scale (VAS) in accordance with previous research findings, which showed that when respondents scored the 3 L first, there could be a tendency to not use the “inbetween” levels 2 and 4 of the 5 L [24]. All statistical analyses were carried out using an R statistic program we developed.
Statistical analysis
The algorithm of the improved RIDIT method
Consider two randomly independent samples Y (healthy control or receiving a treatment in clinical trial) and X (patients with chronic disease or receiving a treatment of reference/placebo) from a study of healthrelated quality of life. The data obtained for each dimension of EQ5D was summarized in a contingency table where rows refer to the Y and X samples and columns refer to the EQ5D levels (3 and 5 columns for EQ5D3 L and EQ5D3 L, respectively) (Table 1). For each dimension of the EQ5D data, the categorical random variable is ordered from the least serious level (“No Problem” is coded “1”) to the most serious level (“Extreme Problem” is coded “5” for the 5 L version).
Let y_{j} (j = 1, …, n_{y}) denote the n_{y} observations of the reference sample (Y) and x_{i} (i = 1, …, n_{x}) denote the n_{x} observations of the comparison sample (X).
Define the function ω as:
From the function ω, we define the following probabilities:
π^{+}, π^{0} and π^{−} indicate that subjects randomly selected from the X sample are in a worse, in a same or in a better health state, respectively than subjects randomly selected from the Y sample.
The Improved RIDIT test is based on the statistic:
With
From Eqs. 1, 2, 3 and 4, we obtain:
Then an estimate of (π^{+} − π^{−}) is:
From (Eq. 1) and Table 1, we rewrite (Eq. 5) as:
The variance of W was given by Flora [20]:
From (eq.3), an estimate of π^{0} is:
From the probability propriety: π^{0} + π^{−} + π^{+} = 1 and (Eq. 9), estimates of π^{+} and π^{−} are respectively:
Therefore, estimation of π^{−} and π^{+} requires calculation of W from Eq. 10, which is not trivial. However, we propose the following simple algorithm:

1.
Transform the contingency table (Table 1) into a square matrix ω_{(k, k)} whose elements are the products of the frequencies \( {f}_X^{(l)} \) and \( {f}_Y^{(m)} \) (l, m = 1, …, k) and k is the level number of the ordinal variable.
$$ \boldsymbol{\omega} =\left(\begin{array}{ccccc}{f}_Y^{(1)}{f}_X^{(1)}& {f}_Y^{(1)}{f}_X^{(2)}& \cdots & {f}_Y^{(1)}{f}_X^{\left(k1\right)}& {f}_Y^{(1)}{f}_X^{(k)}\\ {}{f}_Y^{(2)}{f}_X^{(1)}& {f}_Y^{(2)}{f}_X^{(2)}& \cdots & {f}_Y^{(2)}{f}_X^{\left(k1\right)}& {f}_Y^{(2)}{f}_X^{(k)}\\ {}\vdots & \vdots & \ddots & \vdots & \vdots \\ {}{f}_Y^{\left(k1\right)}{f}_X^{(1)}& {f}_Y^{\left(k1\right)}{f}_X^{(2)}& \cdots & {f}_Y^{\left(k1\right)}{f}_X^{\left(k1\right)}& {f}_Y^{\left(k1\right)}{f}_X^{(k)}\\ {}{f}_Y^{(k)}{f}_X^{(1)}& {f}_Y^{(k)}{f}_X^{(2)}& \cdots & {f}_Y^{(k)}{f}_X^{\left(k1\right)}& {f}_Y^{(k)}{f}_X^{(k)}\end{array}\right) $$
Looking at the ω_{(k, k)} matrix, we notice that:

The sum of all entries in ω_{(k, k)} is equal to n_{x}n_{y}.

The sum of all entries above the main diagonal is equivalent to the left term in (eq. 10). These values \( {f}_Y^{(m)}\ {f}_X^{(l)}\kern0.50em \left(l>m\right) \) indicate that a randomly subject selected from X is in a high level of the ordinal variable (worse heath state in the EQ5D dimension) than a randomly subject selected from Y.

The sum of all entries below the main diagonal is equivalent to the right term in eq. 10. These entries \( {f}_Y^{(m)}\ {f}_X^{(l)}\kern0.50em \left(l<m\right) \) indicate that a randomly selected subject from X is in a low level of the ordinal variable (better heath state in the EQ5D dimension) than a randomly selected subject from Y.

The sum of the diagonal elements (the trace of the ω_{(k, k)} matrix) divided by n_{x}n_{y} is equivalent to (Eq. 12).

2.
Estimate π^{0} by \( {\hat{\pi}}^0= Tr\left(\boldsymbol{\omega} \right)/{n}_x{n}_y \) where Tr is the trace.

3.
Construct the W_{(k, k)} matrix by multiplying the diagonal elements of the ω_{(k, k)} matrix by zero (x_{i} = y_{j} from Eq. 1), the elements above the diagonal of ω_{(k, k)} by 1 (x_{i} > y_{j} from Eq. 1) and the elements below the diagonal of ω_{(k, k)} by − 1 (x_{i} < y_{j} from Eq. 1):
$$ \boldsymbol{W}=\left(\begin{array}{ccccc}0& {f}_Y^{(1)}{f}_X^{(2)}& \cdots & {f}_Y^{(1)}{f}_X^{\left(k1\right)}& {f}_Y^{(1)}{f}_X^{(k)}\\ {}{f}_Y^{(2)}{f}_X^{(1)}& 0& \cdots & {f}_Y^{(2)}{f}_X^{\left(k1\right)}& {f}_Y^{(2)}{f}_X^{(k)}\\ {}\vdots & \vdots & \ddots & \vdots & \vdots \\ {}{f}_Y^{\left(k1\right)}{f}_X^{(1)}& {f}_Y^{\left(k1\right)}{f}_X^{(2)}& \cdots & 0& {f}_Y^{\left(k1\right)}{f}_X^{(k)}\\ {}{f}_Y^{(k)}{f}_X^{(1)}& {f}_Y^{(k)}{f}_X^{(2)}& \cdots & {f}_Y^{(k)}{f}_X^{\left(k1\right)}& 0\end{array}\right) $$ 
4.
Calculate W (Eq. 10) as the sum of all the entries of the W_{(k, k)} matrix.

5.
Estimate \( {\hat{\pi}}^{+} \) as the sum of all the entries above the diagonal of the W_{(k, k)} matrix divided by n_{x}n_{y}.

6.
Estimate \( {\hat{\pi}}^{} \) as the absolute value of the sum of all the entries below the diagonal of the W_{(k, k)} matrix divided by n_{x}n_{y}.
Notice that \( {\hat{\pi}}^{},{\hat{\pi}}^{+},{\hat{\pi}}^0 \) and W estimated in the steps 2 and 4–6 verify the Eqs (13) and (14).
The improved RIDIT statistic
Once the quantities π^{+} and π^{−} are estimated, the Improved RIDIT method leads to study the inference between two randomly samples X and Y under two experimental conditions by three statistics:
The standard distribution test
If the two samples Y and X came from the same population so the null hypothesis H0: π^{+} = π^{−}. Consequently, the expected value of W is zero. Alternatively, if X and Y come from two different populations then the alternative hypothesis H1: π^{+} ≠ π^{−}. For example, let X be a sample of persons suffering of back pain and Y a sample of healthy persons and consider the dimension pain/discomfort in EQ5D. Testing the hypothesis that these two populations are equivalent with respect to pain/discomfort severity is equivalent to testing the null hypothesis H0 (π^{+} = π^{−}) against the alternative one H1 (π^{+} ≠ π^{−}). H1 indicates that persons with back pain will be more seriously in a state of pain/discomfort than healthy persons. Conover [22] demonstrated that for large samples X and Y (sample size greater than 30, which is a general rule of thumb for the large sample condition), W is approximately normally distributed. Therefore, the test statistic for the hypothesis that the two populations are the same regarding the level of severity is constructed as the following:
The null hypothesis is rejected in favor to the alternative one if the Z value is greater than the threshold value of the standard distribution at the significance level α. In the context of EQ5D data, the Improved RIDIT performing the simultaneous comparability of the five EQ5D dimensions. However, these simultaneous tests are not independent since they share the same data. In this case of multiple testing, a common strategy is to adjust the significance level as a function of the number of tests we are running. In order to reduce the possibility of getting false positives, we used the BonferroniHolm procedure [25] to adjusting the pvalues calculated by the Improved RIDIT.
The absolute risk reduction (ARR) and needed number to treat (NNT)
Considering the case of a study of the impact of treatment on healthrelated quality of life where X and Y are two randomly independent samples representing Control/or placebo and Treatment groups, respectively. If the treatment has a positive impact on the quality of life, then π^{−} = P(Y < X) indicates the probability that a random individual from the Treatment group has a better outcome than a random individual from the Control/placebo group and π^{+} = P(Y > X) indicates the probability that a random subject from the Control/placebo group has a preferable outcome than a random individual from the Treatment group. For the ordinal variable, Shepstone [26] showed that P(Y < X)  P(Y > X) is equivalent to the Absolute Risk Reduction (ARR). Walters estimated the ARR in the case of ordinal data by the use of the MWW U statistics. Walters proposed ARR = (U_{YX} − U_{XY})/n_{x}n_{y} where U_{YX} and U_{XY} are the U statistics and n_{x} and n_{y} are the sample sizes of X and Y [27]. Walters discussed the limitation of this estimation. If there are no ties in the data (π^{0} = P(X = Y) = 0) then the ARR can be estimated exactly. But in the presence of a large amount of ties (π^{0} > 0) in the data, which is the case of EQ5D outcomes, the ARR using the MWW approach can only be estimated approximately. Hence, approaches leading to exact estimation of P(Y < X) and P(Y > X) (π^{−} and π^{+}, respectively) are welcome. Applying the Improved RIDIT to the healthrelated data, the ARR can easily be estimated by:
Also, applying the Improved RIDIT to healthrelated quality of life data allows estimating the Needed Number to Treat (NNT). In clinical trials it is very important to estimate the NNT to evaluate the efficacy of a treatment. The NNT is defined as the number of patients who would have to receive a treatment to prevent one additional bad outcome [28]. It is estimated by NNT = 100/ARR (where ARR is %) or NNT = 1/ARR (where ARR is proportion).
In the case of studies concerning the evaluation of the negative effect of pathologies on the healthrelated quality of life, we can define the Absolute Risk Increase (ARI) similarly to the ARR in the clinical trials studies. The ARI is estimated by:
where π^{+} is the probability that a randomly patient from the disease group has a worse outcome than a randomly selected healthy subject from the control group.
Both ARR and ARI are positive estimates of the impact of a treatment or a disease on the quality of life. The confidence interval for the ARR/ARI at the confidence level (1α) is done by:
The impact of the treatment/disease is statistically significant if the confidence interval of ARR/ARI does not include zero.
The odds ordinal
Agresti [29] has proposed a generalization of the odds to two independent ordinal variables X and Y as:
Applying the Improved RIDIT to the EQ5D data obtained from two randomly independent samples X and Y, we can easily estimate the ordinal odds by:
Considering one of the five dimensions of the EQ5D, odds_{ordinal} means that patients are odds_{ordinal} times in worse health dimension than the healthy ones.
Similarly, for clinical trials if X and Y are two independent samples of subjects receiving randomly a reference treatment/placebo and a treatment of interest, respectively the odds is estimated by:
which means that patients are odds_{ordinal} times in better health dimension after receiving treatment.
Simulation
EQ5D data generation
To generate EQ5D outcomes that are likely to occur in practice, we used the exponential distribution f(x) = λe^{−λx}. Different values of the rate λ (Fig. 1) were used to generate different shapes of frequencies for the two EQ5D5 L and EQ5D3 L versions according to the following procedure:

For each value of λ, 10,000 observations have been generated.

We replaced each observation by its integer value plus one.

For the simulation of 5 L data, if an observation is greater than 5 it is replaced by 5. For 3 L data, if an observation is greater than 3 it is replaced by 3.
Figure 1 shows the decrease of the level “1” (no problems) frequency (from Fig. 1a to Fig. 1c) and the increase of the other level frequencies. This illustrates typically what could be occurring in patient outcomes when the disease affects the healthrelated quality of life. Interpreting Fig. 1 in the strand sense (from Fig. 1c to Fig. 1a), we observe improvement of the patient outcomes (increasing of the “no problem” level and decreasing of the other levels). This indicates a positive effect of a medical treatment on the quality of life.
Comparison between the improved RIDIT and MannWhitneyWilcoxon (MWW)
To compare the two methods, the computer simulation was based on the comparison between the false positive rate and the power estimated by each of the two approaches. This algorithm simulation was based on the analysis of false positive rate and power estimation.
False positive rate estimation

1.
For each of the three profiles (Fig. 1a, b and c), draw with replacement two independent random samples X and Y of size n_{x} and n_{y}, respectively. The sample sizes used were 30, 50, 100, 200 and 300. In theses simulations, we consider situations where the distributions have the same shape (the null hypothesis is true).

2.
The test statistic, Improved RIDIT or MWW, was calculated between the X and Y samples. If the calculated test is lower than the critical value of the test statistic for a nominal α, a false positive was recorded. This step was carried out for all the combinations (n_{x}, n_{y}) sizes.

3.
Steps 1 and 2 were repeated S times. We estimated the false positive rate (\( \hat{FP}\Big) \) by the proportion of false positives among the S repetitions (in all simulations, S = 10,000).
Power estimation
Consider situations where the two distributions to be compared do not have the same shape. For this, suppose the distribution profile in Fig. 1a representing a control/treatment group (sample X) and the profiles in Fig. 1b and c representing a disease group (sample Y).

1.
Draw two independent random samples X and Y (X from Fig. 1a and y from Fig. 1b or Fig. 1c) of size n_{x} (30, 50, 100, 200 or 300) and n_{y} (30, 50, 100, 200 or 300), respectively.

2.
For each combination (n_{x}, n_{y}), the Improved RIDIT and MWW tests were calculated. If the test value is greater or equal than the critical value of the test statistic for a nominal α, a success is recorded (the alternative hypothesis is true).

3.
Steps 1 and 2 were repeated S times. We estimated the power (\( \hat{P}\Big) \) by the proportion of successes among the S repetitions (in all simulations, S = 10,000).
On the use of the odds ordinal and the ARR/ARI
The ARR/ARI and the odds_{ordinal} are connected by the relationship:
\( {odds}_{ordinal}=1+\frac{ARR}{\pi^{}} \), in the case of comparing a treatment group to a control group.
\( {odds}_{ordinal}=1+\frac{ARI}{\pi^{+}} \), in the case of comparing a healthy group to a control group.
To analyze deeply this relationship, we simulated values of π^{−} and π^{+} for fixed values of ARR/ARI. At first the ARR/ARI values were generated according to a normal distribution. At a second time, at each simulation 10, 000 values of π^{−} and π^{+} were generated from a uniform distribution. After, we retained only the values π^{−} and π^{+} whose difference was equal to the ARR/ARI values generated in the first step.
Results
Simulation results
Simulation result of the comparison between MWW and improved RIDIT
Table 2 summarizes comparisons of the false positive rates obtained by Improved RIDIT and MWW over the three selected distributions depicted in Fig. 1 for each of the 25 sample size combinations. These results demonstrated the robustness of the Improved RIDIT and MWW tests since all the false positive rate values obtained by the two methods did not greatly exceed the nominal level of significance (α = 0.05) (Table 2). In addition, for the three distributions and all the sample sizes examined, the Improved RIDIT always estimated less false positives than MWW. In contrast, MWW some times gave falsely significant results, i.e., the false positive rate exceeded 5% (Table 2).
Investigation of the power was based on the comparison between two simulated EQ5D dimension distributions like as those presented in Fig. 1. We simulated two randomly independent samples X and Y whose distributions were similar to Fig. 1a and Fig. 1c, respectively. Analysis of the power values estimated by MWW and Improved RIDIT (Fig. 2) demonstrated that the power increases with the sample size. The power of the two tests was greater than 0.8 except in the case where the two samples had a size equal to 30. If there are fewer subjects in one of the two samples (size < 50), we observed a small difference between the power values of the MMW and improved RIDIT tests. This difference vanished for sample sizes greater than 50.
Simulation result on the use of the ordinal odds and ARR/ARI
Figure 3 shows odds_{ordinal} versus π^{−} obtained from two simulations corresponding to ARR = 0.5 and ARR = 0.7. We observed that the lower the risk in the control group (π^{−}), the larger the difference between odds_{ordinal} and ARR. Also, the odds_{ordinal} overestimated the risk when π^{−} is too low (< 10%).
Data set 1
We applied the Improved RIDIT statistics (Table 3) on EQ5D5 L data used to evaluate 1857 individuals (1016 women and 841 men) with diabetes living in Spain [23]. The study results revealed that diabetes affects the health related quality of life outcomes of subjects, affecting the health of women more than that of man by acting on all five dimensions of health. Findings from this study were based on comparisons between two frequencies f1 (obtained for “level 1”) and f2 (obtained as the sum of the four others), which masks information on the severity of each level. By applying the Improved RIDIT method on the same data set, the likelihood that diabetic women are in a worse state than diabetic men was estimated by three statistics (Table 3): the standard statistic test (Wtest), The ARI and the ordinal odds. In term of risk difference, the ARI values showed that diabetic women have an absolute risk increase of 22% in mobility, 16% in self care, 23% in usual activities, 29% in pain discomfort and 20% in anxiety depression than men. The odds ordinal values showed similar result demonstrating that diabetes affect the mobility, self care, usual activities, pain discomfort and anxiety for about 2 times in women than in men (Table 3). Also, the difference between diabetic women and men on healthrelated quality was demonstrated by the Wtest. Indeed, the BonferroniHolm adjusted pvalues calculated for the five dimensions (Table 3) were highly significant.
Data set 2
We also tested the Improved RIDIT method on the healthrelated quality of life data obtained through the EQ5D3 L version administered by Devlin et al. [14]. In contrast to Mateo et al. [23] analyzing the effect of pathology on the quality of life of patients, Devlin et al. evaluated the impact of a treatment (Enzalutamide) on the quality of life of patients aged ≥ 18 with castrationresistant prostate cancer. In total, 1717 patients were randomized to receive Enzalutamide (n = 872) or placebo (n = 845) [14].
Using the chisquare test, they showed that the effect of enzalutamide on healthrelated quality of life was significant in the pain/discomfort and anxiety/depression dimensions at week 13. Also, they find that enzalutamide improved the usual activities and pain/discomfort at week 25. There were no significant betweengroup differences in the mobility and selfcare dimensions (Table 4). All these analyses were based on comparison of the percentages of level 1 “no problems” to the other levels, which does not take into account the ordinal nature of the data. Applying the Improved RIDIT method on Devlin’s data and using the BonferroniHolm adjusting for multiple testing (Table 4), we found that pain/discomfort (P_{BH}value = 0.007) and anxiety/depression (P_{BH}value = 0.013) dimensions were improved at week 13. At week 25, a difference was observed in pain/discomfort after adjusting pvalues (P_{BH}value = 0.033) but not for usual activities (P_{BH}value = 0.124) (Table 4). The ARR statistics showed that enzalutamide decreased significantly the pain/discomfort by 8,9% and by 9,1% at week 13 and week 25, respectively. Also, the treatment reduced the anxiety/depression by 7,5% at week 13. The usual activities dimension was in the limit of signification at week 13 (pvalue = 0.058), which was confirmed by the BonferroniHolm adjusting (P_{BH}value = 0.175). Enzalutamide, however, increased the usual activities at week 25 by ARR = 6,2% but not significantly due to the adjusted pvalue (P_{BH}value = 0.124). In terms of odds ordinal ratio, the positive effect of enzalutamide on pain/discomfort was estimated by 1.42 times and 1.44 times at week 13 and week 25, respectively (Table 4). At week 13, the patients receiving enzalutamide were 1.43 times less anxious.
Data set 3
Table 5 shows the characteristics of the study participants. For the 5 samples, females comprise 55–60% of participants in all groups other than hypertension for which 86% of the group is female. The mean age was 36 years (SD = 10 years) for the control group, 38 years (SD = 10 years) for back pain, 52 years (SD = 16 years) for diabetic patients, 61 years (SD = 11 years) for hypertensive patients and 55 years (SD = 17 years) for people with renal insufficiency. Patients with chronic hypertension or renal insufficiency were the oldest with 58.7 and 51.8% with an age greater than 60 years. The response frequencies across 3 and 5 levels of EQ5D dimensions are presented in Table S1 (Supplementary Material).
Analyses of the 5 L questionnaire data by the Improved RIDIT statistic are presented in Table 6. Considering the control sample as a reference, the absolute risk increase and the ordinal odds values for each disease and each health dimension were estimated. Also, the Wtest and the associated BonferroniHolm pvalues corrections were calculated (Table 6). The ARI values show that hypertension reduced mobility by 65%, selfcare by 34% and usual activities by 32%. This pathology increased pain/discomfort by 49% and anxiety/depression by 66%. Regarding the ordinal odds values, hypertensive patients were 12.71 times less mobile, 14.42 times less autonomous and 4.24 times less active than healthy individuals. They suffer from pain 4 times and anxiety 8 times higher than the healthy ones. These findings were consolidated by the Wtest where all the BonferroniHolm adjusted values were much smaller than the significance level.
For individuals with back pain, the adjusted pvalues for multiple testing demonstrated that the healthrelated quality of life was significantly affected in patients (Table 6). The most affected dimension was pain/discomfort with an absolute risk increase of 69% and ordinal odds of 12.11. The four other dimensions were also affected from 27 to 36% (Table 6).
Renal insufficiency impacts the quality of life of patients by acting on mobility (ARI = 57%, odds_{ordinal} = 9.95 and P_{BH} < 0.0001), usual activities (ARI = 44%, odds_{ordinal} = 6.41 and P_{BH} < 0.0001), selfcare (ARI = 35%, odds_{ordinal} = 6.96 and P_{BH} < 0.0001) and pain/discomfort (ARI = 37%, odds_{ordinal} = 3.12 and P_{BH} < 0.0001). Anxiety/depression was the less affected dimension in patients with renal insufficiency (ARI = 27%, odds_{ordinal} = 2.45 and P_{BH} < 0.001).
Diabetic patients showed values of ARI and odds_{ordinal} lower than those estimated for individuals with hypertension, renal insufficiency or back pain. The only significant values were obtained for anxiety/depression (ARI = 50%, odds_{ordinal} = 4.8 and P_{BH} < 0.0001) and mobility (ARI = 22%, odds_{ordinal} = 3.1 and P_{BH} = 0.0006) dimensions.
Similar results were obtained from the Improved RIDIT analysis of the EQ5D3 L data of the four diseases in comparison to the healthy group. These findings demonstrated the effectiveness of the Improved RIDIT method.
Discussion
Researchers interested in understanding health related quality of life data adopt different strategies of analyzing data from the EQ5D instrument. Most of these strategies are based on the use of the utility scores. The utility scoresbased strategies provide information on the global health changes between two groups (ill/healthy, before/after treatment). Our findings suggest that using the Improved RIDIT approach provides a more detailed understanding of a health intervention or impact of pathology on quality of life outcomes for each dimension separately.
One example of EQ5D analysis is to summarize frequencies corresponding to the responses obtained for each of the ordered severity levels. The chisquare test is often applied for statistical inference but does not take into consideration the ordinal nature of the severity of the questionnaire levels. The MWW test is an adequate test for ordinal data. Comparison between the Improved RIDIT and MWW tests using computer simulation showed that the improved RIDIT had less false positive rate than MWW for all sample sizes studied (Table 2). For sample sizes greater than 50, our simulation found that the two tests had similar power (Fig. 2). In addition, the MWW does not allow estimating the ARR/ARI and the NNT that the improved RIDIT provides.
Devlin et al. developed an approach based on the Pareto principle [14]. In this approach, an EQ5D state of health is considered better than another if it is better in at least in one dimension and remains equivalent or better for the other dimensions. Conversely, a state of health is considered worse than another if it is worse at least in one dimension and it is not better in the remaining dimensions. On the basis of this principle to compare health states between two groups of patients, Devlin et al. found that following a hip replacement, 82% of patients had improved health, less than 5% had no change, less than 5% had worse health, and under 10% had a ‘mixed’ change [14].
Precise results were obtained when we applied the Improved RIDIT analyses on Devlin et al.’s data [14]. Devlin et al.’s approach makes it possible to evaluate the benefit of the treatment by considering the state of health in all five dimensions at the same time. The Improved RIDIT analysis allowed us to analyze the 5 dimensions separately, which makes it possible to analyze the effect of the treatment on each of the five dimensions. Our analyses revealed that the treatment resulted in significant improvement in pain/discomfort, anxiety/depression and usual activities. Understanding a specific impact of a treatment helps clinical teams better articulate potential quality of life outcomes to patients.
The Improved RIDIT statistic allows estimating the ARR/ARI and odds_{Ordinal}, which are not estimated by the Pareto approach used by Devlin et al. Calculations using Pareto are based on the comparisons of the health status. When the EQ5D3 L is used there are 243 × 243 possible comparisons. This number is larger with the EQ5D5 L system (3125 × 3125). Pareto allows information about the global health and the improved RIDIT allows more information per dimension. The two approaches are complementary and can be used synergistically to analyze the EQ5D data for evaluating the impact of treatment or disease on health outcomes.
By applying the Improved RIDIT method on diabetic Moroccan patients and the diabetic sample from Mateo’s study [23], we found similar findings. Diabetic Moroccan women were in a worse health state than diabetic men. The improved RIDIT added the information that significant differences between Moroccan diabetic women and men were observed for pain/discomfort (ARI = 0.33, odds_{Ordinal} = 3.06 and) and depression/anxiety (ARI = 0.43 and odds_{Ordinal} = 3.38) dimensions.
We also studied the relationship between the ordinal odds ratio and the absolute risk reduction/the absolute risk increase (ARR/ARI). The different analyzes of the EQ5D data showed that in some cases the values of odds_{Ordinal} were high compared to their corresponding ARR/ARI. This can be explained, in these cases, by the small values of π^{−} (Fig. 3). We suggest to use only the ARR when π^{−} < 0.1, otherwise both ARR and odds_{ordinal} can be used to evaluate the impact of treatment or disease on the healthrelated quality of life.
Analyzing the EQ5D data by the improved RIDIT allowed also estimating the NNT. For example, application of the approach on the data of Devlin [14], showed that the ARR of the Enzalutamide at week 13 was 7.5% (Table 4) for the anxiety dimension. It means that one have to treat NNT = 13 patients to prevent one additional bad anxiety outcome.
To summarize all the information provided by the improved RIDIT method, we suggest using this statistic procedure to conduct analyses of the EQ5D healthrelated quality of life outcomes from two randomly independent samples following the scheme in Fig. 4. The improved RIDIT method showed advantages of its use in the case of EQ5D data obtained from two independent random samples by (i) calculating the Wtest and taking into account for multiple testing (ii) estimating the ARR/ARI and the associated NNT and (iii) estimating the ordinal odds. It is of interest to extend the improved RIDIT approach to paired samples and for more than two samples.
The ARR/ARI, NNT and ordinal odds estimation are not allowed by either the MWW test or the utilitybased approaches. Providing depth analysis of the EQ5D data, the Improved RIDIT can be used to complete the information derived from utility scores.
Conclusion
We utilized the Improved RIDIT statistic to examine the results of the EQ5D surveys in two previously published studies and a data set we collected from Moroccans with chronic health conditions compared to a group of healthy individuals. This statistic procedure takes into account the severity level, which permits to estimate the absolute risk reduction/absolute risk increase and the ordinal odds ratio. The improved RIDIT analysis permits us to analyze the five dimensions of the EQ5D separately, which gives clinical teams more precision in understanding the effect of a pathology or treatment on the health status.
Availability of data and materials
All data generated or analyzed during this study are included in this article.
Abbreviations
 ARI:

Absolute Risk Increase
 ARR:

Absolute Risk Reduction
 MWW:

MannWhitneyWilcoxon
 NNT:

Number Needed To Treat
References
 1.
Safita N, Islam SM, Chow CK, Niessen L, Lechner A, Holle R, Laxy M. The impact of type 2 diabetes on health related quality of life in Bangladesh: results from a matched study comparing treated cases with nondiabetic controls. Health Qual Life Outcomes. 2016;14(1):129–37.
 2.
Feng Y, Devlin N, Herdman M. Assessing the health of the general population in England: how do the three and fivelevel versions of EQ5D compare? Health Qual Life Outcomes. 2015;13:171–86.
 3.
Scalone L, Ciampichini R, Fagiuoli S, Gardini I, Fusco F, Gaeta L, et al. Comparing the performance of the standard EQ5D 3L with the new version EQ5D 5L in patients with chronic hepatic diseases. Qual Life Res. 2013;22(7):1707–16.
 4.
Tran B, Ohinmaa A, Nguyen L. Quality of life profile and psychometric properties of the EQ5D5 L in HIV/AIDS patients. Health Qual Life Outcomes. 2012;10(1):132–40.
 5.
Jia YX, Cui FQ, Li L, Zhang DL, Zhang GM, Wang FZ, et al. Comparison between the EQ5D5 L and the EQ5D3 L in patients with hepatitis B. Qual Life Res. 2014;23(8):2355–63.
 6.
Peeters Y, Stiggelbout AM. Health state valuations of patients and the general public analytically compared: a metaanalytical comparison of patient and population health state utilities. Value Health. 2010;13(2):306–9.
 7.
Sullivan PW, Lawrence WF, Ghushchyan V. A national catalog of preference based scores for chronic conditions in the United States. Med Care. 2005;43:736–49.
 8.
Devlin N, Hansen P, Herbison P. Variations in selfreported health status: results from a New Zealand survey. N Z Med J. 2000;113:517–20.
 9.
Hinz A, Kohlmann T, StöbelRichter Y, Zenger M, Brähler E. The quality of life questionnaire EQ5D5L: psychometric properties and normative values for the general German population. Qual Life Res. 2014;23:443–7.
 10.
Ferreira LN, Ferreira PL, Pereira LN, Oppe M. EQ5D Portuguese population norms. Qual Life Res. 2014;23(2):425–30.
 11.
Devine A, Taylor SJ, Spencer A, DiazOrdaz K, Eldridge S, Underwood M. The agreement between proxy and selfcompleted EQ5D for care home residents was better for index scores than individual domains. J Clin Epidemiol. 2014;67(9):1035–43.
 12.
Liu H, Cella D, Gershon R, Shen J, Morales LS, Riley W, et al. Representativeness of the patientreported outcomes measurement information system internet panel. J Clin Epidemiol. 2010;63(11):1169–78.
 13.
Hakimi Z, Herdman M, Pavesi M, Devlin N, Nazir J, Hoyle C, Odeyemi IA. Using EQ5D3L and OAB5D to assess changes in the healthrelated quality of life of men with lower urinary tract symptoms associated with benign prostatic hyperplasia. Qual Life Res. 2017;26(5):1187–95.
 14.
Devlin N, Herdman M, Pavesi M, Phung D, Naidoo S, Beer TM, Tombal B, et al. Healthrelated quality of life effects of enzalutamide in patients with metastatic castrationresistant prostate cancer: an indepth post hoc analysis of EQ5D data from the PREVAIL trial. Health Qual Life Outcomes. 2017;15(1):130.
 15.
Otto C, Barthel D, Klasen F, Nolte S, Rose M, Meyrose AK, et al. Predictors of selfreported healthrelated quality of life according to the EQ5DY in chronically ill children and adolescents with asthma, diabetes, and juvenile arthritis: longitudinal results. Qual Life Res. 2018;27(4):879–90.
 16.
Scott D, Ferguson GD, Jelsma J. The use of the EQ5DY health related quality of life outcome measure in children in the Western cape, South Africa: psychometric properties, feasibility and usefulness  a longitudinal, analytical study. Health Qual Life Outcomes. 2017;15(1):12.
 17.
Kildesø J, Wyon D, Skov T, Schneider T. Visual analogue scales for detecting changes in symptoms of the sick building syndrome in an intervention study. Scand J Work Environ Health. 1999;25(4):361–7.
 18.
Wewers ME, Lowe NK. A critical review of visual analogue scales in the measurement of clinical phenomena. Res Nurs Health. 1990;13(4):227–36.
 19.
Bross DJ. How to use Ridit analysis. Biometrics. 1958;14:18–38.
 20.
Flora JD. A Note on Ridit Analysis. Technical report No. UMHSRISA7410. 1974; Highway Safety Research Institute, University of Michigan. 1974.
 21.
Mann HB, Whitney DR. On a test of whether one of two random variables is stochastically larger than the other. Ann Math Stat. 1947;18(1):50–60.
 22.
Conover WJ. Rank tests for one sample, two samples and K samples without the assumptions of a continuous distribution function. Ann Stat. 1973;1:1105–25.
 23.
Mateo DC, Gordillo MG, Olivares PR, Adsuar JC. Normative values of EQ5D5L for diabetes patients from Spain. Nutr Hosp. 2015;32(4):1595–602.
 24.
Janssen MF, Birnie E, Haagsma JA, Bonsel GJ. Comparing the standard EQ5D threelevel system with a five level version. Value Health. 2008;11(2):275–84.
 25.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6(2):65–70.
 26.
Shepstone L. Reconceptualising and generalising the absolute risk difference: a unification of effect sizes, odds ratios and numberneededtotreat. J Epidemiol Community Health. 2001;55:A7.
 27.
Walters SJ. Sample size and power estimation for studies with health related quality of life outcomes: a comparison of four methods using the SF36. Health Qual Life Outcomes. 2004;2:26.
 28.
Wen L, Badgett R, Cornell J. Number needed to treat: a descriptor for weighing therapeutic options. Am J Health Syst Pharm. 2005;62(19):2031–6.
 29.
Agresti A. Generalized odds ratios for ordinal data. Biometrics. 1980;36(1):59–67.
Acknowledgments
We would like to thank all study participants, clinicians and pharmacy staff for their contribution to this study. We would like also to thank the EuroQol group for authorizing us to use the EQ5D questionnaire.
Funding
The study was conducted without funding.
Author information
Affiliations
Contributions
MA, YM I, developed the study concept and aims, and implemented study protocol. MA developed the statistical methodology. SE, EY, NC and HA contributed to the analysis of data. MA, YMI and JT drafted the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Approval provided from the ethics committee of the Hassan 1st University (Settat, Morocco).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1 Table S1.
Sample response frequencies (%) across 3 and 5 levels of EQ5D dimensions from healthrelated quality of life outcomes of healthy subjects and sick subjects
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Marfak, A., YoulyouzMarfak, I., El Achhab, Y. et al. Improved RIDIT statistic approach provides more intuitive and informative interpretation of EQ5D data. Health Qual Life Outcomes 18, 63 (2020). https://doi.org/10.1186/s12955020013133
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12955020013133
Keywords
 Improved RIDIT
 Statistical analysis
 EQ5D5 L
 EQ5D3 L
 Absolute risk reduction
 Ordinal odds
 Number needed to treat