Next Article in Journal
Multichannel Classifier for Recognizing Acoustic Impacts Recorded with a phi-OTDR
Previous Article in Journal
A Quantitative Detection Algorithm for Multi-Test Line Lateral Flow Immunoassay Applied in Smartphones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Autonomic Regulation of Facial Temperature during Stress: A Cross-Mapping Analysis

by
Federica Gioia
1,2,
Mimma Nardelli
1,2,
Enzo Pasquale Scilingo
1,2 and
Alberto Greco
1,2,*
1
Dipartimento di Ingegneria dell’Informazione, University of Pisa, 56122 Pisa, Italy
2
Research Center “E. Piaggio”, University of Pisa, Largo Lucio Lazzarino, 1, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(14), 6403; https://doi.org/10.3390/s23146403
Submission received: 25 May 2023 / Revised: 6 July 2023 / Accepted: 12 July 2023 / Published: 14 July 2023
(This article belongs to the Section Biomedical Sensors)

Abstract

:
Skin temperature reflects the Autonomic Nervous System (ANS)’s response to emotions and mental states and can be remotely measured using InfraRed Thermography. Understanding the physiological mechanisms that affect facial temperature is essential to improve the precision of emotional inference from thermal imaging. To achieve this aim, we recorded thermal images from 30 volunteers, at rest and under acute stress induced by the Stroop test, together with two autonomic correlates, i.e., heart rate variability and electrodermal activity, the former serving as a measure of cardiovascular dynamics, and the latter of the activity of the sweat glands. We used a Cross Mapping (CM) approach to quantify the nonlinear coupling of the temperature from four facial regions with the ANS correlates. CM reveals that facial temperature has a statistically significant correlation with the two autonomic time series, under both conditions, which was not evident in the linear domain. In particular, compared to the other regions, the nose shows a significantly higher link to the electrodermal activity in both conditions, and to the heart rate variability under stress. Moreover, the cardiovascular activity seems to be primarily responsible for the well-known decrease in nose temperature, and its coupling with the thermal signals significantly varies with gender.

1. Introduction

InfraRed Thermal Imaging (IRT) is taking hold in many sectors as a powerful contactless monitoring tool. Thermal cameras convert the IR energy emitted by any object at a higher temperature than absolute zero to a thermal image, also called a thermogram. Hence, through IRT we can visualize, measure and monitor an entire scene in thermal. In particular, in the psychophysiological field, thermal imaging has captured the attention of researchers since it could potentially reveal the subjects’ emotional state [1].
Numerous studies have investigated the level of autonomic emotion specificity by examining the various trends in peripheral physiological reactions to a range of emotional stimuli [2,3,4]. However, researchers seek a non-intrusive and contactless monitoring device to replace the commonly used wearable systems, which rely on electrodes positioned on the skin of the subject and may include cables connecting them to the acquisition unit. In fact, these tools are likely to interfere with the subject, compromising the natural physiological responses under study.
The body’s thermal response to affective stimuli is likely a dynamic process influenced by various physiological and psychological factors. This complexity presents a challenge in the accurate interpretation of thermal imaging data in psychophysiological research. The autonomic nervous system (ANS), which regulates the visceral responses to a change in mental state, seems to be responsible for the thermal modulation of the skin. Indeed, ANS regulates sweat production and influences the vasomotor activity of the subcutaneous vessels and consequently the dynamics of the blood circulation which produces the temperature variation [5,6,7,8]. A higher blood flux would cause a rise in local skin temperature, in contrast to perspiration which would cause a drop in temperature, due to the production of cold sweat. Moreover, respiration, muscle tension, and metabolic processes may also play a role in thermal regulation [6,9,10].
Therefore, thermal imaging has the potential to provide valuable insights into the physiological mechanisms underlying psychological states and behaviours. However, there is limited knowledge about how the interplay of these phenomena leads to the final temperature. Modelling these processes could ease the evaluation of experimental results, leading to higher specificity and interpretability. To achieve this aim, it is important to consider the nonlinear nature of the vast majority of biological systems including the ANS processes [11]. In fact, the analysis of nonlinear dynamics of ANS correlates, such as heart rate variability (HRV) and electrodermal activity (EDA), has shed light on the regulation of the organism’s homeostasis by the ANS in the face of numerous internal and external stimuli [12,13,14]. Furthermore, previous studies have shown that many of the biological processes related to human thermoregulation also show chaotic behaviour [15,16,17].
Another important aspect to take into consideration when investigating thermoregulation is that thermal dynamics differ based on gender. In particular, differences in thermoregulation have been associated with subcutaneous fat content, menstrual cycles, and sweat evaporation efficiency [18,19]. According to Christensen et al., males have a higher average facial skin temperature than females due to differences in blood circulation in the skin or the higher basic metabolic rate in males [20].
In our study, we aim to provide insights into the ANS nonlinear processes underlying facial skin temperature regulation in response to emotional stimuli. Specifically, we analyze how cardiovascular and sweat gland dynamics, in terms of HRV and EDA, influence thermal modulation, focusing also on gender differences. Indeed, HRV dynamics reflects the interplay between the sympathetic and parasympathetic nervous systems, EDA is the result of the sympathetic regulation only, and together they could quantify the two ANS sub-system activities. To quantify the nonlinear correlations between IRT and both HRV and EDA time series, we used a cross-mapping approach under the assumption that two time series generated by the same underlying system should have comparable dynamical characteristics and behavioural patterns [21,22,23]. In fact, crossmapping provides a powerful tool for exploring the similarity between the two time series and determining whether they are produced by the same dynamical system. Particularly, we applied this analysis to thermal signals extracted from four facial regions recorded from healthy subjects during rest and under a stressful stimulus, to address the question both under parasympathetic and sympathetic nervous system predominance experimental conditions, respectively.

2. Materials and Methods

2.1. Study Population

We recruited 30 healthy Caucasian volunteers (20 females, age = 26.6 ± 3.6 years) to participate in this study. We specifically chose individuals whose body mass index (BMI) was in the healthy range and who had no history of mental illnesses. While the BMI was controlled to reduce any potential effects of fat tissue on the thermal behavior of the skin, mental problems can possibly affect a variety of physiological and psychological aspects. The experiments were approved by the “Bioethics Committee of the University of Pisa” (n. 15/2019), and each volunteer signed the informed consent. All subjects were asked not to drink coffee or smoke on the same day of the experiment to avoid the potential vasomotor effects of these substances. Finally, the day of the experiment, subjects were not exposed to extrinsic factors affecting skin temperature (e.g., massages, epilation or physical activity).

2.2. Experimental Protocol

Before the experiment, the subject sat quietly in the experimental room for at least 10 min to acclimatize to the environmental setting. The experimental timeline is shown in Figure 1. The experiment consisted of a 5-min resting session followed by a 5-min stressor task.
The instructions were presented on a laptop as the experimenter left the room to avoid interactions with the participant and equalize conditions for all the volunteers. During the resting session, the subject was required to rest with eyes open. Then, the affective elicitation was performed using a computerized-paced Stroop Test, aimed to induce mental stress [24]. This is a well-known task that requires the resolution of incongruous color/semantic-meaning stimuli. Here, the subject was required to press a button corresponding to the tint of the word presented in the middle of the screen. The stimuli were paced every 2 s. At the top of the screen, a counter kept track of the number of successes, as a motivational stressor. At the start and the end of the experiment, the subject was required to fill in the State-Trait Anxiety Inventory (STAY-Y1), a psychological inventory consisting of 20 self-report items on a four-point Likert scale to measure state anxiety. Before and after the stressor, a self-assessed perceived stress (PS) level was reported on a Likert scale from 0 (not at all) to 10 (very stressed), and it served as a reference for our analysis.
Environmental conditions were kept constant during each experimental session using a thermostat. In particular, room atmospheric temperature and humidity ranged around 24–26 ° C and 40–60%, respectively. Moreover, within the same experimental session, the temperature change was never over 1 ° C and humidity variation was a maximum of 5%. All recordings were obtained in the same room, and the participant was never exposed to any direct light sources or ventilation. A chinrest guaranteed that the camera was always perpendicular to the face of the subject and at a fixed distance of 1 m. The camera’s emissivity was set to 0.98, the recommended value for a human body surface. Twenty minutes before the start of the recording, the camera was turned on to allow for sensor stabilization. All the points of the checklist described in [25] were addressed to guarantee reliable skin temperature recordings with infrared thermography.

2.3. Equipment and Data Acquisition

Facial thermograms were acquired using a FLIR T640 thermal camera with a 24.6 mm lens, 640 × 480 pixels, NETD < 0.04 mK @ +30 ° and spectral range of 7.8–14 μ m (Long-wave-infrared, LWIR). The sampling frequency was 5 Hz. In addition, RGB images in the visible spectrum were recorded using a Logitech HD WebCam C270. Physiological signals (i.e., electrodermal and cardiac activity) were acquired using a BIOPAC MP150 system, at a sampling rate of 250 Hz. Specifically, EDA electrodes were placed on the distal phalanx of the index and ring fingers of the non-dominant hand. Electrocardiogram (ECG) was recorded using three electrodes, respectively below the right and left clavicle, and on the lower left chest.

2.4. Data Pre-Processing

We used thermal imaging to monitor the temperature of several regions of interest (ROIs) of the face (i.e., nose, forehead, and cheeks), as shown in Figure 2 [26,27].
Firstly, IR images were co-registered using a motion correction framework based on optical flow [5]. We used RGB images to localize the position of the ROIs in the first frame of each IR recording. To accomplish this aim, RGB and IR videos were synchronized in time. In addition, the first frames from the two techniques were spatially aligned using fiducial marker-based registration with an affine-2D transformation matrix. Then, we identified 68 facial landmarks in the RGB domain using the Yuval Nirkin algorithm [28]. Based on the position of these landmarks, the centers of the ROIs were automatically detected. The algorithm to locate the ROIs was developed in our previous study [27], and it is based on the intersection of lines connecting specific landmarks. In order to adapt the ROIs to each face shape and size, the dimension of each ROI was calculated to be proportional to the number of face pixels. Table 1 presents the mean and standard deviation of the number of pixels per ROI across all subjects in our study.
Accordingly, the face was separated from the background using a landmark-intensity-based segmentation to count the number of face pixels. Finally, the coordinates of the ROIs in the first RGB frame were referred to the co-registered IR images, in order to obtain the thermal time-series of each region. Each point of the thermal time-series is the median temperature value of the pixels within a ROI of a specific frame. The temperature was quite uniform within each ROI, as indicated by the average Median Absolute Deviation (MAD) of the temperature within each ROI, across frames and subjects, reported in Table 2. As a result, we focused on the median value because it provides a reliable approximation of the temperature inside the ROI and is unaffected by artifacts.
For the analysis, we selected 150-second-long segments of the thermal time-series. In particular, the last half of the resting session was adopted as a baseline, while the first half of the stressor task was chosen as the stimulation phase. On the one hand, this choice insured that during the baseline the subject had completely recovered from the distress caused by the setup process. On the other, it minimized the habituation effect on the stimulation phase. Hereinafter, we will refer to these segments as Rest and Stress, respectively. Finally, we extracted the median temperature during Rest and Stress.
ECG signals were processed using Kubios HRV software [29] to obtain the Heart Rate Variability (HRV) series. The interbeat (RR) signals were obtained using the Pan–Tompkins algorithm. In particular, a pre-processing phase was followed by decision rules to detect the QRS complex. Firstly, the ECG was bandpass filtered to eliminate power line noise and baseline drift. Then, the data samples were squared to highlight peaks, and moving average filtered to soften nearby peaks. The amplitude threshold and comparison to the predicted value between neighbouring R-waves are part of the decision rules. Every time a new R-wave was successfully identified, both of these rules were adaptively modified. The artifacts present in some of the signals were corrected using the threshold-based artifact correction algorithm proposed by the same software. Essentially, artefacts are corrected by comparing each RR interval value to a local average interval, obtained by median filtering the time series. If an RR interval differs from the local average more than a threshold, the interval is identified as an artifact and is marked for correction. The lowest threshold (0.45 s) was chosen for artifact detection. A cubic spline interpolation is used to replace the identified artefacts with interpolated values to make the fix. Afterwards, to obtain uniformly sampled signals, the RR series were interpolated using a piecewise cubic interpolation at 5 Hz, matching the IR sampling frequency. For the analysis, we selected the 150-s segments of HRV time-series corresponding to the ones selected for the thermal signal, for analogous reasons. For each HRV segment, we extracted several features related to the cardiovascular activity, as reported in Table 3.
The EDA signal reflects the activity of the sweat glands, which causes variations in the electrical properties of the skin. Sweat secretion in response to emotional stimuli is regulated by the sympathetic branch of the ANS. Therefore, EDA is a well-established signal to monitor the activity of the SNS in relation to changes in mental state. The EDA signal has a tonic component, henceforth referred to as EDAtonic, whose spectrum is below 0.05 Hz. The EDAtonic reflects the sympathetic tone and is related to changes in the overall state of the subject. Conversely, short-term and event-related responses are reflected in the phasic component. Here, we used the cvxEDA algorithm to extract these two components from the raw EDA signal. This algorithm performs a model-based signal decomposition exploiting Bayesian statistics, convex optimization, and sparsity, without the need for pre- and post-processing steps [30]. Afterwards, the tonic and phasic components were downsampled to 5 Hz, to obtain the same time points as the thermal and HRV series. In addition, we selected the 150-s segments Rest and Stress.
Of note, two female subjects were excluded from the analysis due to the presence of artifacts. Finally, from each segment of both the tonic and phasic components, we extracted well-known features related to the SNS activity, as listed in Table 4.

2.5. Cross-Mapping

In this study, we applied the cross-mapping procedure to investigate the nature of the facial thermal dynamics of healthy subjects in response to an external emotional stimulus. Based on the previous literature, we hypothesize that the main phenomena underlying a change in skin temperature are related to both sweat glands and cardiovascular dynamics [5,31,32]. Since these two phenomena are directly correlated to HRV and EDA, we applied the cross-mapping procedure to each couple of time series composed of a ANS correlate (i.e., HRV, EDAtonic) and a thermal signal, respectively denoted as X = (X 1 , X 2 , …, X n ) and Y = (Y 1 , Y 2 , …, Y n ). Of note, only the tonic component of the EDA signal was selected for this analysis, as it has comparable dynamics with the thermal signals. Specifically, we estimated the thermal signals of each ROI from the reconstructed manifolds of each physiological signal M X . Firstly, we used Takens’ Time-Delay Embedding Theorem to reconstruct the attractor of each time-series [33]. All the time series were simultaneously recorded and uniformly resampled at a frequency of 5 Hz. According to Takens’ theorem, we can trace the trajectory of a dynamic system in the phase space having only one time series and two parameters: the dimension of the phase space m and the time delay τ . Both of these parameters are unknown but can be estimated using well-known methods. Here, the value of τ was selected as the lag corresponding to the first minimum of the automutual information (MI) function of the time series used to reconstruct the phase space, ensuring statistical independence between the values, in both linear and nonlinear terms [34]. The embedding dimension m was selected using the False Nearest Neighbors algorithm (FNN), which is the most popular tool for the selection of the minimal embedding dimension [35]. This method is based on the assumption that if m is too small to unfold the attractor, points that are far from each other can appear as neighbours simply because the actual geometry of the attractor has been projected down onto a smaller space. These points are called false neighbours as they will no longer be neighbours when increasing the embedding dimension. Therefore, the algorithm iteratively increases the embedding dimension and at each step calculates the number of nearest neighbours. The final m is chosen as the dimension for which the change in the number of false neighbours compared to the previous iteration is below a threshold. Here, we used the algorithms proposed in [36,37] to select the values of τ and m, respectively.
The reconstructed attractors, are referred to as shadow manifolds in [21], and they are approximations of the true attractor. Indeed, the shadow manifolds of each variable is diffeomorphic from the original manifold. Takens’ theorem demonstrates that when two different variables are representations of the same dynamical system, their shadow manifolds are diffeomorphic to each other. Hence, if both X and Y dynamics are driven by the same underlying dynamical system, their manifolds are diffeomorphic and we can obtain the estimate of Y from the manifold M X , i.e., Y ^ |M X . In this study, we obtained Y ^ |M X following the procedure described in [38]. More specifically, since M X is diffeomorphic to M Y , points that are nearby in M X correspond to points that are also nearby in M Y . Therefore, we identified the m+1 nearest-neighbours in M X over time, and their corresponding points in M Y . The number of points depends on the maximum value of embedding dimension m between the two time series. The latter were used to map a region of m+1 points around Y t . The estimate of each point of Y starting from the reconstructed manifold of X is found as follows:
Y ^ t | M X = i = 1 m + 1 ω i Y t i
where ω i are the weights obtained as the Euclidean distances between X t and the nearest-neighbours points ( | | · | | indicates the Euclidean distance in R m ):
ω i = u i j = 1 m + 1 u j
u i = e x p | | X t X t i | | | | X t X t 1 | |
In Figure 3, we show, for a single subject in the two sessions, the estimates of the thermal signal of the nose from each physiological time series.
Finally, we computed the Spearman correlation coefficients between each original thermal signal and the corresponding estimate obtained using cross-mapping, for each subject during each condition, and we will refer to them as ρ Y Y ^ .

2.6. Statistical Analysis

We performed statistical analysis to explore at the group level the ANS correlates and the psychometric variation between the two sessions, Rest and Stress. The non-parametric Wilcoxon sign-rank test was used for each thermal, EDA and HRV feature between Rest and Stress ( α = 0.05). Similarly, the same test was used to assess variations between the PS level and the STAY-Y1 total score reported before and after the Stress. We controlled false discovery rate through the Benjamini, Krieger, and Yekutieli (FDR-BKY) correction for multiple hypothesis testing [39].
We assessed the significance level of each correlation coefficient ρ Y ^ Y using a non-parametric permutation test with 1000 repetitions. The original thermal signal was divided into blocks of length L, equal to τ , which were randomly shuffled at each iteration to generate the null distribution. Similarly, we computed the Spearman correlation coefficients between the original thermal signals and each physiological signal, ρ X Y , and their significance levels were assessed using the same permutation test.
For each physiological signal used to obtain the cross-mapped estimate of a thermal signal, we performed Wilcoxon sign-rank test with Bonferroni correction among the four different facial ROIs to test for variations in the correlation coefficient ρ Y ^ Y .
Finally, we performed a split-plot analysis of variance (SPANOVA), to statistically compare group-wise means between ρ Y ^ Y values obtained after the application of the cross-mapping procedure to each couple of time-series X, Y. Specifically, this statistical analysis consisted of testing for significant differences in the ρ Y ^ Y values between female and male participants, while subjecting participants to repeated measures. Thus, there is one between-subjects variable (Female vs. Male) and one within-subjects variable (Rest vs. Stress). This analysis was performed in RStudio, using the package rstatix.

3. Results

The statistical analysis on both psychometric tests indicated that subjects reported being more stressed after the stressor compared to before. In particular, the Wilcoxon test on the PS levels reported a p-value equal to 0.01, and the p-value for the STAY-Y1 total scores was equal to 0.03.
Among the facial regions, the nose was the only ROI to change significantly in the two conditions. In particular, the median temperature of the nose decreased significantly during Stress, as shown in Figure 4.
Regarding the ANS correlates data analysis, except for the LF/HF ratio, all the features extracted from the EDA and HRV confirmed a significant variation between Rest and Stress. Specifically, all the EDA features significantly increased during the stimulation, as shown in Figure 5.
Likewise, the mean of HRV increased significantly during the Stress. Conversely, the remaining HRV features, except the LF/HF ratio, decreased significantly compared to the Rest, as shown in Figure 6.
Regarding the correlation coefficients ρ X Y , over 50% of the subjects showed a significant correlation between the time series of the tonic component of the EDA signal and each thermal signal. The correlation between the HRV time series and the thermal signals was significant for less than 50% of the subjects. However, the median ρ X Y absolute values were always lower than 0.26, in both conditions. On the contrary, the p-values obtained for the correlation coefficients ρ Y ^ Y were lower than 0.05, and the median ρ Y ^ Y absolute values were always higher than 0.62, as shown in Table 5.
The results of the statistical analysis performed on the ρ Y ^ Y between facial ROIs are shown in Figure 7. The ρ Y ^ Y related to the thermal signal of the Nose region significantly differed from the others (i.e., Forehead, RCheek, LCheek) when estimated from the EDAtonic during both Rest and Stress, and when estimated from the HRV during Stress. In particular, the ρ Y ^ Y between the original thermal signal of the Nose and the same thermal signal estimated from each physiological signal, was always higher compared to the corresponding correlation coefficient for the other regions.
Finally, the results of the SPANOVA on the ρ Y ^ Y values are shown in Figure 8. As for the “between subjects” effect, the null hypothesis was that females and males have equal ρ Y ^ Y values, on average. This hypothesis was always rejected at the alpha level equal to 0.05 when comparing the ρ Y ^ Y values obtained after applying the cross-mapping approach to the HRV time series and each thermal signal.
The main effect for the condition, whose levels are Rest and Stress, was significant only for the ρ Y ^ Y values obtained after applying the cross-mapping approach to the HRV time series and the thermal signal of the nose (p-value = 0.03). In particular, the ρ Y ^ Y values increased during Stress. In this case, the interaction between gender and condition was not significant, indicating that the effect of the condition variable did not differ significantly between males and females. Finally, no significant differences were found between condition or gender in the correlation coefficients obtained after applying the cross-mapping approach to the EDA time series and the thermal signals from any ROI. However, the interaction term for the ρ Y ^ Y , where Y ^ is the cross-mapped estimate of the Forehead and Left Cheek’s temperature signal from the EDAtonic time series, was significant (p-value equal to 0.031 and 0.034, respectively). Hence, the effect of the condition variable on the correlation coefficient ρ Y ^ Y varied depending on the level of the gender variable. Specifically, for the Forehead, the effect of the condition variable on the correlation coefficient switched direction for males and females, resulting in a crossover pattern, as visible in Figure 8. In fact, ρ Y ^ Y increases in men during Stress compared to Rest, but it decreases in women.

4. Discussion

This study explores the physiological phenomena driving facial temperature modulation, under both parasympathetic and sympathetic dominance. Many studies have associated the temperature drop in the nose to a sympathetically driven vasoconstriction of the subcutaneous vessel, in response to a stressful stimulus [31,32,40,41]. Similar to this, perspiration brought on by stress or high arousal stimuli was linked to a temperature drop in the upperlip region with a specific dotted pattern [42,43]. Here, we wanted to give a foundation to the hypotheses that have prevailed in the interpretation of the results obtained in thermographic studies in psychophysiology. In particular, we performed a nonlinear analysis based on chaos theory, to investigate the imprint that the cardiovascular and sudomotor activity have on the temperature extracted from different regions of the face. Moreover, we quantified the influence that these two systems have on skin temperature, in two different conditions, namely resting state and stress. In addition, we explored potential gender differences in the physiology behind thermal modulation.
First of all, we confirmed the effectiveness of the protocol in inducing stress during the stimulation by comparing the psychometric scores reported before and after the stressor. Objective physiological measures gave us further evidence of the increased arousal during the stressor. Although the features we extracted are naturally standardized in the physiological ranges, outliers are frequently present in studies of physiological aspects due to the high intrinsic subject variability. To limit the potential impact of these outliers on our results, we have used non-parametric tests that are less influenced by extreme values. All the features extracted from the EDA signal, increased significantly during Stress. Similarly, all the HRV features except one (the LF/HF ratio), both in the time and frequency domain, varied significantly in the expected direction between the two conditions. Finally, the median temperature of the nose decreased during Stress, as widely described in the literature [26].
Nonlinear analysis of the thermal and HRV time series revealed a strong relation between the thermal and cardiovascular dynamics that was not obvious in the linear domain. Indeed, the linear correlation between the original HRV and thermal series was significant in the majority of the subjects only for the nose region during Stress, and the correlation values ρ X Y reported a median value lower than 0.15. On the other hand, the correlation values between the thermal signals and their estimated version obtained using the cross-mapping approach with the HRV time series were significant for all the subjects during both Rest and Stress, with a median value higher than 0.7. Specifically, the nose region presented the highest median correlation value, that increased during Stress compared to Rest.
Conversely, the relation between the thermal and electrodermal dynamics was already visible in the linear domain, as the linear correlation was significant in the majority of the subjects for all the facial regions and in both conditions. However, median correlation values ρ X Y were lower than 0.25. In this case, the correlation values between the thermal signals and their estimated version obtained using the cross-mapping approach with the time series of tonic component of the EDA were always significant and their median value was close to 1. This result suggests that the thermal signals from the four selected regions and the tonic electrodermal component describe the dynamics of the same system, thus perspiration seems to play an important role in the facial thermal modulation, possibly higher than the vascular activity.
The nonlinear correlation values quantify how the knowledge of the shadow manifold obtained from the time series of the physiological time series can be used to estimate values of the thermal time series. Hence, we used the correlation value ρ Y ^ Y as an index of the ‘imprint’ of the EDA and HRV time series on the thermal signals of the most commonly selected regions of the face. A comparison of the nonlinear correlation values obtained for the different facial regions reveals that, although the correlation coefficient is very high for all the regions, the nose has a significantly higher influence of the electrodermal activity compared to the other face areas, during both conditions. On the other hand, the influence of cardiovascular activity seems to be similar in all the regions during Rest, while it is significantly higher in the nose during Stress. This different behaviour of EDA- and HRV-estimated thermal signals between rest and stress is confirmed by the SPANOVA comparison among conditions and genders. In fact, the comparison between conditions and genders was never significant when considering the EDA-estimated thermal signal of the nose, as well as of the other regions. In fact, the EDA appears to have a significant impact on all facial temperatures, at a similar level in both conditions and without obvious gender variations. The imprint of the HRV on the variation of the thermal signal from the four ROIs varies significantly between the rest and stress. In particular, the contribution of the HRV is on average higher during stress in the nose and forehead. Hence, we can speculate that the well-known decrease in temperature of the nose tip during stress is mostly due to the sympathovagal-driven vasoconstriction of the subcutaneous vessel. The median temperature of the forehead did not change significantly during stress, as shown from the first analysis by means of the Wilcoxon signed-rank test. Thus, we cannot associate the increase in the ρ Y ^ Y with a particular physiological phenomenon with the same confidence. Our findings indicate a gender-dependent variation between stress and rest states. Specifically, in contrast to their male counterparts, female subjects demonstrated a significant decrease in correlation values under stress conditions, albeit the average correlation values remained high. This observed difference may potentially be indicative of inherent biological gender disparities within these regions of interest. However, it is worthwhile noting that these are preliminary findings and further comprehensive investigations are required to draw definitive conclusions.
A limitation of this study is gender imbalance. The disparity between males and females in our study population could have, in fact, hidden some other results that would have led to a better interpretation. In general, the sample size we had available was limited. Hence, future work will expand the study to a larger cohort of subjects to confirm the findings and strengthen the conclusions drawn. Moreover, we intent to compare thermal modulation in different populations. Indeed, the sympathetic and parasympathetic nervous systems’ responses to emotional stress can vary depending not only on gender, but also age or clinical conditions. Hence, comparing thermal modulation in various groups, can help to comprehend the underlying mechanisms and potentially lead to customised therapies.

5. Conclusions

In conclusion, our study presented a novel approach to the study of the nonlinear nature of thermal modulation. We have researched and, in a preliminary manner, shown the effect of both the sympathetic or parasympathetic autonomic subsystems on the facial temperature variation in response to emotional stress. Future work will explore the nonlinear coupling between the autonomic nervous system and the temperature of a broader range of facial regions, ideally encompassing the entire face. Advanced approaches like finite element methods may be employed to complete this thorough investigation. Understanding the mechanisms underlying these temperature variations may make it possible to create more precise techniques for stress detection and stress management, lie detection, and the prediction of underlying psychiatric conditions. Finally, the imprint of autonomic correlates in the thermal signal, whose existence we have demonstrated, paves the way for future research aimed at extracting information specifically related to HRV and EDA directly from the thermal signal.

Author Contributions

Conceptualization, all authors; methodology, M.N., A.G. and F.G.; software, F.G. and M.N.; validation, A.G., M.N. and F.G.; formal analysis, F.G.; investigation, F.G., M.N. and A.G.; resources, F.G., M.N. and A.G.; data curation, F.G.; writing—original draft preparation, F.G.; writing—review and editing, all authors; visualization, F.G., M.N. and A.G.; supervision, M.N., E.P.S. and A.G.; project administration, E.P.S. and A.G.; funding acquisition, E.P.S. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results has received partial funding from the Italian Ministry of Education and Research (MIUR) in the framework of the FoReLab project (Departments of Excellence), from European Union Horizon 2020 Programme under grant agreement n 824153 of the project “POTION: Promoting Social Interaction through Emotional Body Odours”, and from PNRR-M4C2-Investimento 1.3, Partenariato Esteso PE00000013-“FAIR- Future Artificial Intelligence Research”-Spoke 1 “Human-centered AI”, funded by the European Commission under the NextGeneration EU programme.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Bioethics Committee of the University of Pisa Review No. 14/2019, 3 May 2019.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data are not publicly available due to privacy restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Al Qudah, M.M.; Mohamed, A.S.; Lutfi, S.L. Affective State Recognition Using Thermal-Based Imaging: A Survey. Comput. Syst. Sci. Eng. 2021, 37, 47–62. [Google Scholar] [CrossRef]
  2. Wang, Z.; Zhou, X.; Wang, W.; Liang, C. Emotion recognition using multimodal deep learning in multiple psychophysiological signals and video. Int. J. Mach. Learn. Cybern. 2020, 11, 923–934. [Google Scholar] [CrossRef]
  3. Smith, T.W.; Deits-Lebehn, C.; Williams, P.G.; Baucom, B.R.; Uchino, B.N. Toward a social psychophysiology of vagally mediated heart rate variability: Concepts and methods in self-regulation, emotion, and interpersonal processes. Soc. Personal. Psychol. Compass 2020, 14, e12516. [Google Scholar] [CrossRef]
  4. Lang, P.J. The cognitive psychophysiology of emotion: Fear and anxiety. In Anxiety and the Anxiety Disorders; Routledge: Abingdon, UK, 2019; pp. 131–170. [Google Scholar]
  5. Sonkusare, S.; Breakspear, M.; Pang, T.; Nguyen, V.T.; Frydman, S.; Guo, C.C.; Aburn, M.J. Data-driven analysis of facial thermal responses and multimodal physiological consistency among subjects. Sci. Rep. 2021, 11, 12059. [Google Scholar] [CrossRef] [PubMed]
  6. Cho, Y.; Julier, S.J.; Bianchi-Berthouze, N. Instant stress: Detection of perceived mental stress through smartphone photoplethysmography and thermal imaging. JMIR Ment. Health 2019, 6, e10140. [Google Scholar] [CrossRef] [Green Version]
  7. Perpetuini, D.; Cardone, D.; Filippini, C.; Chiarelli, A.M.; Merla, A. Modelling impulse response function of functional infrared imaging for general linear model analysis of autonomic activity. Sensors 2019, 19, 849. [Google Scholar] [CrossRef] [Green Version]
  8. Filippini, C.; Cardone, D.; Perpetuini, D.; Chiarelli, A.; Petitto, L.; Merla, A. Assessment of autonomic response in 6–12-month-old babies during the interaction with robot and avatar by means of thermal infrared imaging. Quant. Infrared Thermogr. J. 2022, 20, 78–91. [Google Scholar] [CrossRef]
  9. Kiashari, S.E.H.; Nahvi, A.; Bakhoda, H.; Homayounfard, A.; Tashakori, M. Evaluation of driver drowsiness using respiration analysis by thermal imaging on a driving simulator. Multimed. Tools Appl. 2020, 79, 17793–17815. [Google Scholar] [CrossRef]
  10. He, C.; Mahfouf, M.; Torres-Salomao, L.A. Thermal Imaging for Psychophysiological State Detection in the Human-Machine Interface (HMI) Control System. In Proceedings of the Informatics in Control, Automation and Robotics: 15th International Conference, ICINCO 2018, Porto, Portugal, 29–31 July 2018; Revised Selected Papers 15. Springer: Cham, Switzerland, 2020; pp. 214–229. [Google Scholar]
  11. Goldberger, A.L.; West, B.J. Fractals in physiology and medicine. Yale J. Biol. Med. 1987, 60, 421. [Google Scholar]
  12. Henriques, T.; Ribeiro, M.; Teixeira, A.; Castro, L.; Antunes, L.; Costa-Santos, C. Nonlinear methods most applied to heart-rate time series: A review. Entropy 2020, 22, 309. [Google Scholar] [CrossRef] [Green Version]
  13. Delliaux, S.; Delaforge, A.; Deharo, J.C.; Chaumet, G. Mental workload alters heart rate variability, lowering non-linear dynamics. Front. Physiol. 2019, 10, 565. [Google Scholar] [CrossRef] [Green Version]
  14. Nardelli, M.; Greco, A.; Sebastiani, L.; Scilingo, E.P. ComEDA: A new tool for stress assessment based on electrodermal activity. Comput. Biol. Med. 2022, 150, 106144. [Google Scholar] [CrossRef]
  15. Vuksanović, V.; Sheppard, L.W.; Stefanovska, A. Nonlinear relationship between level of blood flow and skin temperature for different dynamics of temperature change. Biophys. J. 2008, 94, L78–L80. [Google Scholar] [CrossRef] [Green Version]
  16. Kodera, S.; Kamiya, T.; Miyazawa, T.; Hirata, A. Correlation between estimated thermoregulatory responses and pacing in athletes during marathon. IEEE Access 2020, 8, 173079–173091. [Google Scholar] [CrossRef]
  17. Kitney, R. An analysis of the nonlinear behaviour of the human thermal vasomotor control system. J. Theor. Biol. 1975, 52, 231–248. [Google Scholar] [CrossRef]
  18. Binek, M.; Drzazga, Z.; Socha, T.; Pokora, I. Do exist gender differences in skin temperature of lower limbs following exercise test in male and female cross-country skiers? J. Therm. Anal. Calorim. 2022, 147, 7373–7383. [Google Scholar] [CrossRef]
  19. Kaciuba-Uscilko, H.; Grucza, R. Gender differences in thermoregulation. Curr. Opin. Clin. Nutr. Metab. Care 2001, 4, 533–536. [Google Scholar] [CrossRef]
  20. Christensen, J.; Væth, M.; Wenzel, A. Thermographic imaging of facial skin—gender differences and temperature changes over time in healthy subjects. Dentomaxillofacial Radiol. 2012, 41, 662–667. [Google Scholar] [CrossRef] [Green Version]
  21. Sugihara, G.; May, R.; Ye, H.; Hsieh, C.h.; Deyle, E.; Fogarty, M.; Munch, S. Detecting causality in complex ecosystems. Science 2012, 338, 496–500. [Google Scholar] [CrossRef]
  22. Nardelli, M.; Vanello, N.; Galperti, G.; Greco, A.; Scilingo, E.P. Assessing the quality of heart rate variability estimated from wrist and finger ppg: A novel approach based on cross-mapping method. Sensors 2020, 20, 3156. [Google Scholar] [CrossRef]
  23. Nardelli, M.; Greco, A.; Vanello, N.; Scilingo, E.P. Reliability of Pulse Rate Variability in Elderly Men and Women: An Application of Cross-Mapping Approach. In Proceedings of the 2021 43rd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Guadalajara, Mexico, 1–5 November 2021; IEEE: New York, NY, USA, 2021; pp. 492–495. [Google Scholar]
  24. Stroop, J.R. Studies of interference in serial verbal reactions. J. Exp. Psychol. 1935, 18, 643. [Google Scholar] [CrossRef]
  25. Gomes Moreira, D.; Costello, J.T.; Brito, C.J.; Sillero Quintana, M. A checklist for measuring skin temperature with infrared thermography in sports and exercise medicine. Thermol. Int. 2017, 27, 141–143. [Google Scholar]
  26. Cho, Y.; Bianchi-Berthouze, N. Physiological and affective computing through thermal imaging: A survey. arXiv 2019, arXiv:1908.10307. [Google Scholar]
  27. Gioia, F.; Pascali, M.A.; Greco, A.; Colantonio, S.; Scilingo, E.P. Discriminating Stress From Cognitive Load Using Contactless Thermal Imaging Devices. In Proceedings of the 2021 43rd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Guadalajara, Mexico, 1–5 November 2021; IEEE: New York, NY, USA, 2021; pp. 608–611. [Google Scholar]
  28. Nirkin, Y.; Masi, I.; Tuan, A.T.; Hassner, T.; Medioni, G. On face segmentation, face swapping, and face perception. In Proceedings of the 2018 13th IEEE International Conference on Automatic Face & Gesture Recognition (FG 2018), Xi’an, China, 15–19 May 2018; IEEE: New York, NY, USA, 2018; pp. 98–105. [Google Scholar]
  29. Tarvainen, M.P.; Niskanen, J.P.; Lipponen, J.A.; Ranta-Aho, P.O.; Karjalainen, P.A. Kubios HRV–heart rate variability analysis software. Comput. Methods Programs Biomed. 2014, 113, 210–220. [Google Scholar] [CrossRef] [PubMed]
  30. Greco, A.; Valenza, G.; Lanata, A.; Scilingo, E.P.; Citi, L. cvxEDA: A convex optimization approach to electrodermal activity processing. IEEE Trans. Biomed. Eng. 2015, 63, 797–804. [Google Scholar] [CrossRef] [Green Version]
  31. Filippini, C.; Spadolini, E.; Cardone, D.; Merla, A. Thermal imaging based affective computing for educational robot. Proceedings 2019, 27, 27. [Google Scholar]
  32. Cho, Y.; Bianchi-Berthouze, N.; Oliveira, M.; Holloway, C.; Julier, S. Nose heat: Exploring stress-induced nasal thermal variability through mobile thermal imaging. In Proceedings of the 2019 8th International Conference on Affective Computing and Intelligent Interaction (ACII), Cambridge, UK, 3–6 September 2019; IEEE: New York, NY, USA, 2019; pp. 566–572. [Google Scholar]
  33. Takens, F. Detecting strange attractors in turbulence. In Dynamical Systems and Turbulence, Warwick 1980; Springer: Berlin/Heidelberg, Germany, 1981; pp. 366–381. [Google Scholar]
  34. de Pedro-Carracedo, J.; Fuentes-Jimenez, D.; Ugena, A.M.; Gonzalez-Marcos, A.P. Phase space reconstruction from a biological time series: A photoplethysmographic signal case study. Appl. Sci. 2020, 10, 1430. [Google Scholar] [CrossRef] [Green Version]
  35. Krakovská, A.; Mezeiová, K.; Budáčová, H. Use of false nearest neighbours for selecting variables and embedding parameters for state space reconstruction. J. Complex Syst. 2015, 2015, 932750. [Google Scholar] [CrossRef] [Green Version]
  36. Grinsted, A. Average Mutual Information, 2022.
  37. Kennel, M.B.; Brown, R.; Abarbanel, H.D. Determining embedding dimension for phase-space reconstruction using a geometrical construction. Phys. Rev. A 1992, 45, 3403. [Google Scholar] [CrossRef] [Green Version]
  38. Mønster, D.; Fusaroli, R.; Tylén, K.; Roepstorff, A.; Sherson, J.F. Causal inference from noisy time-series data—testing the convergent cross-mapping algorithm in the presence of noise and external influence. Future Gener. Comput. Syst. 2017, 73, 52–62. [Google Scholar] [CrossRef] [Green Version]
  39. Benjamini, Y.; Krieger, A.M.; Yekutieli, D. Adaptive linear step-up procedures that control the false discovery rate. Biometrika 2006, 93, 491–507. [Google Scholar] [CrossRef]
  40. Engert, V.; Merla, A.; Grant, J.A.; Cardone, D.; Tusche, A.; Singer, T. Exploring the use of thermal infrared imaging in human stress research. PLoS ONE 2014, 9, e90782. [Google Scholar] [CrossRef] [Green Version]
  41. Perpetuini, D.; Cardone, D.; Filippini, C.; Spadolini, E.; Mancini, L.; Chiarelli, A.M.; Merla, A. Can Functional Infrared Thermal Imaging Estimate Mental Workload in Drivers as Evaluated by Sample Entropy of the fNIRS Signal? In Proceedings of the 8th European Medical and Biological Engineering Conference: Proceedings of the EMBEC 2020, Portorož, Slovenia, 29 November–3 December 2020; Springer: Cham, Switzerland, 2021; pp. 223–232. [Google Scholar]
  42. Merla, A.; Romani, G.L. Thermal signatures of emotional arousal: A functional infrared imaging study. In Proceedings of the 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Lyon, France, 22–26 August 2007; IEEE: New York, NY, USA, 2007; pp. 247–249. [Google Scholar]
  43. Shastri, D.; Papadakis, M.; Tsiamyrtzis, P.; Bass, B.; Pavlidis, I. Perinasal imaging of physiological stress and its affective potential. IEEE Trans. Affect. Comput. 2012, 3, 366–378. [Google Scholar] [CrossRef]
Figure 1. Protocol Timeline. The dashed line indicates non-timed tasks. PS—Perceived Stress self-assessment.
Figure 1. Protocol Timeline. The dashed line indicates non-timed tasks. PS—Perceived Stress self-assessment.
Sensors 23 06403 g001
Figure 2. Facial ROIs location.
Figure 2. Facial ROIs location.
Sensors 23 06403 g002
Figure 3. Cross-mapped estimates of the thermal signal of the nose (Nose) from each physiological time-series (i.e., EDAtonic, HRV), for a single subject in the two conditions Rest (Left), and Stress (Right). Each cross-mapped estimate is superimposed on the original time series.
Figure 3. Cross-mapped estimates of the thermal signal of the nose (Nose) from each physiological time-series (i.e., EDAtonic, HRV), for a single subject in the two conditions Rest (Left), and Stress (Right). Each cross-mapped estimate is superimposed on the original time series.
Sensors 23 06403 g003
Figure 4. Statistical comparison of the median temperature, for each facial region of interest. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05).
Figure 4. Statistical comparison of the median temperature, for each facial region of interest. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05).
Sensors 23 06403 g004
Figure 5. Statistical comparison for each of the EDA features extracted. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05; *** = p < 0.001, **** = p < 0.0001).
Figure 5. Statistical comparison for each of the EDA features extracted. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05; *** = p < 0.001, **** = p < 0.0001).
Sensors 23 06403 g005
Figure 6. Statistical comparison for each of the HRV features extracted. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05; ** = p < 0.01; *** = p < 0.001).
Figure 6. Statistical comparison for each of the HRV features extracted. Significant differences between Rest = R and Stress = S after the Wilcoxon signed rank test with FDR-BKY correction are highlighted with an asterisk (* = p < 0.05; ** = p < 0.01; *** = p < 0.001).
Sensors 23 06403 g006
Figure 7. Statistical comparison of the ρ Y ^ Y obtained as the correlation between the Y ^ estimated from the same physiological signal and Y corresponding to each facial ROI (Nose = N; Forehead = F; RC = Right Cheek; LC = Left Cheek), during both conditions (R = Rest and S = Stress). Significant differences between facial ROIs after the Wilcoxon signed rank test with Bonferroni correction are highlighted with an asterisk (* = p < 0.05; *** = p < 0.001; **** = p < 0.0001).
Figure 7. Statistical comparison of the ρ Y ^ Y obtained as the correlation between the Y ^ estimated from the same physiological signal and Y corresponding to each facial ROI (Nose = N; Forehead = F; RC = Right Cheek; LC = Left Cheek), during both conditions (R = Rest and S = Stress). Significant differences between facial ROIs after the Wilcoxon signed rank test with Bonferroni correction are highlighted with an asterisk (* = p < 0.05; *** = p < 0.001; **** = p < 0.0001).
Sensors 23 06403 g007
Figure 8. Split-Plot Anova results for each correlation coefficient ρ Y ^ Y obtained as the correlation between the time-series Y ^ estimated from each physiological signal (i.e., HRV, EDAtonic) and Y corresponding to each facial ROI (Nose, Forehead, Right Cheek, Left Cheek); ns = non-significant; * = significant; Rest = R and Stress = S.
Figure 8. Split-Plot Anova results for each correlation coefficient ρ Y ^ Y obtained as the correlation between the time-series Y ^ estimated from each physiological signal (i.e., HRV, EDAtonic) and Y corresponding to each facial ROI (Nose, Forehead, Right Cheek, Left Cheek); ns = non-significant; * = significant; Rest = R and Stress = S.
Sensors 23 06403 g008
Table 1. Number of Pixels per ROI (mean ± standard deviation).
Table 1. Number of Pixels per ROI (mean ± standard deviation).
NoseForeheadLCheekRCheek
275.1 ± 37.52540.2 ± 346.41288 ± 174.41288 ± 174.4
Table 2. Average Temperature within each ROI (Median ± MAD).
Table 2. Average Temperature within each ROI (Median ± MAD).
NoseForeheadLCheekRCheek
34.23 ± 0.1935.98 ± 0.1334.38 ± 0.3934.54 ± 0.4
Table 3. Features extracted from the HRV.
Table 3. Features extracted from the HRV.
FeatureDescription
meanHRVMean value of HRV
stdHRVStandard deviation of HRV
RMSSDRoot mean square of successive RR interval differences
pNN50Percentage of successive RR intervals that differ by more than 50 ms
LFPercentage of the total power in the low-frequency (0.04–0.15 Hz)
HFPercentage of the total power in the high-frequency (0.15–0.40 Hz)
LF/HF ratioRatio of LF to HF power
Table 4. Features extracted from the tonic and phasic components.
Table 4. Features extracted from the tonic and phasic components.
FeatureDescription
TonicMeanMean value of the tonic component *
TonicStdStandard deviation of the tonic component *
PhasicMeanMean value of the phasic component **
PhasicStdStandard deviation of the phasic component **
PksMaxMaximum peak **
NPksNumber of peaks **
PksSumSum of peaks **
EDASympPower of the skin conductance spectrum in the range of 0.04–0.25 Hz
*: value averaged within non-overlapped 30-second-long time windows. **: value averaged within non-overlapped 5-second-long time windows.
Table 5. Percentage of significant correlation coefficients ρ X Y and ρ Y ^ Y (p < 0.05) and median ρ value.
Table 5. Percentage of significant correlation coefficients ρ X Y and ρ Y ^ Y (p < 0.05) and median ρ value.
X
EDAtonicHRV
RESTSTRESSRESTSTRESS
% ρ % ρ % ρ % ρ
YNose ρ X Y 75−0.2685.710.1317.86−0.0260.71−0.17
ρ Y ^ Y 1000.991000.991000.671000.74
Forehead ρ X Y 60.71−0.0257.140.027.140.0239.290.03
ρ Y ^ Y 1000.98100 0.981000.661000.64
Rcheek ρ X Y 53.57−0.1760.710.0514.29−0.0128.57−0.01
ρ Y ^ Y 1000.981000.971000.651000.66
Lcheek ρ X Y 71.43−0.1164.290.0210.710.0232.140.01
ρ Y ^ Y 1000.981000.971000.621000.66
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gioia, F.; Nardelli, M.; Scilingo, E.P.; Greco, A. Autonomic Regulation of Facial Temperature during Stress: A Cross-Mapping Analysis. Sensors 2023, 23, 6403. https://doi.org/10.3390/s23146403

AMA Style

Gioia F, Nardelli M, Scilingo EP, Greco A. Autonomic Regulation of Facial Temperature during Stress: A Cross-Mapping Analysis. Sensors. 2023; 23(14):6403. https://doi.org/10.3390/s23146403

Chicago/Turabian Style

Gioia, Federica, Mimma Nardelli, Enzo Pasquale Scilingo, and Alberto Greco. 2023. "Autonomic Regulation of Facial Temperature during Stress: A Cross-Mapping Analysis" Sensors 23, no. 14: 6403. https://doi.org/10.3390/s23146403

APA Style

Gioia, F., Nardelli, M., Scilingo, E. P., & Greco, A. (2023). Autonomic Regulation of Facial Temperature during Stress: A Cross-Mapping Analysis. Sensors, 23(14), 6403. https://doi.org/10.3390/s23146403

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop