Next Article in Journal
Loss of XBP1 Leads to Early-Onset Retinal Neurodegeneration in a Mouse Model of Type I Diabetes
Next Article in Special Issue
Best-Evidence Rehabilitation for Chronic Pain Part 2: Pain during and after Cancer Treatment
Previous Article in Journal
Shared Decision-Making in Chronic Patients with Polypharmacy: An Interventional Study for Assessing Medication Appropriateness
Previous Article in Special Issue
A Meta-Epidemiological Appraisal of the Effects of Interdisciplinary Multimodal Pain Therapy Dosing for Chronic Low Back Pain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Moderate and Stable Pain Reductions as a Result of Interdisciplinary Pain Rehabilitation—A Cohort Study from the Swedish Quality Registry for Pain Rehabilitation (SQRP)

1
Department of Neurosurgery and Pain Rehabilitation, Skåne University Hospital, SE-221 85 Lund, Sweden
2
Pain and Rehabilitation Centre, and Department of Medical and Health Sciences, Linköping University, SE-581 85 Linköping, Sweden
3
Department of Social and Welfare Studies, Linköping University, SE-602 21 Norrköping, Sweden
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2019, 8(6), 905; https://doi.org/10.3390/jcm8060905
Submission received: 10 June 2019 / Revised: 18 June 2019 / Accepted: 20 June 2019 / Published: 24 June 2019
(This article belongs to the Special Issue Rehabilitation for Persistent Pain Across the Lifespan)

Abstract

:
Few studies have investigated the real-life outcomes of interdisciplinary multimodal pain rehabilitation programs (IMMRP) for chronic pain. This study has four aims: investigate effect sizes (ES); analyse correlation patterns of outcome changes; define a multivariate outcome measure; and investigate whether the clinical self-reported presentation pre-IMMRP predicts the multivariate outcome. To this end, this study analysed chronic pain patients in specialist care included in the Swedish Quality Registry for Pain Rehabilitation for 22 outcomes (pain, psychological distress, participation, and health) on three occasions: pre-IMMRP, post-IMMRP, and 12-month follow-up. Moderate stable ES were demonstrated for pain intensity, interference in daily life, vitality, and health; most other outcomes showed small ES. Using a Multivariate Improvement Score (MIS), we identified three clusters. Cluster 1 had marked positive MIS and was associated with the overall worst situation pre-IMMRP. However, the pre-IMMRP situation could only predict 8% of the variation in MIS. Specialist care IMPRPs showed moderate ES for pain, interference, vitality, and health. Outcomes were best for patients with the worst clinical presentation pre-IMMRP. It was not possible to predict who would clinically benefit most from IMMRP.

1. Introduction

Pain is an unpleasant experience with complex interactions between sensorimotoric, affective, and cognitive brain networks. As such, pain, especially chronic pain, is influenced by and interacts with physical, psychological, social, and contextual factors [1,2,3]. One-fifth of the European population has moderate to severe chronic pain conditions [4]. These conditions are associated with psychological distress, low health, sick leave, and high socioeconomic costs [5]. Therefore, a biopsychosocial (BPS) framework should be considered in clinical practice [6,7,8].
Unlike single/unimodal interventions, interdisciplinary multimodal pain rehabilitation programs (IMMRPs) for chronic pain—an interdisciplinary treatment according to the International Association for the Study of Pain (IASP)—distinguish themselves as well-coordinated complex interventions. Typically, IMMRPs are based on cognitive behavioural therapy (CBT) models (including Acceptance Commitment Therapy, ACT) and are administered over several weeks to months [9,10,11,12]. The Swedish programs generally include group activities such as pain education, supervised physical activity, training in simulated environments, and CBT coordinated by an interdisciplinary team (e.g., physician, occupational therapist, physiotherapist, psychologist, and social worker) based on a BPS framework [9,10,11,12]. The components of IMMRP are most often chosen based on the available evidence for unimodal interventions for chronic pain, for example, with respect to education, exercise, psychological interventions, and interventions for return to work. The core goals of rehabilitation programs in general [13] and especially for patients with chronic pain [14] are broad and multifactorial in combination with the individualised goals of the patient. These include increased ability to participate in valued activities such as work. Hence, IMMRP is a complex intervention [13,15] and, unlike pharmacological intervention, focusses on the whole person rather than just biochemical processes, implying complex patient conditions matched with complex IMMRPs [16,17]. The components of IMMRP can be active independently or interdependently [15], resulting in a combination of effects explained by known and unknown mechanisms. The effects are assumed to be greater than the sum of its components [18].
Systematic reviews (SRs) have generally reported higher efficacy both on a general level and for specific outcomes of IMMRP compared with single-treatment or treatment-as-usual programs [10,12,19,20,21,22,23]. SRs and Randomised Controlled Trials (RCTs) may be associated with risk for bias resulting from, for example, an unrepresentative selection of patients and researcher allegiance [24,25,26]. Thus, it is necessary to investigate whether the evidence obtained from SRs and RCTs can be replicated within a consecutive non-selected flow of patients in practice settings using prospective observational cohort study designs such as practice-based evidence (PBE). PBE has also been applied in the field of rehabilitation research [27]. The importance of such an approach is also emphasised in the real-effectiveness medicine framework [28]. IMMRPs are time consuming and expensive, even when most of the activities are group-based. From an ethical, individual, and socioeconomic perspective, it is indeed remarkable to note the lack of studies investigating effect sizes (ES) in patient populations in real-life practice settings. A recent study from two Swedish university clinics reported effect sizes of 0.51–0.61 (i.e., moderate ES) for two pain intensity variables at 12-month follow-up [29]. These effect sizes should be confirmed in larger studies based not only on patients at university hospitals, but also on specialist units in general. It would be motivating for patients to endure increases in pain, which is often observed in clinical practise during the start-up period of rehabilitation characterised by an increase in activity levels, if it were known that the long-term effects include the reduction of pain levels.
Complex interventions such as IMMRP should have several outcomes measured at multiple levels and strategies for handling multiple outcomes [17,30]. IMMRPs are evaluated using many outcomes. For example, one SR including 46 RCTs reported nine outcomes per RCT (median) [10]. However, outcomes are not usually divided into primary and secondary outcomes [10]. In addition, although it is most likely that changes in several of the selected outcomes are correlated, most SRs of IMMRPs evaluate the outcomes as independent from each other. Patterns of potential correlations (i.e., multivariate correlation patterns) are mainly unknown/uninvestigated, even though they could give valuable information regarding how to optimise IMMRPs. Hence, there is a need to develop clinically applicable ways to evaluate the multiple outcomes of MMPRs both for individual patients and within research studies.
The above knowledge gaps motivated this PBE study of chronic pain patients based on patient reported outcome measures (PROMs) from the Swedish Quality Registry for Pain Rehabilitation (SQRP) [31]. This registry offers an opportunity to investigate clinical outcomes and patterns of change, since all the relevant specialist care units throughout Sweden deliver data to SQRP. Hence, this PBE study has the general aim of investigating the effects of IMMRP in specialist care in Sweden considering the multivariate complexity of outcomes. We hypothesised that IMMRP in special care is associated with small-to-medium ES, that changes in outcomes generally are intercorrelated, and that the baseline situation (pre-IMMRP) can predict the multivariate outcomes. More specifically, we defined the following four aims:
  • To investigate the outcome effect sizes of IMMRP immediately post-IMMRP and at 12-month follow-up.
  • To analyse the multivariate correlation patterns of changes in outcomes of IMMRP: pre-IMMRP versus post IMMRP and pre-IMMRP versus 12-month follow-up.
  • To define a multivariate outcome measure of IMMRP.
  • To investigate if the clinical self-reported presentation pre-IMMRP can predict the multivariate outcome measure.

2. Materials and Methods

2.1. The Swedish Quality Registry for Pain Rehabilitation (SQRP)

The SQRP, recognised by the Swedish Association of Local Authorities and Regions, receives data from all specialist care units in Sweden [31]. The SQRP is based on PROM questionnaires that capture biopsychosocial data such as the patient’s background, pain distribution and intensity, pain-related cognitions, and psychological distress symptoms (e.g., depression and anxiety), as well as activity/participation aspects and health-related quality of life variables. Patients complete the PROM questionnaires on up to three occasions: (1) during assessment at the first visit to the unit (pre-IMMRP); (2) immediately after the IMMRP (post-IMMRP); and (3) at the 12-month follow-up (FU) after IMMRP discharge (12-month FU).

2.2. Subjects

This study included SQRP data from women and men ≥18 years old with complex chronic (≥3 months) non-malignant pain who were referred to specialist pain and rehabilitation units (i.e., specialist care centres) between 2008–2016. These patients can be characterised as complex, as their health profiles included psychiatric comorbidities such as depression and anxiety, low levels of acceptance, high levels of kinesiophobia, decreased working life and participation in social activities, and/or did not respond to routine pharmacological/physiotherapeutic treatments delivered in a monodisciplinary fashion. Strict inclusion and exclusion criteria for inclusion in the registry is not available, since this is a registry study of patients with complex chronic pain conditions referred from mainly the primary care to specialist care in Sweden. A minority of patients were referred from other specialist clinics e.g., orthopedic and rheumatology clinics. The following general inclusion criteria for IMMRP were used: (i) disabling chronic pain (on sick leave or experiencing major interference in daily life due to chronic pain); (ii) age 18 years and above; (iii) no further medical investigations needed; and (iv) written consent to participate and attend IMMRP. General exclusion criteria for IMMRP included severe psychiatric morbidity, abuse of alcohol and/or drugs, diseases that did not allow physical exercise, and specific pain conditions with other treatment options available (i.e., red flags).
The proportions of patients within primary health care with chronic pain conditions are not exactly known, but 10–20% are estimates [32,33]. Furthermore, the proportion of chronic pain patients within primary health care that are referred to specialist clinics is not known.
The study was conducted in accordance with the Helsinki Declaration and Good Clinical Practice and approved by the Ethical Review Board in Linköping (Dnr: 2015/108-31). All the participants received written information about the study and gave their written consent.

2.3. Variables

Background variables that were collected pre-IMMRP and symptom-related self-reported variables that were collected at all three times (pre, post, and 12-month FU) were used in the analyses. The variables and instruments used are mandatory for the units registering their data with the SQRP.

Background Variables

The following background variables were collected: age (years), gender (man or woman), education level, and country of birth. Education level was dichotomised into university and the other alternatives (i.e., upper secondary school, elementary school, or other); this variable was labelled as University. Country of birth was dichotomised as from Europe and outside Europe and labelled as Outside-Europe. In addition, self-reported pain duration (days), persistent pain duration (days), and number of days off work (Days no work) were obtained.
Pain distribution was registered using 36 predefined anatomical areas (18 on the front and 18 on the back of the body) and the patients registered the areas with pain: (1) head/face, (2) neck, (3) shoulder, (4) upper arm, (5) elbow, (6) forearm, (7) hand, (8) anterior aspect of chest, (9) lateral aspect of chest, (10) belly, (11) sexual organs, (12) upper back, (13) low back, (14) hip/gluteal area, (15) thigh, (16) knee, (17) shank, and (18) foot. The number of areas with pain (range: 1–36) were summed, and the obtained variable was denoted as the Pain Region Index (PRI).

2.4. Repeated Self-Reported Measures

For reports of the psychometric aspects of the self-reported measures, the reader is referred to other studies summarising these [7,34,35,36].

2.4.1. Pain Aspects

Pain intensity average during the previous seven days was registered using a 0–10 (0 = no pain and 10 = worst possible pain) numeric rating scale (NRS)—NRS-7days.

2.4.2. The Multidimensional Pain Inventory (MPI)

MPI is a 61-item self-report questionnaire that measures the psychosocial, cognitive, and behavioural effects of chronic pain [37,38]. Part 1 consists of five scales: Pain severity—measuring several aspects of the pain experience (MPI-Pain-severity); Interference—pain-related interference in everyday life (MPI-Pain-interfer); Perceived Life Control (MPI-LifeCon); the level of affective distress (MPI-Distress); and Social Support—perceived support from a spouse or significant others (MPI-SocSupp). Part 2 assesses the perception of responses to displays of pain and suffering from significant others and consists of three scales: Punishing Responses (MPI-Punish); Solicitous Responses (MPI-Solict); and Distracting Responses (MPI-Distract). Part 3 measures to what extent the patients participate in various activities using four scales. These scales can be combined into a composite scale—the General Activity Index (MPI-GAI)—which was used in the present study [39].

2.4.3. Psychological Distress Variables

Symptoms of anxiety and depression were registered using the Hospital Anxiety and Depression Scale (HADS) [40,41]. This instrument comprises seven items in each of two subscales: depression (HADS-D) and anxiety (HADS-A) symptoms. Both subscale scores have a range of 0 to 21. A score of 7 or less in each subscale indicates a non-case, a score of 8–10 indicates a possible case, and a score of 11 or more indicates an almost definite case [40].

2.4.4. The Short Form Health Survey (SF36)

The Short Form Health Survey (SF36) attempts to represent multidimensional health concepts and measurements of the full range of health states, including levels of well-being and personal evaluations of health [42]. Scores are standardised into eight dimensions with a scale from 0 to 100 where higher scores indicate a better perception of health [42]: (1) physical functioning (sf36-pf), physical activity level including activities of daily living; (2) role limitations due to physical functioning (sf36-rp), to what extent physical health limits the performance of work and other regular activities; (3) bodily pain (sf36-bp), pain and related disability; (4) general health (sf36-gh), evaluation of health situation; (5) vitality (sf36-vt), how rested and energetic; (6) social functioning (sf36-sf), disturbances of social life due to physical or mental illness; (7) role limitations due to emotional problems (sf36-re), difficulties in performing work or other regular activities due to emotional problems; and (8) mental health (sf36-mh), anxiety and depressive symptoms. Based on the eight scales, a physical summary component and a mental (psychological) summary component can be calculated, but these two summary component variables were not used in the present study.

2.4.5. The European Quality of Life Instrument (EQ-5D)

The European Quality of Life Instrument (EQ-5D) captures a patient’s perceived state of health [43,44,45]. The first part of the instrument defines five dimensions: mobility, self-care, usual activities, pain/discomfort, and anxiety/depression. Each of these were measured at three levels. An EQ-5D-index is derived by applying a formula that essentially attaches values (weights) to each of the levels in each dimension. The collection of index values (weights) for all the possible EQ-5D states is called a value set. Most EQ-5D value sets have been obtained from a standardised valuation exercise where a representative sample of the general population in a country/region is asked to place a value on EQ-5D health states. The EQ5D also measures the self-estimation of today’s health according to a 100-point scale, which is a thermometer-like scale (EQ-VAS) with defined end points (high values indicate good health and low values indicate bad health).

2.4.6. Estimations of Changes in Pain and in Life Situation

The patients post-IMMRP and at the 12-month FU estimated the degree of positive change in pain (Change-pain) and in their ability to handle life situations in general (Change-life situation). The Change-pain item was rated on a five-point Likert scale from markedly increased pain (0) to markedly decreased pain (4). The Change-life situation item was rated on a five-point Likert scale from markedly worsened (0) to markedly improved (4).

2.5. Statistics

All the statistics were performed using the statistical package IBM SPSS Statistics (version 24.0) and SIMCA-P+ (version 15.0; Umetrics Inc., Umeå, Sweden). A probability of <0.001 (two-tailed) was accepted as the criteria for significance due to the large number of subjects.
The text and tables generally report the mean value ± one standard deviation (±1 SD) together with a median and range of continuous variables. Percentages (%) are reported for categorical variables. The detailed analyses also report 95% confidence intervals (95% CI). SQRP uses predetermined rules when handling single missing items of a scale or a subscale; details about this have been reported elsewhere [29]. To compare groups, we used Student´s t-test for unpaired observations, analysis of variance (ANOVA with post hoc test if significant difference), and Chi square test. Effect sizes (ES; Cohen’s d) for within-group analysis were computed using a calculator when appropriate (https://webpower.psychstat.org/models/means01/effectsize.php). Hedges’ g, which provides a measure of effect size weighted according to the relative size of each sample, was used for between ES using a calculator (https://www.socscistatistics.com/effectsize/default3.aspx). The absolute effect size was considered very large for values ≥ 1.3, large for values between 0.80–1.29, moderate for values between 0.50–0.79, small for values between 0.20–0.49, and insignificant for values < 0.20 [46].
Common methods such as logistic regression (LR) and multiple linear regression (MLR) can quantify the level of relations of individual factors but disregard interrelationships among different factors and thereby ignore system-wide aspects [47]. Moreover, such methods assume variable independence when interpreting results [48], and there are several risks when considering one variable at a time [49]. To account for our aims, the problems related to handling missing data (see below), and the risks associated with multicollinearity problems (see above), we used advanced multivariate analyses (MVDA).
Hence, using SIMCA-P+, we applied advanced Principal Component Analysis (PCA) for the multivariate correlation analyses of all investigated variables and Orthogonal Partial Least Square Regressions (OPLS) for the multivariate regressions. These techniques do not require normal distributions of the included variables [50]. Note that the PCA of SIMCA-P+ differs considerably from the simpler version implemented (e.g., the version used in SPSS).
PCA extracts and displays systematic variation in the data matrix. All the variables were log transformed before the statistical analyses if data were skewed. Using PCA, we analysed the multivariate correlation pattern for the changes in the 22 outcome variables for all the subjects. Note that changes in outcomes are calculated so that a positive value indicates an improvement. A cross-validation technique was used to identify nontrivial components (p). Variables loading on the same component p were correlated, and variables with high loadings but with opposing signs were negatively correlated. Variables with high absolute loadings were considered significant. The obtained components are per definition not correlated and are arranged in decreasing order with respect to explained variation. The loading plot reports the multivariate relationships between variables. A corresponding plot reporting the relationships between subjects (i.e., t-scores) can also be used (score plot), and each subject receives a score (t) for each of the significant components. The t-score was used to calculate a Multivariate Improvement Score (MIS). R2 describes the goodness of fit—the fraction of sum of squares of all the variables explained by a principal component [51]. Q2 describes the goodness of prediction—the fraction of the total variation of the variables that can be predicted using principal component cross-validation methods [51]. Outliers were identified using two methods: score plots in combination with Hotelling’s T2 and distance to model in the X-space. No extreme outliers were detected.
OPLS was used for the longitudinal multivariate regression analyses of the t-scores of the PCA mentioned above using pre-IMMRP data (i.e., baseline data) [51]. The variable influence on projection (VIP) indicates the relevance of each X-variable pooled over all dimensions and Y-variables—the group of variables that best explain Y. VIP ≥ 1.0 was considered significant if VIP had a 95% jack-knife uncertainty confidence interval non-equal to zero. p(corr) was used to note the direction of the relationship (positive or negative). This is the loading of each variable scaled as a correlation coefficient, and thus standardising the range from −1 to +1. [50] p(corr) is stable during iterative variable selection and comparable between models [50]. Thus, a variable/regressor was considered significant when VIP > 1.0. For each regression, we report R2, Q2, and the p-value of a cross-validated analysis of variance (CV-ANOVA). SIMCA-P+ uses the Non-linear Iterative Partial Least Squares (NIPALS) algorithm to handle missing data: maximum 50% missing data for variables/scales and maximum 50% missing data for subjects.
To identify clusters based on the t-scores of the PCA mentioned above, we performed hierarchical clustering analysis (HCA). Based on the identified clusters (subgroups) defined by HCA, we performed partial least squares discriminant analysis (PLS-DA). In addition, we applied a bottom–up HCA to the principal component score vectors using the default Ward linkage criterion to identify relevant subgroups of patients. HCA can find subtle clusters in the multivariate space. In the resulting dendrogram, clusters were identified and, based on these groups, we performed PLS-DA using group belonging as the Y-variable and the psychometric data as predictors (X-variables). The PLS-DA model was computed to identify associations between the X-variables and the subgroups. Based on the HCA defined clusters, traditional inferential statistics (ANOVA including post hoc tests when appropriate) were computed using SPSS.

3. Results

3.1. Background Data

There were 14,666 chronic pain patients registered in the SQRP that fulfilled the inclusion criteria: chronic pain; >18 years of age; and completed the SQRP questionnaire before and on at least one of the two time points after the IMMRP. More than half (60%) of the patients answering the questionnaires pre-IMMRP and post-IMMRP also answered the questionnaires at 12-m FU. Most of the patients (76.3%) were women, 25.2% had studied at university, and 10.4% were born outside of Europe. More men were born outside Europe (men: 13.4% versus women: 9.5%; Chi2 = 43.437, p < 0.001), and fewer men had university education (men: 18.0% versus women: 27.4%; Chi2 = 123.672; p < 0.001). Continuous background variables are shown in Table 1. Women were slightly younger than men (42.9 ± 10.7 versus 44.5 ± 10.7; p < 0.001) and reported more spreading of pain according to PRI (15.4 ± 8.8 versus 10.8 ± 7.0; p < 0.001). The other variables in Table 2 were not affected by gender.

3.2. Pairwise Comparisons of Repeated Measures

The results for pre-IMMRP and post-IMMRP are shown in Table 2. Significant improvements were generally found except for two of the three scales of the second part of the MPI. In addition, the comparisons between pre-IMMRP and the 12-month FU generally revealed significant improvements except for one of the scales on the second part of the MPI (Table 3). Some outcomes were associated with moderate effect sizes. For the pre-IMMRP versus post-IMMRP comparisons, three variables had moderate effects sizes: MPI-pain-severity, sf36-bp, and sf36-vt (Table 2). At the 12-month FU, MPI-pain-severity and sf36-bp were associated with moderate effect sizes; this was also the case for MPI-pain-interference and EQ5d-index (Table 3). However, generally small effect sizes were found for the significant improvements (Table 2 and Table 3). The variables of the second part of the MPI had insignificant effect sizes both post-IMMRP and 12-month FU.

3.3. Patients Not Participating in the 12-Month FU

There were only small differences between those reporting PROM data at the 12-month FU and those not reporting their situation pre-IMMRP (Supplementary Table S1). Although those not reporting had a somewhat worse situation for most of the PROM variables, the differences were of no clinical importance.

3.4. Estimations of Changes in Pain and in Life Situation

At both post-IMMRP and 12-month FU, most patients reported that their pain situation had improved as well as their ability to handle their life situation (Table 4).

3.5. Multivariate Correlation Pattern of Changes in Outcomes

PCAs of the changes (i.e., the difference) were performed for pre-IMMRP versus post-IMMRP and pre-IMMRP versus 12-month FU. Significant models were achieved for both analyses (Table 5). Similar patterns were obtained for the first significant component of the two PCAs (Table 5). The first component of both analyses, reflecting the most important variations, showed that changes in HAD-D, MPI-pain-severity, MPI-pain interference, MPI-control, MPI-distress, sf-36-bp, sf-36-vt, sf-36-sf, and sf36-mh were most important and intercorrelated significantly. Hence, it was obvious that the changes in outcome variables are intercorrelated. That is, rather than representing 22 independent variables, the multivariate analyses show that most changes in these variables are highly intercorrelated.
At 12-month FU, the PCA also identified two additional components (Table 5). The second component mainly reflected the intercorrelation pattern between the social support scale of the MPI and the scales of part 2 of the MPI. A third significant component only explaining 6% of the variation in the dataset was also obtained in the analysis of changes at the 12-month follow-up versus baseline (Table 5).
The loading plot (i.e., the intercorrelations between variables of the two most important components for the changes pre IMMRP versus 12-month FU) is shown in Figure 1 (Figure 1a is a graphic presentation of the first two components reported in Table 5). Figure 1b shows the corresponding score plot (i.e., the relationships between subjects/patients). Each patient can be described with a score (t-score) for each significant component. Patients with high positive t-scores on the first component show prominent changes in the important variables constituting the first component, whereas patients near zero do not benefit, and patients with negative t-scores (located to the left in the score plot) deteriorate in the multivariate context. Hence, the t-score of the first component of both analyses can be considered as a Multivariate Improvement Score, in the following denoted MIS-post-IMMRP and MIS-12-month FU.

3.6. Identification of Subgroups Based on the Multivariate Improvement Scores (MIS)

An HCA based on MIS-post-IMMRP was performed. Three subgroups/clusters were identified. Based on the HCA, a PLS-DA model with two predictive components was obtained with group belonging as the Y-variable (R2 = 0.35; Q2 = 0.35; CV-ANOVA p < 0.001; n = 14,666). Using a similar approach, we performed an HCA based on MIS-12-month FU. This analysis identified three subgroups/clusters. Based on the HCA, a PLS-DA model with two predictive components was obtained with group belonging as the Y-variable (R2 = 0.37; Q2 = 0.37; CV-ANOVA p < 0.001; n = 8851).
The MIS (i.e., t-score) showed clear positive values (i.e., improvements) for cluster 1 and negative scores (i.e., deterioration) for cluster 3 (Table 6 and Table 7). Cluster 2 was an intermediary cluster with overall slightly positive improvements. Prominent effect sizes in the pairwise comparisons were observed post-IMMRP: cluster 1 versus cluster 2 = 3.33; cluster 1 versus cluster 3 = 5.36; and cluster 2 versus cluster 3 = 2.77; 12-month FU: cluster 1 versus cluster 2 = 2.92; cluster 1 versus cluster 3 = 4.99; and cluster 2 versus cluster 3 = 2.34. Thus, distinct differences in improvement levels were detected between the three clusters.
To facilitate the understanding of the identified clusters, the clusters were compared for the variables in each PCA (Table 6 and Table 7). The three clusters differed significantly for all changes according to the ANOVAs performed. The post hoc tests showed that 20 of the 22 changes post-IMMRP differed significantly between all three clusters. The corresponding figure at 12-month FU was 21 of 22 changes.
The estimations of changes in pain (Change-pain) and in management of life (Change-life situation) were not included in the PCAs and thus not included in the calculations of MIS. However, these estimations showed a similar pattern: the most prominent positive changes were in cluster 1, and the least positive changes were in cluster 3 (Table 6 and Table 7).
In the next step, the identified three clusters of both analyses were compared for their pre-IMMRP values (Table 8 and Table 9). For the clusters obtained post-IMMRP (Table 8), small irrelevant cluster differences existed for age. The proportion with university education was significantly highest in cluster 1 and lowest in cluster 3, although the differences were small. Generally, cluster 1 was associated with the worst situation for the most variables followed longitudinally except for social support, two of the scales of Section 2 of the MPI, and sf36-pf. In contrast, cluster 3 had the best situation, and cluster 2 was intermediate (Table 8). A very similar pattern was found when using the clusters obtained from the 12-month FU (Table 9).

3.7. Longitudinal Regression of MIS Using Baseline Data

The outcome data at baseline (pre-IMMRP) together with the background variables were used to regress MIS-post-IMMRP and MIS-12-month FU (Table 10). For both MIS, psychological distress variables were the most important regressors, but life impact variables, pain aspects, and health and vitality aspects contributed significantly. The directions of the correlations revealed that a more severe clinical situation (e.g., psychological distress, lack of control, low vitality and health, pain interference, and high pain intensity) were associated with high MIS (i.e., multivariate improvements). Although the obtained regressions were highly significant according to the CV-ANOVA, the explained variations in MIS were less than 10% (R2 = 0.08 in both analyses). Hence, most of the variations in the two MIS were not possible to predict.
Similar analyses for each of the clusters (Supplementary Tables S2 and S3) revealed that regressions were highly significant, but only explained a minority of the variations in MIS. Although the relative importance of the variables pre-IMMRP differed somewhat between the three clusters, no marked differences existed; that is, psychological distress aspects were the most important post-IMMRP (Supplementary Table S2). For the 12-month FU, somewhat more pronounced differences existed between the clusters: in cluster 2, the pain intensity aspects were the most important for MIS, and in cluster 1 and cluster 3, psychological distress variables together with pain interference were the most important for MIS.

4. Discussions

The major findings of the present large PROM study from SQRP are listed below:
  • Moderate long-term ES were found for pain intensity (MPI Pain severity and SF-36 bodily pain), interference in daily life (MPI Interference), and state of health (EQ-5D-index); most other variables showed small ES. Vitality also displayed moderate effect sizes immediately after IMMRP but fell slightly under cut-off for moderate change at 12-month follow-up. The majority of the 22 investigated outcomes were significantly improved.
  • Significant intercorrelations between changes in pain intensity, interference, control, psychological distress, and mental health were confirmed. The changes in 22 outcomes reflected one (pre-IMMRP versus post-IMMRP) or three (pre-IMMRP versus 12-month follow-up) latent components (groups of variables).
  • The outcomes were best for patients with the worst self-reported clinical presentation pre-IMMRP. Based on a defined multivariate improvement score (MIS), three clusters were identified. Cluster 1—overall, the worst situation pre-IMMRP—showed positive multivariate improvements in outcomes. Cluster 3—deteriorated—showed negative scores. Cluster 2, the intermediate cluster, was associated with overall slightly positive multivariate improvements.
  • Certain variables (especially psychological distress and life impact variables, pain, and health and vitality aspects) pre-IMMRP were associated with improvements according to MIS both post-IMMRP and at 12-month FU. However statistically significant, the pre-IMMRP situation could only explain a small part of the variation in MIS (8%); therefore, for clinical use, it was not possible to predict those who would benefit most from IMMRP.
The outcome variables mandatory in SQRP and presented in the present study are in good agreement with the BSP model of chronic pain and the outcome domains presented by the Initiative on Methods, Measurement, and Pain Assessment in Clinical Trials (IMMPACT) [7,52] and the Validation and Application of a patient-relevant core set of outcome domains to assess multimodal PAIN therapy (VAPAIN) [14] initiatives.
The present study was not primarily designed to evaluate the efficacy of IMMRP, which requires RCTs and SRs/meta-analyses. However, our results for the repeated measurements (Table 2 and Table 3) of chronic pain patients in real-life practice settings are in agreement with the positive evidence for IMMRPs reported in SRs [10,11,12] and in other studies [22,23,53]. As such, the small to moderate ES are noteworthy as these patients, who receive pain rehabilitation in specialist care centres, often have tried other treatments for their chronic pain with no or little effect. That is, these patients have severe problems and relative treatment resistance. Interestingly, the changes in outcomes with moderate ES are broad and not limited to a single outcome domain, and the most stable moderate ES were demonstrated for pain intensity aspects with moderate improvement both immediately after IMMRPs and at 12-month follow-up. Pain interference demonstrated moderate ES improvement at 12-month follow-up, and vitality was moderately improved immediately after IMMRPs. Both objective registrations (e.g., sick-leave registrations and actigraphic recordings [54]) and subjective PROM data may be important for understanding the efficacy of IMMRPs. Very recently, a SQRP study using a subgroup of the same cohort of patients reported that sick leave benefits according to the Swedish Social Insurance Agency decreased as a consequence of IMMRP [55]. Hence, both PROM data and objective sick leave data indicate clinically important positive changes in response to IMMRPs for patients in real-life practice settings. As a comparison, SRs conclude that common pharmacological treatments—e.g., paracetamol, non-steroidal anti-inflammatory drugs, and opioids—for patients with chronic pain have no effects, small effects, and/or lack of long-term follow-up effects [56,57,58].
The present study reported medium ES for two of three pain intensity variables both post-IMMRP and at 12-month FU (i.e., for MPI-pain-severity and sf36-bp); the third pain intensity variable had effect sizes near medium ES. These results contrasted some SRs reporting of no evidence for efficacy with respect to pain intensity [10,11]. However, not all RCTs of IMMRP included pain intensity outcomes, since the interventions are not focused on the pain itself but rather on its consequences [10,11]. Obviously, many pain patients consider pain intensity improvement to be the most important aspect of treatments [59]. However, changing this perspective is considered important in IMMPRs, since focusing on pain reduction in many cases leads to short-sighted attempts to control pain, and this may, when not successful, lead to increased physical and psychological disability and reduced life quality [60,61]. Thus, specialist care IMMRPs in Sweden have largely adopted the idea of introducing acceptance as a cornerstone of the psychological component of IMMRP (i.e., the willingness to have the experiences of pain as it is and to encourage patients to set up activity-related rehabilitation goals and risk initial pain flare-ups). This also means that patients are advised against establishing pain reduction goals. Thus, it could be considered problematic to communicate the present results showing medium effect sizes in real-world practice settings on pain. On the other hand, it may also be ethically problematic if both clinical practice and research ignore the reports and wishes of the patients regarding pain intensity. However, health care providers should not underestimate their patients’ ability to grasp, once explained, the complex pain experience. Therefore, health care providers should emphasise pain education, including descriptions of the affective and cognitive elements of pain as rational for the different components of IMMRPs, and stress the need to experiment with new behaviours and risk short-term pain flare-ups. Since the results are obtained in this context, no change in clinical practise as far as pain communication is called for.
SRs of IMMRP report that it is an effective intervention with small to moderate effects for patients with chronic pain conditions [11,12,62,63]. The present results concerning ES agree with most SRs of IMMRP, but it may also be appropriate to compare with ES results reported in other clinical studies. The moderate ES for two of the pain intensity variables agree with studies in clinical routine care (n = 65–395), and therefore, for long-term follow-up (6–12 months), such studies report small (Cohen’s d: 0.20–0.33 [64,65]) to moderate (Cohen’s d: 0.59–0.70 [26,66,67]) ES for pain intensity. For psychological distress variables, these studies agree with the present results: they generally found small ES for long-term follow-up (i.e., Cohen’s d = 0–0.38 for depressive symptoms [26,64,65] and Cohen’s d = 0.22–0.34 for anxiety) [26,65]. In a recent RCT comparing transdiagnostic emotion-focused exposure treatment (Hybrid) and Internet-delivered pain management treatment (ICBT) for chronic pain patients with comorbid anxiety and depression, we found that within group ES pre versus follow-up for pain interference were reported both for hybrid (ES = 1.17) and for ICBT (ES = 0.65) compared the present effect size of 0.49 [68]. However, the patients were not exclusively recruited from specialist care (i.e., clinical departments of pain rehabilitation); they were also recruited via advertisements in local newspapers and social media [68], and the numbers investigated were considerably smaller. An important observation from the present study is that moderate ES found at 12-month follow-up covered broad aspects (e.g., pain intensity, interference in daily life, and perceived health).
The number of outcomes in IMMRPs in RCTs are generally high, which reasonably reflects the broad goals of the complex intervention. The present study used 22 outcomes that are mandatory in SQRP measured on up to three occasions (i.e., pre-IMMRP, post-IMMRP, and 12-month follow-up). PCA was applied to handle the pattern of changes in potentially intercorrelated outcomes as suggested by the Medical Research Council of the United Kingdom [69]. From these analyses, it can be concluded that changes in pain intensity, pain interference, psychological distress, vitality, etc. were positively intercorrelated (Table 5). In fact, our study showed that the changes in the majority of the 22 outcomes are significantly intercorrelated. Hence, the changes in these variables cannot be considered independent of each other. As a consequence of this observation, the appropriateness to evaluate changes in outcomes separately, as done in a recent SR [70], must be questioned, since the treatment was not designed to target only a single outcome. Moreover, the ES must be seen in this complex context. Thus, small changes in many outcomes may be more important than one prominent change in a single or few outcomes. Furthermore, the Grading of Recommendations Assessment, Development and Evaluation (GRADE) used for evidence ratings in SRs may not adequately describe the evidence base of complex interventions [71]. Different definitions of positive outcomes of IMMRP interventions exist (e.g., the majority of outcomes had to be significantly better than for the control intervention) [10,11]. Another approach was that the authors of the SR predetermined primary and secondary outcomes and what was necessary to classify an intervention as positive before reviewing the RCTs [12].
The presented PCAs also highlight that it may be possible to reduce the number of outcome variables, since several of these appear to measure similar aspects of the chronic pain condition. The fact that 22 outcomes were analysed (Table 2 and Table 3) may raise an issue of multiple comparisons. In such situations, Bonferroni corrections are frequently used [72,73]. This is a conservative approach when the number of tests increases [72,74,75], the chances to detect real treatment effects decrease, and corrections were designed for corrections of independent comparisons [74]. The latter is obviously not present for most changes in outcome variables according to the PCAs performed (Table 5). Hierarchal or ‘gatekeeping’ procedures do not require adjustment for multiplicity [73], but require a natural hierarchy of the outcomes, as such a hierarchy is not obvious for IMMRP, as discussed above. Another approach is that outcomes are combined into a single composite outcome (i.e., a composite outcome consists of two or more component outcomes) [76], but this may be problematic with respect to missing cases and when the components of the composite endpoint are measured on different scales (i.e., non-commensurate outcomes) [76]. However, some multivariate methods such as PCA and OPLS can handle non-commensurate outcomes [76]. We used advanced PCA, including the NIPALS algorithm, to handling missing data and non-commensurate outcomes. We calculated the t-scores for the most relevant latent factor (component). Hence, we defined an objective Multivariate Improvement Score (MIS; the t-score of the first PCA component), which on an individual patient level defines the multivariate improvement; a positive MIS indicates multivariate improvements because of IMMRP.
Three clearly separated clusters based on MIS were identified. On a group level, clusters 1 and 2 were associated with various degrees of improvements, whereas cluster 3 showed negative MIS, indicating deterioration. Although the greater improvement in cluster 1 can be interpreted as a sign of regression to the mean and that these patients did not benefit from IMMRP more than cluster 2, this cluster still improves from IMMRP at least as well as those with e.g., less severe psychological distress symptoms (clusters 2 and 3). This may seem unexpected, but it is important to recognise that addressing psychological symptoms with CBT is an important component of IMMRPs. The patients at post-IMMRP and 12-month follow-up estimated the degree of positive change in pain (i.e., Change-pain) and the ability to handle life situation in general (i.e., Change-life situation). Most patients reported improvements according to both the Change-pain and Change-life situations (Table 4). Relatively small proportions of the patients reported worse situations post-IMMRP and at the 12-month follow-up, which are results that agree with other studies [29,77,78]. These two variables have retrospective elements even though they are not explicitly expressed. There are several problems with such items in general—e.g., desirability and memory aspects, recall time [79,80,81], and in treatment context (e.g., overly optimistic assessments) [82]. However, on a general level, these estimations and the two MIS variables (Table 6 and Table 7) agreed.
We found that cluster 1, which had high MIS values (i.e., prominent improvements), had a more severe clinical picture at baseline/pre-IMMRP than those with lower MIS (i.e., less improvements). These results agree with another SQRP study (N > 35,000) that identified clusters based on the clinical presentation at assessment (decision not taken about participation in IMMRP); the study found that patients with the most severe clinical situation who later participated in IMMRP had the most prominent improvements in six investigated outcomes [34]. Although IMMRP has been commended for its effectiveness (‘of all approaches to the treatment of chronic pain, none has a stronger evidence basis for efficacy, cost-effectiveness, and lack of iatrogenic complications’) [83], both this and our recent study [34] indicate that not all patients show important improvements in several domains of outcome after IMMRP. Both this and our previous study identified a large subgroup of patients that do not seem to significantly benefit from IMMRP. Presumably, these patients—in the present study, those with negative MIS (i.e., cluster 3)—need other interventions. In a relative context, they have a somewhat less complicated self-reported clinical picture pre-IMMRP than those in clusters 1 and 2, even though they are referred to specialist care and hence represent patients with complex needs.
The longitudinal regressions of MIS using background variables and pre-IMMRP data as regressors were significant (Table 10). A blend of variables was important; psychological distress variables were most important, but life impact variables, pain aspects, and health and vitality aspects contributed. Our results appear to be in line with a recent meta-analysis on prognostic factors for IMMRP outcome, demonstrating that both pre-treatment general emotional distress and pain-specific cognitive behavioural factors are related to worse long-term (>6 months) physical functioning [84]. Unfortunately, these regressions cannot be used clinically, since they only explained 8% of variations in MIS. Although the prediction does not work clinically, this and a previous study from our group give clear indications that patients with a severe clinical situation benefit from IMMRP [34].

4.1. Important Clinical Implications

Outcomes of IMMRP in real-life practice settings agree with the conclusions from SRs. Partly in contrast to SRs, this registry study of patients managed within specialist care found that pain intensity was positively affected because of IMMRP. It was also obvious that not all patients benefit from IMMRP. Hence, there is a need to develop better matching between clinical presentation and participation in MMRP in real-life practice settings. Moreover, the intercorrelations of most changes in outcomes also opens up the possibility of reducing the number of outcome variables and hereby reduce the burden upon patients included in the SQRP.

4.2. Strengths and Limitations

This study’s strengths include a large number of patients with complex chronic pain conditions with a nation-wide representation. However, these patients were referred to specialist clinics and thus represent a selection of the most difficult cases, so our results cannot be generalised to other settings. Another strength was the use of MVDA methods such as PCA and OPLS to handle correlation patterns, repeated measures, and regressions when there were obvious risks for multicollinearity. Changes in the social context may have changed and influenced the longitudinal analyses; however, we used validated and well-known instruments. Repeated evaluations using PROM questionnaires in treatment studies may be problematic [85]. Thus, the changes that the patient undergo because of the intervention (i.e., IMMRP) may affect the interpretations of the questions when presented at follow-up. The fact that no control group or treatment-as-usual group was available, which ethically is complicated to arrange for a registry of real-life practice patients, might have influenced our interpretation of changes after IMMRP. Data for the time period 2008–2016 from the SQRP was used in the present study, and changes in the content of IMMRP may have occurred. Unfortunately, no data concerning such changes are available.

5. Conclusions

This large-scale study of IMMRPs in real life practise settings demonstrates significant outcome changes in almost all measures. Most short-term and long-term effect sizes were small, but interestingly, moderate long-term effect sizes were demonstrated for pain, pain interference in daily life, and perceived health. In addition, patients reporting higher levels of perceived disability and suffering displayed greater improvement.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/8/6/905/s1, Table S1: Pre-IMMRP situation for patients reporting their outcomes at 12-m FU and those not reporting their outcomes at 12-m FU, Table S2: OPLS regressions of MIS at post IMMRP in the three clusters, Table S3: OPLS regressions of MIS at 12-month FU in the three clusters.

Author Contributions

Conceptualization, Å.R., E.D., M.B., B.L. and B.G.; Data curation, Å.R. and B.G.; Formal analysis, B.L. and B.G.; Methodology, Å.R., M.B. and B.G.; Validation, E.D.; Writing—original draft, B.G.; Writing—review & editing, Å.R., E.D., M.B. and B.L. All authors commented on different versions of the manuscript and all authors have approved the final version of the manuscript.

Funding

This study was supported by external funding from the Swedish Research Council (2018-02470), County Council of Östergötland (forsknings-ALF; LIO-608021), and AFA insurance (140340). AFA Insurance, a commercial founder, is owned by Sweden’s labour market parties: The Confederation of Swedish Enterprise, the Swedish Trade Union Confederation (LO), and The Council for Negotiation and Co-operation (PTK). They insure employees in the private sector, municipalities, and county councils. AFA Insurance does not seek to generate a profit, which implies that no dividends are paid to shareholders. The sponsors of the study had no role in study design, data collection, data analysis, data interpretation, writing of the report, or the decision to submit for publication. The authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Acknowledgment

The authors are very grateful for valuable comments from Marcelo Rivano Fischer, PhD.

Conflicts of Interest

The authors declare no conflict of interest.

List of Abbreviations

ACT Acceptance Commitment Therapy
ANOVA Analysis of Variance
BPSbiopsychosocial
CBT Cognitive Behavioral Therapy
Change-pain positive change in pain
Change-life situation change in ability to handle life situations in general
CI confidence interval
CV-ANOVA ANOVA of the cross-validated residuals
ES effect size
EQ-5D European Quality of Life instrument
EQ-5D-index index of EQ-5D based on five items
EQ-VAS health scale of EQ-5D
FU follow-up
GRADE Grading of Recommendations Assessment, Development and Evaluation
HADS Hospital Anxiety and Depression Scale
HADS-A Hospital Anxiety and Depression Scale—anxiety subscale
HADS-D Hospital Anxiety and Depression Scale—depression subscale
HCA Hierarchical Clustering Analysis
HR-QoL health-related quality of life
IASP the International Association for the Study of Pain
IMMPACT the Initiative on Methods, Measurement, and Pain Assessment in Clinical Trials
IMMRP Interdisciplinary Multimodal Pain Rehabilitation Program
LR logistic regression
MIS Multivariate Improvement Score
MLR multiple linear regression
MPI Multidimensional Pain Inventory
MPI-Pain-severity MPI subscale concerning pain severity
MPI-Pain-interfer MPI subscale concerning pain-related interference
MPI-Distress MPI-SocSupp affective distress
MPI-LifeCon MPI subscale concerning life control
MPI-SocSupp MPI subscale concerning social support
MPI-Punish MPI subscale concerning punishing responses
MPI-Solict MPI subscale concerning solicitous responses
MPI-Distract MPI subscale concerning distracting responses
MPI-GAI MPI subscale General Activity Index
MVDA advanced multivariate analysis
NIPALS Non-linear Iterative Partial Least Squares
NRS Numeric Rating Scale
NRS-7days average pain intensity the last week
OPLS Orthogonal Partial Least Square Regression
Outside-Europe born outside Europe
PBE practice-based evidence
P(corr) loading scaled as a correlation coefficient between −1.0 and +1.0
PCA Principal Component Analysis
PLS-DA partial least square discriminant analysis
PRI Pain Region Index
PROM patient reported outcome measures
RCT Randomised Controlled Trial
SQRP Swedish Quality Registry for Pain Rehabilitation
sf36 Short Form Health Survey
sf36-pf sf36 subscale concerning physical functioning
sf36-rp sf36 subscale concerning role limitations due to physical functioning
sf36-bp sf36 subscale concerning bodily pain
sf36-gh sf36 subscale concerning general health
sf36-vt sf36 subscale concerning vitality
sf36-sf sf36 subscale concerning social functioning
sf36-re sf36 subscale concerning role limitations due to emotional problems
sf36-mh sf36 subscale concerning mental health
SR systematic review
University University education
VAPAIN Validation and Application of a patient-relevant core set of outcome domains to assess multimodal PAIN therapy
VIP variable influence on projection

References

  1. Linton, S.J.; Bergbom, S. Understanding the link between depression and pain. Scand. J. Pain 2011, 2, 47–54. [Google Scholar] [CrossRef] [PubMed]
  2. Ossipov, M.H.; Dussor, G.O.; Porreca, F. Central modulation of pain. J. Clin. Investig. 2010, 120, 3779–3787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gatchel, R.J.; Peng, Y.B.; Peters, M.L.; Fuchs, P.N.; Turk, D.C. The biopsychosocial approach to chronic pain: Scientific advances and future directions. Psychol. Bull. 2007, 133, 581–624. [Google Scholar] [CrossRef]
  4. Breivik, H.; Collett, B.; Ventafridda, V.; Cohen, R.; Gallacher, D. Survey of chronic pain in Europe: Prevalence, impact on daily life, and treatment. Eur. J. Pain 2006, 10, 287. [Google Scholar] [CrossRef] [PubMed]
  5. Bergman, S. Chronic Musculoskeletal Pain: A Multifactorial Process; Lund University: Lund, Sweden, 2001. [Google Scholar]
  6. World Health Organization (WHO). International Classification of Functioning, Disability and Health (ICF); World Health Organization: Geneva, Switzerland, 2001. [Google Scholar]
  7. Dworkin, R.H.; Turk, D.C.; Farrar, J.T.; Haythornthwaite, J.A.; Jensen, M.P.; Katz, N.P.; Kerns, R.D.; Stucki, G.; Allen, R.R.; Bellamy, N.; et al. Core outcome measures for chronic pain clinical trials: IMMPACT recommendations. Pain 2005, 113, 9–19. [Google Scholar] [CrossRef] [PubMed]
  8. Fillingim, R.B. Individual Differences in Pain: Understanding the Mosaic that Makes Pain Personal. Pain 2017, 158, S11–S18. [Google Scholar] [CrossRef] [PubMed]
  9. Bennett, M.; Closs, S. Methodological issues in nonpharamacological trials for chronic pain. Pain Clin. Updates 2010, 18, 1–6. [Google Scholar]
  10. Swedish Council on Health Technology Assessment (SBU). Methods for Treatment of Chronic Pain a Systematic Review of the Literature; SBU-Rapport; Swedish Council on Health Technology Assessment: Stockholm, Sweden, 2006; Volume 177, (In Swedish: Metoder för behandling av långvarig smärta: En systematisk litteraturöversikt).
  11. Swedish Council on Health Technology Assessment (SBU). Rehabilitation of Chronic Pain; SBU-Rapport; Swedish Council on Health Technology Assessment: Stockholm, Sweden, 2010; Volume 198, (In Swedish: Rehabilitering vid långvarig smärta. En systematisk litteraturöversikt).
  12. Scascighini, L.; Toma, V.; Dober-Spielmann, S.; Sprott, H. Multidisciplinary treatment for chronic pain: A systematic review of interventions and outcomes. Rheumatology 2008, 47, 670–678. [Google Scholar] [CrossRef] [PubMed]
  13. Wade, D.T. Describing rehabilitation interventions. Clin. Rehabil. 2005, 19, 811–818. [Google Scholar] [CrossRef] [Green Version]
  14. Kaiser, U.; Kopkow, C.; Deckert, S.; Neustadt, K.; Jacobi, L.; Cameron, P.; De Angelis, V.; Apfelbacher, C.; Arnold, B.; Birch, J.; et al. Developing a core outcome-domain set to assessing effectiveness of interdisciplinary multimodal pain therapy: The VAPAIN consensus statement on core outcome-domains. Pain 2017, 159, 673–683. [Google Scholar] [CrossRef] [PubMed]
  15. Campbell, M.; Fitzpatrick, R.; Haines, A.; Kinmonth, A.L.; Sandercock, P.; Spiegelhalter, D.; Tyrer, P. Framework for design and evaluation of complex interventions to improve health. BMJ 2000, 321, 694–696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Paterson, C.; Baarts, C.; Launsø, L.; Verhoef, M.J. Evaluating complex health interventions: A critical analysis of the ’outcomes’ concept. BMC Complement. Altern. Med. 2009, 9, 18. [Google Scholar] [CrossRef] [PubMed]
  17. Shiell, A.; Hawe, P.; Gold, L. Complex interventions or complex systems? Implications for health economic evaluation. BMJ 2008, 336, 1281–1283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Hawe, P.; Shiell, A.; Riley, T. Complex interventions: How “out of control” can a randomised controlled trial be? BMJ 2004, 328, 1561–1563. [Google Scholar] [CrossRef] [PubMed]
  19. Weiner, S.S.; Nordin, M. Prevention and management of chronic back pain. Best Pract. Res. Clin. Rheumatol. 2010, 24, 267–279. [Google Scholar] [CrossRef] [PubMed]
  20. Kamper, S.J.; Apeldoorn, A.T.; Chiarotto, A.; Smeets, R.J.E.M.; Ostelo, R.W.J.G.; Guzman, J.; Van Tulder, M.W.; Van Tulder, M. Multidisciplinary biopsychosocial rehabilitation for chronic low back pain: Cochrane systematic review and meta-analysis. BMJ 2015, 350, h444. [Google Scholar] [CrossRef] [PubMed]
  21. Norlund, A.; Ropponen, A.; Alexanderson, K. Multidisciplinary interventions: Review of studies of return to work after rehabilitation for low back pain. J. Rehabil. Med. 2009, 41, 115–121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Busch, H.; Bodin, L.; Bergström, G.; Jensen, I.B. Patterns of sickness absence a decade after pain-related multidisciplinary rehabilitation. Pain 2011, 152, 1727–1733. [Google Scholar] [CrossRef] [PubMed]
  23. Jensen, I.B.; Busch, H.; Bodin, L.; Hagberg, J.; Nygren, A.; Bergström, G. Cost effectiveness of two rehabilitation programmes for neck and back pain patients: A seven year follow-up. Pain 2009, 142, 202–208. [Google Scholar] [CrossRef] [PubMed]
  24. Munder, T.; Brütsch, O.; Leonhart, R.; Gerger, H.; Barth, J. Researcher allegiance in psychotherapy outcome research: An overview of reviews. Clin. Psychol. Rev. 2013, 33, 501–511. [Google Scholar] [CrossRef]
  25. Margison, F.R.; Barkham, M.; Evans, C.; McGrath, G.; Clark, J.M.; Audin, K.; Connell, J. Measurement and psychotherapy. Evidence-based practice and practice-based evidence. Br. J. Psychiatry 2000, 177, 123–130. [Google Scholar] [CrossRef] [PubMed]
  26. Preis, M.A.; Vögtle, E.; Dreyer, N.; Seel, S.; Wagner, R.; Hanshans, K.; Reyersbach, R.; Pieh, C.; Mühlberger, A.; Probst, T. Long-Term Outcomes of a Multimodal Day-Clinic Treatment for Chronic Pain under the Conditions of Routine Care. Pain Res. Manag. 2018. [Google Scholar] [CrossRef] [PubMed]
  27. Whiteneck, G.G.; Gassaway, J. SCIRehab Uses Practice-Based Evidence Methodology to Associate Patient and Treatment Characteristics with Outcomes. Arch. Phys. Med. Rehabil. 2013, 94, S67–S74. [Google Scholar] [CrossRef] [PubMed]
  28. Malmivaara, A. Assessing the effectiveness of rehabilitation and optimizing effectiveness in routine clinical work. J. Rehabil. Med. 2018, 50, 849–851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Gerdle, B.; Molander, P.; Stenberg, G.; Stålnacke, B.-M.; Enthoven, P. Weak outcome predictors of multimodal rehabilitation at one-year follow-up in patients with chronic pain—A practice based evidence study from two SQRP centres. BMC Musculoskelet. Disord. 2016, 17, 287. [Google Scholar] [CrossRef]
  30. Craig, P.; Dieppe, P.; MacIntyre, S.; Michie, S.; Nazareth, I.; Petticrew, M. Medical Research Council Guidance. Developing and evaluating complex interventions: The new Medical Research Council guidance. BMJ 2008, 337, a1655. [Google Scholar] [CrossRef]
  31. Bromley Milton, M.; Borsbo, B.; Rovner, G.; Lundgren-Nilsson, A.; Stibrant-Sunnerhagen, K.; Gerdle, B. Is Pain Intensity Really That Important to Assess in Chronic Pain Patients? A Study Based on the Swedish Quality Registry for Pain Rehabilitation (SQRP). PLoS ONE 2013, 8, e65483. [Google Scholar] [CrossRef]
  32. Hasselström, J.; Liu-Palmgren, J.; Rasjö-Wrååk, G.; Liu-Palmgren, J.; Rasjö-Wrååk, G. Prevalence of pain in general practice. Eur. J. Pain 2002, 6, 375–385. [Google Scholar] [CrossRef]
  33. Mäntyselkä, P.; Kumpusalo, E.; Ahonen, R.; Kumpusalo, A.; Kauhanen, J.; Viinamäki, H.; Halonen, P.; Takala, J. Pain as a reason to visit the doctor: A study in Finnish primary health care. Pain 2001, 89, 175–180. [Google Scholar] [CrossRef]
  34. Gerdle, B.; Åkerblom, S.; Brodda Jansen, G.; Enthoven, P.; Ernberg, M.; Dong, H.-J.; Stålnacke, B.; Äng, B.; Boersma, K. Who benefit from multimodal rehabilitation—An exploration of pain, psychological distress and life impacts in over 35,000 chronic pain patients identified in the Swedish Quality Registry for Pain Rehabilitation (SQRP). J. Pain Res. 2019, 12, 891–908. [Google Scholar] [CrossRef]
  35. Rovner, G.S.; Sunnerhagen, K.S.; Björkdahl, A.; Gerdle, B.; Börsbo, B.; Johansson, F.; Gillanders, D. Chronic pain and sex-differences; women accept and move, while men feel blue. PLoS ONE 2017, 12, 0175737. [Google Scholar] [CrossRef] [PubMed]
  36. Bernfort, L.; Gerdle, B.; Husberg, M.; Levin, L.-A. People in states worse than dead according to the EQ-5D UK value set: Would they rather be dead? Qual. Life Res. 2018, 27, 1827–1833. [Google Scholar] [CrossRef] [PubMed]
  37. Turk, D.C.; Rudy, T.E. Toward an empirically derived taxonomy of chronic pain patients: Integration of psychological assessment data. J. Consult. Clin. Psychol. 1988, 56, 233–238. [Google Scholar] [CrossRef] [PubMed]
  38. Turk, D.C.; Rudy, T.E. Towards a comprehensive assessment of chronic pain patients. Behav. Res. Ther. 1987, 25, 237–249. [Google Scholar] [CrossRef]
  39. Bergström, G.; Jensen, I.B.; Bodin, L.; Linton, S.J.; Nygren, A.L.; Carlsson, S.G. Reliability and factor structure of the Multidimensional Pain Inventory—Swedish Language Version (MPI-S). Pain 1998, 75, 101–110. [Google Scholar] [CrossRef]
  40. Zigmond, A.S.; Snaith, R.P. The Hospital Anxiety and Depression Scale. Acta Psychiatr. Scand. 1983, 67, 361–370. [Google Scholar] [CrossRef] [Green Version]
  41. Bjelland, I.; Dahl, A.A.; Haug, T.T.; Neckelmann, D. The validity of the Hospital Anxiety and Depression Scale. J. Psychosom. Res. 2002, 52, 69–77. [Google Scholar] [CrossRef]
  42. Sullivan, M.; Karlsson, J.; Ware, J. The Swedish 36 Health survey. Evaluation of data quality, scaling assumption, reliability and construct validity across general populations in Sweden. Soc. Sci. Med. 1995, 41, 1349–1358. [Google Scholar] [CrossRef]
  43. EuroQol Group. EuroQol—A new facility for the measurement of health-related quality of life. Health Policy 1990, 16, 199–208. [Google Scholar] [CrossRef]
  44. Brooks, R. EuroQol: The current state of play. Health Policy 1996, 37, 53–72. [Google Scholar] [CrossRef]
  45. Dolan, P.; Sutton, M.; Sutton, M. Mapping visual analogue scale health state valuations onto standard gamble and time trade-off values. Soc. Sci. Med. 1997, 44, 1519–1530. [Google Scholar] [CrossRef]
  46. Bäckryd, E.; Persson, E.B.; Larsson, A.I.; Fischer, M.R.; Gerdle, B. Chronic pain patients can be classified into four groups: Clustering-based discriminant analysis of psychometric data from 4665 patients referred to a multidisciplinary pain centre (a SQRP study). PLoS ONE 2018, 13, e0192623. [Google Scholar] [CrossRef] [PubMed]
  47. Jansen, J.J.; Szymanska, E.; Hoefsloot, H.C.; Jacobs, D.M.; Strassburg, K.; Smilde, A.K. Between Metabolite Relationships: An essential aspect of metabolic change. Metabolomics 2012, 8, 422–432. [Google Scholar] [CrossRef] [PubMed]
  48. Pohjanen, E.; Thysell, E.; Jonsson, P.; Eklund, C.; Silfver, A.; Carlsson, I.-B.; Lundgren, K.; Moritz, T.; Svensson, M.B.; Antti, H. A Multivariate Screening Strategy for Investigating Metabolic Effects of Strenuous Physical Exercise in Human Serum. J. Proteome Res. 2007, 6, 2113–2120. [Google Scholar] [CrossRef] [PubMed]
  49. Eriksson, L.; Byrne, T.; Johansson, E.; Trygg, J.; Vikström, C. Multi—And Megavariate Data Analysis—Basic Principles and Applications, 3rd ed.; Umetrics Academy: Umeå, Sweden, 2013. [Google Scholar]
  50. Wheelock Åsa, M.; Wheelock, C.E. Trials and tribulations of ‘omics data analysis: Assessing quality of SIMCA-based multivariate models using examples from pulmonary medicine. Mol. BioSyst. 2013, 9, 2589. [Google Scholar] [CrossRef] [PubMed]
  51. Eriksson, L.; Johansson, E.; Kettaneh-Wold, N.; Trygg, J.; Wikström, C.; Wold, S. Multi—And Megavariate Data Analysis: Part I and II, 2nd ed.; Umetrics AB: Umeå, Sweden, 2006. [Google Scholar]
  52. Turk, D.C.; Dworkin, R.H.; Allen, R.R.; Bellamy, N.; Brandenburg, N.; Carr, D.B.; Cleeland, C.; Dionne, R.; Farrar, J.T.; Galer, B.S.; et al. Core outcome domains for chronic pain clinical trials: IMMPACT recommendations. Pain 2003, 106, 337–345. [Google Scholar] [CrossRef] [PubMed]
  53. Norrefalk, J.; Ekholm, K.; Linder, J.; Borg, K.; Ekholm, J. Evaluation of a multiprofessional rehabilitation programme for persistent musculoskeletal-related pain: Economic benefits of return to work. Acta Derm. Venereol. 2008, 40, 15–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Matthias, M.S.; Miech, E.J.; Myers, L.J.; Sargent, C.; Bair, M.J. There’s more to this pain than just pain: How patients’ understanding of pain evolved during a randomized controlled trial for chronic pain. J. Pain 2012, 13, 571–578. [Google Scholar] [CrossRef] [PubMed]
  55. Rivano-Fischer, M.; Persson, E.; Stålnacke, B.; Schult, M.; Löfgren, M. Return to work after interdisciplinary pain rehabilitation: One- and two-years follow-up based on the Swedish Quality Registry for pain rehabilitation. J. Rehabil. Med. 2019, 51, 281–289. [Google Scholar] [CrossRef] [PubMed]
  56. Busse, J.W.; Wang, L.; Kamaleldin, M.; Craigie, S.; Riva, J.J.; Montoya, L.; Mulla, S.M.; Lopes, L.C.; Vogel, N.; Chen, E.; et al. Opioids for Chronic Noncancer Pain: A Systematic Review and Meta-analysis. JAMA 2018, 320, 2448–2460. [Google Scholar] [CrossRef] [PubMed]
  57. Enthoven, W.T.; Roelofs, P.D.; Deyo, R.A.; Van Tulder, M.W.; Koes, B.W. Non-steroidal anti-inflammatory drugs for chronic low back pain. Cochrane Database Syst. Rev. 2016, 2, 012087. [Google Scholar] [CrossRef] [PubMed]
  58. Saragiotto, B.T.; Machado, G.C.; Ferreira, M.L.; Pinheiro, M.B.; Abdel Shaheed, C.; Maher, C.G. Paracetamol for low back pain. Cochrane Database Syst. Rev. 2016. [Google Scholar] [CrossRef] [PubMed]
  59. Henry, S.G.; Bell, R.A.; Fenton, J.J.; Kravitz, R.L. Goals of Chronic Pain Management: Do Patients and Primary Care Physicians Agree and Does It Matter? Clin. J. Pain 2017, 33, 955–961. [Google Scholar] [CrossRef] [PubMed]
  60. Thompson, M.; McCracken, L.M. Acceptance and Related Processes in Adjustment to Chronic Pain. Curr. Pain Headache Rep. 2011, 15, 144–151. [Google Scholar] [CrossRef] [PubMed]
  61. McCracken, L.M.; Zhao-O’Brien, J.; Zhao-O’Brien, J. General psychological acceptance and chronic pain: There is more to accept than the pain itself. Eur. J. Pain 2010, 14, 170–175. [Google Scholar] [CrossRef] [PubMed]
  62. Skelly, A.; Chou, R.; Dettori, J.; Turner, J.; Friedly, J.; Rundell, S.; Fu, R.; Brodt, E.; Wasson, N.; Winter, C.; et al. Noninvasive Nonpharmacological Treatment for Chronic Pain: A Systematic Review [Internet]; Agency for Healthcare Research and Quality (US): Rockville, MD, USA, 2018.
  63. Salathé, C.R.; Melloh, M.; Crawford, R.; Scherrer, S.; Boos, N.; Elfering, A. Treatment Efficacy, Clinical Utility, and Cost-Effectiveness of Multidisciplinary Biopsychosocial Rehabilitation Treatments for Persistent Low Back Pain: A Systematic Review. Glob. Spine J. 2018, 8, 872–886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Ruscheweyh, R.; Dany, K.; Marziniak, M.; Gralow, I. Basal Pain Sensitivity does not Predict the Outcome of Multidisciplinary Chronic Pain Treatment. Pain Med. 2015, 16, 1635–1642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Borys, C.; Lutz, J.; Strauss, B.; Altmann, U. Effectiveness of a Multimodal Therapy for Patients with Chronic Low Back Pain Regarding Pre-Admission Healthcare Utilization. PLoS ONE 2015, 10, 0143139. [Google Scholar] [CrossRef] [PubMed]
  66. Letzel, J.; Angst, F.; Weigl, M.B. Multidisciplinary biopsychosocial rehabilitation in chronic neck pain: A naturalistic prospective cohort study with intraindividual control of effects and 12-month follow-up. Eur. J. Phys. Rehabil. Med. 2018, in press. [Google Scholar] [CrossRef]
  67. Moradi, B.; Hagmann, S.; Zahlten-Hinguranage, A.; Caldeira, F.; Putz, C.; Rosshirt, N.; Schonit, E.; Mesrian, A.; Schiltenwolf, M.; Neubauer, E. Efficacy of multidisciplinary treatment for patients with chronic low back pain: A prospective clinical study in 395 patients. J. Clin. Rheumatol. 2012, 18, 76–82. [Google Scholar] [CrossRef]
  68. Boersma, K.; Södermark, M.; Hesser, H.; Flink, I.; Gerdle, B.; Linton, S. The efficacy of a transdiagnostic emotion-focused exposure treatment for chronic pain patients with comorbid anxiety and depression: A randomized controlled trial. Pain 2019, in press. [Google Scholar] [CrossRef] [PubMed]
  69. Craig, P.; Dieppe, P.; Macintyre, S.; Michie, S.; Nazareth, I.; Petticrew, M. Developing and Evaluating Complex Interventions: New Guidance. 2008. Available online: https://www.researchgate.net/publication/32899190_Developing_and_Evaluating_Complex_Interventions_New_Guidance_Online (accessed on 4 June 2019).
  70. Kamper, S.J.; Apeldoorn, A.T.; Chiarotto, A.; Smeets, R.J.; Ostelo, R.W.J.G.; Guzman, J.; Van Tulder, M.W. Multidisciplinary biopsychosocial rehabilitation for chronic low back pain. Cochrane Database Syst. Rev. 2014, 350. [Google Scholar] [CrossRef] [PubMed]
  71. Movsisyan, A.; Melendez-Torres, G.; Montgomery, P. Outcomes in systematic reviews of complex interventions never reached “high” GRADE ratings when compared with those of simple interventions. J. Clin. Epidemiol. 2016, 78, 22–33. [Google Scholar] [CrossRef] [PubMed]
  72. Feise, R.J. Do multiple outcome measures require p-value adjustment? BMC Med. Res. Methodol. 2002, 2, 8. [Google Scholar] [CrossRef]
  73. Turk, D.C.; Dworkin, R.H.; McDermott, M.P.; Bellamy, N.; Burke, L.B.; Chandler, J.M.; Cleeland, C.S.; Cowan, P.; Dimitrova, R.; Farrar, J.T.; et al. Analyzing multiple endpoints in clinical trials of pain treatments: IMMPACT recommendations. Pain 2008, 139, 485–493. [Google Scholar] [CrossRef] [PubMed]
  74. Bagiella, E. Clinical Trials in Rehabilitation: Single or Multiple Outcomes? Arch. Phys. Med. Rehabil. 2009, 90, S17–S21. [Google Scholar] [CrossRef] [PubMed]
  75. Tyler, K.M.; Normand, S.L.; Horton, N.J. The use and abuse of multiple outcomes in randomized controlled depression trials. Contemp. Clin. Trials 2011, 32, 299–304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Teixeira-Pinto, A.; Mauri, L. Msc Statistical Analysis of Noncommensurate Multiple Outcomes. Circ. Cardiovasc. Qual. Outcomes 2011, 4, 650–656. [Google Scholar] [CrossRef]
  77. Boonstra, A.M.; Reneman, M.F.; Waaksma, B.R.; Schiphorst Preuper, H.R.; Stewart, R.E. Predictors of multidisciplinary treatment outcome in patients with chronic musculoskeletal pain. Disabil. Rehabil. 2015, 37, 1242–1250. [Google Scholar] [CrossRef]
  78. Morley, S.; Williams, A.; Hussain, S. Estimating the clinical effectiveness of cognitive behavioural therapy in the clinic: Evaluation of a CBT informed pain management programme. Pain 2008, 137, 670–680. [Google Scholar] [CrossRef] [Green Version]
  79. Pina-Sánchez, J.; Koskinen, J.; Plewis, I. Measurement Error in Retrospective Reports of Unemployment; CCSR Working Paper; The Cathie Marsh Centre for Census and Survey Research, University of Manchester: Manchester, UK, 2012; pp. 1–56. [Google Scholar]
  80. Bernard, H.R.; Killworth, P.; Kronenfeld, D.; Sailer, L. The Problem of Informant Accuracy: The Validity of Retrospective Data. Annu. Rev. Anthropol. 1984, 13, 495–517. [Google Scholar] [CrossRef]
  81. Van Der Vaart, W.; Van Der Zouwen, J.; Dijkstra, W. Retrospective questions: Data quality, task difficulty, and the use of a checklist. Qual. Quant. 1995, 29, 299–315. [Google Scholar] [CrossRef]
  82. Schwartz, N. Retrospective and concurrent self-reports: The rationale for real-time data capture. In The Science of Real-Time Data CAPTURE: Self-Reports in Health Research; Oxford University Press: New York, NY, USA, 2007; pp. 11–26. [Google Scholar]
  83. Schatman, M. Interdisciplinary Chronic Pain Management: International Perspectives. Pain Clin. Updates 2012, 20, 1–6. [Google Scholar]
  84. Tseli, E.; Stalnacke, B.M.; Boersma, K.; Enthoven, P.; Gerdle, B.; Ang, B.O.; Grooten, W.J.A. Prognostic Factors for Physical Functioning After Multidisciplinary Rehabilitation in Patients with Chronic Musculoskeletal Pain: A Systematic Review and Meta-analysis. Clin. J. Pain 2018, 35, 148. [Google Scholar] [CrossRef] [PubMed]
  85. Westlander, G. Refined use of standardized self-reporting in intervention studies (In Swedish: Förfinad användning av standardiserad självrapportering i interventionstudier). Soc. Tidskr. 2004, 2, 168–181. [Google Scholar]
Figure 1. Loading plot of changes (pre-IMMRP vs. 12-month follow-up) in the 22 outcome variables—i.e., the relationships between the changes (a) and score plot ((b); the relationships between the patients). diff = change in a certain variable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; For explanations of the subscale abbreviations, see Methods.
Figure 1. Loading plot of changes (pre-IMMRP vs. 12-month follow-up) in the 22 outcome variables—i.e., the relationships between the changes (a) and score plot ((b); the relationships between the patients). diff = change in a certain variable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; For explanations of the subscale abbreviations, see Methods.
Jcm 08 00905 g001aJcm 08 00905 g001b
Table 1. Continuous background variables; mean ± SD and 95% confidence intervals (95% CI).
Table 1. Continuous background variables; mean ± SD and 95% confidence intervals (95% CI).
VariablesMean ± SD95% CI
Lower Bound
95% CI
Upper Bound
Age (years)43.2 ± 10.743.343.9
Days no work1055 ± 24619681095
Pain duration3057 ± 334129703170
Persistent pain duration 2368 ± 298022392414
PRI15.4 ± 8.615.115.6
Notes: SD = standard deviation; CI = confidence intervals; PRI = Pain Region Index.
Table 2. Outcome variables at baseline (pre-IMMRP) and immediately after IMMRP (post-IMMRP). Statistical comparisons are presented furthest to the right together with effects sizes (i.e., Cohen’s d). Effect sizes in bold were moderate, i.e., Cohen’s d ≥ 0.50. IMMRP: interdisciplinary multimodal pain rehabilitation programs.
Table 2. Outcome variables at baseline (pre-IMMRP) and immediately after IMMRP (post-IMMRP). Statistical comparisons are presented furthest to the right together with effects sizes (i.e., Cohen’s d). Effect sizes in bold were moderate, i.e., Cohen’s d ≥ 0.50. IMMRP: interdisciplinary multimodal pain rehabilitation programs.
Baseline vs. After IMMRP Pre-IMMRP Post-IMMRP
NMeanSDMeanSDp-ValueCohen’s d
NRS-7days14,1466.861.725.952.09<0.0010.45
HADS-A14,7749.004.767.784.55<0.0010.32
HADS-D14,7728.494.446.704.31<0.0010.47
MPI-Pain-severity14,6924.390.933.871.16<0.0010.52
MPI-Pain-interfer14,5524.381.023.941.19<0.0010.49
MPI-LifeCon14,6872.721.103.301.18<0.0010.47
MPI-Distress14,6973.461.262.891.38<0.0010.42
MPI-Socsupp14,6184.161.343.951.35<0.0010.21
MPI-punish13,0541.741.361.721.330.0370.02
MPI-protect12,9992.981.402.851.38<0.0010.12
MPI-distract13,0482.541.192.561.170.0430.02
MPI-GAI14,6762.440.842.630.82<0.0010.26
EQ-5D-index13,9890.260.310.390.33<0.0010.40
EQ-VAS13,77741.2219.0950.9921.38<0.0010.44
sf36-pf14,25352.7620.5857.6721.17<0.0010.30
sf36-rp13,94512.5324.4022.4633.12<0.0010.30
sf36-bp14,26824.3614.4932.9617.41<0.0010.52
sf36-gh13,98841.7020.2246.6921.88<0.0010.29
sf36-vt14,20623.9518.4835.6722.76<0.0010.54
sf36-sf14,22947.2925.1954.9325.91<0.0010.30
sf36-re13,70142.7742.9251.1543.48<0.0010.18
sf36-mh14,19455.0321.3562.5521.55<0.0010.38
NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey. For explanations of the subscale abbreviations, see Methods.
Table 3. Outcome variables at baseline (pre-IMMRP) and at the 12-month follow-up (FU). Statistical comparisons are presented furthest to the right together with effects sizes (i.e., Cohen’s d). Effect sizes in bold were moderate (i.e., Cohen’s d ≥ 0.50).
Table 3. Outcome variables at baseline (pre-IMMRP) and at the 12-month follow-up (FU). Statistical comparisons are presented furthest to the right together with effects sizes (i.e., Cohen’s d). Effect sizes in bold were moderate (i.e., Cohen’s d ≥ 0.50).
Baseline vs.
12-Month Follow-Up
Pre IMMRP 12-Month FU
NMeanSDMeanSDp-ValueCohen’s d
NRS-7days85686.841.725.782.32<0.0010.47
HADS-A88658.734.697.384.70<0.0010.33
HADS-D88658.184.376.744.66<0.0010.35
MPI-Pain-severity89044.360.913.711.33<0.0010.56
MPI-Pain-interfer88294.341.023.731.37<0.0010.54
MPI-LifeCon88712.771.103.281.27<0.0010.40
MPI-Distress88893.421.272.921.45<0.0010.35
MPI-Socsupp88304.171.333.771.42<0.0010.35
MPI-punish78241.691.341.691.350.6760.01
MPI-protect77842.961.392.781.40<0.0010.16
MPI-distract78112.521.172.451.17<0.0010.06
MPI-GAI88592.470.832.640.86<0.0010.20
EQ-5D-index88440.270.310.440.34<0.0010.50
EQ-VAS860741.9019.2952.9622.87<0.0010.46
sf36-pf845953.0720.3059.7322.57<0.0010.36
sf36-rp830113.0724.9127.7436.32<0.0010.39
sf36-bp845824.6014.1135.4120.05<0.0010.56
sf36-gh834242.5920.4947.3523.52<0.0010.25
sf36-vt844124.9618.7934.4123.85<0.0010.41
sf36-sf845948.9525.5057.6627.05<0.0010.32
sf36-re815944.6943.1755.6043.53<0.0010.22
sf36-mh843556.3421.1562.7022.53<0.0010.30
NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36); Health Survey; FU = Follow-up. For explanations of the subscale abbreviations see Methods.
Table 4. Estimations of pain situation (Change-pain) and in the ability to handle life situation in general (Change-life situation) made immediately after IMMRP (post-IMMRP) and at the 12-month FU.
Table 4. Estimations of pain situation (Change-pain) and in the ability to handle life situation in general (Change-life situation) made immediately after IMMRP (post-IMMRP) and at the 12-month FU.
Change-Pain Post-IMMRP 12-Month FU
n%n%
0. Markedly increased pain4473.22252.6
1. Partially increased pain1517115906.9
2. No change400829.1290534
3. Partially diminished pain617844.9366242.8
4. Markedly diminished pain160711.7117413.7
Total13 7571008 556100
Change-Life situation Post-IMMRP 12-Month FU
n%n%
0. Markedly deteriorated740.51081.3
1. Partially deteriorated2481.82823.3
2. No change192313.9161518.8
3. Partially improved841260.9462854
4. Markedly improved316122.9193722.6
Total13 8181008 570100
FU = Follow-up.
Table 5. Principal component analysis (PCA) of changes pre-IMMRP vs. post-IMMRP (left part) and pre-IMMRP vs. 12-month FU (right part). The significant components (p) are shown. Absolute loadings ≥ 0.25 are shown in bold to facilitate interpretation. Changes in outcomes are calculated so that a positive value indicates an improvement.
Table 5. Principal component analysis (PCA) of changes pre-IMMRP vs. post-IMMRP (left part) and pre-IMMRP vs. 12-month FU (right part). The significant components (p) are shown. Absolute loadings ≥ 0.25 are shown in bold to facilitate interpretation. Changes in outcomes are calculated so that a positive value indicates an improvement.
Changes Pre-IMMRP vs. Post-IMMRPChanges Pre-IMMRP vs. 12-Month FU
p[1]p[1]p[2]p[3]
diff-NRS-7days0.230.23−0.150.29
diff-HADS-A0.230.220.19−0.33
diff-HADS-D0.260.250.17−0.24
diff-MPI-Pain-sever0.270.27−0.160.26
diff-MPI-Pain-interfer0.260.28−0.110.10
diff-MPI-LifeCon0.260.250.09−0.05
diff-MPI-distress0.270.260.13−0.21
diff-MPI-SOCsupp−0.03−0.070.410.21
diff-MPI-punish0.070.080.320.11
diff-MPI-protect−0.02−0.020.510.33
diff-MPI-distract0.010.000.450.36
diff-MPI-GAI0.120.150.000.07
diff-EQ-5D-index0.220.22−0.070.12
diff-EQ-VAS0.220.23−0.030.09
diff-sf36-pf0.200.21−0.140.20
diff-sf36-rp0.190.21−0.100.17
diff-sf36-bp0.260.27−0.150.25
diff-sf36-gh0.210.210.020.02
diff-sf36-vt0.270.260.03−0.03
diff-sf36-sf0.250.240.05−0.09
diff-sf36-re0.180.170.13−0.26
diff-sf36-mh0.270.250.21−0.30
R20.310.360.100.06
Q20.250.310.040.02
N14,6668851
diff = change in a certain variable; p = principal component; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; FU = Follow-up. For explanations of the subscale abbreviations see Methods.
Table 6. Clusters from hierarchical clustering analysis (HCA) based on Multivariate Improvement Score (MIS) (t-scores) of the first component of the PCA for the changes in outcomes from pre-IMMRP to post-IMMRP (denoted MIS post-IMMRP). To facilitate understanding, the changes for all the outcomes are shown (mean, SD, and 95% confidence interval). The two bottom rows show the estimations of changes (not included in PCA and the calculation of MIS).
Table 6. Clusters from hierarchical clustering analysis (HCA) based on Multivariate Improvement Score (MIS) (t-scores) of the first component of the PCA for the changes in outcomes from pre-IMMRP to post-IMMRP (denoted MIS post-IMMRP). To facilitate understanding, the changes for all the outcomes are shown (mean, SD, and 95% confidence interval). The two bottom rows show the estimations of changes (not included in PCA and the calculation of MIS).
Cluster 1 (15.0%)Cluster 2 (54.1%)Cluster 3 (30.8%)
NMeanSD95% CI NMeanSD95% CI NMeanSD95% CI ANOVA
Variables LBUB LBUB LBUBp-ValuePost Hoc
MIS post-IMMRP22054.371.614.314.4479380.331.080.300.354523−2.741.16−2.78−2.71<0.001all different
diff-NRS-7days20863.042.052.953.1375111.011.690.971.044267−0.281.58−0.32−0.23<0.001all different
diff-HADS-A21915.123.724.965.2878731.593.081.521.664458−1.303.24−1.39−1.20<0.001all different
diff-HADS-D21885.933.555.786.0878692.242.952.172.304462−0.963.04−1.05−0.87<0.001all different
diff-MPI-Pain-severity21851.771.021.731.8278740.560.740.540.584500−0.170.69−0.19−0.15<0.001all different
diff-MPI-Pain-interfer21731.530.991.491.5778110.500.700.480.514453−0.170.66−0.18−0.15<0.001all different
diff-MPI-LideCon21891.941.051.891.9878610.710.960.690.734491−0.330.99−0.36−0.30<0.001all different
diff-MPI-distress21872.181.172.132.2378690.711.020.690.734487−0.461.05−0.49−0.43<0.001all different
diff-MPI-SOCsupp2180−0.341.12−0.39−0.307828−0.240.98−0.26−0.214462−0.110.98−0.14−0.08<0.001all different
diff-MPI-punish19910.391.200.340.4570160.041.120.020.074008−0.201.15−0.24−0.17<0.001all different
diff-MPI-protect1988−0.181.17−0.23−0.136984−0.141.00−0.17−0.123989−0.051.05−0.09−0.02<0.001cl1 NE cl2, cl2 NE cl3
diff-MPI-distract19920.111.100.060.1670120.000.96−0.030.0240060.011.03−0.020.04<0.001cl1 NE cl2, cl3
diff-MPI-GAI21870.570.820.530.6078660.230.680.210.244489−0.060.69−0.08−0.03<0.001all different
diff-EQ-5D-index21050.440.300.430.4574940.160.280.150.174205−0.070.28−0.08−0.06<0.001all different
diff-EQ-VAS206830.5119.6629.6631.36740911.6118.5811.1912.034126−3.8419.15−4.42−3.25<0.001all different
diff-sf36-pf214718.7416.9718.0219.4676446.0213.635.726.334324−3.8314.22−4.25−3.40<0.001all different
diff-sf36-rp211039.9339.2838.2541.60751310.4629.709.7911.134215−5.9625.56−6.74−5.19<0.001all different
diff-sf36-bp214628.0416.2727.3528.7376639.7212.969.4310.014323−2.9312.38−3.30−2.56<0.001all different
diff-sf36-gh212621.0617.5920.3121.8175296.0514.515.726.374235−4.8914.57−5.33−4.45<0.001all different
diff-sf36-vt213937.4518.7436.6538.24762613.5817.2613.1913.964311−4.2016.16−4.69−3.72<0.001all different
diff-sf36-sf214634.2622.1733.3235.20765210.1220.399.6610.584324−9.8621.10−10.49−9.23<0.001all different
diff-sf36-re208743.5244.5841.6045.43739011.9843.2111.0012.974121−15.8142.02−17.09−14.52<0.001all different
diff-sf36-mh213929.6117.4828.8630.35762010.1214.909.7810.454307−8.0115.89−8.49−7.54<0.001all different
Change-Pain20593.280.693.253.3172802.530.872.512.5540152.080.932.052.11<0.001all different
Change-Life situation20673.510.583.483.5373153.070.633.053.0840322.760.722.742.79<0.001all different
LB = Lower Bound; UB = Upper Bound; diff = change in a certain variable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey. For explanations of the subscale abbreviations, see Methods.
Table 7. Clusters from HCA based on MIS (t-scores) of the first component of the PCA for the changes in outcomes from pre-IMMRP to 12-month FU (denoted as MIS 12-m FU). To facilitate understanding, the changes for all the outcomes are shown (mean, SD, and 95% confidence interval). The two bottom rows show the estimations of changes (not included in PCA and the calculation of MIS).
Table 7. Clusters from HCA based on MIS (t-scores) of the first component of the PCA for the changes in outcomes from pre-IMMRP to 12-month FU (denoted as MIS 12-m FU). To facilitate understanding, the changes for all the outcomes are shown (mean, SD, and 95% confidence interval). The two bottom rows show the estimations of changes (not included in PCA and the calculation of MIS).
Cluster 1 (12.4%)Cluster 2 (46.6%)Cluster 3 (41.0%)
NMeanSD95% CI NMeanSD95% CI NMeanSD95% CI ANOVA
Variables LBUB LBUB LBUBp-ValuePost Hoc
MIS -12-m FU10995.011.784.905.1141230.781.350.740.823629−2.431.39−2.47−2.38<0.001all different
diff-NRS-7days10313.672.153.543.8038761.461.891.401.523435−0.161.66−0.21−0.10<0.001all different
diff-HADS-A10955.723.845.495.9540872.063.341.962.163588−0.803.53−0.91−0.68<0.001all different
diff-HADS-D10966.133.755.906.3540862.393.132.292.483588−1.053.39−1.16−0.94<0.001all different
diff-MPI-Pain-severity10922.301.142.242.3741080.880.900.850.913619−0.090.76−0.12−0.07<0.001all different
diff-MPI-Pain-interfer10922.301.142.242.3741080.880.900.850.913589−0.140.74−0.17−0.12<0.001all different
diff-MPI-LifeCon10862.271.102.212.3440760.830.850.810.863604−0.251.05−0.28−0.21<0.001all different
diff-MPI-distress10912.101.102.032.1640950.751.000.710.783618−0.381.13−0.41−0.34<0.001all different
diff-MPI-SOCsupp1086−0.691.25−0.77−0.624081−0.601.16−0.64−0.573583−0.101.04−0.13−0.07<0.001all different
diff-MPI-punish9790.571.190.500.653633−0.021.22−0.060.023194−0.131.22−0.18−0.09<0.001all different
diff-MPI-protect980−0.231.24−0.31−0.153616−0.351.20−0.39−0.3131730.031.09−0.010.06<0.001all different
diff-MPI-distract9810.021.15−0.050.103627−0.221.09−0.25−0.1831880.081.020.040.11<0.001CL NE cl3, cl2 NE cl3
diff-MPI-GAI10900.760.950.700.8240940.260.730.230.283613−0.110.73−0.14−0.09<0.001all different
diff-EQ-5D-index10480.530.300.510.5538860.250.300.240.263351−0.040.30−0.05−0.03<0.001all different
diff-EQ-VAS102236.0620.4734.8037.32383316.2719.8215.6416.893316−3.3119.66−3.98−2.64<0.001all different
diff-sf36-pf105426.7719.0925.6127.92395310.0015.139.5310.473445−3.3115.24−3.81−2.80<0.001all different
diff-sf36-rp104357.3439.5154.9459.74388919.1534.3018.0720.223365−3.6927.14−4.60−2.77<0.001all different
diff-sf36-bp105437.4419.5636.2638.63395414.2914.7913.8314.753444−1.3313.03−1.76−0.89<0.001all different
diff-sf36-gh104425.2619.0924.1026.4239087.7816.337.268.293388−5.0415.87−5.57−4.50<0.001all different
diff-sf36-vt105440.7119.7239.5241.91395113.3217.8312.7713.883433−4.6016.67−5.16−4.05<0.001all different
diff-sf36-sf105340.7323.2639.3242.14395814.4321.2113.7715.093445−7.6622.38−8.40−6.91<0.001all different
diff-sf36-re103153.1044.2750.4055.81384917.2245.0815.7918.643277−9.7945.5911.35−8.22<0.001all different
diff-sf36-mh105432.7718.8131.6433.91394810.6216.3310.1111.133430−6.6617.25−7.24−6.09<0.001all different
Change-Pain10493.320.743.283.3738872.720.812.692.7533732.200.852.172.22<0.001all different
Change-Life situation10493.520.643.483.5639013.050.703.033.0833752.620.832.592.65<0.001all different
LB = Lower Bound; UB = Upper Bound; diff = change in a certain variable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey. For explanations of the subscale abbreviations, see Methods.
Table 8. Pre-IMMRP values for the three clusters based on MIS obtained post-IMMRP.
Table 8. Pre-IMMRP values for the three clusters based on MIS obtained post-IMMRP.
Cluster 1Cluster 2Cluster 3
BaselineNMeanSD95% CI NMeanSD95% CI NMeanSD95% CI ANOVA
Variables LBUB LBUB LBUBp-ValuePost Hoc
Gender22050.220.410.200.2379380.240.430.230.2545230.250.430.240.260.014NA
Age220542.711.242.243.2793843.510.743.343.8452343.010.642.743.30.001cl1 NE cl3, cl2 NE cl3
Outside-Europe21850.110.310.090.1278770.100.300.090.1044840.110.320.100.120.011NA
University21680.280.450.260.2978260.260.440.250.2744450.230.420.220.240.000cl1 NE cl2, cl2 NE cl3
Days no work716889291267511022976103723119541120177411522502103612690.045NA
PRI220513.88.313.414.1793814.48.314.214.6452314.58.314.314.70.002NA
NRS-7days21587.11.77.07.278016.91.76.86.944406.71.76.76.80.000all different
HADS-A219910.44.710.210.678919.04.78.99.144948.24.68.18.40.000all different
HADS-D21979.64.49.49.878928.64.48.58.744947.74.37.67.90.000all different
MPI-Pain-severity21934.50.94.54.678954.40.94.44.445104.30.94.34.30.000all different
MPI-Pain-interfer21874.61.04.64.678534.41.04.44.444854.31.04.24.30.000all different
MPI-LifeCon21952.41.12.42.578802.71.12.72.745022.91.12.93.00.000all different
MPI-Distress21913.91.23.94.078913.51.33.53.545023.21.33.13.20.000all different
MPI-Socsupp21904.21.44.24.378614.21.34.14.244884.11.44.14.20.005NA
MPI-punish20691.91.51.81.973191.81.41.71.841691.71.31.61.70.000all different
MPI-protect20653.01.52.93.172952.91.42.93.041583.01.42.93.00.136NA
MPI-distract20692.61.22.52.673122.51.22.52.541682.51.22.52.60.075NA
MPI-GAI21932.40.82.42.478852.40.82.42.545002.50.82.52.50.000cl1 NE cl2, cl2 NE cl3
EQ-5D-index21260.20.30.20.275870.30.30.20.342770.30.30.30.30.000all different
EQ-VAS209738.018.237.238.8752940.719.040.341.1422443.919.443.344.50.000all different
sf36-pf215152.121.051.252.9768352.820.652.353.2434953.220.252.653.80.119NA
sf36-rp21399.019.88.29.9762911.823.811.312.4428315.326.814.516.10.000all different
sf36-bp215121.113.820.521.7769424.114.323.824.4434726.414.626.026.80.000all different
sf36-gh213939.820.038.940.6761041.420.141.041.9429943.320.442.743.90.000all different
sf36-vt215119.716.619.020.4767923.318.422.923.7434727.319.026.727.80.000all different
sf36-sf215140.023.739.041.0769046.624.946.147.2434852.125.351.452.90.000all different
sf36-re212630.639.628.932.2754941.642.640.742.6422050.643.649.351.90.000all different
sf36-mh215147.320.746.548.2767254.521.254.054.9434360.120.559.460.70.000all different
LB = Lower Bound; UB = Upper Bound; NA = not applicable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; PRI = Pain Region Index. For explanations of the subscale abbreviations, see Methods.
Table 9. Pre-IMMRP values for the three clusters based on MIS obtained at 12-month FU.
Table 9. Pre-IMMRP values for the three clusters based on MIS obtained at 12-month FU.
Cluster 1Cluster 2Cluster 3
BaselineNMeanSD95% CI NMeanSD95% CI NMeanSD95% CI ANOVA
Variables LBUB LBUB LBUBp-ValuePost Hoc
Gender10990.240.430.210.2641230.210.410.200.2236290.250.430.240.27<0.001all different
Age109941.911.141.242.5412343.711.143.444.0362944.210.443.844.5<0.001all different
Outside-Europe10880.090.290.070.1140910.080.280.080.0936060.100.300.090.110.012NA
University10790.290.450.260.3140600.270.450.260.2935670.230.420.210.24<0.001all different except cl2 vs. cl3
Days no work3586612587392930140910462339923116813231270260411301411<0.001all different
PRI109912.88.012.313.2412314.28.314.014.5362914.88.314.515.1<0.001all different
NRS-7days10717.01.76.97.140496.81.76.86.935746.81.76.86.90.049NA
HADS-A10969.74.79.410.041018.74.78.68.936068.44.68.38.6<0.001all different
HADS-D10968.84.48.69.141028.24.38.18.336068.04.47.98.1<0.001all different
MPI-Pain-severity10944.50.94.44.541144.40.94.34.436254.30.94.34.4<0.001cl1NE cl2, cl2 NE cl1, cl3
MPI-Pain-interfer10934.50.94.54.641014.41.04.34.436094.31.04.24.3<0.001all different
MPI-LifeCon10962.51.12.52.641132.81.12.72.836142.91.12.82.9<0.001all different
MPI-Distress10953.81.23.83.941143.41.33.43.536213.31.33.23.3<0.001all different
MPI-Socsupp10904.31.34.24.440974.21.34.24.336044.11.44.14.1<0.001all different
MPI-punish10361.81.41.71.938331.71.31.61.733551.71.31.71.80.047NA
IMP-protect10383.01.42.93.138242.91.42.93.033412.91.42.93.00.583NA
MPI-distract10392.61.22.52.638312.61.22.52.633522.51.22.42.50.003NA
MPI -GAI10942.40.92.42.541082.50.82.52.536232.50.82.52.50.071NA
EQ-5D-index10580.20.30.20.239350.30.30.30.334070.30.30.30.3<0.001all different
EQ-VAS103739.218.238.140.3390541.319.240.742.0338943.119.642.443.7<0.001all different
sf36-pf105651.920.450.753.1396953.020.452.453.7345953.420.152.754.10.120NA
sf36-rp10538.520.67.29.7393812.424.111.613.1341715.326.914.416.2<0.001all different
sf36-bp105821.513.320.722.3397424.414.223.924.8346125.914.025.426.3<0.001all different
sf36-gh105243.121.041.844.4394442.620.542.043.2342442.320.341.743.00.556NA
sf36-vt105821.417.820.322.5397224.719.124.125.3345226.318.625.726.9<0.001all different
sf36-sf105642.424.240.943.9397648.225.647.449.0345851.925.251.052.7<0.001all different
sf36-re104632.140.329.734.5390644.043.242.745.4336349.143.247.650.5<0.001all different
sf36-mh105849.620.948.350.9397056.021.055.356.7345158.720.958.059.4<0.001all different
LB = Lower Bound; UB = Upper Bound; NA = not applicable; NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory; EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; PRI = Pain Region Index. For explanations of the subscale abbreviations, see Methods.
Table 10. Orthogonal Partial Least Square Regressions (OPLS) regressions of MIS post-IMMRP (left part) and at 12-month FU (right part) using the variables pre-IMMRP as regressors. Variables in bold type are significant regressors.
Table 10. Orthogonal Partial Least Square Regressions (OPLS) regressions of MIS post-IMMRP (left part) and at 12-month FU (right part) using the variables pre-IMMRP as regressors. Variables in bold type are significant regressors.
Post-IMMRPVIPp(corr)12-Month FUVIPp(corr)
Variables Pre-IMMRPVariables Pre-IMMRP
sf36-mh1.80−0.80sf36-mh1.63−0.62
MPI-Distress1.720.76MPI-Distress1.590.61
HADS-D1.560.68sf36-sf1.44−0.53
MPI-LifeCon1.48−0.65sf36-re1.39−0.54
HADS-A1.480.65MPI-LifeCon1.35−0.49
sf36-sf1.47−0.63HADS-D1.340.46
sf36-re1.43−0.64HADS-A1.280.46
sf36-vt1.27−0.54MPI-Pain-interfer1.260.39
MPI-Pain-interfer1.260.45Persistent-Pain-duration1.16−0.46
EQ-5D-index1.12−0.41Pain duration1.15−0.45
sf36-bp1.07−0.35EQ-5D-index1.14−0.37
MPI-Pain-severity1.050.28sf36-bp1.11−0.35
EQ-VAS1.04−0.39sf36-vt1.11−0.36
sf36-gh0.95−0.30MPI-Pain-severity0.990.21
NRS-7days0.880.20EQ-VAS0.97−0.28
sf36-pf0.87−0.07sf36-rp0.96−0.37
sf36-rp0.80−0.33PRI0.95−0.26
PRI0.72−0.09sf36-gh0.87−0.10
MPI-GAI0.68−0.20NRS-7days0.810.11
MPI-punish 0.590.25Days no work0.80-0.31
Outside-Europe0.45−0.01sf36-pf0.76−0.05
Days no work0.44−0.16Age0.75−0.30
University0.430.12MPI-GAI0.72−0.22
Persistent Pain duration0.42−0.16University0.540.18
MPI-protect0.41−0.04MPI-punish0.430.12
Pain duration0.38−0.15Outside-Europe0.39−0.03
MPI-Socsupp0.32−0.05MPI-protect0.310.03
MPI-distract0.29−0.01MPI-distract0.260.06
Age0.23−0.10MPI-Socsupp0.220.07
Gender0.04−0.01Gender0.10−0.04
R20.08 R20.08
Q20.08 Q20.07
n14 657 n7 976
CV-ANOVA p-value<0.001 CV-ANOVA p-value<0.001
VIP (VIP > 1.0 is significant) and p (corr) are reported for each regressor. The sign of p (corr) indicates the direction of the correlation with the dependent variable (+ = positive correlation; − = negative correlation). The four bottom rows of each regression report R2, Q2, and p-value of the CV-ANOVA and number of patients included in the regression (n). NRS-7days = Pain intensity as measured by a numeric rating scale for the previous seven days; HADS = Hospital Anxiety and Depression Scale; MPI = Multidimensional Pain Inventory EQ-5D-index = The index of the European quality of life instrument; EQ-VAS = The European quality of life instrument thermometer-like scale; sf36 = The Short Form (36) Health Survey; PRI = Pain Region Index. For explanations of the subscale abbreviations see Methods.

Share and Cite

MDPI and ACS Style

Ringqvist, Å.; Dragioti, E.; Björk, M.; Larsson, B.; Gerdle, B. Moderate and Stable Pain Reductions as a Result of Interdisciplinary Pain Rehabilitation—A Cohort Study from the Swedish Quality Registry for Pain Rehabilitation (SQRP). J. Clin. Med. 2019, 8, 905. https://doi.org/10.3390/jcm8060905

AMA Style

Ringqvist Å, Dragioti E, Björk M, Larsson B, Gerdle B. Moderate and Stable Pain Reductions as a Result of Interdisciplinary Pain Rehabilitation—A Cohort Study from the Swedish Quality Registry for Pain Rehabilitation (SQRP). Journal of Clinical Medicine. 2019; 8(6):905. https://doi.org/10.3390/jcm8060905

Chicago/Turabian Style

Ringqvist, Åsa, Elena Dragioti, Mathilda Björk, Britt Larsson, and Björn Gerdle. 2019. "Moderate and Stable Pain Reductions as a Result of Interdisciplinary Pain Rehabilitation—A Cohort Study from the Swedish Quality Registry for Pain Rehabilitation (SQRP)" Journal of Clinical Medicine 8, no. 6: 905. https://doi.org/10.3390/jcm8060905

APA Style

Ringqvist, Å., Dragioti, E., Björk, M., Larsson, B., & Gerdle, B. (2019). Moderate and Stable Pain Reductions as a Result of Interdisciplinary Pain Rehabilitation—A Cohort Study from the Swedish Quality Registry for Pain Rehabilitation (SQRP). Journal of Clinical Medicine, 8(6), 905. https://doi.org/10.3390/jcm8060905

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop