 Research
 Open access
 Published:
How is disease severity associated with quality of life in psoriasis patients? Evidence from a longitudinal populationbased study in Sweden
Health and Quality of Life Outcomes volume 15, Article number: 151 (2017)
Abstract
Background
Assessing the impact of disease severity on generic quality of life (QOL) is a critical step in outcomes research and in the development of decisionanalytic models structured around health states defined by clinical measures. While data from routine clinical practice found in healthcare registers are increasingly used for research, more attention should be paid to understanding the relationship between clinical measures of disease severity and QOL. The purpose of this work was therefore to investigate this relationship in psoriasis using a populationbased dataset.
Methods
Severity was measured by the Psoriasis Area and Severity Index (PASI), which combines severity of erythema, induration, and desquamation into a single value ranging from 0 to 72. The generic EQ5D3L utility instrument, under the UK tariff, was used to measure QOL. The association between PASI and EQ5D3L was estimated using a populationbased dataset of 2674 patients with moderate to severe psoriasis enrolled over ten years in the Swedish psoriasis register (PsoReg). Given the repeated measurement of patients in the register data, a longitudinal fixedeffects model was employed to control for unobserved patientlevel heterogeneity.
Results
Marginal changes in PASI are associated with a nonlinear response in EQ5D3L: Moving from PASI 10 to 9 (1 to 0) is associated with an increase of 0.0135 (0.0174) in EQ5D3L. Furthermore, unobserved patientlevel heterogeneity appears to be an important source of confounding when estimating the relationship between QOL and PASI.
Conclusions
Using register data to estimate the impact of disease severity on QOL while controlling for unobserved patientlevel heterogeneity shows that PASI appears to have a larger impact on QOL than previously estimated. Routine collection of generic QOL data in registers should be encouraged to enable similar applications in other disease areas.
Trial registration
Not applicable.
Background
Economic evaluations of health technologies compare the efficiency of a new health intervention with one or more alternatives in costeffectiveness analyses with qualityadjusted life years (QALYs) as outcomes. In order to be relevant for decision making clinical efficacy, generic quality of life (QOL), and cost data are often synthesized in a decisionanalytic modelling framework [1]. Decisionanalytic models used to evaluate costeffectiveness are often based on clinicallydefined health states, each of which has a cost and QOL profile [2]. Although establishing the relationship between measures of clinicallydefined disease severity and QOL measures is an essential part of an economic evaluation, realworld datasets such as registers are underutilised and the methodologies employed often admit avoidable sources of bias.
Register data have increasingly been used to develop inputs for decisionanalytic models, although deriving healthcare resource use and associated costs has likely been the most common research objective [3]. In Sweden, many populationbased registers are available that include both disease severity and QOL measures, while reflecting actual clinical practice. A recent review shows that out of 103 Swedish healthcare registers, all include measures of disease severity, 46 measure QOL with EQ5D, and 14 measure QOL with SF36 [4]. Using unique personal identification numbers, these registers can be linked to other databases containing demographic, medical, and socioeconomic data.
Registers also provide follow up over time, making it possible to exploit the advantages of longitudinal statistical models. A primary advantage of a longitudinal approach is that time trends can be identified, allowing for observation of otherwise hidden relationships.
The objective of this study was to show how register data can be used to provide longitudinal estimates of the association between QOL and disease severity. This was illustrated using a case study in psoriasis, where similar estimates do not currently appear to be available.
Methods
Case study in psoriasis
Psoriasis is a chronic autoimmune disease affecting approximately 2–4% of the population in Western countries [5], and can be a debilitating disease that has negative consequences for patient QOL, due to physical and psychosocial factors. In psoriasis, disease severity is most often measured using the Psoriasis Area and Severity Index [6] (PASI), a multidimension questionnaire completed by the caregiver. This tool divides the body into four areas (head, arms, trunk, legs) where each area has an associated importance weight. The areas are scored individually for severity, measured by levels of erythema, induration, and desquamation, resulting in a score ranging from 0 (minimal severity) to a theoretical 72 (maximal severity). Economic evaluations commonly measure QOL using EQ5D, a multiattribute utility instrument with five dimensions [7]. The version with three levels of response was used in this study (EQ5D3L). Respondents fill in whether they have no, some, or extreme problems in each of the five dimensions (mobility, selfcare, usual activities, pain/discomfort, and anxiety/depression). Patients’ responses were mapped to a single index value using the UK population value set [8], where the values are bounded between −0.59 and 1. A sensitivity analysis was conducted using the Swedish experience based value set, ranging from 0.34 to 0.97 [9].
Previous research on the association between PASI and EQ5D3L is scarce. Existing results are limited in terms of their ability to control for confounding variables [10,11,12], and have not controlled for unobserved heterogeneity such as genetic factors. Since there are unobserved variables that are correlated with both QOL and disease severity, estimates of the impact of disease severity could be biased by omitting these unmeasured factors.
Overall analytic approach and data
Longitudinal fixed effects modelling is one technique to reduce bias arising from certain types of unobserved variables by using the patient as her own control. This method is recommended for certain types of data structures and research questions such as the ones in this study [13] as it exploits the longitudinal nature of the register data, as described in detail in the statistical methods section.
PsoReg is Sweden’s national psoriasis register, including patients with moderate to severe psoriasis who are eligible for systemic treatment from a specialist in dermatology [14]. Patients are registered on a continuous basis, where the current data cut includes patients registered from 2006 to 2014, including 4449 unique patients. The database contains diseaserelated data including severity and QOL, ideally recorded at each examination, and has been estimated to contain 65% of all psoriasis patients in Sweden receiving biologic treatment.
From the initial patient population of 4449 patients (20,048 observations), observations with missing data or input errors were excluded (Fig. 1). Patients with only one observation were also excluded, since the statistical methodology described below requires two or more observations per patient. In total, 2674 patients and 14,044 observations were included in the analysis with an average of 5.3 followup visits (minimum of 2, maximum of 36). The characteristics of patients and their disease in the analysis population (measured at patients’ first observation in PsoReg) and in the full PsoReg population are similar (Table 1).
Statistical methods
Quality of life is a complex phenomenon that is influenced by various factors, many of which are not observable. However, the fact that they are unobservable does not diminish their effect, and in order to reduce bias in the estimation of the impact of PASI on QOL those unobservable factors should be controlled for. As PsoReg data is longitudinal, analyses can address confounding due to timeinvariant, patientlevel, unobserved characteristics. Specifically, there are aspects of each patient, such as genetics, personality, and sex that do not tend to change over time. The first two are typically unobserved. While traditional models can control for observed timeinvariant characteristics such as sex, more sophisticated (longitudinal) models, such as fixed or random effects models, can also account for unobserved timeinvariant effects.
The estimation of the treatment effect is more efficient in a randomeffects model than a fixedeffect model, but the former makes the assumption that the unobserved individual level effects are uncorrelated with the observed variables. Fixed and random effects models are both unbiased when this assumption holds, but only fixed effects are unbiased when it does not. In the context of this analysis, independence between the observed and unobserved effects is unreasonable. For example, since genetics (unobserved fixed effect) predispose certain people to psoriatic arthritis (observed effect), and both are correlated with QOL, the independence assumption does not hold and a fixedeffects model should be estimated. Although fixed effects methods used in this analysis have been well described in statistical literature, different definitions of the term “fixed effects” are sometimes used. For clarity, the theoretical model is written in Eqs. (1), (2), (3).
The index i is the patient and t is the time index. y is QOL, X is the matrix of observed variables, δ is the vector of unobserved individual fixed effects, and e is the random error component. The model in Eq. (3) is actually estimated, and provides equivalent estimates of β as in Eq. (1) according to the FrischWaughLovell theorem [15, 16]. Eq. (1) cannot be estimated because the fixed effects (δ_i) are not observed. The double dot accent (Eq. 3) represents the value of each variable after subtracting the average over all observations for each patient (Eq. 2).
In the analysis, the dependent variable is QOL, measured by EQ5D3L. The quadratic term of each variable was tested for influence in the estimated model by minimising the Akaike information criteria (AIC). The model with the lowest AIC (−12,930) included quadratic PASI and BMI terms and therefore these were included in the model. Interaction terms between PASI and age, PASI and psoriatic arthritis, and age and psoriatic arthritis were also included. The fixed effects estimator used a constraint that the sum of the individual fixed effects equal 0, implying that the intercept is calculated so that the predicted values of EQ5D3L equal the average value of EQ5D3L at the mean of each variable.
To confirm the necessity of using fixed effects instead of random effects as described above, a Hausman test [17] was performed which rejected the null hypothesis that both estimators are consistent, at an alpha level of 0.01. As such, the final model used fixed effects. Robust standard errors were also used following a Wald test for groupwise homoscedasticity [18], which was rejected at an alpha level of 0.01.
Stata v11.2 ( www.stata.com ) was used to analyse data and R version 3.1.1 ( www.rproject.org ) for graphics.
Results
The increase in EQ5D3L due to the decrease in PASI depends on the value of PASI and its square as well as the interactions with age and PSA. The linear and quadratic effects of PASI are both statistically significant (p < 0.01). The interaction of PASI with age is also significant (p < 0.01) while the interaction with PSA is significant at a 10% level (p = 0.08) (Table 2).
Based on the regression results, EQ5D3L was predicted as a function of PASI holding all other variables constant at their respective mean (Fig. 2 and Table 3). The xaxis in Fig. 3 extends to where PASI is equal to 30; the 99th percentile of the PASI values in PsoReg is 25.1.
The model explained 61.2% of variance in EQ5D3L according to the R^{2} measure of goodness of fit. Within the sample, the individual fixed effects can be calculated. In this case, the model’s (Eq. 1) root mean square error (RMSE) of the actual EQ5D3L value compared to the model’s prediction was 0.152 and the mean absolute error (MAE) was 0.104 (Fig. 3).
As discussed in the methods, the fixed effects regression cannot estimate the effects of timeinvariant effect such as sex. However, the effect of PASI on EQ5D3L can be calculated within sex stratifications (Table 4), where it appears that the direct linear effect of PASI is stronger in females than males.
The results described above use the UK EQ5D3L tariff. A sensitivity analysis was conducted using the Swedish EQ5D3L tariff, which has a smaller possible range of values. Using the same analysis as above where the only change was to replace the UKbased EQ5D3L index value with the Swedish EQ5D3L index value, the coefficient of the linear PASI estimate was −0.0086 (95% CI 0.0104, −0.0069; p < 0.001). The estimate of the PASI squared estimate was positive and statistically significant, but very small. The resulting R^{2} using the Swedish tariff was 0.665.
Discussion
Developing QOL inputs for economic evaluations, based on clinical definitions of disease severity, is a core component of an economic assessment. In practice, this is often hindered by a lack of available data and methodology that allows researchers to address selection bias, obtain relevant covariates, and observe patients for a sufficiently long duration. This study illustrates how register data can be used to derive the association between disease severity and QOL by exploiting the unique content and structure of registers. Registers contain many relevant clinical and QOL variables, and may be linked to other databases containing sociodemographic information, allowing for adjustment of many sources of confounding. Even when an important confounding variable is not available, the longitudinal data structure can be exploited to control for factors that do not change within patients over time, such as genetics. Adjusting for all relevant factors, observed or unobserved, is one of the most important steps in reducing bias in a statistical analysis, which can be accomplished by employing the methods and type of data proposed in this study to estimate the association between QOL and disease severity. In turn, this may also reduce bias and uncertainty in economic evaluations where these estimates are used as parameter inputs.
In the case study presented here, the impact of longitudinal modelling methods in reducing bias can be observed by comparing predicted values where fixed effects are included or excluded: At PASI 10, all else equal, the predicted EQ5D3L controlling for fixed effects is 0.673, while excluding fixed effects is 0.708 (a difference of 0.035). In an economic evaluation, this relatively small difference in statespecific QOL may result in large differences when accumulated over time, particularly important when evaluating chronic diseases such as psoriasis, and thus have a large impact on the results of a costeffectiveness analysis.
In terms of performance, the model’s predictive ability compared well to previous research. A metaanalysis of 119 models mapping disease specific measures to generic preferencebased measures shows that previous studies have an R^{2} between 0.17–0.51, an MAE between 0.00–0.19, and an RMSE between 0.08–0.20 [19]. The predicted EQ5D3L values including fixed effects result in MAE and MRSE values that fall in the middle of this range. The R^{2} and adjusted R^{2} values of this study are above or at the top of this range, implying that the variation in the predictions from the estimated model is comparable to other mappingfocused research, and that the model explains more variation than most previous mappings.
In the sensitivity analysis using the Swedish EQ5D3L tariff, the finding that the PASI coefficient estimates using the Swedish tariff are smaller than under the UK tariff are unsurprising given the smaller range of the EQ5D3L index value using the Swedish tariff compared to the UK tariff. Interestingly, the linear PASI coefficient is approximately 2.6 times larger in the analysis using the UK tariff compared to the Swedish tariff, and the range of EQ5D3L values is also approximately 2.6 times larger in the UK tariff compared to the Swedish tariff. This may suggest that in some sense, the magnitude of the association between PASI and EQ5D3L is consistent in both analyses, where the difference may reflect the range of the dependent variable. The goodness of fit estimates suggest that the model may explain somewhat more variation in EQ5D3L under the Swedish tariff than the UK tariff.
In addition to illustrating how register data can generally be used to associate disease severity and QOL, this psoriasis case study also provides the first results of this association that control for unobserved patient characteristics in psoriasis. The results show that there is a statistically significant, negative, nonlinear relationship between disease severity and QOL in psoriasis. Previous estimates of this relationship using PsoReg data provided estimates of the impact of PASI on QOL using EQ5D3L [11]. These results were not fully adjusted and used a less complex functional form, but it is the only study that estimated the association in a comparable way to our study, and the estimates show that each unit reduction in PASI results in a reduction of EQ5D3L by 0.0089 [11]. In the present study, changes in PASI result in nonlinear changes in EQ5D3L. However, the size of the effect is larger than previously estimated for PASI values below 21: Moving from PASI 10 to 9 results in a marginal EQ5D3L change of 0.0135 and a move from PASI 1 to 0 results in an increase in EQ5D3L of 0.0174. At the upper end of the PASI scale, smaller changes in EQ5D3L measured QOL are predicted: moving from PASI 25 to 24 results in an increase in QOL of 0.0071, highlighting the nonlinearity of the relationship. For low PASI scores, the impact of a marginal change in PASI on QOL is therefore greater than previously estimated.
The results can also be put into the context of typical clinical endpoints in psoriasis, which are often based on percentage PASI improvements (i.e. PASI 90). The change in absolute QOL can be calculated based on a patient’s baseline value and percentage improvement. For example, if it is known that a patient achieved 90% PASI improvement and started with PASI equal to 10, we can infer that patient improved to at least PASI equal to 1, a total of 9 units. Using the results of the regression, the value of EQ5D3L at PASI = 10 and PASI = 1 can be calculated (Table 3), where the difference is equal to 0.1372. This value represents the estimated difference in EQ5D3L between a PASI = 9 and PASI = 1 health state after controlling for each observed variable as well as timeinvariant unobserved fixed effects.
An important aspect of decisionanalytic modelling it to account for patient heterogeneity, as costeffectiveness may differ in patients with different characteristics. This implies that input parameters, including QOL, should ideally be estimated for different subgroups of patients. The data and analytic approach employed in this study can be used to achieve this by disaggregating QOL based on applicable patient characteristics. For example, using the information derived in the study, it is possible to calculate the expected EQ5D3L specifically for a 60 yearold patient with PSA, a BMI of 30, and a PASI score of 10 (including the interactions). For a patient with these characteristics, the predicted EQ5D3L is 0.656 compared to the average aggregated patient, where predicted EQ5D3L is 0.673. As illustrated above, a benefit of the model estimated in this study is that it is a classic linear model, as it is linear in parameters while allowing for nonlinear relationships between EQ5D3L and each regressor. This means that the coefficients estimated for each variable have the simple interpretation that a one unit increase in the variable value corresponds to a change in EQ5D3L equal to the size and direction of the coefficient.
Since fixedeffects models cannot estimate timeinvariant parameters such as sex, stratification can be used to estimate the results in different patient groups. In the psoriasis case study, the results show that the effect of reducing PASI by one unit may have a larger impact on QOL for female patients than males (Table 4). Put in a clinical context, overall health gains in psoriasis patients may depend on baseline PASI score [20], treatment effect and the corresponding absolute change from baseline PASI [21, 22], and patient characteristics such as sex.
Limitations and future research
Registers are a rich source of information, and we have illustrated how they can be utilized if they include both disease severity and QOL. However, they can vary in both content and quality, and it is therefore important to evaluate the usefulness and representativeness of the data on a casebycase basis with respect to the percentage of eligible patients, potential for bias, and which measures of severity and QOL are available.
Provided that reliable data are available in registers, we believe that future research should employ the methods outlined in this work in order to estimate the association between QOL and clinical measures of disease severity in disease areas other than psoriasis. Differences in predicted QOL between various data sources, such as registers and trial data, could also be explored further from both a qualitative and quantitative perspective.
Conclusions
Register data collected from routine clinical practice can be a valuable source of information to estimate the relationship between clinical measures of disease severity and QOL, a central component in economic evaluations. Given the data structure of registers, longitudinal modelling techniques provide a means to adjust for unobserved patientlevel heterogeneity. In a case study combining register data with longitudinal modelling techniques, a nonlinear relationship between clinically defined disease severity and QOL was found.
Abbreviations
 AIC:

Akaike information criteria
 EQ5D3L:

3level EuroQoL 5dimension instrument
 PASI:

Psoriasis Area and Severity Index
 QALY:

Qualityadjusted life years
 QOL:

(Generic) quality of life
References
Sculpher MJ, Claxton K, Drummond M, McCabe C. Whither trialbased economic evaluation for health care decision making? Health Econ. 2006;15:677–87.
Siebert U, Alagoz O, Bayoumi AM, Jahn B, Owens DK, Cohen DJ, et al. Statetransition modeling: a report of the ISPORSMDM modeling good research practices task force3. Value Health. 2012;15:812–20.
Asaria M, Grasic K, Walker S. Using linked electronic health records to estimate healthcare costs: key challenges and opportunities. PharmacoEconomics. 2016;34:155–60.
Emilsson L, Lindahl B, Köster M, Lambe M, Ludvigsson JF. Review of 103 Swedish healthcare quality registries. J Intern Med. 2015;277:94–136.
Parisi R, Symmons DP, Griffiths CE, Ashcroft DM. Global epidemiology of psoriasis: a systematic review of incidence and prevalence. J Invest Dermatol. 2013;133:377–85.
Langley RG, Ellis CN. Evaluating psoriasis with psoriasis area and severity index, psoriasis global assessment, and lattice system physician’s global assessment. J Am Acad Dermatol. 2004;51:563–9.
The EuroQol Group. EuroQola new facility for the measurement of healthrelated quality of life. Health policy. 1990;16:199–208.
Dolan P. Modeling valuations for EuroQol health states. Med Care. 1997;35:1095–108.
Burström K, Sun S, Gerdtham UG, Henriksson M, Johannesson M, Levin LÅ, et al. Swedish experiencebased value sets for EQ5D health states. Qual Life Res. 2014;23:431–42.
Herédi E, Rencz F, Balogh O, Gulácsi L, Herszényi K, Holló P, et al. Exploring the relationship between EQ5D, DLQI and PASI, and mapping EQ5D utilities: a crosssectional study in psoriasis from Hungary. European Journal of Health Economics. 2014;15:111–9.
Norlin JM, Steen Carlsson K, Persson U, SchmittEgenolf M. Analysis of three outcome measures in moderate to severe psoriasis: a registrybased study of 2450 patients. Br J Dermatol. 2012;166:797–802.
Weiss SC, Kimball AB, Liewehr DJ, Blauvelt A, Turner ML, Emanuel EJ. Quantifying the harmful effect of psoriasis on healthrelated quality of life. J Am Acad Dermatol. 2002;47:512–8.
Haro JM, Kontodimas S, Negrin MA, Ratcliffe M, Suarez D, Windmeijer F. Methodological aspects in the assessment of treatment effects in observational health outcomes studies. Appl Health Econ Health Policy. 2006;5:11–25.
SchmittEgenolf M. PsoRegthe Swedish registry for systemic psoriasis treatment. The registry’s design and objectives. Dermatology. 2007;214:112–7.
Frisch R, Waugh FV. Partial time regressions as compared with individual trends. Econometrica J Econ Soc. 1933;1(4):387–401.
Lovell MC. Seasonal adjustment of economic time series and multiple regression analysis. J Am Stat Assoc. 1963;58:993–1010.
Hausman JA. Specification tests in econometrics. Econometrica. 1978;36(6):1251–71.
Greene WH. Econometric Analysis. 6th Edition. Upper Saddle River: Pearson Prentice Hall; 2008.
Brazier JE, Yang Y, Tsuchiya A, Rowen DL. A review of studies mapping (or cross walking) nonpreference based measures of health to generic preferencebased measures. Eur J Health Econ. 2010;11:215–25.
Hägg D, Eriksson M, Sundström A, SchmittEgenolf M. The higher proportion of men with psoriasis treated with biologics may be explained by more severe disease in men. PLoS One. 2013;8:e63619.
Hyrich KL, Watson KD, Silman AJ, Symmons DP, BSR BR. Predictors of response to antiTNFalpha therapy among patients with rheumatoid arthritis: results from the British Society for Rheumatology biologics register. Rheumatology. 2006;45:1558–65.
Kvien TK, Uhlig T, Ødegård S, Heiberg MS. Epidemiological aspects of rheumatoid arthritis. Ann N Y Acad Sci. 2006;1069:212–22.
Acknowledgements
The authors would like to thank all patients and healthcare professionals for using and advancing PsoReg.
Funding
PsoReg is supported by the Swedish Board of Health and Welfare and the Swedish Association of Local Authorities and Regions. This specific research was not sponsored and was completed as part of a doctoral program. The authors independently designed the study, interpreted the data, and wrote the manuscript.
Availability of data and materials
The data that support the findings of this study were retrieved from the Swedish Psoriasis Registry (PsoReg). Restrictions apply to the availability of these patient data, which were used under license for the current study.
Author information
Authors and Affiliations
Contributions
All authors contributed to the preparation and review of the manuscript, and to the interpretation of the data and results. KG drafted the manuscript and conducted the statistical analyses. All authors read and approved the final version. MSE is the guarantor.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This research was done in adherence to the Declaration of Helsinki and approved by the Umeå Ethical Review Board.
Consent for publication
Not applicable.
Competing interests
KG was employed at PAREXEL International at the time the research was conducted and the manuscript was prepared. MSE is responsible for dermatology in the project management for the national guidelines for psoriasis at the Swedish Board of Health and Welfare. MH has no disclosures.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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
Geale, K., Henriksson, M. & SchmittEgenolf, M. How is disease severity associated with quality of life in psoriasis patients? Evidence from a longitudinal populationbased study in Sweden. Health Qual Life Outcomes 15, 151 (2017). https://doi.org/10.1186/s129550170721x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s129550170721x