- Research
- Open access
- Published:
Advancing the integration of biosignal-based automated pain assessment methods into a comprehensive model for addressing cancer pain
BMC Palliative Care volume 23, Article number: 198 (2024)
Abstract
Background
Tailoring effective strategies for cancer pain management requires a careful analysis of multiple factors that influence pain phenomena and, ultimately, guide the therapy. While there is a wealth of research on automatic pain assessment (APA), its integration with clinical data remains inadequately explored. This study aimed to address the potential correlations between subjective and APA-derived objectives variables in a cohort of cancer patients.
Methods
A multidimensional statistical approach was employed. Demographic, clinical, and pain-related variables were examined. Objective measures included electrodermal activity (EDA) and electrocardiogram (ECG) signals. Sensitivity analysis, multiple factorial analysis (MFA), hierarchical clustering on principal components (HCPC), and multivariable regression were used for data analysis.
Results
The study analyzed data from 64 cancer patients. MFA revealed correlations between pain intensity, type, Eastern Cooperative Oncology Group Performance status (ECOG), opioids, and metastases. Clustering identified three distinct patient groups based on pain characteristics, treatments, and ECOG. Multivariable regression analysis showed associations between pain intensity, ECOG, type of breakthrough cancer pain, and opioid dosages. The analyses failed to find a correlation between subjective and objective pain variables.
Conclusions
The reported pain perception is unrelated to the objective variables of APA. An in-depth investigation of APA is required to understand the variables to be studied, the operational modalities, and above all, strategies for appropriate integration with data obtained from self-reporting.
Trial registration
This study is registered with ClinicalTrials.gov, number (NCT04726228), registered 27 January 2021, https://classic.clinicaltrials.gov/ct2/show/NCT04726228?term=nct04726228&draw=2&rank=1
Background
Cancer pain is a prevalent concern among cancer patients. Significantly, according to the World Health Organization (WHO), 55% of individuals undergoing anti-cancer treatment suffer from moderate to severe pain. Furthermore, an even larger percentage, potentially reaching 70%, of those facing advanced or terminal illnesses, experience pain [1, 2].
The challenge of pain management lies in tailoring effective strategies that address the unique nature of each individual's pain experience while considering factors such as pain intensity, underlying causes, and personal responses to treatments [3,4,5]. The pain experience often reflects an intricate interplay of multiple factors, such as psychosocial elements, levels of distress, tumor attributes, disease advancement, and the underlying pain pathophysiology[4]. In this complex scenario, pain assessment is the cornerstone for planning effective pain management. Nevertheless, this crucial step is often challenging. The subjective nature of pain presents one of the most intricate aspects. The description of pain can vary significantly among individuals, highlighting the importance of employing methods that can accurately capture this diversity. Additionally, cultural influences play a role in how pain is expressed and communicated, leading to potential misunderstandings or underreporting [6]. Furthermore, certain populations, such as children [7, 8] or those with cognitive impairments [9], present unique challenges, requiring specialized approaches to accurately assess their pain experiences [10, 11].
Self-report measurement instruments, including the visual analog scale (VAS) and the numeric rating scale (NRS), are frequently employed to evaluate the intensity of pain. These scales heavily rely on the individual's ability to accurately recall and convey their pain intensity, which can be affected by memory biases, emotional states, and cognitive impairments. Hence, due to the susceptibility of subjective methods to individual interpretations and biases, they fall short of fully encompassing the intricate complexity of pain [12].
Automatic Pain Assessment (APA) is focused on objective measures for the assessment of pain intensity, offering a more unbiased substitute to subjective pain scales [13]. APA methodologies encompass a variety of behavioral and physiological approaches. Behavioral approaches include the analysis of facial expressions [14], linguistic patterns using qualitative and quantitative analyses [15], and the observation of non-verbal physical behaviors such as body movements and gestures [16]. On the other hand, physiological indicators concern the study of biosignals such as electrocardiogram (ECG) and electrodermal activity (EDA), as well as advanced brain imaging techniques such as functional magnetic resonance imaging (fMRI) and electroencephalography (EEG), for addressing pain-related brain activity [17]. In particular, biosignals like EDA and ECG could also offer the potential of being used in continuously monitoring settings and easily available as wearable devices, thus allowing implementing systems that can register objective physiological variables reflecting both the sympathetic and parasympathetic activity in daily situations and conditions. This would enable the long-term monitoring of subjective autonomous nervous system responses, as reflected by specific sympathetic and parasympathetic indicators that can be derived from galvanic skin responses and variability of cardiac rhythm.
Nevertheless, while APA methods have the potential to offer more objective insights into pain intensity, it is crucial to acknowledge their limitations. These mostly include the lack of high-quality validation studies, uncertainty about which parameters should be implemented in different settings, and technical issues such as timing of utilization [18]. Consequently, a comprehensive pain assessment approach should ideally combine both subjective self-reporting and objective measures to provide a more holistic understanding of pain experiences and enable better-informed pain management strategies.
Given these premises, the objective of this study was to explore potential correlations between subjective and objectively derived variables using APA modalities in cancer patients, as an initial step towards establishing a comprehensive pain assessment pathway.
Methods
Study design and population
The investigation was granted approval by the local Ethics Committee of Istituto Nazionale Tumori, Fondazione Pascale, Naples, Italy (protocol code 41/20 Oss; approval date: 26 November 2020), and all participants provided written informed consent. All methods were performed following the ethical standards as laid down in the Declaration of Helsinki and its later amendments or comparable ethical standards and it was registered with ClinicalTrials.gov, number (NCT04726228). The study employed a prospective observational design, utilizing medical records as its primary data source. The cohort for this study comprised adult patients undergoing treatment for cancer-related pain at the Istituto Nazionale Tumori, Fondazione Pascale, Italy. Clinical-instrumental assessments and collection of anamnestic data were conducted during a single outpatient visit at the pain therapy clinic. The process of recording biosignals was executed spanning a timeframe of approximately 5 min.
Variables considered
The analysis involved the examination of numerous variables, including demographic factors such as age, and body mass index (BMI). Clinical measurements included the number of comorbidities (categorized as “None” and at least one), Eastern Cooperative Oncology Group (ECOG) performance status (categorized as low ≤ 2, high > 2), metastases (yes/no), bone metastases (yes/no). Other clinical data concerned the type of anticancer therapy (chemotherapy, immunotherapy, surgery, and radiotherapy). Pain-related variables were pain type (nociceptive, neuropathic) and intensity (0–10 NRS), breakthrough cancer pain (BTcP) [19], and its features (type: nociceptive or neuropathic; 0–10 NRS intensity). Pain therapy included opioids, calculated as morphine equivalent dose (MED) for background pain, rapid onset opioids for BTcP, and pain adjuvants (corticosteroids, antidepressants, anticonvulsants, muscle relaxants, and anti-anxiety medications).
Biosignals
The objective variables were obtained from the analysis of the EDA and ECG biosignals. EDA reflects the sympathetic nervous system activity and serves as a valuable indicator for evaluating pain-induced neurocognitive stress. It detects alterations in the skin's electrical properties, triggered by sweat gland activation, ultimately resulting in increased skin conductance [20]. The continuous variations in skin conductance are termed the Skin Conductance Level (SCL), while the rapid responses occurring within seconds are identified as the Galvanic Skin Response (GSR). Both SCL and GSR play roles in the tonic and phasic components (Fig. 1).
For acquiring and analyzing biosignals, we employed a previously validated approach [21]. Regarding EDA biosignals, two EDA-derived analyses have been performed including the continuous decomposition analysis (CDA), and the trough-to-peak (TTP) analysis, both conducted in accordance to [22,23,24]. For the analysis, the mean amplitude of all above-threshold GSRs was calculated for both CDA and TTP analyses and adopted as a representative parameter of the EDA signal since it provides an overall synthetic measure of the average amplitude of registered skin conductance responses.
Concerning the ECG signal, temporal changes in inter-beat intervals offer a gauge of heart rate variability (HRV), closely associated with the autonomic nervous system (ANS) activity. These fluctuations can also signal physiological reactions to stressful or painful situations. Following our previous study [21], we focused on the R-R series of interbeat intervals, defined as the time between successive R waves of the QRS complex on the ECG waveform, [25], and then used it to extract time-domain parameters of the HRV. In particular, as reported in [21], R peaks were detected using a modified version of the Pan-Tompkins technique based on the ECG envelope and subsequent flattening of the signal to enhance the QRS-complexes as proposed in [26,27,28], and the raw series of differences between consecutive peaks was obtained (R-R). Then, using a recursive filtering process, the ECG-derived R-R time series were filtered to remove outliers and intervals that deviated most from the mean of the nearby RR intervals. Subsequently, from the ECG-derived R-R time series (Fig. 2), the heart rate (HR) time series has been obtained, and the standard deviation (SD) of both R-R and HR series has been calculated and used as a reference time-domain indicator of the variability of the heart rhythm.
As we described in a previous study [21], EDA and ECG signals have been acquired by using a BITalino device equipped with sensors for the recording of ECG and EDA signals. The BITalino platform is a hard-ware-affordable and open-source biosignals platform developed for physiological computing. Signals have been collected at a 1000 Hz sampling rate.
Statistical analyses
The sensitivity analysis tests were conducted using the Li and Yu [29] and Raikov test [30] on the continuous variables for assessing the potential Missing Not At Random (MNAR) process presence. Therefore, based on the information obtainable from the data and the context, the best conclusion on the mechanism generating the missing data will be reached. A Multiple Imputation analysis was finally performed on data to preserve the sample size.
A Multiple Factorial Analysis (MFA) was performed on data to detect multi-correlations and associations between numerical and categorical variables within the context of subjective and objective pain features. It is a statistical technique used for analyzing databases with large sets of variables that are suspected of being statistically related and describe the same set of observations or individuals [31]. This strategy is often implemented for addressing different types of variables (e.g., numerical, categorical) and for analyzing their relationships in a single analysis by creating principal components (PCs) which are assigned as a part of data behavior and variability. MFA allows analysts to explore patterns, similarities, and differences between the different sets of variables and observations [32]. In our analysis, a focus was applied to the main interest variables namely EDA and ECG derivates and BTcP presence (groups: No BTcP, BTcP), by reading 2D-projected modalities in the factorial plane respect with the BTcP intensity groups.
Within the realm of MFA, when dealing with complex multidimensional datasets featuring numerous categorical variables, the utilization of Multiple Correspondence Analysis (MCA) proves to be a valuable strategy. It enables the condensation of dataset dimensions into a selected set of categorical variables, effectively capturing the most essential information. Therefore, we adopted a data-driven approach known as Hierarchical Clustering on Principal Components (HCPC). This approach harmoniously integrates three established methods including MCA, or PC analysis (PCA) for numerical variables, hierarchical clustering, and the k-means algorithm, to derive a refined and enhanced cluster solution [33]. Through this amalgamation of techniques, we aimed to extract a more robust and accurate representation of patterns and relationships present within the data. Operatively, the HCPC method provides a PCA and, subsequently, a hierarchical clustering by implementing an agglomerative hierarchical tree; the ideal number of groups is obtained by a pruning approach; then k-means was performed on data by setting the number of groups, obtained as described above. Thus, different variables were inserted in the routine to find the BTcP-related risk groups. They included BTcP type, ECOG, metastases (No, Yes, bone metastases), chemotherapy, radiotherapy, and surgery, as well as drug therapy including MED (≤ 60 mg; > 60 mg) and adjuvants. Finally, a multivariable linear regression analysis was conducted to assess the main associations between MED (treated as a number and binary variable) and potentially related variables.
The data were analyzed using the R software version 4.2.3 (R Core Teams, R Foundation for Statistical Computing, Vienna, Austria). The toolkit included car, purr, boot, snow, misty, and naniar. The Mice package was adopted for the imputation of the missing data. FactoMineR was used for the implementation of factor analysis methods. Moreover, graphical packages were adopted for the visualization of the plots. The graphics packages included ggplot2 and factoextra. Biosignal processing and analysis were conducted using MatLab v. R2021b from The Math-Works Inc.
Results
Among 153 patients deemed eligible for the study, biosignals were recorded from 118 individuals. Data from 54 patients were excluded from the analysis due to issues with signal quality or incomplete information. Ultimately, the exploratory statistical analysis incorporated clinical, demographic, and biosignal variables from 64 patients (Fig. 3).
Demographic data, clinical variables, as well as subjective and objective data are reported in Table 1.
Multiple factorial analysis for BTcP analysis
The MFA focused on the main two components or dimensions of variables correlation, expressed as Factorial Axes (FA1 and FA2). The total explained variance was 30.3% (17.7% on FA1, 12.6% on FA2) for the complete set of variables.
The first dimension (FA1) was mainly described by BTcP presence (23.6%) and Type of BTcP (22.8%), followed by ECOG, Opioids, and Metastases (13.1%, 11.2%, and 10.3%, respectively) (Fig. 4a). The second dimension (FA2) was strongly characterized by the objective variables SDHR and SDRR (42.9% and 42.0%, respectively), followed by metastases (7.1%) (Fig. 4b).
Moreover, poor correlations were found for Metastases with FA1 (cos2 = 0.05); concerning FA2, strong indicators were found for both SDHR and SDRR (both contributed > 42% and > 0.8 cos2’s, and almost null values were calculated for other variables); mean TTP and CDA were mainly collected into the third dimension (absolute contributes: 26.2% and 24.5%, and similar cos2’s).
The factorial plane highlighted the characteristics of the BTcP groups. Patients without BTcP tended to be treated with lesser opioid dosages (MED) for background pain (67% vs 93%, p = 0.03) and, although not significantly, they were less prone to undergo chemotherapy (33% vs 57%, p = 0.1), and radiotherapy (22% vs 32%, p = 0.54) compared to BTcP affected cancer patients. Patients with BTcP were more likely to have metastases (64% vs. 22%, p = 0.18, not shown in table), especially bone metastases (46% vs. 22%, p = 0.12) compared with no BTcP group. Moreover, the BTcP group was characterized by a higher ECOG (57% vs. 11%, p < 0.01). MFA detected no potential correlations between subjective variables (BTcP presence) and objective ones (Fig. 5).
Multiple Factorial Analysis. Factorial planes for modalities and individuals. Legend: The first two dimensions were considered. Strong components were given by the type of BTcP, the intensity of BTcP, ECOG, Opioids, and metastases. Very poor contributions were finally detected for objective pain variables (red spots and labels): TTP, CDA SDHR, SDNN. The projection from the center to the variable's point revealed two key aspects including an orthogonal alignment with the existing variables' cloud pattern and a predominantly centralized contribution to FA2. Therefore, null correlations were confirmed between subjective variables (BTcP Intensity as NRS measure) and objective ones. The red arrow (zoomed mean of coordinates from the objective variables) points out the variable cloud shape, stating that poor correlations were detected. Each black point represents a variable modality which was plotted as “name of the variable”-underscore-modality: “N” = “No”, “Y” = “Yes”, “neuro” = “neuropathic”, “nocic” = “nociceptive”. Abbreviations: BMI, body mass index; ECOG, Eastern Cooperative Oncology Group performance status; BTcP, Breakthrough cancer pain; CDA, Continuous Decomposition Analysis; NRS, Numeric Rating Scale; TTP, Trough-To-Peak; SDRR, standard deviation R-R (interbeat intervals); SDHR, standard deviation heart rate; MED, Morphine Equivalent Dose
The results of the univariable group analysis for BTcP are presented in Table 2.
Clustering on principal components for BTcP
The optimal partition based on principal components yielded 3 distinct clusters. MCA reported a global variability explanation of 35.2% on 10 categorical variables and 22 modalities. None of the cancer patients belonging to Cluster 1 (27 cancer patients) had MED (p < 0.01), almost 93% of them had a low ECOG (p < 0.01), more than 88% did not have chemotherapy (p < 0.01), 96.3% radiotherapy (p < 0.01), did not present metastases (74.1%, p < 0.01). This cluster was mainly characterized by no pain nor nociceptive pain (0% and 18.5%, respectively, p < 0.01); 59.4% of cancer patients with no pain belong to cluster 1. Cluster 2 was composed of 14 cancer patients, all of whom had chemotherapy (p < 0.01) and was characterized by a low ECOG (p < 0.01). An amount of 57% of them suffered from neuropathic pain (p = 0.04), only one had nociceptive pain, and another one did not assume opioids. Cluster 3 (23 individuals) was characterized by opioids assumption (100%, p < 0.01), high ECOG (90% of them belonged to cluster 3, p < 0.01, 78.3% of cluster 3), had MED > 60 (65%; corresponding to the 75% of the whole set of patients with MED > 60, p < 0.01); moreover, cancer patients of cluster 3 were affected by bone metastases BTcP (81% of all bone metastases belonged to cluster 3, of them the 73.9% was bone metastases, p < 0.01) (Fig. 6).
Cluster analysis for opioid therapy. Legend: The study identified three distinct clusters through principal component analysis. Multiple Correspondence Analysis (MCA) explained 35.2% of global variability across 10 categorical variables and 22 modalities. Cluster 1 (27 patients) had no medication, mostly low Eastern Cooperative Oncology Group (ECOG) performance status scores, no chemotherapy, and no metastases, with no or minimal pain reported. Cluster 2 (14 patients) consisted of patients all on chemotherapy, with low ECOG scores, and some experiencing neuropathic pain. Cluster 3 (23 patients) showed 100% opioid usage, high ECOG scores, MED > 60, and a high prevalence of bone metastases breakthrough cancer pain (BTcP)
The principal component analysis identified three distinct clusters. Multiple Correspondence Analysis (MCA) explained 35.2% of global variability across 10 categorical variables and 22 modalities. Cluster 1 (27 patients) had no MED, predominantly low ECOG, no chemotherapy, and no metastases. It was characterized by the absence of pain, with 59.4% of pain-free patients belonging to this cluster. Cluster 2 (14 patients) had all patients on chemotherapy, low ECOG, and some experiencing neuropathic pain. Cluster 3 (23 patients) had 100% opioid usage, high ECOG, MED > 60, and a high prevalence of bone metastases BTcP.
Results from the k-mean were reported in Table 3. Comparing clusters with BTcP measures 87% of cancer patients, in cluster 1 were not affected by BTcP while cluster 3 was characterized by higher NRS pain (85% vs. 17%, p < 0.01). Mean CDA, mean TTP, SDNN (defined as the normal R-R intervals standard deviation, after outliers removal from R-R series), and SDHR were not found as significantly different among clusters.
Multivariable regression analyses for opioid therapy
In multivariable analysis, age and BMI were not associated with opioid therapy (see Model 1, Table 4). Similarly, no significant differences based on SDNN, SDHR, and adjuvant therapy were detected; the opioid amount was significantly and positively associated with ECOG as cancer patients with higher ECOG meanly received 72.4 MED more (95% CI: [12.1,133.0], p = 0.02). A positive correlation was found with the type of BTcP (p < 0.01): neuropathic pain likely received more opioids (147 units of MED) compared with No pain. Mean TTP showed a non-significant trend with the MED therapy (p = 0.09). In the logistic analysis (Model 2, Table 4) a significant association with higher MED assumption was given by the Intensity of BTcP, as expected: in particular, cancer patients with high pain (> 8 NRS) were 8 times more likely to receive higher opioid dosages (MED > 60) (OR = 8.3, 95%CI = [1.8, 48.4], p = 0.02). A borderline association was detected for SDHR (OR = 1.5, 95%CI = [0.98, 2.39], p = 0.06).
Discussion
To effectively manage pain, a thorough assessment of the various components that characterize different pain manifestations is essential [34]. This is particularly crucial for oncological pain, which exhibits complex and distinctive features [35]. For instance, BTcP is a sudden, intense flare-up of pain that occurs in cancer patients who are already receiving pain medication for their chronic pain. It is a pain phenomenon with a complex pathophysiology that has yet to be fully characterized [36]. In this convoluted situation, integrating objective data in a multidimensional model and analyzing clinical and therapeutic elements can yield significant insights for tailoring the therapy [17].
While research on APA processes is a rapidly advancing field [17, 18, 37], to our knowledge, this study represents one of the first attempts to integrate data derived from biosignal analysis and subjective, patient-reported, elements. For this aim, we followed a multidimensional exploratory statistical pathway. In the MFA, variability is condensed into a few components (dimensions) by creating FAs. This approach is employed to clarify the general pattern of data when there is a large number of statistical relations among a set of numerous variables [32]. The obtained latent variables (i.e., FA1 and FA2) are linear combinations of the data variables. They collect the sample's variability in a descending manner with the first being the most significant, followed by the second, and so on. These latent variables are uncorrelated with each other and reflect distinct aspects of the whole set of variables. The analysis showed that the presence or not of BTcP, and BTcP type are mainly projected onto the first dimension (FA1). Notably, they represent “subjective" variables. This result suggests that clinical and subjective factors, including patients’ overall health, strongly describe BTcP features. On the other hand, statistical measures (contributions and squared cosines) demonstrated that the objective APA-extracted variables are not correlated with the other clinical factors. Therefore, although as expected, patients with BTcP had worse clinical conditions, for example in terms of ECOG and metastases [38], the analysis failed to find a correlation between subjective and objective pain variables to explain the BTcP phenomenon. These findings were confirmed by the clustering analysis as well as in the multivariate analysis conducted for therapeutic profile characterization (MED). Remarkably, the results indicated that performance status and the type of pain are correlated with a higher demand for opioid analgesics, objective variables alone did not provide a clear or comprehensive explanation for the overall opioid dosage.
The lack of statistical correlation between APA variables and pain features and management can be attributed to the inherent complexity of pain, the individual nature of pain perception, and the limitations of measurement methods. Objective pain measures are based on physiological or clinical indicators, such as HR, blood pressure, or the presence of specific physical symptoms. Subjective pain assessments, on the other hand, rely on patients' self-reported experiences, including pain intensity, quality, and emotional aspects. These two approaches may capture different aspects of pain, making them less likely to correlate. Besides the barrier between objective and subjective measures, the lack of significance of quantitative physiological parameters based on EDA and ECG should be investigated from a broader and more complex perspective as it can be affected by other measuring and technical factors, such as the duration of the recordings and the type of features considered, which is limited to the most widespread and simple indicators GSRs and HR variations and does not still consider frequency-domain parameters that could better enlighten the role of sympathovagal balance as well as more complex and sophisticated measures of nonlinear and chaotic dynamics (e.g., entropy indices, fractal analysis) that may be involved in physiological control systems in pathological states. As mentioned in the Introduction section, by elucidating the role of biosignals in cancer pain assessment, the potential of such sources of information to be used in continuously monitoring settings could then be fully exploited, such as through the integration of pain indicators within wearable platforms for long-term patients’ monitoring. Moreover, pain is a complex and highly individual experience [39]. The variability in pain perception can make it challenging to establish a direct relationship between objective physiological changes and subjective pain reports. Furthermore, pain perception is influenced by psychological factors, such as anxiety, depression, and coping strategies [40]. These aspects can significantly impact a person's subjective pain experience but may not be reflected in objective measures. Another key element is the clinical setting adopted for APA recording. Objective measures may be more reliable in acute pain situations, such as post-surgery, where there is a clear physiological response. In chronic pain conditions, the relationship between objective and subjective measures may be less straightforward due to long-term adaptations and psychosocial factors [41]. Consequently, serious challenges arise for the analysis of acute cancer pain phenomena such as BTcP. More importantly, for an effective physiological characterization, objective measures should capture the phenomenon in real time. It would require almost continuous recording of physiological data. Likely, the availability of detection systems with wearable platforms and more complex feature extraction approaches for biosignal analysis may provide us with more reliable data. In this initial investigation, the purpose was to understand the potential for integrating data to subsequently enhance both data collection and analysis strategies.
Study limitations
While the study aims to offer insights into the multimodal assessment of cancer pain, its limitations underscore the need for cautious interpretation of its findings. For example, the study’s sample size might not be representative of the broader population of cancer patients. Moreover, the investigation was conducted at a single medical institution, which could limit the generalizability of the findings to other healthcare settings and diverse populations. Biosignal measurement represents the main limitation of the study. For example, the timing and methods of recording are crucial aspects to be defined to calibrate a model. The accuracy and reliability of these signals in assessing pain can be influenced by factors such as noise, signal artifacts, and individual variability in physiological responses. Furthermore, additional biosignals-based features derived from time- and frequency-domain approaches as well as time–frequency analysis and newly proposed deep learning-based features for both ECG and EDA have not been considered in this study but will be adopted in future works to widen the range of quantitative physiological parameters that could provide new descriptors and predictors of pain, thus yielding higher sensitivity in quantifying cancer pain. Finally, the study lacked extensive validation of its biosignal-based pain assessment methods against established pain measurement techniques. Notably, the absence of high-quality validation studies could limit the reliability of these methods.
Conclusions
Recognizing the significance of thorough pain assessment is a pivotal step in enhancing pain management strategies. Although subjective pain assessment modalities are widely used to understand and manage pain, they are not without limitations. Striking a balance between patient reports and objectivity in pain assessment can offer incredible opportunities to provide optimal care. Nevertheless, achieving a cohesive integration between the objective physiological parameters assessed through APA and the subjective insights gathered from self-reporting is a complex endeavor. It is crucial to recognize the significance of each considered feature and thoughtfully harmonize them to derive comprehensive conclusions that truly capture the essence of pain phenomena. Future research with larger and more diverse samples, rigorous validation of biosignal-based methods, and controlled experimental designs could provide a more comprehensive understanding of pain assessment and management strategies.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from Marco Cascella on reasonable request.
Abbreviations
- ANS:
-
Autonomic nervous system
- APA:
-
Automatic pain assessment
- BMI:
-
Body mass index
- BTcP:
-
Breakthrough cancer pain
- CDA:
-
Continuous decomposition analysis
- ECG:
-
Electrocardiogram
- ECOG:
-
Eastern Cooperative Oncology Group
- EDA:
-
Electrodermal activity
- EEG:
-
Electroencephalography
- FA:
-
Factorial Axis
- fMRI:
-
Functional magnetic resonance imaging
- GSR:
-
Galvanic Skin Response
- HCPC:
-
Hierarchical clustering on principal components
- HR:
-
Heart rate
- HRV:
-
Heart rate variability
- MCA:
-
Multiple Correspondence Analysis
- MED:
-
Morphine Equivalent Dose
- MFA:
-
Multiple factorial analysis
- MNAR:
-
Missing Not At Random
- NRS:
-
Numeric rating scale
- PCA:
-
PC analysis
- SCL:
-
Skin Conductance Level
- SDHR:
-
Heart rate standard deviation
- SDNN:
-
Normal RR intervals standard deviation
- SDRR:
-
RR intervals standard deviation
- SDs:
-
Standard deviations
- VAS:
-
Visual Analogic Scale
- WHO:
-
World Health Organization
References
WHO Guidelines for the pharmacological and radiotherapeutic management of cancer pain in adults and adolescents. https://www.who.int/publications-detail-redirect/9789241550390. Accessed 2 Sept 2023.
Fallon M, Giusti R, Aielli F, Hoskin P, Rolke R, Sharma M, et al. Management of cancer pain in adult patients: ESMO Clinical Practice Guidelines. Ann Oncol. 2018;29(Suppl 4):iv166-91.
Brant JM. The assessment and management of acute and chronic cancer pain syndromes. Semin Oncol Nurs. 2022;38:151248.
Fink RM, Gallagher E. Cancer pain assessment and measurement. Semin Oncol Nurs. 2019;35:229–34.
Cascella M, Vittori A, Petrucci E, Marinangeli F, Giarratano A, Cacciagrano C, et al. Strengths and Weaknesses of cancer pain management in Italy: findings from a nationwide SIAARTI Survey. Healthcare (Basel). 2022;10:441.
Mustajoki M, Forsén T, Kauppila T. Pain assessment in native and non-native language: difficulties in reporting the affective dimensions of pain. Scand J Pain. 2018;18:575–80.
Cascella M, Bimonte S, Saettini F, Muzio MR. The challenge of pain assessment in children with cognitive disabilities: Features and clinical applicability of different observational tools. J Paediatr Child Health. 2019;55:129–35.
Marchetti G, Vittori A, Cascella M, Mascilini I, Piga S, Petrucci E, et al. Pain prevalence and pain management in children and adolescents in an italian third level pediatric hospital: a cross-sectional study. Ital J Pediatr. 2023;49:41.
Fry M, Arendts G, Chenoweth L. Emergency nurses’ evaluation of observational pain assessment tools for older people with cognitive impairment. J Clin Nurs. 2017;26:1281–90.
Vittori A, Cascella M, Petrucci E, Cortegiani A, Bignami EG, Innamorato MA, et al. Strategies to build and maintain competence in pain management: Insights from a SIAARTI survey on educational needs among Italian anesthesiologists. Pain Pract. 2023. https://doi.org/10.1111/papr.13207.
Can G, Mushani T, Rajhi BHA, Brant JM. The global burden of cancer pain. Semin Oncol Nurs. 2019;35:315–21.
Caraceni A, Shkodra M. Cancer pain assessment and classification. Cancers (Basel). 2019;11:510.
Hammal Z, Cohn JF. Automatic detection of pain intensity. Proc ACM Int Conf Multimodal Interact. 2012;2012:47–52.
Cascella M, Vitale VN, Mariani F, Iuorio M, Cutugno F. Development of a binary classifier model from extended facial codes toward video-based pain recognition in cancer patients. Scand J Pain. 2023;23:638–45.
Vandenbussche N, Van Hee C, Hoste V, Paemeleire K. Using natural language processing to automatically classify written self-reported narratives by patients with migraine or cluster headache. J Headache Pain. 2022;23:129.
Szczapa B, Daoudi M, Berretti S, Pala P, Del Bimbo A, Hammal Z. Automatic estimation of self-reported pain by interpretable representations of motion dynamics. Proc IAPR Int Conf Pattern Recogn. 2021;2020.
Cascella M, Schiavo D, Cuomo A, Ottaiano A, Perri F, Patrone R, et al. Artificial intelligence for automatic pain assessment: research methods and perspectives. Pain Res Manag. 2023;2023:6018736.
Gkikas S, Tsiknakis M. Automatic assessment of pain based on deep learning methods: A systematic review. Comput Methods Programs Biomed. 2023;231:107365.
Davies AN. Breakthrough cancer pain. Curr Pain Headache Rep. 2014;18:420.
Bari DS, Aldosky HYY, Tronstad C, Kalvøy H, Martinsen ØG. Electrodermal activity responses for quantitative assessment of felt pain. J Electr Bioimpedance. 2018;9:52–8.
Cascella M, Vitale VN, D’Antò M, Cuomo A, Amato F, Romano M, et al. Exploring biosignals for quantitative pain assessment in cancer patients: a proof of concept. Electronics. 2023;12:3716.
Daviaux Y, Bonhomme E, Ivers H, de Sevin É, Micoulaud-Franchi J-A, Bioulac S, et al. Event-related electrodermal response to stress: results from a realistic driving simulator scenario. Hum Factors. 2020;62:138–51.
Benedek M, Kaernbach C. omposition of skin conductance data by means of nonnegative deconvolution. Psychophysiology. 2010;47:647–58.
Benedek M, Kaernbach C. A continuous measure of phasic electrodermal activity. J Neurosci Methods. 2010;190:80–91.
Pan J, Tompkins WJ. A real-time QRS detection algorithm. IEEE Trans Biomed Eng. 1985;32:230–6.
Moeyersons J, Amoni M, Van Huffel S, Willems R, Varon C. R-DECO: An open-source Matlab based graphical user interface for the detection and correction of R-peaks.
Moeyersons J, Amoni M, Huffel SV, Willems R, Varon C. R-DECO: an open-source Matlab based graphical user interface for the detection and correction of R-peaks. PeerJ Comput Sci. 2019;5:e226.
Goldberger AL, Amaral LA, Glass L, Hausdorff JM, Ivanov PC, Mark RG, et al. PhysioBank, PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals. Circulation. 2000;101:E215-220.
Li J, Yu Y. A nonparametric test of missing completely at random for incomplete multivariate data. Psychometrika. 2015;80:707–26.
Raykov T, Schneider BC, Marcoulides GA, Lichtenberg PA. Examining measure correlations with incomplete data sets. Struct Equ Modeling. 2014;21:318–24.
EUDML | Analyse factorielle multiple appliquée aux variables qualitatives et aux données mixtes. https://eudml.org/doc/106525. Accessed 13 Feb 2024.
Sundin P, Crespi CM. Power analysis for stepped wedge trials with multiple interventions. Stat Med. 2022;41:1498–512.
Lever J, Krzywinski M, Altman N. Principal component analysis. Nat Methods. 2017;14:641–2.
Zhang H. Cancer pain management-new therapies. Curr Oncol Rep. 2022;24:223–6.
Cuomo A, Cascella M, Forte CA, Bimonte S, Esposito G, De Santis S, et al. Careful breakthrough cancer pain treatment through rapid-onset transmucosal fentanyl improves the quality of life in cancer patients: results from the BEST multicenter study. J Clin Med. 2020;9:1003.
Mercadante S. Treating breakthrough pain in oncology. Expert Rev Anticancer Ther. 2018;18:445–9.
Nagireddi JN, Vyas AK, Sanapati MR, Soin A, Manchikanti L. The analysis of pain research through the lens of artificial intelligence and machine learning. Pain Physician. 2022;25:E211–43.
Camps Herrero C, Batista N, Díaz Fernández N, Escobar Álvarez Y, Gonzalo Gómez A, Isla Casado D, et al. Breakthrough cancer pain: review and calls to action to improve its management. Clin Transl Oncol. 2020;22:1216–26.
Treede R-D, Rief W, Barke A, Aziz Q, Bennett MI, Benoliel R, et al. Chronic pain as a symptom or a disease: the IASP classification of chronic pain for the International Classification of Diseases (ICD-11). Pain. 2019;160:19–27.
Fillingim RB. Individual differences in pain: understanding the mosaic that makes pain personal. Pain. 2017;158(Suppl 1 Suppl 1):S11-8.
Moscato S, Orlandi S, Giannelli A, Ostan R, Chiari L. Automatic pain assessment on cancer patients using physiological signals recorded in real-world contexts. Annu Int Conf IEEE Eng Med Biol Soc. 2022;2022:1931–4.
Acknowledgements
The authors would like to thank Giacomo Sollecchia for his cooperation.
Funding
This work was supported also by the University of L’Aquila.
Author information
Authors and Affiliations
Contributions
Conceptualization, M.C., P.D.G., A.V. and S.C.; methodology, M.C., P.D.G., A.C., A.V., E.P., S.S., F.M., A.M.P., M.R., C.O., A.O., F. S., F.P., O.P. and S.C.; software, M.C., P.D.G., A.C., A.V., A.M.P., M.R., C.O., O.P. and S.C; validation, M.C., P.D.G., A.C., AV., A.M.P., M.R., C.O., O.P. and S.C.; formal analysis, M.C., P.D.G., A.C., AV., A.M.P., M.R., C.O., O.P. and S.C.; investigation, M.C., P.D.G., A.C., AV., A.M.P., M.R., C.O. and S.C.; resources, M.C., A.V. and F.M.; data curation, M.C., P.D.G., A.C., AV., A.M.P., M.R., C.O., O.P. and S.C.; writing-original draft preparation, M.C., P.D.G., A.C., A.V. and S.C.; writing-review and editing, M.C., P.D.G., A.C., A.V., E.P., S.S., F.M., A.M.P., M.R., C.O., A.O., F. S., F.P.,F.S., O.P. and S.C.; visualization, M.C., A.V., O.P. and S.C.; supervision, M.C., O.P. and S.C.; project administration, M.C., A.V., F.M., and O.P. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The investigation was granted approval by the local Ethics Committee of Istituto Nazionale Tumori, Fondazione Pascale, Naples, Italy (protocol code 41/20 Oss; approval date: 26 November 2020), and all participants provided written informed consent. All methods were performed in accordance with the ethical standards as laid down in the Declaration of Helsinki and its later amendments or comparable ethical standards.
Consent for publication
Not competent.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Cascella, M., Di Gennaro, P., Crispo, A. et al. Advancing the integration of biosignal-based automated pain assessment methods into a comprehensive model for addressing cancer pain. BMC Palliat Care 23, 198 (2024). https://doi.org/10.1186/s12904-024-01526-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12904-024-01526-z