Ò
PAIN 155 (2014) 179–189
www.elsevier.com/locate/pain
Dose or content? Effectiveness of pain rehabilitation programs for patients with chronic low back pain: A systematic review Franka P.C. Waterschoot a,⇑, Pieter U. Dijkstra a,b, Niek Hollak c, Haitze J. de Vries a, Jan H.B. Geertzen a, Michiel F. Reneman a a b c
Department of Rehabilitation Medicine, Center for Rehabilitation, University Medical Center Groningen, University of Groningen, The Netherlands Department of Oral and Maxillofacial Surgery, University Medical Centre Groningen, University of Groningen, The Netherlands Center for Human Movement Sciences, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
Sponsorships or competing interests that may be relevant to content are disclosed at the end of this article.
a r t i c l e
i n f o
Article history: Received 9 February 2013 Received in revised form 2 October 2013 Accepted 4 October 2013
Keywords: Chronic low back pain Dose Effectiveness Pain rehabilitation programs
a b s t r a c t We sought to systematically analyze the influence of dose of pain rehabilitation programs (PRPs) for patients with chronic low back pain (CLBP) on disability, work participation, and quality of life (QoL). Literature searches were performed in PubMed, Cochrane Library, Cinahl, and Embase up to October 2012, using MeSH terms, other relevant terms and free-text words. Randomized controlled trials in English, Dutch, and German, analyzing the effect of PRPs, were included. One of the analyzed interventions had to be a PRP. Outcomes should be reported regarding disability, work participation, or QoL. To analyze dose, the number of contact hours should be reported. Two reviewers independently selected titles, abstracts, and full-text articles on the basis of inclusion and exclusion criteria. Data were extracted and risk of bias was assessed. Effect sizes (ES) were calculated for each intervention, and influence of dose variables was analyzed by a mixed model analysis. Eighteen studies were identified, reporting a wide variety of dose variables and contents of PRPs. Analyses showed that evaluation moment, number of disciplines, type of intervention, duration of intervention in weeks, percentage of women, and age influenced the outcomes of PRPs. The independent effect of dose variables could not be distinguished from content because these variables were strongly associated. Because dose variables were never studied separately or reported independently, we were not able to disentangle the relationship between dose, content, and effects of PRPs on disability, work participation, and QoL. Ó 2013 International Association for the Study of Pain. Published by Elsevier B.V. All rights reserved.
1. Introduction Multidisciplinary pain rehabilitation programs (PRPs) are effective in improving daily functioning of patients with chronic low back pain (CLBP) [12,14,30]. Most studies investigating the effects of PRPs focused on the relationship between therapy content and effect. Guidelines for managing CLBP are based on evidence obtained from these studies [1,27]. However, this evidence might be biased. In a literature search, we could not identify any study analyzing the relationship between dose and effect as a primary objective. Two systematic reviews [14,30] were identified that analyzed dose and effect as a secondary objective. These reviews presented conflicting conclusions on dose effects. ⇑ Corresponding author. Address: Department of Rehabilitation Medicine, Center for Rehabilitation, University Medical Center Groningen, University of Groningen, PO Box 30.002, 9750 RA Haren, The Netherlands. Fax: +31 505338258. E-mail address:
[email protected] (F.P.C. Waterschoot).
In a review by Guzman et al. [14], 10 randomized controlled trials (RCTs) were included reporting on 12 PRPs. PRPs were divided into 2 categories: daily intensive programs with more than 100 h of therapy, and once- or twice-weekly programs with less than 30 h of therapy. On the basis of that distinction, it appeared that multidisciplinary PRPs of more than 100 h were superior to monodisciplinary treatment, and PRPs of less than 30 h were not. The authors concluded that intensive multidisciplinary PRPs are superior to less intensive multidisciplinary PRPs [14]. These results have been used in guidelines and clinical practice [1,31]. However, it may be debated whether the conclusion regarding dose is valid to support its clinical implications. First, there is an absence of knowledge about PRPs with doses between 30 to 100 h. Second, no distinction was made within intensive and less intensive PRPs, although there was a wide variety within the groups. Less intensive PRPs ranged from 17.5 to 30 h, while intensive PRP ranged from 100 to 280 h. Third, the review [14] was designed to assess the effect of different PRPs on clinically relevant outcomes. All PRPs differed not only in dose but also in content, setting, and number
0304-3959/$36.00 Ó 2013 International Association for the Study of Pain. Published by Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.pain.2013.10.006
180
Ò
F.P.C. Waterschoot et al. / PAIN 155 (2014) 179–189
of disciplines involved, which may have confounded the dose aspect of the conclusion. Additionally, the authors discussed whether the improvements gained with intensive PRP were worth the expense. The review of van Geen et al. [30] also had some limitations. It assessed the long-term effect of PRPs with different contents and doses. Ten studies were included, and a distinction was made between 30 h of training a week or more (intensive therapy) vs less than 30 h of training a week (low-intensive therapy). Both the intensive and low-intensive PRPs showed positive effects. The conclusion in that review was that the dose of the intervention had no substantial influence on the effectiveness of the intervention. In summary, while many studies provide evidence for effectiveness of multidisciplinary PRP, dose of multidisciplinary PRP to achieve these effects remains unclear. The objective of the current study was to analyze the influence of dose variables on the outcome of PRPs for patients with CLBP.
Table 1 Criteria risk of bias analysis. 1 2
3 4 5 6 7 8
9 10 11 12
Was the method of randomization adequate? Was the treatment allocation concealed? Was knowledge of the allocated interventions adequately prevented during the study? Was the patient blinded to the intervention? Was the care provider blinded to the intervention? Was the outcome assessor blinded to the intervention? Were incomplete outcome data adequately addressed? Was the drop-out rate described and acceptable? Were all randomized participants analyzed in the group to which they were allocated? Are reports of the study free of suggestion of selective outcome reporting? Other sources of potential bias: Were the groups similar at baseline regarding the most important prognostic indicators? Were cointerventions avoided or similar? Was the compliance acceptable in all groups? Was the timing of the outcome assessment similar in all groups?
2. Methods Publications were retrieved by computer-aided search on PubMed, the Cochrane Library, Cinahl, and Embase up to October 2012. A specific search was developed using MeSH terms and other relevant terms for each database. The PubMed search is described in the Supplement. Refworks was used to store the results of the searches and to remove duplicates. 2.1. Selection of studies Selection criteria were applied independently by 2 reviewers (FW and NH). The retrieved studies were first selected by title and abstract. Doubtful cases were discussed by the reviewers and included or excluded for full-text analysis by consensus. Full-text reports of studies eligible for inclusion were analyzed. Disagreements were resolved by consensus or when necessary by a third reviewer (HdV). Studies were selected on the basis of the following inclusion criteria: (1) RCTs written in English, Dutch, or German; (2) the objective was to assess effectiveness of a multidisciplinary PRP for patients with CLBP; PRP was defined as a rehabilitation program on the basis of the biopsychosocial model [33] with 3 or more disciplines providing the program (with or without a medical doctor); (3) total number of contact hours of PRP was described; (4) participants were between 18 and 65 years with disabling nonspecific CLBP for at least 3 months; and (5) outcome variables were described in the domain of disability, work participation, or quality of life (QoL). Studies were excluded if: (1) the multidisciplinary PRP was given in primary care; (2) the objective was to assess effectiveness of only a biomedical intervention; (3) the study included participants who were diagnosed with specific disorders or severe comorbidities interfering with PRP, such as heart failure, rheumatoid arthritis, or psychiatric disorders; and (4) only total duration of PRP was reported. 2.2. Data management Risk of bias of included studies was assessed according to the Cochrane Back Review Group [5] by 2 reviewers independently (FW and HdV). The criteria are presented in Table 1. Each criterion was scored as positive (Y), negative (N), or unclear (U). The total score was computed by counting the number of criteria scored as positive. Studies with a score of 6 or higher were defined as low risk of bias; a score lower than 6 was defined as high risk of bias. In case of unclear scores, corresponding authors of the studies were contacted by e-mail.
Performance bias was analyzed with different items of blinding. The item regarding blinding of care providers for intervention is frequently impossible in nondrug trials. To analyze the influence of lack of blinding of care providers on the judgment of methodological quality, we performed a sensitivity analysis by excluding item 4 of the risk of bias analysis. A data extraction form was developed and piloted before data extraction. Data were extracted by 2 reviewers independently (FW and HdV). Disagreement was resolved by consensus or, if necessary, by a third reviewer (MR). Data were extracted on general study, participant, and dose characteristics (including total number of contact hours and total duration of the treatment in weeks) and treatment content (including description of treatment, treatment components, and number of disciplines). Reported outcome measures were categorized into disability, work participation, and QoL. The following interventions were distinguished: PRP, no treatment, care as usual (CAU), individual physical treatment, individual psychological treatment, surgery, and multidisciplinary treatment not defined as PRP. For each intervention, the effect size (ES) was calculated by subtracting the posttreatment mean from the pretreatment mean (for each evaluation moment separately), divided by the pretreatment standard deviation (SD). If the mean change (pretreatment posttreatment) was reported, ES was calculated by dividing the mean change by the SD of the change. If outcome was reported as a proportion (eg, work participation), ES was calculated according to Hojat et al. [17]. If medians and ranges were reported [3–6], means and SDs were estimated [18]. Because of the small sample sizes of the different interventions, the ES were corrected with Hedges J as described by Borenstein et al. [7]. To analyze the influence of dose variables on ES, a linear mixed effect model was applied by SPSS version 18.0 software. Analyses were applied for all interventions on each outcome category separately. The study of Hellum et al. [15] was the only one comparing effects of PRP with those of surgery. In that study, patients with CLBP were included, but it was not clear whether they had a specific or nonspecific diagnosis. Because of the unique character of that study, the analyses were done twice: with and without this study for outcome categories disability and QoL. ES per intervention was the response variable. Predictor variables were PRP (yes/no), number of disciplines involved in the program, number of contact hours, treatment duration in weeks, type of intervention (no treatment, CAU, surgery, individual physical therapy [exercise], psychological treatment, other multidisciplinary treatment, PRP), evaluation moment (in months), age, percentage of women, and
Ò
F.P.C. Waterschoot et al. / PAIN 155 (2014) 179–189
risk of bias (high/low). These variables were entered one by one into the model. The variables remained in the model if the regression coefficient was significant (Y % >Y
X X
X X
X X
13 0
Disabilityi Disabilityi
Pre Pre
Post Post
0.473 0.006
Mangels [24]
2009
Traditional orthopedic rehabilitation treatment
131
48.7
14.7
78.6
X
X
X
3
W
196.8
Disabilityg
Pre
Post
0.301
QoLu
Pre
Disabilityg
Pre
12 Post 12 Post
0.333 0.419 0.268 0.487
QoLu
Pre
Disabilityg
Pre
12 Post 12 Post 12
0.344 0.550 0.396 0.370 0.306
QoLu
Pre
Post 12
0.490 0.511
Disabilityb,r
Pre
Post 6
0.602 0.515
12 24 Post 24
0.455 0.549 0.231 0.391
Multidisciplinary rehabilitation
Multidisciplinary rehabilitation + booster sessions
Dufour [10]
2010
Group based multidisciplinary biopsychosocial rehabilitation program
113
119
129
49.5
48.3
46.6 53.4
9.0
15.8
78.8
75.6
56.6
X
X
X
X
X
X
X
X
X
4
4 + 12 FU
12
W
W + FU M
W
223.2
223.2
87
Work
m
Pre
185
(continued on next page)
Ò
n
F.P.C. Waterschoot et al. / PAIN 155 (2014) 179–189
Intervention
186
Table 3 (continued) Study
Year
Intervention
n
Age, y
Intensive individually therapist-assisted back muscle strengthening exercises
143
49.7 50.3
SD
Months (M)/ years (Y)
SD
Total duration
Week (W) or month (M)
Total h
Outcome domain
Evaluation moment, mo
ES
55.9
X
X
X
12
W
46
Disabilityb,r
Pre
Post 6
0.359 0.289
12 24 Post 24
0.211 0.236 0.041 0.362
1.5 3 6 12 24 1.5 3 6 12 24 1.5 3 6 12 24 1.5 3 6 12 24
1.341 1.144 0.726 1.379 1.527 0.894 0.671 0.766 0.862 1.054 1.121 1.787 1.870 2.076 2.131 0.957 1.320 1.254 1.221 1.254
Work Hellum [15]
2011
Rehabilitation
87
40.8
7.1
58.6
X
X
X
3 to 5
W
60
m
Disabilityb,e
Pre Pre
QoLs
Surgery
86
41.1
7.1
46.5
X
X
X
X
X
X
Disabilityb,e
Pre
QoLs
Meng [25]
Roche-Leboucher [26]
2011
2011
New biopsychosocial back school program
197
50.2
7.6
65.2
51.5
% >5 Y
X
24.3 (4.8)
Mean days
6.4
Disabilityj
Pre
6
0.702
Traditional back school program (CAU)
185
49.5
7.7
63
39.1
% >5 Y
X
24.9 (5.2)
Mean days
5.5
Disabilityj
Pre
12 6
0.377 0.569
12
0.420
Functional restoration program
68
40.8
7.4
32.4
X
X
X
5
W
150
Disabilityf
Pre
12
1.747
Pre Pre Pre
12 12 12
1.037 0.652 0.751
o
Active individual therapy
64
38.7
6.1
37.5
X
X
X
5
W
15
Work Disabilityf Work
X, no data available; ES, effect size; FU, follow-up; QoL, quality of life; TAU, treatment as usual. Median. Measurement instrument used for operationalization of outcome domain: a, Function; b, SF-36 Physical component; c, SF-36 Physical Functioning + SF-36 Role Physical; d, RMDQ; e, Oswestry; f, Pain impact; g, PDI; h, PCS; i, QBPDS; j, Physical activity (Freiburger); k, Million Index; l, Funktionsfahigkeit; m, work readiness; n, Return to work; o, percentage full time at work; p, work participation; q, subjective work capacity; r, Arbeidsfähigkeit; s, Euroqol; t, general well-being; u, Life Satisfaction Questionnaire; v, Somatische Gesundheit. *
Ò
Duration pain
F.P.C. Waterschoot et al. / PAIN 155 (2014) 179–189
% Women
Ò
187
F.P.C. Waterschoot et al. / PAIN 155 (2014) 179–189 Table 4 Ranges of ES per outcome category for PRPs and control interventions. Outcome
Disability Work participation Quality of life
Including study of Hellum et al. [15]
Excluding study of Hellum et al. [15]
PRP
Controls
PRP
Controls
0.306 to 1.960 0.231 to 1.535 0.147 to 1.054
0.330 to 2.131 0.035 to 1.411 0.114 to 1.320
0.306 to 1.960 0.231 to 1.535 0.147 to 0.610
0.330 to 1.307 0.035 to 1.411 0.114 to 0.574
ES, effect size; PRP, pain rehabilitation program.
Table 5 Results of mixed effect model analysis for the different outcome variables, in and excluding data from the study of Hellum et al. [15]. Outcome
Variable
SE b
P
Disabilitya
Interceptr No treatment CAU Surgery Individual therapy/exercise Psychological treatment Other multidisciplinary treatment 1.5–6 mo FU 12 mo FU 12–24 mo FU Percentage women Interceptr Number of disciplines PRPr CAU Surgery Individual therapy/exercise Other multidisciplinary treatment
0.501 0.617 0.719 0.748 0.142 0.496 0.294 0.173 0.182 0.407 0.002 0.475 0.088 0.538 0.241 0.372 0.011 0.147
0.167 0.383 0.176 0.218 0.108 0.176 0.198 0.116 0.142 0.131 0.001 0.116 0.025 0.104 0.112 0.071 0.096 0.099
.004 .109