Elementary Integrate-and-Fire Process Underlies Pulse Amplitudes in Electrodermal Activity
Authors: Subramanian Sandya , Purdon Patrick L. , Barbieri Riccardo , Brown Emery N. 
Published: January 2021

Electrodermal activity (EDA) is a direct read-out of sweat-induced changes in the skin’s electrical conductance. Sympathetically-mediated pulsatile changes in skin sweat measured as EDA resemble an integrate-and-fire process, which yields an inverse Gaussian model as the inter-pulse interval distribution. We have previously showed that the inter-pulse intervals in EDA follow an inverse Gaussian distribution. However, the statistical structure of EDA pulse amplitudes has not yet been characterized based on the physiology. Expanding upon the integrate-and-fire nature of sweat glands, we hypothesized that the amplitude of an EDA pulse is proportional to the excess volume of sweat produced compared to what is required to just reach the surface of the skin. We modeled this as the difference of two inverse Gaussian models for each pulse, one which represents the time required to produce just enough sweat to rise to the surface of the skin and one which represents the time requires to produce the actual volume of sweat. We proposed and tested a series of four simplifications of our hypothesis, ranging from a single difference of inverse Gaussians to a single simple inverse Gaussian. We also tested four additional models for comparison, including the lognormal and gamma distributions. All models were tested on EDA data from two subject cohorts, 11 healthy volunteers during 1 hour of quiet wakefulness and a different set of 11 healthy volunteers during approximately 3 hours of controlled propofol sedation. All four models represent simplifications of our hypothesis outperformed other models across all 22 subjects, as measured by Akaike’s Information Criterion (AIC), as well as mean and maximum distance from the diagonal on a quantile-quantile plot. Our broader model set of four simplifications offered a useful framework to enhance further statistical descriptions of EDA pulse amplitudes. Some of the simplifications prioritize fit near the mode of the distribution, while others prioritize fit near the tail. With this new insight, we can summarize the physiologically-relevant amplitude information in EDA with at most four parameters. Our findings establish that physiologically based probability models provide parsimonious and accurate description of temporal and amplitude characteristics in EDA.

AUTHOR SUMMARY

Electrodermal activity (EDA) is an indirect read-out of the body’s sympathetic nervous system, or fight-or-flight response, measured as sweat-induced changes in the electrical conductance properties of the skin. Interest is growing in using EDA to track physiological conditions such as stress levels, sleep quality, and emotional states. Our previous worked showed that the times in between EDA pulses obeyed a specific statistical distribution, the inverse Gaussian, that arises from the physiology of EDA production. In this work, we build on that insight to analyze the amplitudes of EDA pulses. In an analysis of EDA data recorded in 11 healthy volunteers during quiet wakefulness and 11 different healthy volunteers during controlled propofol sedation, we establish that the amplitudes of EDA pulses also have specific statistical structure, as the differences of inverse Gaussians, that arises from the physiology of sweat production. We capture that structure using a series of progressively simpler models that each prioritize different parts of the pulse amplitude distribution. Our findings show a physiologically-based statistical model provides a parsimonious and accurate description EDA. This enables increased reliability and robustness in analyzing EDA data collected under any circumstance.