Characterizing psychological dimensions in non-pathological subjects through autonomic nervous system dynamics

Abstract

The objective assessment of psychological traits of healthy subjects and psychiatric patients has been growing interest in clinical and bioengineering research fields during the last decade. Several experimental evidences strongly suggest that a link between Autonomic Nervous System (ANS) dynamics and specific dimensions such as anxiety, social phobia, stress and emotional regulation might exist. Nevertheless, an extensive investigation on a wide range of psycho-cognitive scales and ANS non-invasive markers gathered from standard and nonlinear analysis still needs to be addressed. In this study, we analyzed the discerning and correlation capabilities of a comprehensive set of ANS features and psycho-cognitive scales in 29 non-pathological subjects monitored during resting conditions. In particular, the state of the art of standard and nonlinear analysis was performed on Heart Rate Variability, InterBreath Interval series, and Inter-Beat Respiration series, which were considered as monovariate and multivariate measurements. Experimental results show that each ANS feature is linked to specific psychological traits. Moreover, nonlinear analysis outperforms the psychological assessment with respect to standard analysis. Considering that the current clinical practice relies only on subjective scores from interviews and questionnaires, this study provides objective tools for the assessment of psychological dimensions.

1. Introduction

Psychological assessment refers to the practice of standardized evaluation of performance or impairment in different domains of thinking, learning and behavior. Accordingly, such an assessment can be used to characterize and quantify different behaviors in healthy subjects or to reveal the presence of behavioral disorders such as anxiety and social phobia. Depending on the factors under observation, psychological assessment can be achieved via different routes: behavioral tasks, questionnaires, or interviews. The evaluation is done by a professional (i.e., certified psychologist) in order to obtain a standardized and quantifiable information of the subject under study (Cohen et al., 1992). These approaches are useful in performing an individual assessment for which the performance of one person can be interpreted through pre-existing norms, as well as in group assessment which allows for different comparisons (within a single group or between groups) (Kenny et al., 2008). It is worthwhile noting that self-report questionnaires and interviews currently represent the standard clinical practice in diagnosing psychiatric disorders (Cohen et al., 1992; Valenza et al., 2013a, 2014c).

Nevertheless, several issues in these kinds of approaches still need to be addressed. First, the scores are obtained with subjective procedures which might be biased by possible social desirability thoughts of the subject and possible recent emotional events. Moreover, professionals need to choose the appropriate test for each psychological dimension and subject, and verify that it has good psychometric properties in order to adhere to the evidence-based paradigm (i.e., reliability and validity) (Groth-Marnat, 2003; Hunsley and Mash, 2010). To overcome these problems, several efforts have been made in psycho-physiological and bioengineering research fields to objectify the psychological assessment. In particular, physiological correlates of the central and autonomic nervous systems (CNS and ANS, respectively) have been extensively studied and taken into account (Taillard et al., 1990, 1993; Carney et al., 1995; Glassman, 1998; Stampfer, 1998; Iverson et al., 2002, 2005; Watkins et al., 2002; Calvo and D'Mello, 2010; Lin et al., 2010; Petrantonakis and Hadjileontiadis, 2011; Valenza et al., 2012a,b, 2013a,b, 2014c).

To give some significant examples, physiological correlates of mood disorders such as bipolar disorders have been found on sleep (Stampfer, 1998; Iverson et al., 2002, 2005), hormonal system (Carney et al., 1995; Glassman, 1998; Watkins et al., 2002), and ANS dynamics through heartbeat and respiratory dynamics (Taillard et al., 1990, 1993; Valenza et al., 2013a, 2014c). Moreover, as the psychological dimensions can be related to variations of emotional states, several computational methods for automatic emotion recognition have been developed using electroencephalogram (EEG) and ANS signal analysis (Taillard et al., 1990, 1993; Calvo and D'Mello, 2010; Lin et al., 2010; Petrantonakis and Hadjileontiadis, 2011; Valenza et al., 2012a,b, 2013a,b,2014c).

Here we focus on the link between ANS dynamics and psychological dimensions. This choice is justified by the fact that ANS dynamics cannot be straightforwardly changed by the subject intention and is under direct control of CNS pathways such as the prefrontal cortex, amygdala, and brainstem (Ruiz-Padial et al., 2011). Of note, dysfunctions on these CNS recruitment circuits lead to pathological effects (Heller et al., 2009) such as anhedonia, i.e., the loss of pleasure or interest in previously rewarding stimuli, which is a core feature of major depression and other serious mood disorders. Moreover, ANS monitoring is widely available, cost-effective, and can be easily performed through wearable systems such as sensorized t-shirts (Valenza et al., 2008, 2014c) or gloves (Lanatà et al., 2012), and its dynamics is thought to be less sensitive to artifact events than in the EEG case.

ANS dynamics has been demonstrated to provide effective markers of typical psychological processes. As a matter of fact, previous studies (Freeman and Nixon, 1985; Yeragani et al., 1999; Virtanen et al., 2003; Cohen and Benjamin, 2006; Shinba et al., 2008; Licht et al., 2009; Thayer et al., 2010; 2012) suggest that patients with anxiety are at increased risk for heart disease (e.g., the association between phobic anxiety or panic disorder and somatic morbidity as coronary heart disease, coronary spasm and ventricular arrhythmia). ANS markers of anxiety and panic disorders can be found through the analysis of the Heart Rate Variability (HRV), revealing an increased heart rate and decreased power in low-frequency (LF) and high frequency (HF) bands. A decreased HF spectral power of HRV was also found in patients affected by generalized anxiety disorder (Thayer et al., 1996), whereas a decreased heart rate was also found in autism spectrum disorders (Jansen et al., 2006) in response to stress. This change could be related to abnormal high basal (nor)epinephrine levels. On the contrary, increased mean heart rate associated to a reduced variability has been observed in depressed patients (Carney et al., 2005). Moreover, it has been shown how subjects reporting excessive and persistent fear of social situations are characterized by atypical ANS dynamics which is evident in variables as HRV mean, respiration rate, tidal volume, and blood pressure (Grossman et al., 2001). ANS markers gathered from non-linear analysis were related to phycological dimensions as anxiety (Cohen and Benjamin, 2006) and panic disorder through symbolic analysis (Yeragani et al., 2000). Despite the elevated number of previous studies, none of these researches have reached an acceptable level of accuracy to effectively, reliably, and objectively characterize the psychological dimensions of healthy subjects and psychiatric patients, and to forecast a clinical course. A possible reason can be related to the limited amount of ANS features and specific psychological traits that were taken into account.

Therefore, here we present a detailed study on psychological assessments through an extensive analysis of the ANS dynamics. Psychological dimensions were quantified by means of the 6 psycho-cognitive scales (see details on the series definition, estimation, and parameter extraction in Section 2.3).

In order to perform a comprehensive study, the ANS non-linear dynamics has to be taken into account. Although the detailed physiology behind such complex dynamics has not been completely clarified, it is worthwhile noting that ANS non-linear dynamics plays a crucial role in most of the underlying biological processes, as they have been proven to be of prognostic value in aging and diseases, showing robust and effective discerning and characterizing properties (Poon and Merrill, 1997; Glass, 2001; Goldberger et al., 2002; Stiedl and Meyer, 2003; Tulppo et al., 2005; Atyabi et al., 2006; Glass, 2009; Wu et al., 2009; Citi et al., 2012; Valenza et al., 2014a). Indeed, physiological systems are intrinsically non-linear systems characterized by multi-feedback interactions associated to long-range correlations (Marmarelis, 2004), likely due to the enormous amount of structural units inside them and to the various non-linear neural interactions and integrations occurring at the neuron and receptor levels. The study of the complexity of physiological signals, in particular, has led to important results in recent decades in understanding the mechanisms underlying mental illness (Yang and Tsai, 2012). Several measures of complexity have also been proposed and applied to the study of mental illness based on various biomedical signals, from EEG (Hu et al., 2006; Takahashi et al., 2010; Gao et al., 2011), to MEG (Fernandez et al., 2010), through HRV (Mujica-Parodi et al., 2005; Hu et al., 2009, 2010; Gao et al., 2013; Valenza et al., 2014b). Accordingly, in this study we investigate the role of ANS non-linear dynamics in performing the psychological assessment, with respect to the standard analysis, i.e., analysis in the time and frequency domain.

2. Materials and Methods

2.1. Subjects Recruitment, Experimental Protocol, and Acquisition Set-up

A group of 29 non-pathological subjects (5 males), i.e., not suffering from both cardiovascular and evident mental pathologies, was recruited to participate in the experiment. Subjects were students recruited from the Babes-Bolyai University, via an online screening questionnaire assessing their intention to take part in the study. Participation was voluntary and each subjects signed a written informed consent after the study procedure had been explained. No compensation for participation was offered. Subjects underwent a medical screening interview to assess the presence of any medical condition or medication that might have interfered with their cardiovascular data. Their age ranged from 21 to 35 and were naive to the purpose of the experiment. The group was as heterogeneous as possible in order to have a wide range of psycho-cognitive-behavioral dimensions. The experimental protocol was structured in the following two phases: (1) submission of self report psycho-behavioral tests; (2) recording of the physiological signs. More in detail, all participants were screened by 6 self-report questionnaires (see details below), which were comprised of a total of 25 sub-scales. Then, physiological signals such as ElectroCardioGram (ECG), Respiration (RSP) were simultaneously acquired during resting state condition for 25 min through the BIOPAC MP150 device. The sampling rate was 1000 Hz for all signals. We used the ECG100C Electrocardiogram Amplifier from BIOPAC inc., connected with pregelled Ag/AgCl electrodes placed following Einthoven triangle configuration. The dedicated module of BIOPAC MP150 used to record the respiration activity is RSP100C Respiration Amplifier with the TSD201 sensor, which is a piezo-resistive sensor with the output resistance within the range 5–125 KOhm and bandwidth of 0.05–10 Hz. This piezoresistive sensor changed its electrical resistance if stretched or shortened, and it was sensitive to the thoracic circumference variations occurring during respiration.

The ECG signal was used to extract the HRV series, which refer to the variation of the time intervals between consecutive heartbeats identified with R-waves (RR intervals). Two different time series were extracted from the respiration activity: InterBreath Interval time series (IBI) and InterBeat Respiration (IBR). The IBI series was obtained detecting the local maxima of each respiratory act, whereas IBR consists of the amplitude of the respiration activity signal when sampled at the R-peak times.

2.2. Scales for the Assessment of Psychological Dimensions

In this work, we used a total of 6 self-report questionnaires in which, for most of them, different sub-scales are considered. The total number of sub-scales used in this experiment was 25. A Cronbach's α measure (Bland and Altman, 1997) is assigned to each scale and represents the consistency of the test. Such an α index depends on the number and average inter-correlation among the test questions. Details on each scale and related sub-scales are as follows:

• Positive and Negative Affect Schedule (PANAS, Watson et al., 1988). The PANAS contains 2 sub-scales—positive affect (PA) and negative affect (NA)—of 10 items describing emotions each. The scale has good reliability (Cronbach's α = 0.88 for the PA and 0.87 for the NA sub-scale, respectively) and good construct validity. Cronbach's α for this sample was 0.87 for PA and 0.90 for NA, supporting good internal consistency.

• Liebowitz Social Anxiety Scale (LSAS, Liebowitz, 1987). The LSAS is a self-assessment social phobia questionnaire containing 24 items describing actions done in social situations, grouped at first in 2 sub-scales (social interaction and performance). Subjects rate these situations in terms of fear/anxiety and avoidance, allowing for a total of 4 separate sub-scales. The scale presents a very good internal consistency (Cronbach's α = 0.96) as well as good convergent and divergent validity (Heimberg et al., 1999). Cronbach's α for this sample was 0.92, showing a very good internal consistency.

• Difficulties in Emotion Regulation (DERS, Gratz and Roemer, 2004). The DERS is a 36-item self-report scale measuring emotion dysregulation. The scale offers an overall score as well as scores for each of the 6 sub-scales related to DERS (Non-acceptance of Emotional Responses, Difficulties Engaging in Goal-Directed Behavior, Impulse Control Difficulties, Lack of Emotional Awareness, Limited Access to Emotion Regulation Strategies, and Lack of Emotional Clarity). Internal consistency for this scale is excellent (Cronbach's α = 0.93) and construct and predictive validity are considered adequate.

• Interpersonal Reactivity Index (IRI, Davis, 1980). The IRI is a 28-item questionnaire measuring empathy. The scale provides scores for 4 sub-scales (Fantasy, Perspective-taking, Empathic Concern, and Personal Distress), as well as a general score of empathy. Internal consistency of the four sub-scales is acceptable (ranging from α = 0.70–0.78).

• Behavioral Inhibition/Behavioral activation Scales (BIS/BAS, Carver and White, 1994). The BIS/BAS scale is composed of 20 items comprised in 4 sub-scales (Inhibition, Reward Responsiveness, Drive, and Fun Seeking), measuring behavioral inhibition and activation sensitivity. The scale has been adapted on the Romanian population, showing good construct validity and acceptable internal consistency (ranging from α = 0.62–0.81) (Sava and Sperneac, 2006).

• Zuckerman Kuhlman Personality Questionnaire (ZKPQ, Zuckerman et al., 1993). The ZKPQ represents a five-factor (Impulsive Sensation Seeking, Neuroticism-Anxiety, Aggresion-Hostility, Sociability, and Activity) personality inventory containing 99 true-false items, therefore we used 5 sub-scales. The Romanian adaptation of this scale presents adequate internal consistency (α ranging from 0.69 to 0.88) and good convergent validity (Sârbescu and Neguţ, 2012).

2.3. Methodology of Signal Processing

In this section, the methodology of signal processing applied to the Heart Rate Variability (HRV), InterBreath Interval (IBI), and InterBeat Respiration (IBR) series is reported in detail. HRV refers to the variability of the series comprised of the distances between two consecutive R-waves detected from the Electrocardiogram, i.e., the R-R intervals. IBI is the series comprised of the distances between two consecutive local maxima of the respiration activity (the two maxima within two respiratory acts), whereas IBR series is the respiratory activity sampled at times corresponding to the R-peaks. Standard and non-linear monovariate and multivariate measure are extracted from each series in order to investigate a wide set of parameters characterizing the ANS linear and non-linear dynamics acting of the cardio-respiratory control.

2.3.1. Standard Measures

Standard analysis was performed on HRV series in order to extract parameters defined in the time and frequency domain (Camm et al., 1996; Acharya et al., 2006; Valenza et al., 2012b). Time domain features include statistical parameters and morphological indexes. More specifically, concerning the time domain analysis, in addition to the first (meanRR) and second order moment (SDNN) of the RR intervals, so-called normal-to-normal (NN) intervals, the square root of the mean of the sum of the squares of differences between subsequent NN intervals (RMSSD=1N−1∑j=1N−1(RRj+1−RRj2) and the number of successive differences of intervals which differ by more than 50 ms, expressed as a percentage of the total number of heartbeats analyzed (pNN50=NN50N−1100%) were calculated. Moreover, the triangular index (TINN) was estimated as the base of a triangle which better approximated the NN interval distribution (the minimum square difference is used to find such a triangle).

Concerning the frequency domain analysis, several features were calculated from the Power Spectral Density (PSD) analysis. In this work, PSD was estimated by using the Welch's periodogram, which uses the FFT (Fast Fourier Transform) algorithm. Window's width and overlap were chosen as a best compromise between the frequency resolution and variance of the estimated spectrum. Given the PSD, three spectral bands are defined as follows: VLF (very low frequency) with spectral components below 0.04 Hz; LF (low frequency), ranging between 0.04 and 0.15 Hz; HF (high frequency), comprising frequencies between 0.15 and 0.4 Hz. For each of the three frequency bands, the frequency having maximum magnitude (VLF peak, LF peak, and HF peak), the power expressed as percentage of the total power (VLF power %, LF power %, and HF power %), and the power normalized to the sum of the LF and HF power (LF power nu and HF power nu) were also evaluated. Moreover, the LF/HF power ratio was calculated.

2.3.2. Non-Linear Analysis

From the HRV, IBI, and IBR series, several non-linear measures were calculated. Such indices refer to the estimation and characterization of the phase space (or state space) of the physiological system generating the series. The phase space estimation involved the Takens method (Takens, 1981; Casdagli et al., 1991) and three parameters: m, the embedding dimension, which is a positive integer, τ, the time delay, and r, which is a positive real number and represents the margin of tolerance of the trajectories within the space. Takens theory allows for the reconstruction of the dynamic systems of different nature from time series through the method of “delayed outputs.” Starting from a time series

X=[u(T),u(2T),… u(NT)]

the attractors of the discrete dynamical system are rebuilt in a m-dimensional space, operating a delay τ on the signal. This allows achieving N − (m − 1) signals of length m starting from only one:

{X1=[u(T),u(2T),…u(mT)]X2=[u(2T),u(2T+2τ),…u(2T+(m−1)τ)]…XN−(m−1)=[u(N−(m−1))T),…u(N−(m−1))T +(m−1)τ)]

The various vectors Xj are the “delayed coordinates” and the derived m-dimensional space is called “reconstructed space.” From the state space theory, several ANS non-linear parameters can be derived using the following analyses:

• Poincarí Plot

• Recurrence Plot

• Correlation dimension, Approximate, and Sample Entropy

• Detrended Fluctuation Analysis

• Multiscale Entropy and Multivariate Multiscale Entropy Analysis

2.3.2.1. Poincarí Plot

This technique quantifies the fluctuations of the dynamics of the time series through a map of each point RR(n) of the RR series vs. the previous one. The quantitative analysis from the graph can be made by calculating the standard deviation of the points by the straight line RRj + 1 = RRj. The first standard deviation, SD1, is related to the points that are perpendicular to the line-of-identity and describes the short-term variability, whereas the second, SD2, describes the long-term variability.

2.3.2.2. Recurrence Plot

RP is a graphical method to investigate and quantify the time series complexity. The estimation starts from vectors

uj=(RRj, RRj+τ,…,RRj+(m−1)τ) j=1,2,…,N−(m−1)τ.

RP is a symmetrical square matrix of zeros and ones, whose dimensions are N − (m − 1)τ, and each element is given by

RP(j,k)= {1 if d(uj−uk≤r)0otherwise

where d is the Euclidean distance.

Several features can be extracted from the RP by means of the Recurrence Quantification Analysis (RQA). In particular, in this study the following RQA indices were taken into account: longest diagonal line (RP Lmax) and average diagonal line length (RP Lmean), divergence (RP DIV), the percentage of recurrence points which form diagonal lines recurrence rate, determinism (RP DET), trend (RP REC) and entropy (RP ShanEn) (Zbilut et al., 1990; Marwan et al., 2002, 2007).

2.3.2.3. Correlation dimension, Approximate, and Sample Entropy Measures

Starting from the vectors X1, X2, …, XN − m + 1 in ???^m, the distance between two vectors Xi and Xj, according to the definition of Takens applied to high dimensional deterministic systems is given by Takens (1981) and Schouten et al. (1994):

d[Xi,Xj]=maxk=1,2,…,m|u(i+k−1)−u(j+k−1)| (1)

For each i, with 1 ≤ i ≤ N − m + 1, we measured a parameter C^mi(r):

Cim(r)=Number of j such that(d[Xi,Xj] ≤ r)N−m+1 (2)

and we defined

Cm(r)=∑i=1N−m+1logCim(r)N−m+1 (3)

The correlation dimension (CD) is given by Theiler (1987)

CD=limr→0limN→inflogCm(r)logr

The calculation of ApEn used in this study refers to the expression (Pincus, 1991; Fusheng et al., 2001):

ApEn(m,r,N)=[Cm(r)−Cm+1(r)] (4)

SampEn is a remake of ApEn and measures the number of pairs of vectors of length m considered “neighbors,” i.e., whose distance is less than r, even if the dimension of pattern increases from m to m + 1. Unlike ApEn(m, r, N), SampEn does not include the distance of vectors with themselves, i.e., self-matches, as suggested in the later work of Grassberger and co-workers (Grassberger and Procaccia, 1983; Grassberger, 1988) and it has the advantage of being less dependent on time series length, showing relative consistency over a broader range of possible r-, m-, and N-values. By renaming C^m(r) parameters without self-matches with the notation U^m(r), SampEn is calculated by the following expression (Richman and Moorman, 2000):

SampEn(m,r,N)=−lnUm+1Um (5)

2.3.2.4. Detrended Fluctuation Analysis

The detrended fluctuation analysis features (DFA1 and DFA2) (Peng et al., 1995; Penzel et al., 2003) were evaluated to study short- and long-term autocorrelation of the HRV series. The algorithm foresaw the estimation of the series

y(k)=∑j=1k(RRj−RR¯)

k = 1, …, N. This series was divided into segments of equal length n and for each segment the linear approximation (least square fit, yn) was computed. Then root-mean-square fluctuation was calculated

F(n)=1N∑k=1N(y(k)−yn(k))2

Making a double log graph between log(F(n)) and different values of n, the slope of the regression line is the α scaling exponent. DFA1 and DFA2 features represent this slope between the ranges 4 ≤ n ≤ 16 and 16 ≤ n ≤ 64.

2.3.2.5. Multiscale Entropy and Multivariate Multiscale Entropy Analysis

Multiscale Entropy Analysis (MSE) is a powerful methodology based on the SampEn estimation. MSE was applied in several fields such as study of human gait dynamics (Costa et al., 2003), enhancement of postural complexity (Costa et al., 2007), and synthetic RR time series (Costa et al., 2002). MSE can be an effective non-linear method to collect information about physiological systems whose dynamics is associated to multiple different scales. This method is based on the application of sample entropy method to course-grained time series constructed from a one-dimensional discrete time series by averaging the data points within non-overlapping windows of increasing length, σ. Given a time series {x1, …, xi, …, xN} and a scale factor σ, each element of a course-grained series{y(σ)}is calculated using the equation

yj(σ)=1σ∑i=(j−1)σ+1jσxi,1≤j≤N/σ (6)

The length of each coarse-grained time series is equal to the length of the original time series divided by σ. The second step consists in the computation of SampEn (Richman and Moorman, 2000; Lake et al., 2002) algorithm on these series. Previous studies in which MSE algorithm was applied to physiological data use the standard value m = 2 for the pattern dimension (Costa et al., 2003; Leistedt et al., 2011). In this work the choice of the right r was performed by a method already used in the liter SampEn values were calculated for scale factors σ which were in a range from 1 to 20 and the same process was carried out on HRV, IBI, and IBR series. The complexity index (CI) was measured as the area under the curve of MSE graph and it can be calculated for short time scales, from 1 to 8 (short CI), and for higher time scales, up to 20 (long CI) (Leistedt et al., 2011).

Besides MSE analysis, we performed the Multivariate Multiscale Entropy (MMSE) (Ahmed and Mandic, 2011, 2012) analysis. This algorithm allows performing MSE analysis using multivariate time series. In this work, MMSE was used to quantify the complexity of the series derived from the electrocardiogram and breath. In particular, MMSE results were obtained on the bivariate series HRV-IBI, and HRV-IBR through the estimation of the CI indices (as described above on MSE). Before the MMSE calculation, the involved time series are scaled in the range between 0 and 1 to prevent that the different amplitudes may influence the complexity complexity (Ahmed and Mandic, 2011).

2.3.2.6. Symbolic Analysis

Symbolic analysis (Yeragani et al., 2000; Porta et al., 2001; Baumert et al., 2002; Bella and Montano, 2005; Tobaldini et al., 2009; Caminal et al., 2010) is another powerful non-linear method which was applied on HRV data series. For each HRV series gathered from each subject, 6 levels were constructed evenly dividing the amplitude range of the samples, and a symbol (from 0 to 5) was assigned to each data sample according to the level of belonging. Then, a window of three consecutive points moves along the HRV series, and three possible configurations are identified when running all the signal: the three points belong to the same level, i.e., no variation (0V), two consecutive points belong to the same level and one to another, i.e., one variation (1V), and the remaining cases, i.e., two variations (2V). The number of patterns falling into each group (0V, 1V, 2V) and the percentage of the total (0V%, 1V%, 2V%) were calculated and used as features. Previous studies support the hypothesis that an increase of 0V patterns is related to an activation of the sympathetic activity, an increase of 2V patterns is related to an increase of the parasympathetic activity, and increases of 1V patterns is associated to a simultaneous increase of both parasympathetic and sympathetic activities.

3. Experimental Results

Experimental results are expressed in terms of statistical and correlation analysis. In the literature it can be found the threshold score of each questionnaire above which the behavior of the subject results to show altered psycho-cognitive-behavioral traits. Among all the sub-scales we only considered those where the subjects spread out over a wide range of scores in order to identify two groups, one below and the other above the threshold. For each scale we identified two groups of subjects separated by the median. In order to have two groups numerically equivalent, we selected and investigated only these scales where the median was congruent with the threshold reported in the literature. In addition, for each of the 16 scales we verified that maximum and minimum scores of each group were in the tails of the population distribution reported by the literature. In other words, for each psychological subscale, the median value of the subjects score is calculated to identify two groups: one comprised of the subjects having scores below the median, and one comprised of the subjects having scores above the median. Only 16 out of 25 sub-scales divided the subjects in two groups numerically comparable, therefore we performed the statistical analysis on the scores obtained in these 16 sub-scales. The reference values from the literature about these sub-scales are evaluated on the control groups used in several previous works. For example we considered a sample of 103 subjects (age = 27.00 ± 8.80) for IRI Empathic Concern and IRI Personal Distress sub-scales, referring to a study which explored the relationship among psychological mindedness and several aspects of awareness which comprended this indices of empathy (Beitel et al., 2005) and a sample of 582 subjects for IRI fantasy sub-scale taking this data from a guide study on the empathy scales (Davis, 1980). For the two PANAS sub-scales, a group of 537 volunteers aged 18–91 was in a work that tried to evaluate the reliability and validity of the PANAS (Crawford and Henry, 2004), and 53 participants (age = 34.32 ± 10.50) were asked to answer to the LSAS questionnaires to demonstrate that this method may be employed in the assessment of social anxiety disorder (Fresco et al., 2001; Rytwinski et al., 2009). As a reference for the values of BIS and BAS sub-scales we chose a previous study where the answers of of 2725 individuals aged 18–79 were observed to validate the application of this scale to measure the behavioral inhibition and activation and its correlation with depression and anxiety (Jorm et al., 1998). The threshold value of the answers of a group of 639 participants in a study of the shortened form of the questionnaire, was taken in account for ZKPQ Impulsive Sensation Seeking and Activity sub-scales (age = 22.31 ± 5.08) (Aluja et al., 2003). At last, as regards DERS subscales, a study on 260 subjects in order to explore the factor structure and psychometric properties of DERS measures (age = 23.10 ± 5.67) was used as reference for DERS Awareness (Gratz and Roemer, 2004) and a reference sample of 42 individuals (age = 24.24 ± 4.38) was considered for the other DERS sub-scales, extracted from a research which compared the values of the this psychological tests on depressed patients and healthy subjects (Ehring et al., 2008).

In the statistical analysis, for each psychological sub-scale and for each ANS feature, we applied the Mann-Whitney test in order to evaluate whether the two groups were statistically different. Moreover, the non-parametric Spearman correlation coefficient was calculated between each psychological sub-scale and ANS feature.

3.1. Statistical Analysis

As mentioned above, for each ANS feature, Mann-Whitney non-parametric U-tests were used to test the null hypothesis of having no statistical difference between two groups. The use of such a non-parametric test is justified by having non-gaussian distribution of the samples (p < 0.05 of the null hypothesis of having gaussian samples of the Kolmogorov-Smirnov test).

Concerning features from HRV standard analysis, 8 sub-scales (LSAS Anxiety of Performance, DERS Non-Acceptance, DERS Awareness, IRI fantasy, IRI Empathic Concern, ZKPQ Activity, ZKPQ Impulsive Seeking Sensation, BAS) showed significant discerning capability mostly through frequency domain parameters (see details in Table 1). Concerning ANS features coming from non-linear analysis, 9 sub-scales (PANAS Positive Affect, DERS non-Acceptance, DERS Impulse, DERS Awareness, DERS Strategies, IRI Empathic Concern, BIS, BAS, ZKPQ Activity) showed significant differences considering monovariate and multivariate measures (see details in Table 2). An exemplary plot showing the discerning capability of MMSE analysis on DERS Non-Accept sub-scale is shown in Figure 1.

TABLE 1

Table 1. Statistical results related to standard HRV features (U-test).

TABLE 2

Table 2. Statistical results related to non-linear features (U-test).

FIGURE 1

Figure 1. Exemplary plot of Multivariate Multiscale Entropy analysis applied to HRV-IBI series in discerning the two groups (under the median-lower scores: group 1; over the median-higher scores: group 2) according to scores gathered from the DERS Non-Accept sub-scale.

To summarize the results, all the extracted features were able to discern the two groups in 12 out of 16 sub-scales. More specifically, standard HRV analysis provided exclusive information, i.e., not overlapped with that coming from the non-linear analysis, on the psychological assessment in only 2 sub-scales, whereas features from ANS non-linear dynamics exclusively discriminated the two groups in 4 sub-scales (see details in Figure 2).

FIGURE 2

Figure 2. Number of features with significant p-values (p < 0.05) given by the Mann-Whitney tests. For each phycological scale, the number of significant parameter of standard analysis, non-linear analysis, and the total are shown.

3.2. Correlation Analysis

The Spearman correlation coefficient was used to show the relationship between the values of each features through all the subjects and the relative score for each sub-scale. Accordingly, the coefficient ρ and p−value expressing the probability that no correlation between the two variables exist, were assigned for each sub-scale and each feature. Results are shown in Tables 3, 4.

TABLE 3

Table 3. Spearman correlation test results related to standard HRV features.

TABLE 4

Table 4. Spearman correlation test results related to non-linear HRV, IBI, IBR features.

We found that ANS features related to the linear HRV dynamics are significantly correlated with 5 sub-scales, reaching absolute values of ρ up to 0.52 (BAS and ZKPQ Impulsive Sensation Seeking). Moreover, 10 sub-scales are significantly correlated with markers of ANS non-linear dynamics, reaching absolute values of ρ up to 0.55 (DERS Non-Acceptance).

Although the correlation coefficient is not very high, it is, however, a very interesting result to be further validated and confirmed.

The number of features with significant p-values (p < 0.05) given by such a correlation coefficient is shown in Figure 3 for each phycological dimension.

FIGURE 3

Figure 3. Number of features with significant p-values (p < 0.05) given by the Spearman non-parametric correlation coefficient. For each phycological scale, the number of significant parameter of standard analysis, non-linear analysis, and the total are shown.

4. Discussion and Conclusion

In conclusion, we found several ANS biomarkers of psychological dimensions in non-pathological subjects. Such biomarkers are derived from the standard and complexity analysis of ANS measures such as HRV, IBI, and IBR series. We found that dimensions related to difficulties in emotion regulation (DERS), interpersonal reactivity (IRI), behavioral activation or inhibition (BIS/BAS), sensation-seeking and activity (ZKPQ), and anxiety performance (LSAS) are always associated to changes in the HRV dynamics, quantified using time and frequency domain indices (see Table 1). As all the scale define different psychological dimensions, it is very difficulty to give a common interpretation of features through them. The LF/HF ratio decrease, associated to increased questionnaires scores, characterizes the ZKPQ activity and IRI empathic concern, whereas an opposite trend is found for the awareness of difficulties in emotion regulation (DERS). HRV time domain indices such as TINN, Pnn50, and RMSSD are effective only to characterize the empathic concern and emotion regulation. These results, gathered from statistical analyses of standard HRV parameters, are further confirmed by the correlation analyses whose details are shown in Table 3.

It is worthwhile noting that the HF power decreases with the DERS score. According to the literature (Porges, 1991, 1992), vagal tone is associated to the ability of emotional self-regulation and high flexibility and adaptability to environmental changes. According to our results, when an emotion dysregulation occurs, the sympathetic activity increases.

Other evidences supporting our results can be found in the current literature (Freeman and Nixon, 1985; Yeragani et al., 1999; Virtanen et al., 2003; Cohen and Benjamin, 2006; Shinba et al., 2008; Licht et al., 2009; Thayer et al., 2010; 2012) which suggest that patients with anxiety disorders revealed a decreased power in the HRV-LF bands.

Concerning the ANS non-linear dynamics, several biomarkers of psychological dimensions were found in complexity measures such as sample entropy, monovariate and multivatiate multiscale entropy, short- an long-term correlations, correlation dimension, recurrence and symbolic analysis in characterizing dimensions as positive and negative affect (PANAS), social phobia (Liebowitz Social Anxiety Scale, LSAS), difficulties in emotion regulation (DERS), Interpersonal reactivity (IRI), behavioral inhibition or activation (BIS/BAS), and sensation-seeking and activity (ZKPQ). Our results on non-linear ANS markers for psychological dimensions confirm the previous findings (Yeragani et al., 2000; Cohen and Benjamin, 2006) and provide a wider portrait of the complexity modulation associated with behavioral characters.

Figures 2, 3 report the number of statistically significant features given by Mann-Whitney and Spearman non-parametric correlation, respectively. It is worthy to note that the non-linear features are overall more than those extracted from standard analysis, confirming that complexity dynamics measures play a relevant role in assessing the psycho-physiological dimensions.

Finally, some prudential considerations should be made. The patterns of physiological signals are acquired in rest conditions right after performing the test and the assumption behind the experiment is that the psychological assessment acted as an affective elicitation. Results have to be considered as preliminary to future experiments where subjects experience an actual affective dimension while they are monitored. Nevertheless, it is worthwhile pointing out that complexity measures can be considered promising markers to assess the psychological traits. Is important to underline that such interest is not diminished by the difficulty in giving a physiological meaning to complexity measurements. In this sense and more in generally, we underscored how our data suggested the possibility of an ANV fingerprinting of psychological dimensions. Therefore, beyond their precise physiological meaning, our results have interesting consequences for the psychometric and clinical fields. Our approach may be promising in describing the psychological dimensions as a combination of different features, providing a full classification of psychological characteristics through a baseline ECG acquisition. However, more studies with a much higher number of subjects are needed to test the reliability and the feasibility of these potential clinical implications. Furthermore, to test if our methodology could also be extended to the extremes of the psychological dimensions, these studies should also include pathological samples (e.g., diagnosed subjects). Should that prove to be the case, this approach might hold promise as a tool for providing an external validation to psychological diagnosis.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The research leading to these results has received partial funding from the European Union Seventh Framework Programme FP7/2007–2013 under grant agreement n 601165 of the project “WEARHAP.”

References

Acharya, U. R., Joseph, K. P., Kannathal, N., Lim, C. M., and Suri, J. S. (2006). Heart rate variability: a review. Med. Biol. Eng. Comput. 44, 1031–1051. doi: 10.1007/s11517-006-0119-0

Ahmed, M. U., and Mandic, D. P. (2011). Multivariate multiscale entropy: a tool for complexity analysis of multichannel data. Phys. Rev. E Stat. Nonlin. Soft. Matter Phys. 84:061918. doi: 10.1103/PhysRevE.84.061918

Ahmed, M. U., and Mandic, D. P. (2012). Multivariate multiscale entropy analysis. IEEE Signal Process. Lett. 19, 91–94. doi: 10.1109/LSP.2011.2180713

Aluja, A., Garcia, O., and Garcia, L. F. (2003). Psychometric properties of the zuckerman–kuhlman personality questionnaire (zkpq-iii-r): a study of a shortened form. Pers. Individ. Dif. 34, 1083–1097. doi: 10.1016/S0191-8869(02)00097-1

Atyabi, F., Livari, M., Kaviani, K., and Tabar, M. (2006). Two statistical methods for resolving healthy individuals and those with congestive heart failure based on extended self-similarity and a recursive method. J. Biol. Phys. 32, 489–495. doi: 10.1007/s10867-006-9031-y

Baumert, M., Walther, T., Hopfe, J., Stepan, H., Faber, R., and Voss, A. (2002). Joint symbolic dynamic analysis of beat-to-beat interactions of heart rate and systolic blood pressure in normal pregnancy. Med. Biol. Eng. Comput. 40, 241–245. doi: 10.1007/BF02348131

Beitel, M., Ferrer, E., and Cecero, J. J. (2005). Psychological mindedness and awareness of self and others. J. Clin. Psychol. 61, 739–750. doi: 10.1002/jclp.20095

Bland, J. M., and Altman, D. G. (1997). Statistics notes: Cronbach's alpha. BMJ 314:572. doi: 10.1136/bmj.314.7080.572

Calvo, R. A., and D'Mello, S. (2010). Affect detection: an interdisciplinary review of models, methods, and their applications. IEEE Trans. Affect. Comput. 1, 18–37. doi: 10.1109/T-AFFC.2010.1

Caminal, P., Giraldo, B., Vallverdú, M., Benito, S., Schroeder, R., and Voss, A. (2010). Symbolic dynamic analysis of relations between cardiac and breathing cycles in patients on weaning trials. Ann. Biomed. Eng. 38, 2542–2552. doi: 10.1007/s10439-010-0027-1

Camm, A., Malik, M., Bigger, J., Breithardt, G., Cerutti, S., Cohen, R., et al. (1996). Heart rate variability: standards of measurement, physiological interpretation, and clinical use. Circulation 93, 1043–1065. doi: 10.1161/01.CIR.93.5.1043

Carney, R., Freedland, K., Rich, M., and Jaffe, A. (1995). Depression as a risk factor for cardiac events in established coronary heart disease: a review of possible mechanisms. Ann. Behav. Med. 17, 142–149. doi: 10.1007/BF02895063

Carney, R. M., Freedland, K. E., and Veith, R. C. (2005). Depression, the autonomic nervous system, and coronary heart disease. Psychosom. Med. 67(Suppl. 1), S29–S33. doi: 10.1097/01.psy.0000162254.61556.d5

Carver, C. S., and White, T. L. (1994). Behavioral inhibition, behavioral activation, and affective responses to impending reward and punishment: the bis/bas scales. J. Pers. Soc. Psychol. 67:319. doi: 10.1037/0022-3514.67.2.319

Casdagli, M., Eubank, S., Farmer, J. D., and Gibson, J. (1991). State space reconstruction in the presence of noise. Phys. D 51, 52–98. doi: 10.1016/0167-2789(91)90222-U

Citi, L., Valenza, G., and Barbieri, R. (2012). “Instantaneous estimation of high-order nonlinear heartbeat dynamics by lyapunov exponents,” in Engineering in Medicine and Biology Society (EMBC), 2012 Annual International Conference of the IEEE (San Diego, CA: IEEE), 13–16.

Cohen, H., and Benjamin, J. (2006). Power spectrum analysis and cardiovascular morbidity in anxiety disorders. Auton. Neurosci. 128, 1–8. doi: 10.1016/j.autneu.2005.06.007

Cohen, R., Swerdlik, M., and Smith, D. (1992). Psychological Testing and Assessment: An Introduction to Tests and Measurement. Houston, TX: Mayfield Publishing Co.

Costa, M., Goldberger, A., and Peng, C. (2002). Multiscale entropy analysis of complex physiologic time series. Phys. Rev. Lett. 89:68102. doi: 10.1103/PhysRevLett.89.068102

Costa, M., Peng, C.-K., Goldberger, A.L., and Hausdorff, J. M. (2003). Multiscale entropy analysis of human gait dynamics. Phys. A 330, 53–60. doi: 10.1016/j.physa.2003.08.022

Costa, M., Priplata, A., Lipsitz, L., Wu, Z., Huang, N., Goldberger, A., et al. (2007). Noise and poise: enhancement of postural complexity in the elderly with a stochastic-resonance–based therapy. Europhys. Lett. 77:68008. doi: 10.1209/0295-5075/77/68008

Crawford, J. R., and Henry, J. D. (2004). The positive and negative affect schedule (panas): construct validity, measurement properties and normative data in a large non-clinical sample. Br. J. Clin. Psychol. 43, 245–265. doi: 10.1348/0144665031752934

Davis, M. H. (1980). A multidimensional approach to individual differences in empathy. JSAS Catal. Sel. Doc. Psychol. 10:85.

Ehring, T., Fischer, S., Schnülle, J., Bösterling, A., and Tuschen-Caffier, B. (2008). Characteristics of emotion regulation in recovered depressed versus never depressed individuals. Pers. Individ. Dif. 44, 1574–1584. doi: 10.1016/j.paid.2008.01.013

Fernandez, A., Hornero, R., Gomez, C., Turrero, A., Gil-Gregorio, P., Matías-Santos, J., et al. (2010). Complexity analysis of spontaneous brain activity in alzheimer disease and mild cognitive impairment: an meg study. Alzheimer Dis. Assoc. Disord. 24, 182–189. doi: 10.1097/WAD.0b013e3181c727f7

Freeman, L. J., and Nixon, P. (1985). Are coronary artery spasm and progressive damage to the heart associated with the hyperventilation syndrome? Br. Med. J. (Clin. Res. Ed.) 291:851. doi: 10.1136/bmj.291.6499.851

Fresco, D. M., Coles, M. E., Heimberg, R. G., Liebowitz, M. R., Hami, S., and Stein, M. B. (2001). The liebowitz social anxiety scale: a comparison of the psychometric properties of self-report and clinician-administered formats. Psychol. Med. 31, 1025–1035. doi: 10.1017/S0033291701004056

Fusheng, Y., Bo, H., and Qingyu, T. (2001). “Approximate entropy and its application to biosignal analysis,” in Nonlinear Biomedical Signal Processing: Dynamic Analysis and Modeling, Vol. 2, ed M. Akay (Hoboken, NJ: JohnWiley & Sons Inc), 72–91.

Gao, J., Gurbaxani, B. M., Hu, J., Heilman, K. J., Emanuele, II. V. A., Lewis, G. F., et al. (2013). Multiscale analysis of heart rate variability in non-stationary environments. Front. Physiol. 4:119. doi: 10.3389/fphys.2013.00119

Gao, J., Hu, J., and Tung, W.-w. (2011). Complexity measures of brain wave dynamics. Cogn. Neurodynam. 5, 171–182. doi: 10.1007/s11571-011-9151-3

Glass, L. (2001). Synchronization and rhythmic processes in physiology. Nature 410, 277–284. doi: 10.1038/35065745

Glass, L. (2009). Introduction to controversial topics in nonlinear science: is the normal heart rate chaotic? Chaos 19, 028501. doi: 10.1063/1.3156832

Glassman, A. (1998). Depression, cardiac death, and the central nervous system. Neuropsychobiology 37, 80–83. doi: 10.1159/000026482

Goldberger, A., Peng, C., and Lipsitz, L. (2002). What is physiologic complexity and how does it change with aging and disease? Neurobiol. Aging 23, 23–26. doi: 10.1016/S0197-4580(01)00266-4

Grassberger, P. (1988). Finite sample corrections to entropy and dimension estimates. Phys. Lett. A 128, 369–373. doi: 10.1016/0375-9601(88)90193-4

Grassberger, P., and Procaccia, I. (1983). Estimation of the kolmogorov entropy from a chaotic signal. Phys. Rev. A 28, 2591–2593. doi: 10.1103/PhysRevA.28.2591

Gratz, K. L., and Roemer, L. (2004). Multidimensional assessment of emotion regulation and dysregulation: development, factor structure, and initial validation of the difficulties in emotion regulation scale. J. Psychopathol. Behav. Assess. 26, 41–54. doi: 10.1023/B:JOBA.0000007455.08539.94

Grossman, P., Wilhelm, F. H., Kawachi, I., and Sparrow, D. (2001). Gender differences in psychophysiological responses to speech stress among older social phobics congruence and incongruence between self-evaluative and cardiovascular reactions. Psychosom. Med. 63, 765–777. doi: 10.1097/00006842-200109000-00010

Groth-Marnat, G. (2003). Handbook of Psychological Assessment. Hoboken; New York: John Wiley & Sons Inc.

Guzzetti, S., Borroni, E., Garbelli, P. E., Ceriani, E., Della Bella, P., Montano, N., et al. (2005). Symbolic dynamics of heart rate variability a probe to investigate cardiac autonomic modulation. Circulation 112, 465–470. doi: 10.1161/CIRCULATIONAHA.104.518449

Heimberg, R. G., Horner, K., Juster, H., Safren, S., Brown, E., Schneier, F., et al. (1999). Psychometric properties of the liebowitz social anxiety scale. Psychol. Med. 29, 199–212. doi: 10.1017/S0033291798007879

Heller, A. S., Johnstone, T., Shackman, A. J., Light, S. N., Peterson, M. J., Kolden, G. G., et al. (2009). Reduced capacity to sustain positive emotion in major depression reflects diminished maintenance of fronto-striatal brain activation. Proc. Natl. Acad. Sci. U.S.A. 106, 22445–22450. doi: 10.1073/pnas.0910651106

Hu, J., Gao, J., and Príncipe, J. C. (2006). Analysis of biomedical signals by the lempel-ziv complexity: the effect of finite data size. IEEE Trans. Biomed. Eng. 53, 2606–2609. doi: 10.1109/TBME.2006.883825

Hu, J., Gao, J., and Tung, W.-W. (2009). Characterizing heart rate variability by scale-dependent lyapunov exponent. Chaos 19, 028506. doi: 10.1063/1.3152007

Hu, J., Gao, J., Tung, W.-W., and Cao, Y. (2010). Multiscale analysis of heart rate variability: a comparison of different complexity measures. Ann. Biomed. Eng. 38, 854–864. doi: 10.1007/s10439-009-9863-2

Hunsley, J., and Mash, E. J. (2010). “The role of assessment in evidence-based practice,” in Handbook of Assessment and Treatment Planning for Psychological Disorders, 2nd Edn., eds M. M. Antony and D. H. Barlow (New York, NY: Guilford Press), 3–22.

Iverson, G., Gaetz, M., Rzempoluck, E., McLean, P., Linden, W., and Remick, R. (2005). A new potential marker for abnormal cardiac physiology in depression. J. Behav. Med. 28, 507–511. doi: 10.1007/s10439-009-9863-2

Iverson, G., Stampfer, H., and Gaetz, M. (2002). Reliability of circadian heart pattern analysis in psychiatry. Psychiatr. Q. 73, 195–203. doi: 10.1023/A:1016036704524

Jansen, L. M., Gispen-de Wied, C. C., Wiegant, V. M., Westenberg, H. G., Lahuis, B. E., and van Engeland, H. (2006). Autonomic and neuroendocrine responses to a psychosocial stressor in adults with autistic spectrum disorder. J. Autism Dev. Disord. 36, 891–899. doi: 10.1007/s10803-006-0124-z

Jorm, A. F., Christensen, H., Henderson, A. S., Jacomb, P. A., Korten, A. E., and Rodgers, B. (1998). Using the bis/bas scales to measure behavioural inhibition and behavioural activation: factor structure, validity and norms in a large community sample. Pers. Individ. Dif. 26, 49–58. doi: 10.1016/S0191-8869(98)00143-3

Kenny, M. C., Alvarez, K., Donohue, B. C., and Winick, C. B. (2008). “Overview of behavioral assessment with adults,” in Handbook of Psychological Assessment, Case Conceptualization, and Treatment, Adults, Vol. 1, eds M. Hersen and J. Rosqvist (Hoboken, NJ: John Wiley & Sons Inc).

Lake, D. E., Richman, J. S., Griffin, M. P., and Moorman, J. R. (2002). Sample entropy analysis of neonatal heart rate variability. Am. J. Physiol. Regul. Integr. Comp. Physiol. 283, R789–R797. doi: 10.1152/ajpregu.00069.2002

Lanatà, A., Valenza, G., and Scilingo, E. P. (2012). A novel eda glove based on textile-integrated electrodes for affective computing. Med. Biol. Eng. Comput. 50, 1163–1172. doi: 10.1007/s11517-012-0921-9

Leistedt, S. J., Linkowski, P., Lanquart, J., Mietus, J., Davis, R. B., Goldberger, A. L., et al. (2011). Decreased neuroautonomic complexity in men during an acute major depressive episode: analysis of heart rate dynamics. Transl. Psychiatry 1, e27. doi: 10.1038/tp.2011.23

Licht, C. M., de Geus, E. J., van Dyck, R., and Penninx, B. W. (2009). Association between anxiety disorders and heart rate variability in the netherlands study of depression and anxiety (nesda). Psychosom. Med. 71, 508–518. doi: 10.1097/PSY.0b013e3181a292a6

Liebowitz, M. R. (1987). Social phobia. Mod. Probl. Pharmacopsychiatry 22, 141–173.

Lin, Y., Wang, C., Jung, T., Wu, T., Jeng, S., Duann, J., et al. (2010). EEG-based emotion recognition in music listening. IEEE Trans. Biomed. Eng. 57, 1798–1806. doi: 10.1109/TBME.2010.2048568

Marmarelis, V. Z. (2004). Nonlinear Dynamic Modeling of Physiological Systems. New York, NY: Wiley-Interscience.

Marwan, N., Carmen Romano, M., Thiel, M., and Kurths, J. (2007). Recurrence plots for the analysis of complex systems. Phys. Rep. 438, 237–329. doi: 10.1016/j.physrep.2006.11.001

Marwan, N., Wessel, N., Meyerfeldt, U., Schirdewan, A., and Kurths, J. (2002). Recurrence-plot-based measures of complexity and their application to heart-rate-variability data. Phys. Rev. E 66:026702. doi: 10.1103/PhysRevE.66.026702

Mujica-Parodi, L., Yeragani, V., and Malaspina, D. (2005). Nonlinear complexity and spectral analyses of heart rate variability in medicated and unmedicated patients with schizophrenia1. Neuropsychobiology 51, 10–15. doi: 10.1159/000082850

Peng, C., Havlin, S., Stanley, H., and Goldberger, A. (1995). Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. Chaos 5:82. doi: 10.1063/1.166141

Penzel, T., Kantelhardt, J. W., Grote, L., Peter, J.-H., and Bunde, A. (2003). Comparison of detrended fluctuation analysis and spectral analysis for heart rate variability in sleep and sleep apnea. IEEE Trans. Biomed. Eng. 50, 1143–1151. doi: 10.1109/TBME.2003.817636

Petrantonakis, P. C., and Hadjileontiadis, L. J. (2011). A novel emotion elicitation index using frontal brain asymmetry for enhanced eeg-based emotion recognition. IEEE Trans. Inf. Technol. Biomed. 15, 737–746. doi: 10.1109/TITB.2011.2157933

Pincus, S. (1991). Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. U.S.A. 88:2297. doi: 10.1073/pnas.88.6.2297

Poon, C., and Merrill, C. (1997). Decrease of cardiac chaos in congestive heart failure. Nature 389, 492–495. doi: 10.1038/39043

Porges, S. W. (1991). Vagal Tone: An Autonomic Mediator of Affect. Cambridge, UK: Cambridge University Press.

Porges, S. W. (1992). “Autonomic regulation and attention,” in Attention and Information Processing in Infants and Adults: Perspectives from Human and Animal Research, eds B. A. Campbell, H. Hayne, and R. Richardson (Hillsdale, NJ: Erlbaum), 201–223.

Porta, A., Guzzetti, S., Montano, N., Furlan, R., Pagani, M., Malliani, A., et al. (2001). Entropy, entropy rate, and pattern classification as tools to typify complexity in short heart period variability series. IEEE Trans. Biomed. Eng. 48, 1282–1291. doi: 10.1109/10.959324

Richman, J., and Moorman, J. (2000). Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 278, H2039–H2049.

Ruiz-Padial, E., Vila, J., and Thayer, J. F. (2011). The effect of conscious and non-conscious presentation of biologically relevant emotion pictures on emotion modulated startle and phasic heart rate. Int. J. Psychophysiol 79, 341–346. doi: 10.1016/j.ijpsycho.2010.12.001

Rytwinski, N. K., Fresco, D. M., Heimberg, R. G., Coles, M. E., Liebowitz, M. R., Cissell, S., et al. (2009). Screening for social anxiety disorder with the self-report version of the liebowitz social anxiety scale. Depress. Anxiety 26, 34–38. doi: 10.1002/da.20503

Sârbescu, P., and Neguţ, A. (2012). Psychometric properties of the romanian version of the Zuckerman-Kuhlman personality questionnaire. Eur. J. Psychol. Assess. 29, 241–252. doi: 10.1027/1015-5759/a000152

Sava, F. A., and Sperneac, A.-M. (2006). Sensitivity to reward and sensitivity to punishment rating scales: a validation study on the romanian population. Pers. Individ. Dif. 41, 1445–1456. doi: 10.1016/j.paid.2006.04.024

Schouten, J. C., Takens, F., and van den Bleek, C. M. (1994). Estimation of the dimension of a noisy attractor. Phys. Rev. E 50:1851. doi: 10.1103/PhysRevE.50.1851

Shinba, T., Kariya, N., Matsui, Y., Ozawa, N., Matsuda, Y., and Yamamoto, K.-I. (2008). Decrease in heart rate variability response to task is related to anxiety and depressiveness in normal subjects. Psychiatry Clin. Neurosci. 62, 603–609. doi: 10.1111/j.1440-1819.2008.01855.x

Stampfer, H. (1998). The relationship between psychiatric illness and the circadian pattern of heart rate. Aust. Psychiatry 32, 187–198. doi: 10.3109/00048679809062728

Stiedl, O., and Meyer, M. (2003). Fractal dynamics in circadian cardiac time series of corticotropin-releasing factor receptor subtype-2 deficient mice. J. Math. Biol 47, 169–197. doi: 10.1007/s00285-003-0197-7

Taillard, J., Lemoine, P., Boule, P., Drogue, M., and Mouret, J. (1993). Sleep and heart rate circadian rhythm in depression: the necessity to separate. Chronobiol. Int. 10, 63–72. doi: 10.3109/07420529309064483

Taillard, J., Sanchez, P., Lemoine, P., and Mouret, J. (1990). Heart rate circadian rhythm as a biological marker of desynchronization in major depression: a methodological and preliminary report. Chronobiol. Int. 7, 305–316. doi: 10.3109/07420529009064636

Takahashi, T., Cho, R. Y., Mizuno, T., Kikuchi, M., Murata, T., Takahashi, K., et al. (2010). Antipsychotics reverse abnormal eeg complexity in drug-naive schizophrenia: a multiscale entropy analysis. Neuroimage 51, 173–182. doi: 10.1016/j.neuroimage.2010.02.009

Takens, F. (1981). “Detecting strange attractors in turbulence,” in Dynamical Systems and Turbulence, Warwick 1980, eds D. A. Rand and L.-S. Young (Warwick: Springer), 366–381.

Thayer, J. F., Åhs, F., Fredrikson, M., Sollers, J. J. III., and Wager, T. D. (2012). A meta-analysis of heart rate variability and neuroimaging studies: implications for heart rate variability as a marker of stress and health. Neurosci. Biobehav. Rev. 36, 747–756. doi: 10.1016/j.neubiorev.2011.11.009

Thayer, J. F., Friedman, B. H., and Borkovec, T. D. (1996). Autonomic characteristics of generalized anxiety disorder and worry. Biol. Psychiatry 39, 255–266. doi: 10.1016/0006-3223(95)00136-0

Thayer, J. F., Yamamoto, S. S., and Brosschot, J. F. (2010). The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int. J. Cardiol. 141, 122–131. doi: 10.1016/j.ijcard.2009.09.543

Theiler, J. (1987). Efficient algorithm for estimating the correlation dimension from a set of discrete points. Phys. Rev. A 36:4456. doi: 10.1103/PhysRevA.36.4456

Tobaldini, E., Porta, A., Wei, S.-G., Zhang, Z.-H., Francis, J., Casali, K. R., et al. (2009). Symbolic analysis detects alterations of cardiac autonomic modulation in congestive heart failure rats. Auton. Neurosci. 150, 21–26. doi: 10.1016/j.autneu.2009.03.009

Tulppo, M., Kiviniemi, A., Hautala, A., Kallio, M., Seppanen, T., Makikallio, T., et al. (2005). Physiological background of the loss of fractal heart rate dynamics. Circulation 112, 314. doi: 10.1161/CIRCULATIONAHA.104.523712

Valenza, G., Allegrini, P., Lanatà, A., and Scilingo, E. P. (2012a). Dominant lyapunov exponent and approximate entropy in heart rate variability during emotional visual elicitation. Front. Neuroeng. 5:3. doi: 10.3389/fneng.2012.00003

Valenza, G., Citi, L., Gentili, C., Lanatá, A., Scilingo, E., and Barbieri, R. (2014a). Point-process nonlinear autonomic assessment of depressive states in bipolar patients. Methods Inf. Med. 53, 296–302. doi: 10.3414/ME13-02-0036

Valenza, G., Gentili, C., Lanatà, A., and Scilingo, E. P. (2013a). Mood recognition in bipolar patients through the psyche platform: preliminary evaluations and perspectives. Artif. Intell. Med. 57, 49–58. doi: 10.1016/j.artmed.2012.12.001

Valenza, G., Lanatà, A., Ferro, M., and Scilingo, E. P. (2008). “Real-time discrimination of multiple cardiac arrhythmias for wearable systems based on neural networks,” in Computers in Cardiology, 2008 (Bologna: IEEE), 1053–1056.

Valenza, G., Lanata, A., and Scilingo, E. P. (2012b). The role of nonlinear dynamics in affective valence and arousal recognition. IEEE Trans. Affect. Comput. 3, 237–249. doi: 10.1109/T-AFFC.2011.30

Valenza, G., Lanatá, A., and Scilingo, E. P. (2013b). Improving emotion recognition systems by embedding cardiorespiratory coupling. Physiol. Meas. 34:449. doi: 10.1088/0967-3334/34/4/449

Valenza, G., Nardelli, M., Bertschy, G., Lanata, A., and Scilingo, E. (2014b). Mood states modulate complexity in heartbeat dynamics: a multiscale entropy analysis. Europhys. Lett. 107:18003. doi: 10.1209/0295-5075/107/18003

Valenza, G., Nardelli, M., Lanata, A., Gentili, C., Bertschy, G., Paradiso, R., et al. (2014c). Wearable monitoring for mood recognition in bipolar disorder based on history-dependent long-term heart rate variability analysis. IEEE J. Biomed. Health Inform. 18, 1625–1635. doi: 10.1109/JBHI.2013.2290382

Virtanen, R., Jula, A., Salminen, J. K., Voipio-Pulkki, L.-M., Helenius, H., Kuusela, T., et al. (2003). Anxiety and hostility are associated with reduced baroreflex sensitivity and increased beat-to-beat blood pressure variability. Psychosom. Med. 65, 751–756. doi: 10.1097/01.PSY.0000088760.65046.CF

Watkins, L., Blumenthal, J., and Carney, R. (2002). Association of anxiety with reduced baroreflex cardiac control in patients after acute myocardial infarction. Am. Heart J. 143, 460–466. doi: 10.1067/mhj.2002.120404

Watson, D., Clark, L. A., and Tellegen, A. (1988). Development and validation of brief measures of positive and negative affect: the panas scales. J. Pers. Soc. Psychol. 54:1063. doi: 10.1037/0022-3514.54.6.1063

Wu, G., Arzeno, N., Shen, L., Tang, D., Zheng, D., Zhao, N., et al. (2009). Chaotic signatures of heart rate variability and its power spectrum in health, aging and heart failure. PLoS ONE 4:e4323. doi: 10.1371/journal.pone.0004323

Yang, A. C., and Tsai, S.-J. (2012). Is mental illness complex? from behavior to brain. Prog. Neuropsychopharmacol. Biol. Psychiatry 45, 253–257. doi: 10.1016/j.pnpbp.2012.09.015

Yeragani, V. K., Jampala, V., Sobelewski, E., Kay, J., and Igel, G. (1999). Effects of paroxetine on heart period variability in patients with panic disorder: a study of holter ecg records. Neuropsychobiology 40, 124–128. doi: 10.1159/000026608

Yeragani, V. K., Nadella, R., Hinze, B., Yeragani, S., and Jampala, V. (2000). Nonlinear measures of heart period variability: decreased measures of symbolic dynamics in patients with panic disorder. Depress. Anxiety 12, 67–77. doi: 10.1002/1520-6394(2000)12:2<67::AID-DA2>3.0.CO;2-C

Zbilut, J. P., Koebbe, M., Loeb, H., and Mayer-Kress, G. (1990). “Use of recurrence plots in the analysis of heart beat intervals,” in Proceedings of Computers in Cardiology 1990 (Chicago, IL: IEEE), 263–266.

Zuckerman, M., Kuhlman, D. M., Joireman, J., Teta, P., and Kraft, M. (1993). A comparison of three structural models for personality: the big three, the big five, and the alternative five. J. Pers. Soc. Psychol. 65:757. doi: 10.1037/0022-3514.65.4.757