2. Materials and Methods
Study 1
A group of six participants with lower back pain, pelvic pain, or fibromyalgia were recruited through referrals, posters, and emails. Inclusion criteria included a minimum age of eighteen, internet and app store access, and pain persisting for at least six months. Participants utilized a hand-held biofeedback device for ten minutes at least twice daily and reported their anxiety and pain levels in a mobile app before and after each use over the course of four weeks. Hand temperature data was also recorded. Along with the biofeedback device, participants attended a weekly educational Zoom meeting and were provided with the book “The Way Out” by Alan Gordon. Pain, anxiety, and satisfaction and recovery changes were measured at the end of the four weeks [
19].
Study 2
A group of nine participants with chronic pain were recruited through referrals, posters, and emails. This study followed a similar protocol to Study 1, however, the inclusion criteria expanded to include all types of chronic pain.
An in-depth description of materials and methods for Studies 1 and 2 has been previously published [
19].
Study 3
A group of six participants with chronic musculoskeletal pain lasting longer than 6 months was recruited through clinician referrals, email lists, and posters. Inclusion criteria included a minimum age of eighteen, musculoskeletal pain or fibromyalgia persisting longer than 6 months, and a minimum average pain level of 4 measured by the Visual Analogue Scale [
20,
21]. Participant demographics including age, sex, race, and pain history were collected (
Table 1). Once participants consented to the study, they were assigned a weekly group meeting time.
The study lasted 5 weeks, beginning with an initial in-person interview followed by Zoom meetings once a week for the entire study duration. Each Zoom meeting consisted of educational PowerPoint presentations followed by group discussions on the presented content. Participants were also assigned educational videos and worksheets to complete before each weekly Zoom meeting to reinforce their understanding of the material. Participants were provided with a hand-temperature biofeedback device for use at least twice daily. The biofeedback device was used with the book “The Way Out” by Alan Gordon to learn concepts important to chronic pain recovery. Measurements of pain, anxiety, and satisfaction and recovery were reported in the initial applicant survey, at the start and end of the trial, and a 6-month follow-up. Pain levels were measured using the Short-Form McGill Pain Questionnaire (SF-MPQ) [
22,
23,
24]. Anxiety was measured on the Short-Form State-Trait Anxiety Inventory (STAI) [
25,
26] and health-related satisfaction and recovery was evaluated on the Satisfaction and Recovery Index (SRI) [
27,
28,
29].
The thermal biofeedback device used in this study monitored changes in hand temperature. The device was designed for portability and was small enough to fit in the user’s palm for use in any setting. Its display comprised six lights which collectively represented a range of temperatures. The lights changed color one by one in a clockwise direction as hand temperature increased, indicating greater blood flow as a result of vasodilation that occurs during relaxation. Participants used the device for at least two ten-minute sessions per day for 5 weeks.
Study 4
A group of four participants with chronic pain and headaches persisting longer than 6 months were recruited through clinician referrals, emails, and posters. Participants were included if they met the inclusion criteria for having pain lasting longer than 6 months, a minimum age of 18, and a minimum average pain of 4 measured on the VAS. Participant demographics including age, sex, race, and pain history were collected.
This 5-week-long study followed the same protocol as Study 3 and utilized the same biofeedback device. Pain, anxiety, and satisfaction and recovery measurements were collected before and after the trial, and 6 months later.
Figure 1.
Flow diagram of participant selection and study timeline for Studies 3-5 from initial screening to 6-month follow-up.
Figure 1.
Flow diagram of participant selection and study timeline for Studies 3-5 from initial screening to 6-month follow-up.
Table 1.
Demographics and pain history for participants across all five studies. The studies listed as “Study 1 and Study 2” in this paper were previously listed as “Study 2 and Study 3,” respectively, in the previously published paper. We have renumbered the studies because data from Study 1 was excluded from the research we have done in this paper.
Table 1.
Demographics and pain history for participants across all five studies. The studies listed as “Study 1 and Study 2” in this paper were previously listed as “Study 2 and Study 3,” respectively, in the previously published paper. We have renumbered the studies because data from Study 1 was excluded from the research we have done in this paper.
Characteristic |
Total (n=39) |
Study 1 (n-6) |
Study 2 (n=9) |
Study 3 (n=6) |
Study 4 (n=4) |
Study 5 (n=14) |
Age, mean (SD) |
53.22 (16) |
38.7 (16.4) |
45.9 (13.2) |
63.3 (13.6) |
50.5 (4.8) |
49.9 (18.8) |
Sex |
|
|
|
|
|
|
Male |
13 |
4 |
1 |
2 |
1 |
5 |
Female |
26 |
2 |
8 |
4 |
3 |
9 |
Race |
|
|
|
|
|
|
Caucasian |
33 |
5 |
7 |
6 |
4 |
11 |
Asian/Pacific Islander |
3 |
1 |
1 |
0 |
0 |
1 |
African American |
2 |
0 |
1 |
0 |
0 |
1 |
American Indian/Alaska Native |
0 |
0 |
0 |
0 |
0 |
0 |
Other |
1 |
0 |
0 |
0 |
0 |
1 |
Number of Pain Conditions |
|
|
|
|
|
1 |
14 |
2 |
3 |
1 |
3 |
5 |
2 |
9 |
2 |
3 |
1 |
0 |
3 |
3+ |
16 |
2 |
3 |
4 |
1 |
6 |
Duration of Pain |
|
|
|
|
|
6-12 months |
3 |
1 |
1 |
0 |
0 |
1 |
1-5 years |
15 |
2 |
4 |
1 |
2 |
6 |
>5 years |
21 |
3 |
4 |
5 |
2 |
7 |
Medications taken |
|
|
|
|
|
Opioids |
3 |
0 |
0 |
2 |
0 |
1 |
NSAID/Acetaminophen |
9 |
2 |
3 |
3 |
0 |
1 |
Adjunctive Pain Medication |
3 |
1 |
0 |
0 |
0 |
2 |
None/Did not specify |
24 |
1 |
5 |
0 |
2 |
6 |
Study 5
A group of 14 participants with chronic pain and headaches were recruited through clinician referrals, emails, and posters. Inclusion criteria included pain lasting longer than 6 months, a minimum age of eighteen, and a minimum average pain of 4 measured on the VAS. Participant demographics including age, sex, race, and pain history were collected.
This study followed the same protocol as the previous two studies but used the book “The Explain Pain Handbook: Protectometer” by Lorimer Mosely rather than “The Way Out” by Allen Gordon. Pain, anxiety, and satisfaction and recovery measurements were collected before and after the trial, and 6 months later.
Statistical Analysis
To assess differences in pain outcomes between studies which frequently assessed pain levels vs. those that did not, we fit hierarchical Bayesian models, allowing us to pool data across studies while accounting for intra-/inter-study variability in outcomes (cf. meta-analytic or meta-regressive modeling approaches). Models were specified and estimated using Stan (
https://mc-stan.org/) and brms [
30] in R. Each model followed the same generalized linear form, where differences in pain outcomes (Final - Baseline) were regressed against a categorical predictor indicating whether pain levels were frequently vs. infrequently assessed during the course of the biofeedback program (formally, coded using a sum-to-zero deviation contrast, capturing the mean difference in outcomes between study types). Baseline scores were also included as a covariate to control for the fact that the magnitude of the Final - Baseline difference is likely correlated with initial pain levels. Finally, we further included a random intercept term conditioned on study ID (broadly accounting for inter-study variability in outcomes).
We specified a series of weakly-informative priors over each parameter, robustly-defined in part by the median absolute deviation (MAD) of each outcome measure. For example, the prior densities over the slope parameters were defined by a normal (Gaussian) distribution with mean = 0 and standard deviation = 2.5 x MAD(dy), where dy gives the vector of Final - Baseline difference scores. The model intercept was defined using a Student t-distribution with 3 degrees of freedom and location/scale parameters given by median(dy) and MAD(dy), respectively. Residual variances and random effects variances were similarly defined using a half-Student t-distribution, with 3 degrees of freedom and location/scale parameters given by 0 and MAD(dy), respectively (i.e., constraining the variances to be positive).
Robust exploration of the posterior parameter space for all models was achieved using Hamiltonian Monte Carlo (HMC), with four independent chains comprised of 15000 sampling iterations each (5000 of which were used as warm-up samples). We ensured that all chains properly converged to equilibrium for all parameters using classical benchmarks, including: the effective sample size (considering the autocorrelation between successive HMC transitions); the Monte Carlo standard error (relative to the posterior SD); R-hat (the variance ratio between each chain relative to all chains); and no divergent HMC transitions after warm-up.
For statistical inference, we consider 95% credibility intervals around the posterior mean parameter estimates (using the highest posterior density interval, HDI) along with the ‘probability of direction’ (
pd), capturing the proportion of the posterior density above or below zero. Here,
pd can range from 0.50 (indicating that the posterior density is perfectly-centered at zero) to 1 (in the most extreme case, that
no proportion of the posterior density intersects zero)—thus,
pd values closer to 1 indicate more posterior evidence for the presence of a nonzero effect. We also provide a fully-probabilistic Bayesian analogue of Cohen’s
d to further quantify standardized mean differences (i.e.,
effect sizes): following conventions for hierarchical linear models [
31], here we took the posterior draws corresponding to the mean difference between study types, or measurement occasions, and divided them by the square root of the summed residual variances and random effects variances.
3. Results
Biofeedback-related changes in self-reported pain have previously been detailed for Studies 1 and 2 [
19]. In brief, we found that hand temperature biofeedback programs were associated with reduced pain and anxiety on average (comparing reported pain and anxiety levels before and after each biofeedback session)—and in some individuals, the magnitude of the change in anxiety predicted the magnitude of the change in pain, such that greater reductions in anxiety after each session were associated with greater reductions in pain. We also showed modest, overall reductions in pain and anxiety metrics when comparing ratings given at baseline (prior to beginning the biofeedback program) to those provided after finishing the program.
Our first goal was, therefore, to determine whether we could replicate our prior results [
19] now in Studies 3-5, comparing Final vs. Baseline scores on a number of self-assessed metrics, including: average, minimum, and maximum pain levels (rated on a 0-10 VAS); the McGill Pain Inventory; the Satisfaction and Recovery Index; and the State-Trait Anxiety Inventory. A full set of descriptive statistics for each of these scales is shown in
Table 2. With respect to VAS ratings (
Figure 2A), we observed the largest reductions in self-assessed average pain levels (
d = -0.92, 95CI = [-1.47, -0.34];
b = -1.77,
SD = 0.48, 95CI = [-2.72, -0.85],
pd = 1), followed by maximum pain (
d = -0.70, 95CI = [-1.18, -0.22];
b = -1.73,
SD = 0.52, 95CI = [-2.77, -0.69],
pd = 0.999) and minimum pain levels (
d = -0.65, 95CI = [-1.05, -0.22];
b = -1.65,
SD = 0.44, 95CI = [-2.52, -0.78],
pd = 1). These findings align with our prior work showing consistent reductions in self-assessed pain at multiple levels.
With respect to our additional survey measures (
Figure 2B), we both replicated and exceeded the effects previously reported in Studies 1-2. Where we initially observed a small-to-moderate reduction in anxiety (
d = -0.39), we now see evidence for a much larger effect (
d = -1.00, 95CI = [-1.50, -0.51];
b = -6.46,
SD = 1.23, 95CI = [-8.90, -4.02],
pd = 1). Using the Satisfaction and Recovery Index, our last pair of studies failed to identify any evidence for a positive effect (
d = 0.04); however, in the current group of studies, we observed moderate increases in participants’ health-related satisfaction (
d = -0.46, 95CI = [0.02, 0.93];
b = 9.27,
SD = 4.23, 95CI = [0.90, 17.64],
pd = 0.984)—although we note the posterior variance is quite large here, indicating greater uncertainty. Finally, the reduction in pain assessed by the McGill Pain Inventory exceeded Studies 1-2 (
d = 0.30) nearly fourfold (
d = -1.15, 95CI = [-1.78, -0.55];
b = -8.74,
SD = 1.70, 95CI = [-12.14, -5.35],
pd = 1).
Taken together, these results anecdotally suggest that frequent pain assessments (i.e., before and after each biofeedback session, as in Studies 1-2) may partially hinder recovery, producing weaker effects than those observed when participants are not regularly asked to assess their pain (as in Studies 3-5). We therefore sought to explicitly test this hypothesis by contrasting average reductions in pain between studies using different types of experimental designs. Here we focus primarily on pain assessments derived via the McGill Pain Inventory, which showed the largest effect in Studies 3-5 and appeared to drastically outperform what we previously observed in Studies 1-2 [
19].
When considering the variability both within and between study groups, we were unable to identify a reliable difference between these two design approaches. While biofeedback studies that eschewed frequent pain assessment generally showed stronger reductions in pain, there was nevertheless considerable overlap with studies that
did frequently assess pain (
Figure 3A). Consistent with this, under the model, we observed a moderately large difference between studies, such that infrequent pain assessment was associated with greater changes from Baseline to Final ratings (
d = 0.71, 95CI = [-0.27, 1.74];
b = 4.73,
SD = 3.64, 95CI = [-2.73, 11.97],
pd = 0.921)—while this trend is indeed in the direction we anticipated, the 95% credibility interval around the effect intersects zero, so we cannot conclude that the difference in pain reduction between studies is statistically meaningful (
Figure 3B).
Notably, the trend described above was specific to pain-related outcomes. When assessing changes in trait anxiety levels before and after the biofeedback program (
Figure 4A), we saw no difference between study types—both groups demonstrated similar degrees of anxiety reduction, despite the fact that Studies 1-2 also asked participants to assess their anxiety levels before and after each biofeedback session (
d = 0.04, 95CI = [-1.05, 1.07];
b = 0.16,
SD = 3.38, 95CI = [-7.11, 6,70],
pd = 0.543;
Figure 4B).
While we are ultimately unable to precisely determine the extent to which frequent vs. infrequent pain assessments affect pain outcomes, it is worth noting that the magnitude of the difference is still quite large (
Figure 5). These initial results appear promising despite being a purely retrospective analysis across a number of smaller pilot studies with varying degrees of heterogeneity. It will be essential for future studies to more systematically test these differences in a more controlled fashion.
4. Discussion
This study investigated the impact of frequent pain assessment on chronic pain recovery across five pilot studies. In the studies without frequent pain assessment, we found the reduction in pain measured on the SF-MPI to be almost twice as great compared to those with frequent pain assessment. While this trend was not statistically significant due to large error variance and small sample size, the substantial difference in pain decrease between the two groups strongly suggests that frequent pain assessment may hinder pain reduction.
In the studies without frequent pain assessment, participants were also instructed to avoid discussing their pain with others to further minimize pain attentiveness. We believe this allowed for more effective pain reduction by preventing additional reinforcement of negative thought patterns involved in chronic pain. Frequent pain assessment may catalyze RNT and thus catastrophization, which delays pain recovery by maintaining chronic pain-causing networks in the brain [
32,
33]. Our results are consistent with prior studies that have demonstrated the role of maladaptive thought patterns in the persistence and exacerbation of chronic pain [
34,
35,
36]. Each time participants reassess their pain and are reminded of it, sensitization of neural patterns that create the experience of pain is reinforced [
10,
11,
12,
13].
Additionally, much research has demonstrated the role of operant conditioning in maintaining behaviors and neural patterns involved in chronic pain [
37,
38,
39]. While further research is necessary to determine statistical significance, our results supplement the findings in these studies [
37,
38,
39] as it’s likely that consistent pain assessment, four times daily, over multiple weeks can become a source of negative reinforcement.
There are several limitations to this study. The small sample size and the need for a control group across all five pilot studies limits the generalizability of our findings. These studies also do not address the implications of additional comorbidities that may affect chronic pain recovery. Furthermore, the inclusion criteria in Studies 3-5 was expanded to accept participants with additional types of chronic pain, whereas Studies 1 and 2 accepted exclusively musculoskeletal pain and fibromyalgia participants. The variance in pain inclusion criteria likely caused the high standard deviation in our data and decreased the effect size. While our findings suggest a considerably greater pain reduction in studies without frequent pain assessment, further research with a larger sample population and consistent inclusion criteria is required to establish statistical significance.
The trend we have observed in this research, though not statistically significant, suggests that eliminating frequent pain assessments will improve the efficacy of chronic pain recovery programs.