Pain Intensity Recognition Rates via Biopotential Feature Patterns with Support Vector Machines

Sascha Gruss; Roi Treister; Philipp Werner; Harald C. Traue; Stephen Crawcour; Adriano Andrade; Steffen Walter

doi:10.1371/journal.pone.0140330

Abstract

Background

The clinically used methods of pain diagnosis do not allow for objective and robust measurement, and physicians must rely on the patient’s report on the pain sensation. Verbal scales, visual analog scales (VAS) or numeric rating scales (NRS) count among the most common tools, which are restricted to patients with normal mental abilities. There also exist instruments for pain assessment in people with verbal and / or cognitive impairments and instruments for pain assessment in people who are sedated and automated ventilated. However, all these diagnostic methods either have limited reliability and validity or are very time-consuming. In contrast, biopotentials can be automatically analyzed with machine learning algorithms to provide a surrogate measure of pain intensity.

Methods

In this context, we created a database of biopotentials to advance an automated pain recognition system, determine its theoretical testing quality, and optimize its performance. Eighty-five participants were subjected to painful heat stimuli (baseline, pain threshold, two intermediate thresholds, and pain tolerance threshold) under controlled conditions and the signals of electromyography, skin conductance level, and electrocardiography were collected. A total of 159 features were extracted from the mathematical groupings of amplitude, frequency, stationarity, entropy, linearity, variability, and similarity.

Results

We achieved classification rates of 90.94% for baseline vs. pain tolerance threshold and 79.29% for baseline vs. pain threshold. The most selected pain features stemmed from the amplitude and similarity group and were derived from facial electromyography.

Conclusion

The machine learning measurement of pain in patients could provide valuable information for a clinical team and thus support the treatment assessment.

Citation: Gruss S, Treister R, Werner P, Traue HC, Crawcour S, Andrade A, et al. (2015) Pain Intensity Recognition Rates via Biopotential Feature Patterns with Support Vector Machines. PLoS ONE 10(10): e0140330. https://doi.org/10.1371/journal.pone.0140330

Editor: Daoqiang Zhang, Nanjing University of Aeronautic and Astronautics, CHINA

Received: May 5, 2015; Accepted: September 24, 2015; Published: October 16, 2015

Copyright: © 2015 Gruss et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: Data are available at Dryad - DOI: doi:10.5061/dryad.2b09s.

Funding: Funding Source: This research was part of the DFG/TR233/12 (http://www.dfg.de/) “Advancement and Systematic Validation of an Automated Pain Recognition System on the Basis of Facial Expression and Psychobiological Parameters” project, funded by the German Research Foundation, Fundação de Amparo à Pesquisa do Estado de Minas Gerais, Conselho Nacional de Desenvolvimento Científico e Tecnológico, Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior, and the Brazilian government funded travel, research stay and conference costs.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Quantifying pain is possible with the aid of the Visual Analog Scale or Numeric Rating Scale. However, these methods only work when the patient is sufficiently alert and cooperative, i.e., under conditions not always given in the medical field (e.g., post-surgery phases). There also exist instruments for pain assessment in people with verbal and / or cognitive impairments and instruments for pain assessment in people who are sedated and automated ventilated [1]. Overall, these methods are still in development or in need of validation. If conditions do not allow for a sufficiently valid measurement of pain, then cardiac stress in at-risk patients, under-perfusion of the operating field, or the development of chronic pain may follow in consequence. On the other hand, opiates may alleviate pain sensation, but can also lead to severe addiction, constipation etc. [2]. Hence, the measurement of biopotentials via the autonomic nervous system may be a solution that would permit an objective, reliable, and variable surrogate measurement of pain.

Some studies examined the correlation between a single biopotential and pain [3], [4], [5], [6], [7], [8], [9], [10]. In conventional pain diagnostics it is known that a measurement on a single parameter feature is insufficient for a valid diagnosis. Instead, a combination of multi-parameter features is required [1]. To the best of our knowledge, the study by Treister et al. [11] was the first to take a multi-parameter biopotential approach. Tonic heat was applied to elicit pain for a duration of 1 min, with intensities of “no pain,” “low pain,” “medium pain,” and “high pain.” While all of the features differed significantly in “no pain” and the other thresholds via Friedman Test, only the linear combination of parameters significantly differentiated between pain vs. no pain, as well as between all other pain categories. Bustan et al. also employed multimodal parameters to investigate the relationship between pain intensity, unpleasantness, and suffering [12]. High-intensity stimulation elicited higher skin conductance level compared with low-intensity stimulation under conditions of tonic noxious stimulation; heart rate was higher for short than for long stimulation, while corrugator electromyography showed of no significant effect regarding the response.

In our research we aim at the advancement of pain diagnosis and monitoring of pain states. For the purpose we developed an extensive multimodal dataset in which several levels of pain are induced. In a high density feature space a machine learning model (based on SVM) could be a solution. “An SVM Model could be trained on one set of individuals, and used to accurately classify pain in different individuals”, p. 2 [13].

In Walter et al. [14] preliminary results were presented. A total of 135 features were extracted from the mathematical groupings of amplitude, frequency, stationarity, entropy, linearity, and variability from the facial and trapezius electromyography, skin conductance level, and electrocardiography signals. The following features were statistically chosen as the most selective: 1. electromyography_corrugator_amplitude_peak_to_peak, 2. electromyography_corrugator_entropy_ shannon, and 3. heart_rate_variability_slope_RR. We received a classification rate (based on SVM) for the two class problem baseline vs. pain tolerance threshold of 77.05%. In Werner et al. [15] we received a classification rate for the two class problem baseline vs. pain tolerance threshold of 75.6% (without facial electromyography).

An extension of the feature space and the use of automated feature selection methods could improve classification rates, as compared to the study of Walter et al. [14].

The aim of the present study (see Fig 1d) was:

to select features (“pain pattern”) with a support vector machine learning design
to extract systematically a high-density feature space, for which we included signal similarity features (see Table 1, Equation: 37–42)—Similarity is mentioned in the signal processing literature as a powerful feature [16], [17]
to contribute to the highest recognition rate for pain quantification in a two and multi class problem.

Download:

Fig 1. Experimental procedure, (1a) Thermode on the right arm, (1b) Heat signal with baseline, (1c) Labor setting, (1d) Study procedure.

https://doi.org/10.1371/journal.pone.0140330.g001

Download:

Table 1. Feature information.

https://doi.org/10.1371/journal.pone.0140330.t001

Methods

Subjects

A total of 90 subjects participated in the experiment, recruited from the following age groups [18]: 1. 18–35 years (N = 30 years; 15 men, 15 women), 2. 36–50 years (N = 30; 15 men, 15 women), and 3. 51–65 years (N = 30; 15 men, 15 women). Only 85 subjects were included in the final analysis because four subjects had to be excluded due to limited data quality with regard to the EMG. Recruitment was performed through notices posted at the university for the 18- to 35-year-old age group and through the press for the 36-50- and 51-65-year-old age groups. Only healthy subjects were recruited. The subjects received an expense allowance. The study was conducted in accordance with the ethical guidelines set out in the WMA Declaration of Helsinki (ethical committee approval was granted: 196/10-UBB/bal). The study was approved according the ethics committee of the University of Ulm (Helmholtzstraße 20, 89081 Ulm, Germany). All participants provided a written informed consent to participate in this study. An official written document of the ethics committee approved this consent procedure.

Exclusion criteria

Prior to the experiment, the case history of each applicant was assessed in order to identify persons who met the exclusion criteria. Pre-existing neurological conditions, chronic pain, cardiovascular diseases, regular use of pain medication, and use of pain medication immediately before the experiment were applied as exclusion criteria.

Individual calibration of pain and tolerance threshold

To induce heat pain [18], a Medoc Pathway thermal stimulator was employed. The ATS thermode was attached to the right forearm of the subject (see Fig 1a). Before data recording commenced, we determined each participant’s individual pain threshold and pain tolerance, i.e., the temperatures at which the participant’s perception changed from heat to pain and the level at which the pain became unacceptable. We used these temperature thresholds for the lowest and highest pain levels and added two additional intermediate levels to obtain a ranged set of four equally distributed temperatures (see Fig 1b). The instruction given to subjects for determining the pain threshold was as follows: “Please press the stop button immediately when you experience a burning, stinging, piercing or pulling sensation in addition to the feeling of heat.” In order to determine the tolerance threshold, the following instruction was given: “Please press the stop button immediately when you can no longer tolerate the heat, taking into account the burning, stinging, piercing or pulling sensation.”

Experimental pain stimulation

In the experiment, each of the four temperature levels was applied 20 times in randomized order, resulting in a total of 80 stimuli. The baseline (no pain) was 32°C. For each stimulus, the temperature was maintained for 4 s. The pauses between the stimuli were randomized between 8–12 s (see Fig 1b). The subjects had the option to abort the experiment immediately by pressing an emergency stop button. After the experiment, we asked the subject to apply a cold pack to the spot of the heat stimulation for at least 5 minutes.

Measurement of biopotentials

Biopotentials.

A Nexus-32 amplifier (http://www.mindmedia.nl; accessed May 23, 2014) was used to record biopotential data (see Fig 1c) during the experiment. Biopotential and event data were recorded using Biotrace software. The following parameters were included in the classification process [18].

Electromyography (EMG).

Electrical muscle activity is also an indicator of general psychophysiological stimulation, in which increased muscle tone is associated with increasing activity of the sympathetic nervous system. A decrease in somatomotor activity reflects predominantly parasympathetic stimulation. We used two-channel EMGs for zygomaticus, corrugator, and trapezius muscles. In the area of affective computing, the activity of zygomaticus combined with happiness and the corrugator with negative affectivity [19]. The activity of the trapezius is an indication of a high stress level, which is also to be expected when pain is being experienced.

Skin conductance level (SCL).

To measure the skin conductance level, two electrodes connected to the sensor were positioned on the index and ring fingers. Because the sweat glands are innervated exclusively sympathetically (i.e., without the influence of the parasympathetic nervous system), electrodermal activity is considered a good indicator of the “inner tension" of a person. This phenomenon can be reproduced in particular by the observation of a rapid increase in skin conductance within 1–3 s due to a simple stress stimulus (e.g., deep breathing, emotional excitement, or mental activity).

Electrocardiogram (ECG).

We measured the average action potential of the heart on the skin using two electrodes, one on the upper right and one on the lower left of the body. Common features of the ECG signal are heart rate, interbeat interval, and heart rate variability (HRV). Heart rate variability refers to the oscillation of the interval between consecutive heartbeats and has been used as an indication of mental effort and stress in adults [20].

Preprocessing

We performed the following biopotential preprocessing:

We visualized all biopotentials to check the intensity of the noise and activity with regard to pain stimulation.
We applied a Butterworth filter to the EMG (20–250 Hz) and ECG (0.1–250 Hz) signals.
For the EMG, we also applied an additional filter using the Empirical Mode Decomposition technique developed by [21].
We quantified the pain level caused by the heat applied using four pain thresholds during the “pain window” (5.5 s) and with regard to the baseline during the “non-pain window” (see Fig 2).
We detected bursts of EMG activity using the Hilbert Spectrum [22].

Download:

Fig 2. Pain quantification.

https://doi.org/10.1371/journal.pone.0140330.g002

Feature extraction

We systematically extracted features [23], [24], [25], [26] from the mathematical groups of 1. amplitude (∑ = 40), 2. frequency (∑ = 24), 3. stationarity (∑ = 24), 4. entropy (∑ = 20), 5. linearity (∑ = 8), 6. variability (∑ = 19) and 7. similarity (Fig 3) (∑ = 24) (in total: ∑ = 159). Table 1 provides a detailed information overview of all features. The similarity features of a sample are calculated with regard to the associated mean baseline signal of the person. All features were normalized (z transformed) per person. The dataset of the study, including the raw and preprocessed signals, as well as the extracted features, is available at: https://www-e.uni-magdeburg.de/biovid/.

Download:

Fig 3. Similarity feature.

https://doi.org/10.1371/journal.pone.0140330.g003

Machine learning with Support Vector Machine (SVM)

Machine learning systems are systems that learn from known data and attempt to recognize characteristic patterns. After a 'learning phase' (also referred to 'training phase'), they return a model that can be used to map (i.e., classify) unknown input data into a category [27]. For these classification tasks, there are several machine learners (classifiers), all of which work using different decision algorithms, such as Neural Networks, Decision Trees, K-Nearest Neighbor, and SVM.

For the classification of different pain intensities, we considered using Neural Networks and k-Nearest Neighbor, but finally we chose SVMs (see Fig 4), as these have proven to be highly effective in other studies of affective computing [28] and are capable of maintaining sufficient flexibility with regard to their internal main parameter optimization [29].

Download:

Fig 4. Support Vector Machine hyperplane (H).

https://doi.org/10.1371/journal.pone.0140330.g004

The goal of an SVM is to develop a predictive model based on the given training samples (x_i, y_i) with x_i being a feature vector and y_i its associated class label. This model can subsequently be applied to an unlabeled test dataset to assign a particular class to each sample. With the aid of the feature vectors x_i, the SVM [30] searches for an optimal hyperplane with maximum margin in the feature space that separates the feature vectors of one class from feature vectors of the other. The hyperplane thus serves as the decision function. If the linear separation is not possible in the original feature space, all training vectors can be transformed to a higher dimensional space until the SVM finds a dividing hyperplane. This is done by means of a kernel function. In the present case we used a radial basis function kernel (RBF kernel), because it is able to handle non-linear dependencies between class labels and input attributes. Furthermore, the RBF-kernel has the advantage that the complexity of the model is limited to only two main parameters (C, γ). C controls the cost of misclassification of the training vectors [31], while γ controls the radius of influence of the support vectors [32]. In order to obtain optimal SVM parameters, a systematic grid search was performed for C and γ on 75% of data with a resulting 3-fold-cross-validation. When choosing the parameter values, we followed the procedure by Hsu et al. [29], who recommend exponentially growing sequences for C and γ: C = 2⁻⁵, 2⁻³,…, 2¹⁵ and γ = 2⁻¹⁵, 2⁻¹³,…, 2³. After testing of all combinations, the pair with the highest accuracy was finally selected as the optimal parameter set.

The performance of a classifier in general in a given learning task is measured by its classification rate (accuracy). Simply put, a set of feature vectors with known class labels is divided into two randomized mixed subsets. One subset is used for training the model and the other one is then applied for testing purposes. By putting the remaining feature vectors into the model, comparing their class labels with the classifier’s predictions, and finally counting the correctly predicted testing vectors, one receives the classification rate (accuracy) defined as (1)

We also calculated sensitivity and specificity, which give further information about the performance of the classification task. Both statistical measures are derived from the confusion matrix of the task. Put the case we have class “positive” and class “negative” with corresponding testing vectors for “positive” and “negative”. Then sensitivity is defined as (2) and specificity as (3)

Feature selection

Automatic pattern selection methods are used to further optimize recognition rates. Feature selection is a “method for selecting a subset of features providing optimal classification accuracy of the classification model” [33]. This is accomplished by means of a variety of feature selection (pattern configuration) methods, in combination with a classification procedure.

Before we conducted an automatic feature selection algorithm, we performed a manual pre-selection of the extracted 159 features based on statistical analyses and validation checks. As a first step, we eliminated features (feature groups) containing either a zero or a static number for all conditions. Such values may be attributed to a compromised signal or to a special feature extraction algorithm that appeared as uninterpretable on a particular bio signal. Second, we deleted all features which correlated positively or negatively with other features with at least 0.95 respectively—0.95. By doing so, we attempted to prevent classification of noise and redundant information.

Finally, we applied a feature selection algorithm combined with pain classification tasks to the remaining features. For this step, we chose the forward selection as it has shorter run-times (in most of the cases) as the backward elimination and brute force selection. The forward selection algorithm starts with an empty feature set and adds a new feature in each round. At each round, every new feature is tested for inclusion in the set by calculating a classification accuracy. The feature with the most increased accuracy is then added to the set before the beginning of the next round. The algorithm runs until there is no increase anymore [34].

We conducted the Support Vector Machine classification tasks in conjunction with the previously mentioned feature selection method on varying data sets and different biopotentials. SVMs were trained on 75% of the data (total number of training samples = 6375) tested on the remaining 25% of data (total number of test samples = 2125) (see Fig 5). The number of hereby-obtained optimal features ranges from 5 to 22 (depending on learning task).

Download:

Fig 5. Support Vector Machine learning architecture.

https://doi.org/10.1371/journal.pone.0140330.g005

Validity of classification rates

In many papers that address classification tasks like emotion or pain detection, machine learning classification results are presented without a measure of the validity—like in conventional statistic the p-level or effect size. However, we think this should be included in future analyses. In our case, we chose Cramér’s V [35] as an appropriate measure as it measures the strength of association between nominal variables. The value of Cramér’s V is ranging from 0 to +1 and is interpreted as follows: V < 0.1: negligible association, 0.1 < V < 0.2: weak association, 0.2 < V < 0.4: moderate association, 0.4 < V < 0.6: relatively strong association, 0.6 < V < 0.8: strong association, and 0.8 < V < = 1: very strong association [36].

Results

The classification results are summarized in Fig 6. Significance was tested against chance level. The classification results for a two class problem are between 79.29% - 90.94%, the Cramer’s V results and statistical measures sensitivity and specificity for B vs. T₁ are V = 0.59 (sensitivity = 76.00%, specificity = 82.59%), B vs. T₂ are V = .63 (sensitivity: 80.00%, specificity: 82.59%), B vs. T₃ are 84.94% V = .7 (sensitivity = 84.71%, specificity = 85.18%) and B vs. T₄ are 90.94% V = 0.82 (sensitivity = 92.24%, specificity = 89.65%).

Download:

Fig 6. Comparison between accuracy via support vector machine of study Walter et al. [14], without similarity signal feature [red] vs. support vector machine with similarity feature and automated selected feature [green].

https://doi.org/10.1371/journal.pone.0140330.g006

The Cramer’s V results and statistical measures sensitivity and specificity for the three and five class problem are between 43.05% - 73.09%, the classification results for B vs. T₁ vs. T₄ are V = .62 (B vs. others: sensitivity = 78.64%, specificity = 70.18%; T₁ vs. others: sensitivity = 69.96%, specificity = 74.69%; T₄ vs. others: sensitivity = 70.71%, specificity = 74.34%) and for B vs. T₁ vs. T₂ vs. T₃ vs. T₄ are V = .38 (B vs. others: sensitivity = 76.32%, specificity = 34.50%; T₁ vs. others: sensitivity = 34.89%, specificity = 45.25%; T₂ vs. others: sensitivity = 15.95%, specificity = 49.74%; T₃ vs. others: sensitivity = 28.18%, specificity = 46.52%; T₄ vs. others: sensitivity = 58.71%, specificity = 39.22%). Table 2 points out the most common features in a ranked order.

Download:

Table 2. Top ten importance ranking of selected features.

https://doi.org/10.1371/journal.pone.0140330.t002

Discussion and Conclusions

Our goal was to find significant classification results with high accuracies and an automated selected feature pattern of biopotentials that represents 'pain' and 'no-pain', respectively. We extracted a highly complex and structured mathematical feature space. We could show that (1) similarity features and an (2) automatic feature selection outperform the accuracy results from [14]. We found recognition rates with high accuracies based on a configuration of selected features. Especially the recognition rate for B vs. T₄ (pain tolerance) showed a relatively high quality with regard to Cramer’s V. To the best of our knowledge, this is the first study in the area of automated pain recognition to employ biopotential data with accuracy over 90%. Highly relevant are features based on similarity derived from zygomaticus and corrugator. Similarly to other studies of automated pain recognition via video recording, we pointed out that facial expressions are highly relevant regarding the pain intensity. The importance of features derived from similarity thus needs to be tested systematically.

Study weaknesses

We added two additional intermediate levels to obtain a ranged set of four equally distributed temperatures without a nonlinear correction.

Outlook

Our intention is to optimize our classification algorithm. In reference to that, we are currently planning future experimental procedures.

Generalizability: We will test and improve the generalizability with a complex pain model (phasic and tonic, heat and pressure).
Response Specificity: The specificity of the recognition system will be evaluated with the task to distinguish pain from psychosocial stress.
Assessment modalities: In addition to psychobiological and facial parameters, we will assess paralinguistic properties, skin temperature, body movement, and other modalities for pain recognition.
Reliability of pain recognition will be tested by repeating the experiment after (at least) one week.
All algorithms will be adapted for online processing to advance towards a pain monitoring system.

Further, the classification algorithm requires testing and optimization within a clinical environment. Finally, the goal of the project is the advancement of pain diagnosis and monitoring of pain states. With the use of multimodal sensor technology and highly effective data classification systems, reliable and valid automated pain recognition will be possible. The surrogate measurement of pain with machine learning algorithms will provide valuable information with high temporal resolution for a clinical team, which may help to objectively assess the evolution of treatments (e.g., effect of drugs for pain reduction, information of surgical indication, the quality of care provided to patients).

Author Contributions

Conceived and designed the experiments: SG SW PW HT. Performed the experiments: SG SW. Analyzed the data: SW SG AA. Contributed reagents/materials/analysis tools: AA SG SW RT SC. Wrote the paper: RT SC SW.

References

1. Handel E. Praxishandbuch ZOPA. Bern: Verlag Hans Huber; 2010.
2. Camilleri M. Opioid-induced constipation: challenges and therapeutic opportunities. Am J Gastroenterol [Internet]. American College of Gastroenterology; 2011 May [cited 2015 Feb 17];106(5):835–42; quiz 843. Available from: http://dx.doi.org/10.1038/ajg.2011.30
- View Article
- Google Scholar
3. Colloca L, Benedetti F, Pollo A. Repeatability of autonomic responses to pain anticipation and pain stimulation. Eur J Pain [Internet]. 2006;10(7):659–65. Available from: http://www.ncbi.nlm.nih.gov/pubmed/16337150
- View Article
- Google Scholar
4. Cortelli P, Pierangeli G. Chronic pain-autonomic interactions. Neurol Sci [Internet]. 2003;24 Suppl 2:S68–70. Available from: http://www.ncbi.nlm.nih.gov/pubmed/12811596
- View Article
- Google Scholar
5. Jeanne M, Logier R, De Jonckheere J, Tavernier B. Validation of a graphic measurement of heart rate variability to assess analgesia/nociception balance during general anesthesia. Conf Proc IEEE Eng Med Biol Soc [Internet]. 2009;2009:1840–3. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19963520
- View Article
- Google Scholar
6. Korhonen I, Yli-Hankala A. Photoplethysmography and nociception. Acta Anaesthesiol Scand [Internet]. 2009;53(8):975–85. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19572939
- View Article
- Google Scholar
7. Ledowski T, Ang B, Schmarbeck T, Rhodes J. Monitoring of sympathetic tone to assess postoperative pain: skin conductance vs surgical stress index. Anaesthesia [Internet]. 2009 Jul [cited 2015 Jan 22];64(7):727–31. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19183409
- View Article
- Google Scholar
8. Loggia ML, Juneau M, Bushnell MC. Autonomic responses to heat pain: Heart rate, skin conductance, and their relation to verbal ratings and stimulus intensity. Pain [Internet]. 2011;152(3):592–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21215519
- View Article
- Google Scholar
9. Schlereth T, Birklein F. The sympathetic nervous system and pain. Neuromolecular Med [Internet]. 2008 Jan [cited 2015 Feb 2];10(3):141–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17990126
- View Article
- Google Scholar
10. Storm H. Changes in skin conductance as a tool to monitor nociceptive stimulation and pain. Curr Opin Anaesthesiol [Internet]. 2008 Dec [cited 2015 Jul 25];21(6):796–804. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18997532
- View Article
- Google Scholar
11. Treister R, Kliger M, Zuckerman G, Goor Aryeh I, Eisenberg E. Differentiating between heat pain intensities: the combined effect of multiple autonomic parameters. Pain [Internet]. 2012/06/01 ed. 2012;153(9):1807–14. Available from: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=22647429
- View Article
- Google Scholar
12. Bustan S, Gonzalez-Roldan AM, Kamping S, Brunner M, Löffler M, Flor H, et al. Suffering as an independent component of the experience of pain. Eur J Pain [Internet]. 2015 Aug [cited 2015 Jul 25];19(7):1035–48. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25857478
- View Article
- Google Scholar
13. Brown JE, Chatterjee N, Younger J, Mackey S. Towards a physiology-based measure of pain: patterns of human brain activity distinguish painful from non-painful thermal stimulation. PLoS One [Internet]. 2011 Jan [cited 2015 Apr 27];6(9):e24124. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3172232&tool=pmcentrez&rendertype=abstract
- View Article
- Google Scholar
14. Walter S, Gruss S, Limbrecht-Ecklundt K, Traue HC, Werner P, Al-Hamadi A, et al. Automatic pain quantification using autonomic parameters. Psychol Neurosci [Internet]. 2014 [cited 2015 Feb 17];7(3):363–80. Available from: http://www.psycneuro.org/index.php/path/article/view/372
- View Article
- Google Scholar
15. Werner P, Al-Hamadi A, Niese R, Walter S, Gruss S, Traue HC. Automatic Pain Recognition from Video and Biomedical Signals. 22nd International Conference on Pattern Recognition (ICPR). Stockholm; 2014. p. 4582–7.
16. Hua-Mei C, Varshney PK, Arora MK. Performance of mutual information similarity measure for registration of multitemporal remote sensing images. Geosci Remote Sensing, IEEE Trans. 2003;41(11):2445–54.
- View Article
- Google Scholar
17. Kennedy HL. A New Statistical Measure of Signal Similarity. 2007 Information, Decision and Control [Internet]. IEEE; 2007 [cited 2015 Apr 27]. p. 112–7. Available from: http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=4252487
18. Walter S, Gruss S, Ehleiter H, Tan J, Traue HC, Crawcour S, et al. The biovid heat pain database data for the advancement and systematic validation of an automated pain recognition system. 2013 IEEE International Conference on Cybernetics (CYBCO) [Internet]. IEEE; 2013 [cited 2015 Jan 16]. p. 128–31. Available from: http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=6617456
19. Tan J- W, Walter S, Scheck A, Hrabal D, Hoffmann H, Kessler H, et al. Repeatability of facial electromyography (EMG) activity over corrugator supercilii and zygomaticus major on differentiating various emotions. J Ambient Intell Humaniz Comput [Internet]. Springer-Verlag; 2012;3(1):3–10. Available from: http://dx.doi.org/10.1007/s12652-011-0084-9
- View Article
- Google Scholar
20. Kim J, Andre E. Emotion Recognition Based on Physiological Changes in Music Listening. IEEE Trans Pattern Anal Mach Intell [Internet]. 2008;30(12):2067–83. Available from: <Go to ISI>://000260033900001
- View Article
- Google Scholar
21. Andrade AO, Kyberd P, Nasuto SJ. The application of the Hilbert spectrum to the analysis of electromyographic signals. Inf Sci (Ny) [Internet]. 2008;178(9):2176–93. Available from: http://www.sciencedirect.com/science/article/pii/S0020025508000029
- View Article
- Google Scholar
22. Andrade AO, Nasuto SJ, Kyberd P. Extraction of motor unit action potentials from electromyographic signals through generative topographic mapping. J Franklin Inst [Internet]. 2007;344(3–4):154–79. Available from: http://www.sciencedirect.com/science/article/pii/S0016003206001487
- View Article
- Google Scholar
23. Nakano K, Ota Y, Ukai H. Frequence Detection Method based on Recursive DFT Alhorithm. 14 th PSCC [Internet]. Sevilla; 2002. p. 24–8. Available from: http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=1208376
24. Andrade A. Decomposition and analysis of electromyographic signals. [Internet]. University of Reading. 2005 [cited 2015 Feb 17]. Available from: http://scholar.google.com/citations?view_op=view_citation&hl=de&user=DeBBjAwAAAAJ&citation_for_view=DeBBjAwAAAAJ:qjMakFHDy7sC
- View Article
- Google Scholar
25. Cao C, Slobounov S. Application of a novel measure of EEG non-stationarity as ‘Shannon- entropy of the peak frequency shifting’ for detecting residual abnormalities in concussed individuals. Clin Neurophysiol [Internet]. 2011;122(7):1314–21. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21216191
- View Article
- Google Scholar
26. Chen W, Zhuang J, Yu W, Wang Z. Measuring complexity using FuzzyEn, ApEn, and SampEn. Med Eng Phys [Internet]. Elsevier; 2009;31(1):61–8. Available from: http://www.medengphys.com/article/S1350-4533(08)00072-6/abstract
- View Article
- Google Scholar
27. Mitchell R, Agle B, Wood D. Toward a Theory of Stakeholder Identification and Salience: Defining the Principle of Who and What Really Counts on JSTOR. Acad Manag Rev [Internet]. 1997 [cited 2015 Feb 17];22(4):853–86. Available from: https://www.jstor.org/stable/259247?seq=1#page_scan_tab_contents
- View Article
- Google Scholar
28. Kapoor A, Burleson W, Picard RW. Automatic prediction of frustration. Int J Hum Comput Stud. 2007;65:724–36.
- View Article
- Google Scholar
29. Hsu C, Chang C, Lin C. A Practical Guide to Support Vector Classication [Internet]. 2003. Available from: http://www.csie.ntu.edu.tw/~cjlin/papers/guide/guide.pdf
- View Article
- Google Scholar
30. Boser B, Guyon I, Vapnik V. A Training Algorithm for Optimal Margin Classifiers. Proceedings of the 5th Annual ACM Workshop on Computational Learning Theory [Internet]. 1992 [cited 2015 Feb 17]. p. 1399–44. Available from: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.21.3818
- View Article
- Google Scholar
31. Rychetsky M. Algorithms and Architectures for Machine Learning based on Regularized Neural Networks and Support Vector Approaches. Shaker Verlag GmbH; 2001. 259 p.
- View Article
- Google Scholar
32. scikit-learn: Machine Learning in Python [Internet]. Available from: http://scikit-learn.org/stable/index.html
33. Kolodyazhniy V, Kreibig SD, Gross JJ, Roth WT, Wilhelm FH. An affective computing approach to physiological emotion specificity: Toward subject-independent and stimulus-independent classification of film-induced emotions. Psychophysiology [Internet]. 2011;48(7):908–22. Available from: <Go to ISI>://000291255500004
- View Article
- Google Scholar
34. Fareed A, Caroline H. RapidMiner 5 [Internet]. Operator Reference. 2012 [cited 2015 Feb 17]. Available from: https://rapidminer.com/wp-content/uploads/2013/10/RapidMiner_OperatorReference_en.pdf
35. Cramer H. Mathematical Methods of Statistics (PMS-9). Princeton: Princeton University Press; 1999 [cited 2015 Feb 17]. p. 282. Available from: http://press.princeton.edu/titles/391.html
36. Rea LM, Parker R. Designing and Conducting Survey Research: A Comprehensive Guide [Internet]. 4th ed. Jossey-Bass; 2014 [cited 2015 Feb 17]. 87 p. Available from: http://eu.wiley.com/WileyCDA/WileyTitle/productCd-1118767039,subjectCd-EV00.html

[ref1] 1. Handel E. Praxishandbuch ZOPA. Bern: Verlag Hans Huber; 2010.

[ref2] 2. Camilleri M. Opioid-induced constipation: challenges and therapeutic opportunities. Am J Gastroenterol [Internet]. American College of Gastroenterology; 2011 May [cited 2015 Feb 17];106(5):835–42; quiz 843. Available from: http://dx.doi.org/10.1038/ajg.2011.30
View Article
Google Scholar

[3] View Article

[4] Google Scholar

[ref3] 3. Colloca L, Benedetti F, Pollo A. Repeatability of autonomic responses to pain anticipation and pain stimulation. Eur J Pain [Internet]. 2006;10(7):659–65. Available from: http://www.ncbi.nlm.nih.gov/pubmed/16337150
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref4] 4. Cortelli P, Pierangeli G. Chronic pain-autonomic interactions. Neurol Sci [Internet]. 2003;24 Suppl 2:S68–70. Available from: http://www.ncbi.nlm.nih.gov/pubmed/12811596
View Article
Google Scholar

[9] View Article

[10] Google Scholar

[ref5] 5. Jeanne M, Logier R, De Jonckheere J, Tavernier B. Validation of a graphic measurement of heart rate variability to assess analgesia/nociception balance during general anesthesia. Conf Proc IEEE Eng Med Biol Soc [Internet]. 2009;2009:1840–3. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19963520
View Article
Google Scholar

[12] View Article

[13] Google Scholar

[ref6] 6. Korhonen I, Yli-Hankala A. Photoplethysmography and nociception. Acta Anaesthesiol Scand [Internet]. 2009;53(8):975–85. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19572939
View Article
Google Scholar

[15] View Article

[16] Google Scholar

[ref7] 7. Ledowski T, Ang B, Schmarbeck T, Rhodes J. Monitoring of sympathetic tone to assess postoperative pain: skin conductance vs surgical stress index. Anaesthesia [Internet]. 2009 Jul [cited 2015 Jan 22];64(7):727–31. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19183409
View Article
Google Scholar

[18] View Article

[19] Google Scholar

[ref8] 8. Loggia ML, Juneau M, Bushnell MC. Autonomic responses to heat pain: Heart rate, skin conductance, and their relation to verbal ratings and stimulus intensity. Pain [Internet]. 2011;152(3):592–8. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21215519
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref9] 9. Schlereth T, Birklein F. The sympathetic nervous system and pain. Neuromolecular Med [Internet]. 2008 Jan [cited 2015 Feb 2];10(3):141–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17990126
View Article
Google Scholar

[24] View Article

[25] Google Scholar

[ref10] 10. Storm H. Changes in skin conductance as a tool to monitor nociceptive stimulation and pain. Curr Opin Anaesthesiol [Internet]. 2008 Dec [cited 2015 Jul 25];21(6):796–804. Available from: http://www.ncbi.nlm.nih.gov/pubmed/18997532
View Article
Google Scholar

[27] View Article

[28] Google Scholar

[ref11] 11. Treister R, Kliger M, Zuckerman G, Goor Aryeh I, Eisenberg E. Differentiating between heat pain intensities: the combined effect of multiple autonomic parameters. Pain [Internet]. 2012/06/01 ed. 2012;153(9):1807–14. Available from: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=22647429
View Article
Google Scholar

[30] View Article

[31] Google Scholar

[ref12] 12. Bustan S, Gonzalez-Roldan AM, Kamping S, Brunner M, Löffler M, Flor H, et al. Suffering as an independent component of the experience of pain. Eur J Pain [Internet]. 2015 Aug [cited 2015 Jul 25];19(7):1035–48. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25857478
View Article
Google Scholar

[33] View Article

[34] Google Scholar

[ref13] 13. Brown JE, Chatterjee N, Younger J, Mackey S. Towards a physiology-based measure of pain: patterns of human brain activity distinguish painful from non-painful thermal stimulation. PLoS One [Internet]. 2011 Jan [cited 2015 Apr 27];6(9):e24124. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3172232&tool=pmcentrez&rendertype=abstract
View Article
Google Scholar

[36] View Article

[37] Google Scholar

[ref14] 14. Walter S, Gruss S, Limbrecht-Ecklundt K, Traue HC, Werner P, Al-Hamadi A, et al. Automatic pain quantification using autonomic parameters. Psychol Neurosci [Internet]. 2014 [cited 2015 Feb 17];7(3):363–80. Available from: http://www.psycneuro.org/index.php/path/article/view/372
View Article
Google Scholar

[39] View Article

[40] Google Scholar

[ref15] 15. Werner P, Al-Hamadi A, Niese R, Walter S, Gruss S, Traue HC. Automatic Pain Recognition from Video and Biomedical Signals. 22nd International Conference on Pattern Recognition (ICPR). Stockholm; 2014. p. 4582–7.

[ref16] 16. Hua-Mei C, Varshney PK, Arora MK. Performance of mutual information similarity measure for registration of multitemporal remote sensing images. Geosci Remote Sensing, IEEE Trans. 2003;41(11):2445–54.
View Article
Google Scholar

[43] View Article

[44] Google Scholar

[ref17] 17. Kennedy HL. A New Statistical Measure of Signal Similarity. 2007 Information, Decision and Control [Internet]. IEEE; 2007 [cited 2015 Apr 27]. p. 112–7. Available from: http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=4252487

[ref18] 18. Walter S, Gruss S, Ehleiter H, Tan J, Traue HC, Crawcour S, et al. The biovid heat pain database data for the advancement and systematic validation of an automated pain recognition system. 2013 IEEE International Conference on Cybernetics (CYBCO) [Internet]. IEEE; 2013 [cited 2015 Jan 16]. p. 128–31. Available from: http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=6617456

[ref19] 19. Tan J- W, Walter S, Scheck A, Hrabal D, Hoffmann H, Kessler H, et al. Repeatability of facial electromyography (EMG) activity over corrugator supercilii and zygomaticus major on differentiating various emotions. J Ambient Intell Humaniz Comput [Internet]. Springer-Verlag; 2012;3(1):3–10. Available from: http://dx.doi.org/10.1007/s12652-011-0084-9
View Article
Google Scholar

[48] View Article

[49] Google Scholar

[ref20] 20. Kim J, Andre E. Emotion Recognition Based on Physiological Changes in Music Listening. IEEE Trans Pattern Anal Mach Intell [Internet]. 2008;30(12):2067–83. Available from: <Go to ISI>://000260033900001
View Article
Google Scholar

[51] View Article

[52] Google Scholar

[ref21] 21. Andrade AO, Kyberd P, Nasuto SJ. The application of the Hilbert spectrum to the analysis of electromyographic signals. Inf Sci (Ny) [Internet]. 2008;178(9):2176–93. Available from: http://www.sciencedirect.com/science/article/pii/S0020025508000029
View Article
Google Scholar

[54] View Article

[55] Google Scholar

[ref22] 22. Andrade AO, Nasuto SJ, Kyberd P. Extraction of motor unit action potentials from electromyographic signals through generative topographic mapping. J Franklin Inst [Internet]. 2007;344(3–4):154–79. Available from: http://www.sciencedirect.com/science/article/pii/S0016003206001487
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref23] 23. Nakano K, Ota Y, Ukai H. Frequence Detection Method based on Recursive DFT Alhorithm. 14 th PSCC [Internet]. Sevilla; 2002. p. 24–8. Available from: http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=1208376

[ref24] 24. Andrade A. Decomposition and analysis of electromyographic signals. [Internet]. University of Reading. 2005 [cited 2015 Feb 17]. Available from: http://scholar.google.com/citations?view_op=view_citation&hl=de&user=DeBBjAwAAAAJ&citation_for_view=DeBBjAwAAAAJ:qjMakFHDy7sC
View Article
Google Scholar

[61] View Article

[62] Google Scholar

[ref25] 25. Cao C, Slobounov S. Application of a novel measure of EEG non-stationarity as ‘Shannon- entropy of the peak frequency shifting’ for detecting residual abnormalities in concussed individuals. Clin Neurophysiol [Internet]. 2011;122(7):1314–21. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21216191
View Article
Google Scholar

[64] View Article

[65] Google Scholar

[ref26] 26. Chen W, Zhuang J, Yu W, Wang Z. Measuring complexity using FuzzyEn, ApEn, and SampEn. Med Eng Phys [Internet]. Elsevier; 2009;31(1):61–8. Available from: http://www.medengphys.com/article/S1350-4533(08)00072-6/abstract
View Article
Google Scholar

[67] View Article

[68] Google Scholar

[ref27] 27. Mitchell R, Agle B, Wood D. Toward a Theory of Stakeholder Identification and Salience: Defining the Principle of Who and What Really Counts on JSTOR. Acad Manag Rev [Internet]. 1997 [cited 2015 Feb 17];22(4):853–86. Available from: https://www.jstor.org/stable/259247?seq=1#page_scan_tab_contents
View Article
Google Scholar

[70] View Article

[71] Google Scholar

[ref28] 28. Kapoor A, Burleson W, Picard RW. Automatic prediction of frustration. Int J Hum Comput Stud. 2007;65:724–36.
View Article
Google Scholar

[73] View Article

[74] Google Scholar

[ref29] 29. Hsu C, Chang C, Lin C. A Practical Guide to Support Vector Classication [Internet]. 2003. Available from: http://www.csie.ntu.edu.tw/~cjlin/papers/guide/guide.pdf
View Article
Google Scholar

[76] View Article

[77] Google Scholar

[ref30] 30. Boser B, Guyon I, Vapnik V. A Training Algorithm for Optimal Margin Classifiers. Proceedings of the 5th Annual ACM Workshop on Computational Learning Theory [Internet]. 1992 [cited 2015 Feb 17]. p. 1399–44. Available from: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.21.3818
View Article
Google Scholar

[79] View Article

[80] Google Scholar

[ref31] 31. Rychetsky M. Algorithms and Architectures for Machine Learning based on Regularized Neural Networks and Support Vector Approaches. Shaker Verlag GmbH; 2001. 259 p.
View Article
Google Scholar

[82] View Article

[83] Google Scholar

[ref32] 32. scikit-learn: Machine Learning in Python [Internet]. Available from: http://scikit-learn.org/stable/index.html

[ref33] 33. Kolodyazhniy V, Kreibig SD, Gross JJ, Roth WT, Wilhelm FH. An affective computing approach to physiological emotion specificity: Toward subject-independent and stimulus-independent classification of film-induced emotions. Psychophysiology [Internet]. 2011;48(7):908–22. Available from: <Go to ISI>://000291255500004
View Article
Google Scholar

[86] View Article

[87] Google Scholar

[ref34] 34. Fareed A, Caroline H. RapidMiner 5 [Internet]. Operator Reference. 2012 [cited 2015 Feb 17]. Available from: https://rapidminer.com/wp-content/uploads/2013/10/RapidMiner_OperatorReference_en.pdf

[ref35] 35. Cramer H. Mathematical Methods of Statistics (PMS-9). Princeton: Princeton University Press; 1999 [cited 2015 Feb 17]. p. 282. Available from: http://press.princeton.edu/titles/391.html

[ref36] 36. Rea LM, Parker R. Designing and Conducting Survey Research: A Comprehensive Guide [Internet]. 4th ed. Jossey-Bass; 2014 [cited 2015 Feb 17]. 87 p. Available from: http://eu.wiley.com/WileyCDA/WileyTitle/productCd-1118767039,subjectCd-EV00.html

Figures

Abstract

Background

Methods

Results

Conclusion

Introduction

Methods

Subjects

Exclusion criteria

Individual calibration of pain and tolerance threshold

Experimental pain stimulation

Measurement of biopotentials

Biopotentials.

Electromyography (EMG).

Skin conductance level (SCL).

Electrocardiogram (ECG).

Preprocessing

Feature extraction

Machine learning with Support Vector Machine (SVM)

Feature selection

Validity of classification rates

Results

Discussion and Conclusions

Study weaknesses

Outlook

Author Contributions

References