 Research article
 Open Access
 Open Peer Review
 Published:
A preferencebased item response theory model to measure health: concept and mathematics of the multiattribute preference response model
BMC Medical Research Methodology volume 18, Article number: 62 (2018)
Abstract
Background
A new patientreported health measurement model has been developed to quantify descriptions of health states. Known as the multiattribute preference response (MAPR) model, it is based on item response theory. The response task in the MAPR is for a patient to judge whether hypothetical healthstate descriptions are better or worse than his/her own health status.
Methods
In its most simple form MAPR is a Rasch model where for each respondent on the same unidimensional health scale values are estimated of their own health status and values of the hypothetical comparator health states. These values reflect the quality or severity of the health states. Alternatively, the respondents are offered healthstate descriptions that are based on a classification system (e.g., multiattribute) with a fixed number of health attributes, each with a limited number of levels. In the latter variant, the weights of the levels of the attributes in the descriptive system, which represents the range of the health states, are estimated. The results of a small empirical study are presented to illustrate the procedures of the MAPR model and possible extensions of the model are discussed.
Results
The small study that we conducted to illustrate the procedure and results of our proposed method to measure the quality of health states and patients’ own health status showed confirming results.
Conclusions
This paper introduces the typical MAPR model and shows how it extends the basic Rasch model with a regression function for the attributes of the healthstate classification system.
Background
Health is a sociocultural construct encompassing a wide range of phenomena, so it is not surprising that various actors define it differently. Traditionally, physicians have been guided by a biomedical model and have thus understood health predominantly as a condition that falls within acceptable biological norms. Nowadays, there is an increased awareness of the impact of health and health care on the quality of human life. The conventional clinical healthstatus construct is now often extended to psychological and even social factors, thereby making subjective measures such as (perceived) health status or ‘quality of life’ necessary — and rightly so, because the ultimate goal of all health interventions is to improve a patient’s perceived health condition. The use of these subjective measures has proliferated ever since the World Health Organization published its definition of health in 1946 [1].
There are several ways to express health. We can compile a ‘snapshot’ of a patient’s current health condition from an ‘image bank’ comprised of health states. These health states consist of discrete health attributes (e.g., domains, dimensions, items) each with a number of levels. When combined, they represent a description of a person’s health status or healthrelated quality of life (HRQoL) [2]. Subsequently, such healthstate descriptions can be measured (valued) by assigning meaningful numbers (values) to an individual’s health state. ‘Meaningful’ is here defined as values that reflect the patients’ health status in relationship to other health states. This is different from subjective measures (e.g., visual analogue scale) that reflect the perception of how individuals experience their health status in relationship to their own internal standards. It is convenient to express individuals’ health in single metric values, as these can be used in health outcomes research, for clinical monitoring of the health status of patient groups, and in particular, in disease modeling studies and economic evaluations.
To obtain healthstate values (variously called preferences, utilities, index, or weights), the healthstate descriptions must be quantified in terms of seriousness or quality. Differences between health states values are assumed to correspond to increments of quality differences between these states, which implies that the values are on an intervallevel scale [2]. Most conventional methods of measurement (or valuation) stem from health economics (e.g., standard gamble, time tradeoff) and are susceptible to many disturbing factors such as adaptation, time preference, context, reference point, and other biases [3,4,5]. To control for adaptation, which occurs in most of these conventional methods (especially for chronically ill patients), all economic valuation methods use hypothetical health states that are assessed by a sample of (unaffected) members of the general population. However, it is reasonable to assume that healthy people are not adequately informed or lack the imagination to appropriately judge the impact of health states, particularly severe ones [6, 7].
A new way to quantify health states was recently introduced. This measurement method, the multiattribute preference response (MAPR) model, is based on the Rasch model (an item response theory model) [8]. The MAPR model more or less mimics the situation of a patient with a certain health condition lying in a wardroom where the other occupants have (related) complaints and symptoms. This patient is asked to compare his own health state to that of his roommates by indicating whether his own state is better or worse. The conventional preferencebased measurement methods usually yield an opinion on health states from healthy controls, while the result of the MAPR is an internal positioning of a patient’s health status with respect to other health states. The response mechanism of the MAPR model is less susceptible to various biases that conventional methods are prone to. Moreover, the MAPR is the first generic health preference model that is fully based on patient perception and reporting; as such it is a genuine patientreported outcome measure. Apart from being grounded in a renowned measurement theory, the MAPR response tasks are attractive and easy to perform in a selfcompletion setting.
This article introduces and explains the MAPR model conceptually and mathematically. The first section looks into the background of its measurement mechanism, namely the Rasch model, and expands on its operation in a health setting. The second section describes the MAPR model; the third works through its estimation procedures. Finally, the results of a small empirical study are presented to illustrate the procedures of the MAPR model and possible extensions of the model are discussed.
Methods
Measurement mechanism
A probabilistic measurement model was invented by the Danish mathematician Georg Rasch. While primarily employed to assess educational attainment, it is increasingly used for other purposes [9]. Its original setting was the field of reading skills, where it was intended for use with dichotomous response data (e.g., correct/wrong). Nowadays, the Rasch model or the closely related oneparameter logistic model (OPLM) is considered a variant of the class of item response theory (IRT) models [9, 10]. The Rasch model is built around the idea that the probability of a correct response to an item is modeled as a logistic function of the difference between the difficulty of an item (parameterized by β) and the characteristics of a person (e.g., ability parameter θ):
The Rasch model poses three stringent requirements. The first is unidimensionality: a unique onedimensional latent variable explains the response to the items. The second is monotonicity: the probability of a positive response to an item is a nondecreasing function of the latent variable. And the third is local independence: for any given individual, the item responses are independent conditional on where individuals are on the underlying latent scale.
Under the Rasch model, a Guttman scale is the most likely response pattern for a person when items are ordered from least difficult to most difficult [11]. This means that if someone responds correctly to an item, then that person should succeed on all easier items; conversely, if one responds incorrectly, then he/she should fail on all items that are more difficult (Fig. 1). Unlike the Guttman scale, the Rasch model is a probabilistic model. In the latter, the probability that any person will succeed on an easier item will always be greater than the probability of success on a more difficult item. The Guttman scale is the deterministic limiting case of the Rasch model.
Health context
In the context of health measurement assuming the Rasch model implies that the more positive the difference between the value (the perceived quality) of the health status of a patient (θ) and the value of another health state (β) to be judged, the higher the probability that the patient will indicate that his/her current health status is better than the presented health state. Or the other way around, patients in very poor health will consider many other health states as better than their own. Using the Rasch model, one can estimate the health status of individual patients (i.e., their ability, in Rasch terminology) and the value of the hypothetical health states (i.e., difficulty of the parameters of items) on the same latent scale. In short, patients are asked to respond to hypothetical health states by comparing these with their own health status. For example, “Is this health state better than your own health state?”
In the Rasch model, patients compare their own health status with a few prescribed hypothetical health states. These comparator health states can span the whole continuum from bad to mild (as done in this article in our small empirical study), but they can also denote health states that are closely positioned on the latent scale to the actual health status of the individual patient. Such comparator health states may be based on holistic descriptions or objects. Holistic refers to unstructured verbal descriptions or objects such as people’s faces or skin in photos. In general, holistic objects are often extremely easy to compare and judge. However, features (attributes) by which to describe the object specifically are often absent. Alternatively, health descriptions may be derived from a classification system with multiple attributes, whereby each attribute has a limited number of levels (Fig. 2). The latter approach enables the investigator to predict values for health states that are not part of the empirical study (see below).
Let θ_{ p } be the (unknown) health status of person p (p = 1, ..., P). Suppose that β_{ i } is the (unknown) value of health state i (i = 1, …, I) as measured on a latent scale. Imposing the logistic function of the difference between a person’s health status and the values of the comparator health states on the probability that a person prefers his/her own health state over the comparator description leads to the Rasch model. More formally, let Y_{ ip } be a random variable with a value of one if a person prefers his/her health status over the hypothetical health state and zero otherwise. In this way it is assumed that the health status θ_{ p } of person p is on the same latent scale as the health states i with β_{ i } and that a person will most likely prefer his/her own health status over health state i if θ_{ p } exceeds β_{ i }. Under the Rasch model we assume that
or equivalently in the logit form
where η = log (π/(1 − π)) is the logit link function. This means that if a person’s health status θ_{ p } is equal in quality to the hypothetical comparator health state, so there is no preference difference for either state, the probability of choosing the one over the other is fiftyfifty. Also, the further apart the person’s health state is from the comparator, the larger the probability that the better state is preferred and chosen. In the following, model (1) will be called the holistic MAPR model.
The holistic MAPR model, like the original Rasch, is a descriptive model. It describes the individual patient’s health state (e.g., by localizing patients on the health scale) and the value of the judged comparator states without explaining either of these by characteristics of the patients or the health states. The holistic MAPR is both feasible and attractive in many clinical situations where characteristics cannot be easily discerned, such as body and skin deformations. Comparing and assessing pictures or movies may then be more appropriate. A crucial requirement is that the respondents should be located along the whole health scale from very severe to almost perfect health; otherwise, the model cannot be sufficiently estimated. The typical response task of the MAPR model precludes responses from a sample of the general population. The latter are predominantly healthy and therefore do not provide the information needed to estimate the model.
MAPR model
Several simple classification systems have been developed to capture the major features of health such that they can be used to describe health states. Each system transposes those features into a certain number of health attributes. The health state can then be measured with a discrete response scale for each attribute at a certain number of levels. For example, the classification system of the preferencebased EQ5D3 L instrument consists of five health attributes: mobility, selfcare, usual activities, pain/discomfort, and anxiety/depression, with a value of 1 (best), 2, or 3 (worst) for each attribute [12]. In this way an EQ5D3 L health state can be represented by five digits, with 11111 denoting perfect health and 33333 the worst possible condition. The threelevel version of the EQ5D3 L system defines 3^{5} = 243 possible partially ordered different health states. The SF6D healthstate classification contains six attributes, namely physical functioning, role limitation, social functioning, pain, mental health, and vitality. With response categories ranging from four to six levels, the SF6D can describe 18,000 different health states. Some other classification systems are the Health Utilities Index version 3 (HUI3), 15D, Assessment of Quality of Life (AQoL), and the Quality of WellBeing scale (QWB) [13,14,15,16].
Formal representation
Assume that we now have a classification system wherein a health state is represented as a vector x_{ i } = (x_{i1}, ..., x_{ iJ }) with discrete levels on each of the J attributes. The number of levels in the jth attribute is denoted by n_{ j }, so on attribute j the possible values are 1,2, …, n_{ j }. In this way the vector (1,1,…,1) represents perfect health and (n_{ 1 }, n_{ 2 }, …, n_{ J }) the worst state. Suppose that the value β_{ i } of health state x_{ i } can be described as a function β_{ ι } = f (x_{ i }) to reflect the partial ordering of the health states. In the literature several functions have been proposed to model the value of health states as a function of a set of health attributes. For instance, the simple additive linear model assumes that linearity is present in each attribute and that the value drops by the same amount, for example when moving from level 1 to 2 or from level 2 to 3. A less restrictive and more realistic model can be obtained by taking each attribute as a categorical variable in the regression model:
where d_{ jk }(x_{ ij }) is a dummy variable with d_{ jk }(x_{ ij }) = 1 if x_{ ij } = k and zero otherwise. The contribution to the value β_{ i } of health state x_{ i } of a change in attribute j from level 1 to k is parameterized by α_{ jk }. Notice that the regression equation in (3) has no intercept as this parameter is redundant. Furthermore, additional restrictions on α_{ jks } are required for enforcing the partial ordering on the β s and for identifying the parameters.^{Footnote 1} Substituting linear expression (3) for β_{ i } in the logistic expression (2) gives
The parameterization of the value of health states is not limited to the main effects of the health attributes, as interactions between health attributes can be incorporated in (4) by adding products of (dummies of) health attributes. For identification purposes, the number of parameters should be less than the number of health states that the respondents are asked to compare. In the IRT literature this type of (item/health explanatory) model is called the linear logistic test model (LLTM). It was originally proposed by Scheiblechner [17] and later formalized by Fisher [10, 18,19,20]. LLTM differs from the Rasch model in that the influence of the quality/severity of the comparator health states is reduced to a linear combination of a fixed number of healthstate attributes or interactions between those attributes, with fewer parameters than hypothetical health states. The effects of the attributes and their levels on the health states are estimated instead of the holistic healthstate parameters themselves (Formula 1). Being more restrictive than the Rasch model, it enables one to predict values for the complete set of health states that can be constructed for a specific classification system, so predictions can also be made for health states that are not evaluated in the response study.
Suppose we have a sample of n patients who compared the same m health states β_{ i }(i = 1,...,m) with their own health status. By substituting the parameterization of the items in terms of their attributes as formulated in (3) into formula (1), the holistic MAPR model, we can write the probability of response of patient p on health state i as:
Model (5) will be denoted as the MAPR model. Estimation of the health state parameters of the MAPR model now boils down to estimation of the parameters α_{ jk }. In that way, the value of a health state is reflected in the characteristics of the health states as parameterized with the variables d_{ jk }(x_{ i }).
Adaptive MAPR model
A more adaptive approach is possible. Patients are thereby asked to complete a multiattribute classification (e.g., EQ5D3 L) in advance to classify their own health status, denoted \( {\overset{\sim }{\boldsymbol{x}}}_p=\left({\overset{\sim }{x}}_{p1},..,{\overset{\sim }{x}}_{pJ}\right) \). Then, to perform the MAPR response task, they are confronted with a set of (individualized) hypothetical (comparator) health states that were selected in light of the classification of the patients’ own health state from the first task (Fig. 3). So, in this case patients are shown different subsets of health states, depending on\( {\overset{\sim }{\boldsymbol{x}}}_p \). In principle, this approach allows more precise estimation of the position of the patients’ health status. It also precludes selecting a restricted set of predetermined comparator states to be judged. However, it complicates the analysis of the data, as the subset of presented health states differs between the respondents and depends on the person’s own health state, which is restricted to \( {\theta}_p=f\left({\overset{\sim }{\boldsymbol{x}}}_p\right) \). This adaptive operation of the MAPR model is almost similar to computerized adaptive testing (CAT) that is used for standard IRT models. The difference is that for standard IRT models a routine on a central server determines, from a large item bank of candidate items, the next item offered to an individual respondent. For the MAPR model a simple routine as part of a mobile application (www.healthsnapp.info) determines the comparator states (comprising multiple attributes/items) to be assessed by individual patients.
Estimation of the Rasch model
When assessing health states holistically (i.e., no parameters for the levels of the attributes) as in traditional item response theory, it is assumed that the responses to health states are independent of one another, which gives rise to the following likelihood:
The parameters of the standard Rasch model can easily be estimated by several methods, e.g. full maximum likelihood estimation, conditional maximum likelihood and marginal maximum likelihood. All of these are based on maximum likelihood estimation or Bayesian estimation, and several procedures have been described in the literature [21].We will describe now the conditional maximum likelihood (CML) estimation.
Let \( {R}_p=\sum \limits_{i=1}^I{Y}_{ip} \) be the number of health states that a patient p has compared to his own and were considered worse. This number is a sufficient statistic for estimating the patient’s own health state θ_{ p }. Thus, the conditional likelihood of the responses is independent of θ_{ p } if we condition on R_{ p }. This leads to the (CML) estimation equations, after maximizing the likelihood:
P (Y_{ ip } = 1R_{ p } = r_{ p }, β_{ i }) is the probability that the patient’s health status is better than health state i, given the number of health states found to be worse than the patient’s health state. These I1 equations can be solved using a NewtonRaphson procedure leading to consistent point estimates for the healthstate parameters β_{ i }.
An estimate of the patient’s own health state θ_{ p } can be obtained with a maximum likelihood estimation procedure. In this second step, the conditional maximum likelihood estimates of β_{ i } are assumed to be fixed and are substituted in the estimation Eq. (6). In this way the uncertainty associated to these estimates is not accounted for. One way to incorporate this uncertainty could be to use a Bayesian estimation method. In that case a sample from the posterior distributions of the item person parameters can be used instead of imputing only the estimates itself [22].
The variance of the ML estimates equals:
Note that the estimated βs are not incorporated in this variance. The maximum of the function f (x) = x (1x) is 0.25 for x = 0.5. One can thus see that individual health status can be estimated more precisely when patients have to compare health states that are close to their own state.
For the one parameter logistic model (OPLM), the parameter estimates obtained using CML and marginal maximum likelihood (MML) are usually close. The advantage of CML over the MML procedure is that no a priori assumptions have to be made about a person’s healthstate distribution. When this a priori distribution is misspecified, the MML estimates may be biased. It is expected that the distribution of person’s health states is not normally distributed but typically skewed to the right [23]. On the other hand, it has to be underlined that CML estimation also has some pitfalls, such as the fact that individuals with perfect or zero scores do not provide any information and, missing observations can lead to biases in case of missing not completely at random.
Whether a Rasch model fits the data, thereby yielding a unidimensional scale, can be tested with Andersen’s likelihood ratio test [24]. Note that obeying a Rasch model is a sufficient but not a necessary condition of unidimensionality.
Estimation of the MAPR model
Estimation of the LLTM model is similar to estimation of the Rasch model. Both procedures are based on the fact that the number of worse health states per person is a minimal sufficient statistic for θ_{ p }. As a consequence, the parameters α_{ jk } can be estimated without knowledge of the patient’s health status (known as personfree item assessment). Instead, finding the values for β_{ I, } that maximize the (conditional) likelihood estimation of the LLTM model, boils down to finding the values for α_{ jk }.
Both the existence and uniqueness of the CML estimates depend on whether the data matrix is well conditioned. A response matrix is said to be well conditioned if in every possible partition of the health states into two nonempty subsets some patients have given a response of one on some health state in the first set and a response of zero on some health state in the second set [25, 26].
The fit of the MAPR model (LLTM) can be compared with the fit of the Rasch model by using a likelihood ratio test. The deviance of − 2 loglikelihood of the two nested models is approximately Χ^{2} distributed with df = difference between the number of parameters in the two models [18, 27]. When this test is significant, there is evidence that health states are not sufficiently described by the characteristics of the health states as parameterized with the variables d_{ jk }(x_{ i }). In case there is no statistically significant difference between the Rasch model and the MAPR model, the latter can be used to describe the values of health states. Different formulations of LLTM models can also be compared by performing a likelihood test, as long as these models are nested.
Empirical study
Respondents
The aim of the empirical study was to show first exploratory results in testing the MAPR model. In order to do so we used data from a previously published study that aimed to explore discrepancies in values for health states between the general population and patients that experience specific illness [7]. For this study we used only the data of the patients (n = 75).
Two patient groups from the Radboud University Nijmegen Medical Center (Netherlands) participated in that study, which was approved by the Central Committee on Research Involving Human Subjects (region ArnhemNijmegen) [7]. One group included patients who were diagnosed with cancer within a time frame of 4–6 weeks before they participated in the study. The other group consisted of chronically ill patients living with the symptoms of rheumatoid arthritis (RA) for at least 3 years. The study protocol was administered facetoface by a trained interviewer at the homes of the patients.
This initial sample was extended by including patients with a cerebrovascular accident (CVA) or inflammatory bowel disease (IBD) from the hospital Medisch Spectrum Twente (n = 35) and patients with liver disease or paraplegia from the University Medical Center Groningen (n = 53). The Medical Ethics Review Committee Twente (METC/14124) and Medical Ethics Review Committee UMCG (METC 2015/496) declared that this latter part of the research did not fall under the Medical Research Involving Human Subjects Act.
Study design
In the initial study (Radboud) the judgmental task consisted of ranking 17 EQ5D3 L health states, supplemented with the patient’s own EQ5D3 L description, ‘dead’, and state ‘11111’. Each patient ranked the same 20 health states by putting the card with the ‘best’ health state on top and the ‘worst’ at the bottom. Additionally, the patients unknowingly assessed their own health status in the judgmental task, as their own EQ5D3 L healthstate description had been incorporated in the set, but they did not assess the health states of the other participants. The task in the other two studies was slightly different (patients did not assess their own health status), but is not likely affecting the results in the empirical study as described in this article. Respondents in the latter two studies were asked to compare the same 17 EQ5D3 L health states from the Radboud study with their own health (not explicitly represented in terms of the EQ5D3 L description) and express if the EQ5D3 L health states was worse or better than their own health status. In all three studies, the EQ5D3 L health states were presented in random order to control for potential biases due to presentation order or respondent fatigue.
Analysis of the empirical study
First, we fit the Rasch model to the ranking data. Next, we analyze the following (MAPR) model for the value of health state β_{ i }:
a model with only main effects for all attributes (with dummy variables). To ensure identification of the parameters α_{ jk }, an additional restriction has to be put on these parameters; in this case we choose α_{ jk } = 0 for k = 1, j = 1, . . , 5.
Goodness of fit for the holistic MAPR (i.e., the Rasch) model is tested with the Andersen LR test [28]. Then, MAPR model (8) and the Rasch (i.e., the holistic MAPR) model are compared (LR test, correlation coefficient). Next, for every patient the predicted value of its health state following from the estimated MAPR model (8) are calculated based on the patient’s own EQ5D3 L description. For every health state shown to the patient, it is determined whether the patient’s estimated healthstatus value outperforms (i.e., is preferred by the patient) the estimated value of the shown comparator health state. These predicted preferences will then be compared with the observed preferences using a kappa coefficient as measure of agreement. A kappa larger than 0.75 is considered excellent and between 0.4 and 0.75 fair to good [29]. The eRM package in R was used to estimate the MAPR models (LLTM model) and the Rasch model [30].
Results
In total 163 patients were interviewed for this study. Of these, 48 were cancer patients (34 colorectal cancer, 14 breast cancer) and 42 had a liverrelated disease or transplant. The number of participating RA patients was 27. The mean age differs across the participating hospitals, with the oldest patients coming from the Radboud Medical Center (Table 1). Overall, some or major problems were reported for pain (60.1%) and the least problems were reported for selfcare (17.8%). Major or severe problems for selfcare and mood were reported only by patients with liverrelated disease or transplant or by paraplegic patients (Table 2). As the distribution of the health states across the study sites shows, the UMCG had more patients with a severe health condition. But there was a reasonable spread over the whole HRQoL continuum for the three hospitals (Additional file 1).
The Guttman scalogram reveals that not all health states and persons are perfectly ordered (Fig. 4), this can be seen from the green dots between the red ones that indicate misfit. Given the small number of health states in relation to the small number of patients, this study showed that the Rasch (holistic MAPR) model does not hold on statistical grounds. However, after deleting health states in the analysis that were rather severe and therefore overly judged as worse than the own health conditions of the patients (states: 32211, 33323, 32223, 11133, 32313, 22222, 33333, and 23232) the holistic model does hold. An Andersen LRtest showed a log likelihood value of 7.21 with 8 dfs (p = 0.514). The order of the health states based on their sum score is similar to the order based on the estimates of the Rasch model. This result is as expected since the sum score is a sufficient statistic for the Rasch model. The PersonItem Map shows the distribution of patients’ own health status (the above histogram) compared to the assessed health states, see the histogram below (Fig. 5). This figure shows that more than half of the judged comparator health states were assessed as worse than the patient’s health status.
The estimated regression coefficients for MAPR model (8) reveal logical differences at all levels (Table 3). Some problems with selfcare have the highest impact, followed by some problems with mood, pain, mobility, and usual activities. Severe problems with mood and pain have more impact than the other attributes. Estimates of the health states under the MAPR model (8) give almost the same order as for the Rasch model (Table 4). For the MAPR model, the pairs (11211, 21111), (11131, 11113), and (23232, 32223) have a different order and the estimated value of health state 33323 is much smaller than in the Rasch model.
When comparing the conditional likelihood for the Rasch model and the MAPR model, we found a statistical difference (LR statistic = 87.9; df = 6; p < 0.001). This means that the goodness of fit of the MAPR model is lower than for the Rasch (the holistic MAPR) model. However, the correlation between the item parameters as estimated with the Rasch model and the item parameters of MAPR model is 0.93, so even the elaborated MAPR model performs well in explaining the item parameters. In 88.2% of the comparisons, the observed preferences agree with the predicted preferences based on the MAPR model. The kappa coefficient equals 0.71 (CI: 0.68–0.74), which is considered fair to good.
Discussion
This article presents a novel approach to measuring health: the multiattribute preference response model (MAPR). It was developed to quantify health states and patients’ own health status on the same unidimensional scale. The response mechanism of this model is insensitive to various biases (e.g., time preference, risk aversion, indifference procedure) that arise with conventional methods (i.e., standard gamble, time tradeoff) to derive values for health states. Moreover, this is the first generic health preferencebased model that fully reflects perception and reporting by patients. Besides being grounded in measurement theory, the response tasks are attractive and easy to perform in a selfcompletion setting.
The small study that we conducted to illustrate the procedure and results of our proposed method to measure the quality of health states and patients’ own health status showed confirming results. Although the sample size of our empirical study was very modest for performing an item response theory analysis as done here, the estimated regression weights showed a clear logical structure. Values for the small set of health states included in this study could be computed and showed a valid order of the health states as well as interpretable distances between the health states, compared with results from previous large studies based on conventional measurement methods. Note that in this small study a fixed set of only 17 health states were used and therefore we could not include interaction terms in the regression equation. In the full operational MAPR model patients will not be confronted with a fixed set of states, but with a smaller set of healthstate descriptions that are closely similar to their own health status. This will lead to more efficient and robust estimation of the parameters.
The MAPR model largely eliminates unwanted mechanisms affecting valuations of health states. Prominent among these is adaptation. Healthstate values derived by conventional methods are typically higher when elicited from patients, particularly those with chronic illness or disability, than from nonpatients who only imagine themselves in hypothetical health conditions. Adaptation is manifest in almost all standard methods of health measurement, particularly in multidomain instruments, often based on Likert scales as developed within the setting of classical test theory [31]. Moreover, conventional methods for valuing health states stemming from economics (e.g., standard gamble, time tradeoff) are also complex and require abstract reasoning skills. These drawbacks can now be averted. Measurement with the MAPR model is based on a discrimination principle: a patient’s own health status serves as a comparator state against other (comparator) states. This indirect approach to derive values for health states is different from the conventional valuation techniques used by health economists. These valuation techniques (e.g., standard gamble, time tradeoff) request a direct and absolute score (monadic measurement). Because the response task in the MAPR model is simply a preference (rank order) between a patient’s own health status (that serves as a reference standard) and a (closely) related hypothetical health state, the assessment is less likely affected by ‘subjective’ motives and easier to accomplish. Patients don’t quantity their own health status, they only compare it and rank it. This mode of measurement largely prevents biases such as adaptation and coping. From a theoretical and a practical point of view, the MAPR is more attractive than the existing valuation methods, particularly because both the judgmental task and the analysis are executed within one unifying framework.
A downside of the MAPR model is that it produces relative positions of health states. For application in DALYs and QALYs, however, MAPRderived values need to be rescaled around the position where states are considered to become worse than dead (position of dead = 0) [32, 33]. In conventional valuation methods, ‘dead’ is not only an element of the task itself but also introduces many methodological and practical problems. Separate exercises are needed to localize the position of ‘dead’ for the MAPR model. However, recent studies suggest promising solutions; separate studies can be conducted to localize the juncture where health states are considered worse than dead [34, 35]. Such additional studies probably are better be worked out based on the input from a sample of the general population, instead of patients.
When applying any conventional Rasch model to derive metric measures, it is assumed that the underlying phenomena can be represented on a unidimensional scale. However, this crucial assumption may be questionable when quantifying a subjective phenomenon such as health, a construct with a rather broad scope. Nevertheless, the overall assumption is that health outcomes such as health status, healthrelated quality of life, and wellbeing are unidimensional concepts. Of course, this is true of all data to some extent. As many researchers have convincingly argued, unidimensionality does not imply only one factor or dimension. Rather, it implies the presence of a dominant dimension and possibly of minor dimensions that do not affect the dominant one and the unidimensionality of the model is therefore a reflection of the assumed unidimensionality of the majority of assessments we use [36]. Our health model is comparable to widely accepted models of intelligence. Typically, cognitive abilities are represented as a threelevel hierarchy with numerous narrow factors at the bottom, a handful of broad, more general factors at the intermediate level, and at the apex a single factor, the g factor, which stands for the variance common to all cognitive tasks [37].
The MAPR model can even be extended to offer respondents a large set of candidate attributes (far more than the traditional four to nine attributes in existing instruments). An individual patient could then select those most relevant to his or her assessment. By breaking the fixedset mold, this MAPR variant leads to a truly patientcentered preferencebased health measurement approach. An extended MAPR model would most likely require thousands and thousands of responses from patients. Some solutions for this problem have already been introduced [22, 38, 39].
An alternative method to quantify health states may be the discrete choice (DC) model [40, 41]. It is based on (paired) comparisons of two or more hypothetical health states (and not a person’s health status itself). In that sense, the difference between the MAPR model and the DC model seems only minor, but in fact it is significant. The DC model only scales health states, not respondents (patients).
Several elements related to the MAPR model must be investigated empirically to confirm the assumptions underpinning it and explore its potential limitations. In particular, the data should show a Guttman structure. As data for the MAPR model is collected in patient groups, suboptimal response data may result. This may be due to problems with interpreting the health attributes and their levels, from taking cognitive shortcuts and other factors.
Conclusions
This new patientreported health IRT model can be used as a coherent measurement method and has a profound connection to measurement theories. Apart from developing instruments that can be used in medical settings, the MAPR model can also be used to develop healthoutcome instruments to measure in health care. Operated by dedicated data collection technology with interactive routines, data capturing based on this new measurement method becomes simple and even certain distinct patient populations can be easily approached (e.g., children, elderly). In principle, the MAPR model may be suitable to measure other unidimensional phenomena such as wellbeing, capabilities, and other subjective attributes that are essentially based on quality [42].
Notes
 1.
For example, suppose that two attributes take values 1 to 3. Thus x_{ i } = (x_{i1}, x_{i2}) ∈ {1, 2, 3} × {1, 2, 3}. Since x_{2}=(1,2) is worse than x_{1}=(1,1), we obtain an ordering in the alphas. Due to the ordering we have β_{1} < β_{2}, which implies β_{1} = α_{1, 1} + α_{2, 1} < α_{1, 1} + α_{2, 2} = β_{2} and thus \( {\alpha}_{j,{k}_1}<{\alpha}_{j,{k}_2} \)when k_{1} < k_{2}.
Abbreviations
 CML:

Conditional Maximum Likelihood
 CVA:

Cerebrovascular Accident
 DC:

Discrete Choice
 HrQoL:

Health related Quality of Life
 IBD:

Inflammatory Bowel Disease
 IRT:

Item Response Theory
 LLTM:

Linear Logistic Test Model
 MAPR:

MultiAttribute Preference Response
 MML:

Marginal Maximum Likelihood.
 OPLM:

OneParameter Logistic Model.
 RA:

Rheumatoid Arthritis.
 WHO:

World Health Organization.
References
 1.
WHO LUWHO. Preamble to the Constitution of the World Health Organization as adopted by the International Health Conference, New York, 1922 June, 1946. Geneva, Switzerland. http://whqlibdoc.who.int/hist/official_records/constitution.pdf.
 2.
Krabbe PFM. The Measurement of health and Health Status: Concepts, Methods and Applications from a Multidisciplinary Perspective. San Diego: Elsevier/Academic Press; 2016.
 3.
Nord E, Enge AU, Gundersen V. QALYs: is the value of treatment proportional to the size of the health gain? Health Econ. 2010;19:596–607.
 4.
Salomon JA. Techniques for valuing health states. Encyclopedia of Health Economics. 2014;2:454–8.
 5.
Attema AE, EdelaarPeters Y, Versteegh MM, Stolk EA. Time tradeoff: one methodology, different methods. Eur J Health Econ. 2013;14:53–64.
 6.
Brazier JE, Dixon S, Ratcliffe J. The role of patient preferences in costeffectiveness analysis: a conflict of values. PharmacoEconomics. 2009;27 https://0doiorg.brum.beds.ac.uk/10.2165/1131484000000000000000 LB  Brazier2009.
 7.
Krabbe PFM, Tromp N, Ruers TJM, van Riel PLCM. Are patients’ judgments of health status really different from the general population. Heal Qual Life Outcomes. 2011;9:31.
 8.
Krabbe PFM. A generalized measurement model to quantify health: the multiattribute preference response model. PLoS One. 2013;8
 9.
Rasch G. Probabilistic models for some intelligence and attainment tests expanded edition with foreword and afterword by B.D. Wright. Chicago: University of Chicago Press; 1980.
 10.
Fischer GH. Einführung in die Theorie Psychologischer Tests. In: Bern: Verlag Hans Huber; 1974.
 11.
Andrich D. An elaboration of Guttman scaling with Rasch models for measurement. Sociol Methodol. 1985;15:33–80.
 12.
Dolan P. Modeling valuation for EuroQol health states. Med Care. 1997;35:1095–108.
 13.
Sintonen H. The 15D instrument of healthrelated quality of life: properties and applications. Ann Intern Med. 2001;33:328–36.
 14.
Richardson J, Sinha K, Iezzi A, khan MANV63 MA. Modelling the utility of health states with the assessment of quality of life (AQoL) 8D instrument: overview and utility scoring algorithm. Centre for Health Econ. 2011;
 15.
Anderson JP, Kaplan RM, Berry CC, Bush JW, Rumbaut RG. Interday reliability of function assessment for a health status measure. The Quality of Wellbeing Scale Med Care. 1989;27:1076–84.
 16.
Feeny D, Furlong W, Torrance GW, Goldsmith CH, Zhu Z, DePauw S, et al. Multiattribute and singleattribute utility functions for the health utilities index mark 3 system. Med Care. 2002;40:113–28.
 17.
Scheiblechner H. Das Lernen und Lösen komplexer Denkaufgaben. Zeitschrift für Exp und Angew Psychol. 1972;19:476–506.
 18.
Fischer GH. Linear logistic test model as an instrument in educational research. Acta Psychol. 1973;37:359–74.
 19.
Fischer GH. Logistic latent trait models with linear constraints. Psychometrika. 1983;48:3–26.
 20.
de Boeck P, Wilson M. Explanatory item response models: a generalized linear and nonlinear approach. New York: SpringerVerlag.
 21.
Scheerens J, Glas C, Thomas SM. Educational evaluation, assessment and monitoring. A systematic approach. In: New York: Taylor & Francis; 2003.
 22.
Fox JP. Bayesian item response modeling: theory and applications. New York: Springer; 2010.
 23.
Parkin D, Devlin NJ, Foo Y. What determines the shape of an EQ5D index distribution. Med Decis Mak. 2016;36:941–51.
 24.
Andersen EB. Sufficient statistics and latent trait models. Psychometrika Psychometrica. 1977;24:69–81.
 25.
Fischer GH. On the existence and uniqueness of maximumlikelihood estimates in the Rasch model. Psychometrica. 1981;46:59–77.
 26.
Fischer GH. Unidimensional linear logistic Rasch models. In: Hambleton RK, Van der WJ LI, editors. Handbook of modern item response theory. New York: SpringerVerlag; 1997. p. 225–43.
 27.
Fischer GH. Linear logistic test models. In: Encyclopedia of social measurement; 2005. p. 505–14.
 28.
Andersen EB. A goodness of fit test for the Rasch model. Psychometrica. 1973;38:123–40.
 29.
Fleiss JL. Statistical methods for rates and proportions. 2nd ed. New York: John Wiley; 1981.
 30.
Mair P, Hatzinger R. Mair MJ. eRm: extended Rasch modeling [computer software]. R package version. 0:15–4. http://cran.rproject.org/package=eRm
 31.
Streiner DL, Norman GR, Cairney J. Health measurement scales: a practical guide to their development and use. 5th ed. Oxford: Oxford University Press; 2015.
 32.
Weinstein M, Torrance G, McGuire A. QALYs: the basics. Value Heal 2009;12 Supplement:55–59.
 33.
Fanshel S, Bush J. A healthstatus index and its application to healthservices outcomes. Oper Res. 1970;18:1021–66.
 34.
Scalone L, Stalmeier PFM, Milani S, Krabbe PFM. Values for health states with different life durations. Eur J Health Econ. 2015;16:917–25.
 35.
van Hoorn RA, Donders ART, Oppe M, Stalmeier PM. The better than dead method: feasibility and interpretation of a valuation study. PharmacoEconomics. 2014;32:789–99.
 36.
Panayides P, Robinson C, Tymms P. Rasch measurement: a response to Goldstein. Br Educ Res J. 2015;41:180–2.
 37.
Neisser UBG, Boodoo GBTJ, Bouchard TJ, Boykin AW, Brody N, Ceci SJ, et al. Intelligence: knowns and unknowns. Am Psychol. 1996;51:77–101.
 38.
Albert JH. Bayesian estimation of normal ogive item response curves using Gibbs sampling. J Educ Stat. 1992;17:251–69.
 39.
Chrzan K. Using partial profile choice experiments to handle large numbers of attributes. Int J Mark Res. 2010;52:827–40.
 40.
Bansback N, Brazier JE, Tsuchiya A, Anis A. Using a discrete choice experiment to estimate health state utility values. J Health Econ. 2012;31:306–18.
 41.
Krabbe PFM, Devlin NJ, Stolk EA, Shah KK, Oppe M, van Hout B, et al. Multinational evidence on the feasibility and consistency of a discrete choice model in quantifying EQ5D5L health states. Med Care. 2014;52:935–43.
 42.
Barofsky I. Quality: its definition and measurement as applied to the medically ill. New York: Springer; 2012.
Acknowledgments
We want to thank Jiske Verburg, Martine de Zoete, Kees Brandsema, Katja Morsink, Colin Koel, Gerda Drent, and Noor Tromp for their dedicated effort in collecting data from the patients.
Availability of data and materials
The anonymized dataset is available from the corresponding author.
Author information
Affiliations
Contributions
CGO participated in the design of the study, performed the statistical analysis, interpreted the results of the analysis and drafted the manuscript; EH advised on the mathematics and the statistical analysis, interpreted the results of the analysis, drafted the manuscript; PK participated in the design of the study, interpreted the results of the analysis and drafted the manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Catharina G. M. GroothuisOudshoorn.
Ethics declarations
Ethics approval and consent to participate
The data collection under cancer and RA patients was approved by the Central Committee on Research Involving Human Subjects (region ArnhemNijmegen), the Medical Ethics Review Committee Twente (METC/14124), and the Medical Ethics Review Committee UMCG (METC 2015/496). They issued a waiver for this study, indicating that the pertinent Dutch Legislation (the Medical Research Involving Human Subjects Act) did not apply for this noninterventional study. Formal informed consent was therefore not mandatory. After sending an information letter about the study and the research task, consent to participate was obtained from all patients by a research nurse. Consent to participate was given verbally.
Consent for publication
Not applicable.
Competing interests
All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare that they have no financial relationships with any organizations that might have an interest in the submitted work in the previous 3 years; no competing interest or other relationships or activities that could appear to have influenced the submitted work. As an extension of the IRTbased measurement model presented in this paper, additional measurement models, tools and instruments are developed by PFMK as part of academic/commercial activities.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1:
Classification of the patients in the three studies of their own health status on the EQ5D3 L. (DOCX 21 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
GroothuisOudshoorn, C.G.M., van den Heuvel, E.R. & Krabbe, P.F.M. A preferencebased item response theory model to measure health: concept and mathematics of the multiattribute preference response model. BMC Med Res Methodol 18, 62 (2018) doi:10.1186/s1287401805168
Received
Accepted
Published
DOI
Keywords
 Healthrelated quality of life
 Health status
 Latent logistic test model
 Patientreported measurement
 Rasch model