Abstract
Emotion research typically searches for consistency and specificity in physiological activity across instances of an emotion category, such as anger or fear, yet studies to date have observed more variation than expected. In the present study, we adopt an alternative approach, searching inductively for structure within variation, both within and across participants. Following a novel, physiologically-triggered experience sampling procedure, participantsâ self-reports and peripheral physiological activity were recorded when substantial changes in cardiac activity occurred in the absence of movement. Unsupervised clustering analyses revealed variability in the number and nature of patterns of physiological activity that recurred within individuals, as well as in the affect ratings and emotion labels associated with each pattern. There were also broad patterns that recurred across individuals. These findings support a constructionist account of emotion which, drawing on Darwin, proposes that emotion categories are populations of variable instances tied to situation-specific needs.
Similar content being viewed by others
Introduction
Scientists who study emotion acknowledge that there is variability in the expression of emotion across individuals and cultures, but differ in whether they view this variability as an intrinsic or epiphenomenal property of emotion (see supplementary Figure S1). Classical views of emotion, for example, propose that there is a single pattern of autonomic nervous system (ANS) activity that is diagnostic of each emotion category (e.g., anger, fear), and consider variability to be epiphenomenal to emotion. In these views, variability is assumed to derive from regulatory processes, display rules, an inability to sample intense emotional events in the lab, and/or stochastic noise1,2,3,4,5,6. Prototype views of emotion propose that ANS activity varies across instances of the same emotion category, but it is still assumed that each category has a diagnostic (i.e., typical or frequent) pattern that can be used to consistently and specifically identify its instances (consonant with several theoristsâ work7,8,9). We refer to the idea of an emotion-specific diagnostic pattern as a fingerprints hypothesis10. This hypothesis is often assumed by basic emotion and causal appraisal approaches to emotion1,3,6,7,9. Real fingerprints vary from one instance to the next depending on a host of variables (such as the amount of sweat, pressure or surfaces touched) but are presumed to be sufficiently unique to identify a person across instances. In the same way, if the âfearâ fingerprint were an increase in heart rate and electrodermal (i.e., sweat gland) activity and a decrease in skin temperature11, then this fingerprint may not be present in full for every instance of fear. However, it would be present in many instances of fear both within people across situations and across persons, and it would not be observed at levels above chance in instances of any other emotion category (such as anger or happiness).
Similarly, multiple scientists have proposed that emotion categories are uniquely associated with the affective features of valence and arousal, implying that emotion categories also have affective fingerprints7,9,12,13. From these perspectives, for example, the experience of fear would be consistently and specifically associated with negative, high arousal affect, and would not be associated with pleasant or low arousal affect at levels above chance. Furthermore, by these accounts, affective features should have diagnostic biological correlates, such that experienced pleasantness and activation should be consistently and specifically associated with patterns of biological activity.
Alternatively, other theoretical views in the science of emotion employ a populations hypothesis, proposing that each emotion category is a population of context-dependent instances whose features, including ANS features, vary in a situated fashion. This hypothesis has been proposed by constructionist, behavioral ecological, and functionalist approaches to emotion. These views hypothesize, a priori, that there is substantial within-category variation as well as between-category similarity that far exceeds what is accounted for by a classical or prototype approach14,15,16,17,18,19,20. A populations hypothesis draws on Darwinâs discovery of population thinking21, in which a biological category, such as a species, consists of heterogeneous individuals. Biological categories do not represent discrete physical types (e.g., consistent and specific patterns of features), and in fact, âidealâ instances of categories (i.e., their prototypes) are statistical summaries that need not exist in nature22. Accordingly, these categories are conceptual: they are collections of variable instances that are treated as similar for some function or purpose16,17,18. A populations view of emotions thus posits that emotions, as biological categories, consist of heterogeneous instances which vary in their physiological and psychological features, but which are treated as similar for some function or purpose.
With regard to affect, a populations hypothesis proposes that valence and arousal are features of consciousness that can vary within an emotion category23,24. For example, the experience of fear can be affectively pleasant, such as when riding a roller coaster. Critically, a populations hypothesis proposesâand past research has demonstratedâthat both valence25,26 and arousal27 can vary in their biological correlates.
Consistent with a populations hypothesis, a recent meta-analysis found considerable variability in patterns of peripheral physiological activity within emotion categories as well as considerable similarities across categories10, replicating prior meta-analyses28,29 and demonstrating that ANS patterns for emotion categories are neither consistent nor specific. These findings parallel results showing a lack of consistency and specificity in other biological measures taken during instances of emotion, including facial muscle movements30 and neural correlatesâwhether measured at the level of individual neurons31, as activity in specific brain regions23,32,33, as activity or connectivity in brain networks34,35, or as distributed patterns of activity36. Despite the accumulating evidence in support of a populations hypothesis, however, contemporary research continues to be designed with the assumptions of a fingerprints hypothesis in mind37,38.
Meta-analytic findings, while consistent with the populations hypothesis, cannot rigorously test it; meta-analytic approaches rely on merely summarizing data across people and situations rather than designing studies to specifically assess intra-individual and inter-individual variability across instances. Moreover, standard lab-based emotion induction methods limit the intensity of the emotional instances sampled5, and thus may under-estimate context-specific variation in the physiological features of emotional instances. In the present study, we used a more idiographically-sensitive method that provided us with the opportunity to observe and characterize the extent of physiological variability within and across emotion categories in real-world situations39. Previous studies that have attempted to observe within-category variation using experience-sampling in daily life have found evidence supporting the populations hypothesis. For example, an analysis of momentary reports of experienced affect found that the relationship between valence and arousal varied widely across individuals and situations40. Experience sampling approaches are often employed to assess self-report features of experience in natural settings41,42, yet this work has only just begun to integrate ambulatory monitoring of peripheral physiology39,43,44. Research that collects both self-report and physiological features of experience in everyday life is therefore a critical next step in testing a populations hypothesis.
Our methodâa novel, physiologically-triggered experience-sampling approachâwas specifically designed to sample evocative moments defined in a biological way. Participants were prompted to report on their affective experience upon substantial, sustained changes in their cardiac activity in the absence of physical movement or posture change. By triggering experience samples in this way, we were able to efficiently obtain self-reports when individuals evidenced greater physiological change, and when those changes were most likely due to greater psychological salience rather than gross bodily movement. We then employed an unsupervised machine learning approach to discover the number of distinct, frequently repeated patterns of physiological activity within each individual, and we assessed whether any of these patterns recurred across participants. We also assessed the extent to which any of these patterns were associated with participant-provided labels for emotional experience and/or affective features. This approach allowed us to investigate the relationship between physiological reactivity and emotion categories within individuals in a data-driven way.
As part of a larger study on affective experience and decision making in daily life, 52 participants completed approximately 14, eight-hour days of ambulatory peripheral physiological monitoring. Each day, participants visited the lab and were outfitted with sensors and portable equipment to measure their electrocardiogram (ECG) and impedance cardiogram (ICG) as well as bodily movement and posture (via accelerometers). A custom-built smartphone application initiated an experience sampling prompt any time the interbeat interval (IBI; also called heart period) changed by more thanâ±â167 ms over an eight-second period, with these thresholds adjusted per participant to ensure they received a comparable number of prompts per day. Prompts were not generated if participants had moved substantially or shifted posture within the proceeding 30 s. In response to each experience sampling prompt, participants first rated their felt valence and arousal on 100-point continuous slider scales and then freely labeled their current affective experience. Confirming that our approach helped us to target evocative affective experiences in everyday life, we found that experience sampling events included ratings for valence (Mâ=â10.16, SDâ=â17.43) and arousal (Mâ=â6.40, SDâ=â19.42) that each covered the entire possible range from -50 toâ+â50. To ensure comparability of results across the sample, we excluded participants who had fewer than 70 events devoid of major physiological artifact. This resulted in the removal of six participants, resulting in a final sample size of 46 (see âMethodsâ for participant demographics).
Using the ambulatory physiological data, we derived six cardiovascular features for analysis (Table 1): interbeat interval (IBI), respiratory sinus arrhythmia (RSA), pre-ejection period (PEP), left ventricular ejection time (LVET), stroke volume (SV), and cardiac output (CO). These cardiovascular features are frequently used in the motivated performance literature, making them highly relevant to sampling evocative experiences in-the-world45,46,47,48. Change scores for each feature were computed for each experience sampling event, as the difference in physiological activity between the 30Â s preceding the IBI change that initiated the experience sampling prompt and the 30Â s following (see supplementary Figure S2). Change scores for each participant were submitted to a within-person clustering analysis using Dirichlet Process-Gaussian Mixture Modeling (DP-GMM)49,50,51. Each data point in the model represented a single experience sampling event and was characterized by a six-dimensional vector of within-person standardized change scores for IBI, RSA, PEP, LVET, SV, and CO. The model for each participant yielded a solution comprising clusters of increasingly smaller size, until remaining data represented clusters of single events. Rather than assign each event to only one cluster, DP-GMM assigns and weights cluster membership probabilistically. Thus, each event could be associated with more than one cluster, and contributed proportionally to calculations of weighted mean change scores for the six physiological features across multiple clusters.
Consistent with a populations hypothesis, we predicted that participants would vary in the number of clusters discovered in their data, and that these clusters would represent diverse patterns of change in physiological activity both within and across participants. However, we also predicted that we would find structure within this variation, such that we would find some patterns of change in physiological activity that broadly recurred across participants. We further predicted a many-to-many relationship between clusters (i.e., recurring patterns of change in physiological activity) and the words and affect ratings used to label events in those clusters. Conversely, support for a fingerprints hypothesis would be evidenced by similar patterns of change in physiological activity within and across participants, with each cluster consistently and specifically associated with particular emotion words or valence and/or arousal ratings.
Results
Variation in patterns of change in physiological activity within participants
We examined the number of clusters of physiological activity that accounted for at least five percent of the within-person events submitted for analysis, along with the percentage of events for which physiological data were not included in clusters of this size (i.e., those left unclustered). In total, 219 clusters of physiological activity were identified across participants, with a mean of 4.76 clusters per person (SDâ=â1.25; see additional summary statistics in supplementary Table S4). The discovered clustering solutions fit the data well; eventsâ average probability of cluster membership was high (Mâ=â0.87; SDâ=â0.06). See âMethodsâ and pages 28â30 of the supplementary information for details regarding model and parameter validation. As predicted, we observed variability in the clusters of physiological change across participants, both in their number (minimumâ=â3; maximumâ=â8) and in their nature (supplementary Table S5). Clusters associated with a greater number of events typically evidenced small physiological changes. This makes sense because in everyday life we would expect to see commonly-occurring events associated with smaller physiological changes.
Figure 1 displays the findings for two participants with an average number of clusters. Participant #1â²s data was described by four clusters of physiological activity, in contrast to participant #2â²s data, which was described by five. For participant #1, there was limited change in sympathetic nervous system (SNS) function (i.e., change in PEP) and modest decreases or no change in parasympathetic (PNS) function across the clusters (i.e., change in RSA), whereas for participant #2, physiological changes were primarily driven by changes in SNS function (especially in clusters D and E), with little overall change in PNS function (except in cluster D where there is coactivation of the SNS and PNS).
Mental features associated with patterns of change in physiological activity within participants
Because each event was probabilistically associated with multiple identified clusters per person, the freely labeled emotion words and valence and arousal ratings for each event were assigned weights describing their association with each identified cluster. As predicted, there was a many-to-many relationship between clusters of physiological activity and the words used to label events in those clusters. Within participants, the same emotion word was used to label instances associated with multiple clusters, and, correspondingly, multiple words were used to label events in each cluster. Table 2 lists the number of unique words used by participants over the course of experience sampling (minimumâ=â9, maximumâ=â147), the average number of per-person clusters associated with each emotion word, and the average number of words associated with each cluster per person.
This lack of consistency and specificity in emotion word usage is well illustrated by our example participants. Participant #1 (Fig. 2; left panel) used âtiredâ, âcalmâ, âuncomfortableâ, and âboredâ (among other words) to refer to a substantial number of events in all four clusters (i.e., no specificity), and events in each cluster were labeled with a large number of words (i.e., low consistency). The number and identity of words used to label events varies by cluster, with cluster C associated with the fewest words. Even though most of the words in this cluster are normatively low arousal experiences (e.g., âtiredâ, âboredâ), the normative valence52 of these words does vary (e.g., âcalmâ vs. âsadâ). Participant #2 (Fig. 2; right panel) likewise displays a lack of consistency and specificity in word use. This participant used far more words than participant #1, with each word used relatively infrequently, and words like âhappyâ, ârelaxedâ, âcalmâ, and âexcitedâ being used comparably across most clusters. Taken together, we see virtually no evidence of a one-to-one correspondence between physiological changes and categories of emotional experience. (Bubble plots for all participants are in supplementary Figure S4; these also show virtually no evidence of consistency or specificity.)
A similar pattern of results was observed for ratings of valence and arousal, such that there was a many-to-many relationship between clusters of physiological activity and the affective features of the events in those clusters. Across participants, grand means were large for the within-cluster ranges in valence (Mâ=â45.86, SDâ=â15.06) and arousal (Mâ=â47.16, SDâ=â13.90), suggesting that these affective features varied substantially within clusters. Within participants, separate tests for valence and arousal ratings confirmed that affective features overlapped across clusters. Only two participants evidenced significant differences in valence between one or more clusters (4.35% of KruskalâWallis one-way analysis of variance H values were significant at αâ=â0.05), and only four participants evidenced between-cluster differences in arousal (8.70% of H values were significant at αâ=â0.05). These results were not corrected for multiple comparisons, and do not exceed, or barely exceed, the number of significant differences we would expect to find by chance alone.
The lack of consistency and specificity in affective features again can be seen in our example participants. For participant #1 (Fig. 3, left panel), experience sampling events from all clusters are densely overlapped (i.e., low specificity) and associated with both positive and negative valence ratings, as well as both high and low arousal ratings (i.e., low consistency). These same observations hold for participant #2 (Fig. 3, right panel).
For comparison, we also assessed between-cluster differences in contextual features capturing type of activity, time of day, and posture. Based on previous literature, we predicted that posture (e.g., the likelihood of being seated during an event) would be more strongly related to differences in peripheral physiology53. As expected, we observed a greater number of significant differences in posture across clusters, and found that average effect sizes for posture comparisons were significantly larger than for those for other comparisons, including affect. See page 18 of the supplementary information and supplementary Table S6 for details.
Common patterns of change in physiological activity across participants
To directly compare patterns of physiological activity across participants, we derived within-person effect size estimates as the mean change relative to the standard deviation for each physiological feature. Following published recommendations54, we then classified these effect sizes as indicating ânegligible changeâ, âdecreaseâ, or âincreaseâ in each feature. Using these classifications, we observed 44 broad patterns of change in physiological activity that occurred at least twice across participants, accounting for 145 or 66.21% of the 219 total clusters (see supplementary Table S7; occurrences of these patterns are also identified in supplementary Table S5â²s cluster-level results). These findings suggest that clusters were not completely idiosyncratic, as is expected given the influence of physical context (i.e., most experience sampling prompts occurred while participants were sitting) as well as the constraints of human physiology (e.g., measures of cardiovascular function are not completely independent55). A subset of 10 patterns occurred at least five times across participants, accounting for 65 or 29.68% of the 219 total clusters, as listed in Table 3.
We also examined the emotion words used to describe events in the 65 clusters representing these common patterns. Comparing emotion labels across participants in this way allowed us to assess whether mental features might be more consistently and specifically associated with more commonly recurring patterns of change in physiological activity. We found that the same words were used to label clusters associated with different physiological patterns (low specificity) and that a wide range of emotion words was used to label clusters associated with the same pattern (low consistency). Moreover, the words associated with a given pattern varied in terms of both normative valence (Fig. 4, left panel) and arousal (Fig. 4, right panel). The most frequent pattern (Pattern 1), for example, is associated with âhappyâ, âcalmâ, âannoyedâ, and âboredâ. These findings replicate and extend the within-participant analysis of emotion words, and illustrate the scale of variation in affective experience across participants, even when the patterns of change in physiological activity for those experiences are largely similar.
We assessed all identified clusters according to whether they represented changes in each branch of the ANS, using RSA change scores as an estimate of PNS function57, and PEP change scores as an inverse estimate of SNS function53,55,58 (though it should be noted that RSA is a chronotropic measure and PEP is an inotropic measure of cardiac activity). Changes in PNS function, whether withdrawal or increased activation, were present in about 74% of clusters (withdrawal: 45.21%; activation: 29.38%), whereas changes in SNS function were present in about 55% of clusters (withdrawal: 20.55%; activation: 34.25%). A chi-square test of independence confirmed that the proportion of changes was significantly different between the two branches, X2(2)â=â33.54, pâ<â0.001, suggesting that changes in PNS activity may be a larger contributor to cardiovascular function in daily life. These findings are consistent with prior research on the autonomic space model of cardiovascular function, which found that changes in PNS function alter heart rate over a much larger dynamic range than do changes in SNS function59. The larger proportion of PNS changes also might be explained by the fact that all events occurred in the absence of recent gross movement or posture change, and the majority of events in each cluster occurred when participants were sitting (73.88%) versus standing (12.47%) or reclining (8.27%). Postural changes strongly influence autonomic activity53, and SNS contributions to changes in heart rate are greater during physical activity60,61.
Figure 5 plots clusters according to changes in PNS and SNS activity, to illustrate the potential mode of autonomic cardiac control represented by each59. Across all identified clusters, 50.23% represented an uncoupled control mode, such that there was a change in activity of only one branch of the ANS. ANS modes that were reciprocal (i.e., activity in one branch increased while activity in the other decreased) or coactivational (i.e., activity in both branches increased or decreased) occurred in 19.63% and 20.09% of the clusters, respectively. Negligible change in the activity of both branches accounted for 10.05% of clusters across participants (note that, because these classifications were determined using within-person effect size estimates, ânegligible changeâ merely indicates that changes in PNS or SNS activity were small relative to other changes for that participant). This distribution of clusters across different modes of ANS control is particularly interesting in light of the pervasive assumption that the PNS and SNS are inherently reciprocally organized and functionally antagonistic59. The occurrence of equivalent proportions of clusters associated with coactivational and reciprocal modes may be reflective of the high ecological validity of the present study and its ability to capture changes in everyday life.
Discussion
Many of the ongoing debates about emotion center on questions about the amount of variation in the physical and mental features of emotional events, as well as the nature of that variation. To date, studies have typically been designed to observe consistency and specificity in physiological activity across instances of an emotion categoryâyet have still observed more variation than expected10,23,28,29,30,31,32,33,34,35,36. By using a novel biologically-triggered experience sampling method, the present study builds on a growing body of work to show that empirical approach matters: when studies are designed to observe variation, they find heterogeneity within emotion categories, in addition to similarity across emotion categories23,40,62. Here, participants reported affective experience while their peripheral physiological activity was recorded, specifically during instances of substantial change in cardiac activity in the absence of movement. This method enhanced the likelihood of sampling instances having greater psychological salience, and thus enhanced the chances of observing heterogeneity in emotional experience. Unsupervised clustering analyses revealed variation in the number and nature of patterns of physiological activity that recurred within individuals over time, as well as heterogeneity in the affect ratings and emotion labels associated with each pattern. These findings suggest that physiological and mental features are non-redundant measures of emotional episodes, and may help explain the persistent difficulty in observing robust, stable correlations between self-reports and physiological measures63. Further, the degree of variation that we observed is inconsistent with classical and prototype views of emotion1,3,7. We did not find evidence to support a fingerprints hypothesis. No emotion-specific diagnostic patterns of physiology emerged. Rather, patterns of physiological activity varied in context-dependent ways. These findings support a populations hypothesis positing structure within variation, both within and across participants10,17.
Although a number of theoretical views on emotion are more or less consistent with a populations hypothesis14,15,16,17,18,19,20, it has been most explicitly developed within the theory of constructed emotion17,18. The theory of constructed emotion begins with the observation that the primary purpose of the brain is to predictively regulate physiological resources to most efficiently manage the brainâs and bodyâs energy needsâa process referred to as allostasis64,65. To do this, the brain uses prior experience to build and maintain an internal model of the world that includes its own body66,67,68,69. In a given instance, the brain is hypothesized to construct, ad hoc, a population of temporally dynamic representations that are similar to the current sensory array. Relying on similarity (as opposed to exact matching) is functional, because the representations can vary in their physical and mental features across instances. For example, different patterns of physiological activity are required to support shouting and flailing ones arms versus silently brooding in anger, and only one of these is likely to be functional within a particular anger-relevant context (e.g., being reprimanded by a supervisor). These ad-hoc, in-the-moment representations function as prediction signals that begin as visceromotor66,67,68,69 and motor control signals70,71. Efferent copies of these signals give rise to the predicted sensations that will result from those movements as well as hypotheses about the causes of those sensations17,18. In turn, this entire cascade of sensory and (viscero)motor detail is compressed into a low-dimensional, multi-modal abstraction, which the brain can use as needed to issue future situation-specific predictions. As this abstraction process operates over time, a wide range of highly variable affective and physical experiences become functionally associated with each other, forming a population of instances within the same emotion category (e.g., anger). Even within a single individual, these instances will vary considerably in their physiological features because they are coupled to the energetic needs of each specific situation. In this way, variation in patterns of physiological activity across emotion instances is inherently meaningful. Further, the theory of constructed emotion posits that a particular pattern of physiological activity can be associated with mental features that vary substantially (e.g., valence and arousal, evaluations of novelty or controllability), such that these features are neither consistent with nor specific to any particular emotion category.
Critically, the variability predicted by the theory of constructed emotion is thought to be situated, not random15,16,17,18. This is because peripheral physiological changes are yoked to the energy needs of current or anticipated action64,72,73,74. For example, actions such as freezing and laughing in fear may each be supported by a distinct pattern of physiological activity (resulting in within-category and within-participant variation), whereas smiling in anger and happiness, for example, may share more physiological features (resulting in across-category and across-participant similarity). In the present study, we observed across many participants that a larger proportion of experience sampling events were associated with smaller physiological changes. We also observed some consistency in patterns of change in physiological activity across individuals. These observations suggest that more consistent patterns can occur, but that revealing the sources of these consistencies will require future research to assess a number of situational features, including social and motivational context. Such factors have been shown to influence physiological activity, even for two experiences labeled with the same emotion word75,76. Likewise, future research should further explore how physical bodily contexts, such as activity and posture, play a role in constructing affective experience in everyday life77,78,79.
To more fully understand the situation-specific construction of affective experience, scientists must conduct research that extends beyond the narrow confines of the psychology laboratory to sample a broad range of human activities and experiences39,80. In the present study, we used physiologically-triggered experience sampling that allowed us to observe a much broader range of ecologically-valid affective experiences and their accompanying physiological activity. Using these data, we identified patterns of change in physiological activity that represented multiple modes of autonomic control (reciprocal, coactivation, coinhibition, and uncoupled modes58,59), illustrating the utility of this method for documenting patterns of physiological activity that may be less commonly observed in a controlled laboratory setting. This proof-of-concept study can be expanded in future research to include other measured features, such as more complex motion sensing, self-reported bodily sensations, geolocation data, or additional ambulatory physiological measures. In doing so, future research should also undertake time-dynamic analyses that can account for diurnal trends in physiological activity.
The present study also utilized innovative analytical methods for detecting patterns of change in physiological activity and their association with mental features of affective experience. To date, most studies have used supervised machine learning (i.e., classification) algorithms to deductively assess the relationship between patterns of physiological activity and categories of emotional experience81,82,83. Those using these classification-based approaches often interpret the patterns they identify for each emotion category as a fingerprint of that category81,84,85,86, rather than as the statistical summary that such patterns actually represent36. In contrast, we employed person-specific unsupervised machine learning (i.e., clustering) algorithms to discover the number and nature of patterns of change in physiological activity, and then inductively examined the relationships between these patterns and emotional experience. Our group has recently used similar approaches to understand the relationships between change in physiological activity and subjective experience in the more constrained context of a motivated performance laboratory task to assess physiological patterns of challenge or threat48,87. The current work demonstrates that this method not only can be used on data collected outside the laboratory, but affords an inductive, data-driven approach that enables the discovery of variability, rather than stipulating experimenter-labelled categories as ground truth.
As with any study, there are limitations and constraints on generalization. One limitation is that we triggered experience sampling prompts solely on the basis of IBI change, which itself can be constrained, particularly during coactivation or coinhibition modes of ANS control59. As a result, we may have excluded from sampling some moments of significant psychological salience or of change in physiological features other than IBI. Notwithstanding, the approach here provides a better chance of sampling experience at times salient to the individual than does traditional random sampling. Likewise, by sampling only in the absence of movement, we potentially reduced the number of evocative emotional events available for study. Although this decision was necessary to prevent triggering experience sampling prompts that could have been solely due to movement, future research should consider means of addressing this limitation. For example, future studies could utilize separate triggers for movement- and non-movement-related prompts and assess which, if any, movement-related prompts are also psychologically salient. The homogeneity of our sample and sampling conditions can also be considered a constraint on generalization: most of our participants were university students, whose physiology and self-reports were recorded during weekdays, often while they were in class or working at a desk. At the same time, this limitation gives the advantage to a fingerprints hypothesis, as more constrained sampling should result, if anything, in greater consistency. Further, the relatively small sample size makes it difficult to perform robust between-participants analyses, which are needed to more fully assess the extent of similarities and differences across individuals and to assess the role of theoretically and practically meaningful individual differences (e.g., age, sex, emotional granularity, interoceptive sensitivity) in explaining variability across individuals.
Despite these limitations, however, the current approach opens the empirical door for innovations in emotion research. For over a hundred years, scientists have induced instances of emotion and measured the changes in physiological activity to those inductions88 in an attempt to identify the defining biological features of emotion categories. With the empirical approach adopted in the present study, we are finally heeding the advice of William James, who wrote âthere is no limit to the number of possible different emotions which may exist, and why the emotions of different individuals may vary indefinitely, both as to their constitution and as to objects which call them forthâ (p. 454)89. We now have the requisite methodological and analytical tools to better understand the individual- and situation-specific physiological and mental constituents of affective experience in everyday life.
Methods
All experimental protocols described below were approved by the Northeastern University Institutional Review Board. These methods were carried out in accordance with the relevant guidelines and regulations for research with human subjects.
Participants
Sixty-seven participants ranging in age from 18â36Â years (55% female; 38.8% White, 3.0% Black, 29.8% Asian, 28.4% other; Mageâ=â22.8Â years, SDageâ=â4.4Â years) were recruited from the greater Boston area through posted advertisements, as well as Northeastern University classrooms and online portals. Eligible participants were non-smoking, fluent English-speakers, and excluded if they had a history of cardiovascular illness or stroke, chronic medical conditions, mental illness, asthma, skin allergies, or very sensitive skin. Eligible participants also confirmed they were not taking any medications known to influence physiological arousal including medications for ADHD, insomnia, anxiety, hypertension, rheumatoid arthritis, epilepsy/seizures, cold/flu, or fever/allergies. Informed consent was obtained from all participants before beginning the study. Participants received $490 as compensation for completing all parts of the study, plus up to an additional $55 in compliance and task incentives, as detailed on page 4 of the supplementary information.
Of the 67 recruited participants, six withdrew and an additional nine were dismissed due to poor compliance. A total of 52 participants completed the full protocol (56% female; 38.5% White, 1.9% Black, 38.5% Asian, 21.2% other; Mageâ=â22.5Â years, SDageâ=â4.4Â years). Six of these participants were excluded from the main analyses due to a lack of sufficient usable physiological data, bringing the final sample size to 46 (50% female; 43.5% White, 0% Black, 41.3% Asian, 15.2% other; Mageâ=â22.4Â years, SDageâ=â4.4Â years).
Procedure Overview
Each participant completed approximately 14Â days (Mâ=â14.4, SDâ=â0.6) of context-aware experience sampling distributed across a three- to four-week period (Mâ=â24.9Â days, SDâ=â5.5Â days). On each day of experience sampling, participants came into the lab and were instrumented for peripheral physiological recordings. Participants could not begin the experience sampling portion of the study without functioning ambulatory equipment, and so did not complete the daily protocol if they did not attend the morning session to be instrumented for the recordings. In the event a participant was unable to make a scheduled session or the equipment was not functioning properly during instrumentation, they were rescheduled to complete the day of sampling on another day. Participants generally started each experience sampling day between 8 and 9 am, although this varied between 7:30am and 2:30Â pm according to participantsâ schedules. Participants were provided an Android smartphone with a custom application, MESA (MindWare Technologies LTD, Westerville, OH), which was used to determine when a heart rate response of sufficient duration and amplitude occurred which then triggered an experience sampling prompt. Participants were instructed to continue physiological recordings for eight hours each day, after which they were able to remove and recharge all equipment. Upon completing experience sampling each day, participants automatically received an end-of-day survey via SurveyMonkey (San Mateo, CA), which they used to provide additional details about the prompts they completed throughout the day. Before and after the two-week experience sampling protocol, participants also completed two in-lab sessions. In each session, participants completed tasks and questionnaires that are not reported here (see page 5 of supplementary information and supplementary Table S1 for an overview).
Physiological measurement
All ambulatory peripheral physiological measures were recorded at 500Â Hz on a mobile impedance cardiograph from MindWare Technologies LTD (Model # 50-2,303-02, Westerville, OH), which participants wore clipped onto their clothing on the hip.
ECG and ICG were obtained using pre-gelled ConMed (Westborough, MA) Cleartrace Ag/AgCl sensors, connected via wires to the cardiograph. Sensor sites were cleaned with alcohol and abraded lightly with gauze. ECG was obtained using a modified lead II configuration, with recording electrodes placed on the distal right collarbone and an inferior left rib, respectively, and a reference electrode placed on an inferior right rib. The ECG signal was acquired using a low cutoff of 0.5Â Hz and a high cutoff of 45Â Hz.
ICG was obtained using a four-spot electrode configuration90. Two inner recording electrodes were placed on the front of the torso: one at the base of the neck at the top of the sternum, and a second at the bottom of the sternum over the xiphisternal junction. Two outer source electrodes were placed on the back along the midline approximately 4Â cm above and below the inner recording electrodes, respectively. The source electrodes passed a 4Â mA, 100Â kHz alternating current across the thorax. The distance between the inner recording electrodes (cm) was measured for each day of experience sampling. Basal impedance (Z0) was acquired using a low cutoff of 10Â Hz. The first derivative, dZ/dt, was acquired using a low cutoff of 0.5Â Hz and a high cutoff of 45Â Hz.
Electrodermal activity (EDA) was also recorded over an unobtrusive neck placement, but data were not used in the present analysis due to difficulties in detecting and removing artifacts, as well as generally low mean skin conductance levels across participants.
The mobile impedance cardiograph collected continuous three-axis accelerometry data used to assess movement. Additionally, participants wore two inertial measurement units (IMUs) purchased from LP-Research (Minato-ky, Tokyo, Japan) to derive measures of posture and changes in posture. One IMU was placed medially on the sternum beneath the top inner impedance recording electrode and affixed to the skin using a double-sided adhesive patch. The other IMU was placed on the front of the thigh using either a cloth holder attached to the clothing, or a second adhesive patch affixed to the skin. Participants did not remove sensors until the end of each experience sampling day, unless instructed by the experimenters (e.g., due to synchronizing issues).
Context-aware experience sampling
Peripheral physiological data and accelerometric data were recorded continuously throughout the 8-h sampling period and communicated via Bluetooth to a Motorola Moto G4 smartphone. The smartphone application, MESA, processed the continuous ECG and accelerometer data in real time, and initiated an experience sampling prompt anytime a substantial, sustained change in heart period was detected in the absence of movement or posture change, with an imposed minimum interval of five minutes between prompts. A substantial, sustained change in heart period was operationalized on the first day of sampling as occurring when the interbeat interval (IBI) changed by more thanâ±â167 ms over an eight-second period (at a typical resting heart rate of 60 bpm or IBI of 1,000 ms, this is equivalent to a decrease of about 9 bpm or an increase of about 12 bpm). On subsequent days, this IBI parameter was manually adjusted up or down to ensure each participant received approximately 20 prompts per day. The average adjustment wasâ±â26.17 ms (SDâ=â25.91 ms), with the final parameter setting ranging betweenâ±â65 ms andâ±â266 ms across participants. For background on threshold selection and more information on threshold adjustment, see page 4 of the supplementary information.
Movement was determined from the continuous accelerometer data from the mobile impedance cardiograph. Minimal movement was operationalized as any time none of the three accelerometry channels (alone or in aggregate) exceeded a threshold of 10Â cm/s2 within the preceding 30Â s. Posture (standing, sitting, reclining) was determined by comparing the relative orientation of the two IMUs on participantsâ torso and thigh using their continuous accelerometry data. Absence of posture change was operationalized as any time when the relative orientation of the two IMUs did not change within the preceding 30Â s.
On average, participants received 21.70 (SDâ=â6.90) prompts per day. We observed that prompts were relatively evenly distributed throughout the day within experience sampling days. On average, participants completed 34% of prompts in the morning (before 12Â pm), 52% of prompts in the afternoon (12â5Â pm), and 14% of prompts in the evening (after 5Â pm). This makes sense given the typical start times for participants in our sample. Participants also received on average two ârandomâ prompts per experience sampling day, which occurred in the absence of movement or posture change, but which were not contingent on a change in IBI. Random prompts were spread throughout the experience sampling day, such that they could only receive one in the first four hours and one in the second four hours of experience sampling. Participants were informed that some experience sampling prompts would be generated randomly while others would be generated based on changes in their cardiac activity. By conveying this information, we were able to instruct participants to avoid responding to prompts following specific physiological events (e.g., coughing, sneezing) and minimize the extent to which they paid special attention to their cardiac activity.
To remain in the study, participants were required to complete at least three prompts each day. In addition, as detailed on page 4 of the supplementary information, for the purposes of incentivizing participation and limiting attrition, the experience sampling protocol was broken into three pay periods (days 1â5, 6â10, and 11â14). Participants were required to respond to an average of at least six prompts per day during each pay period to remain in the study, and received a bonus payment for each pay period where they completed an average of eight prompts per day. The average number of completed prompts was 8.80 per day (SDâ=â1.22), which is in line with previous experience sampling studies that have asked participants to complete 10 prompts per day91,92.
At each sampling event, participants were prompted to respond to a series of questions presented in the MESA application. Participants first provided a brief free-text description of what was happening at the time they received the prompt. Next, participants rated their current valence and arousal, each on a 100-point continuous slider scale ranging from -50 (very unpleasant or very deactivated) toâ+â50 (very pleasant or very activated). Participants then self-generated words to label their current affective experience. Specifically, participants were asked to âlist any emotion(s) you were feeling when you received the promptâ. Participants were able to provide as many words as they felt necessary to describe their affective experience but were required to input at least one word. For each self-generated word, participants were asked to provide an intensity rating on a five-point scale: ânot at allâ (1), âa littleâ (2), âmoderatelyâ (3), âa lotâ (4), âvery muchâ (5). We asked participants to self-generate emotion words as this allowed us to capture naturally-occurring variation in how participants categorized their affective experience, as opposed to having them select words from an experimenter-stipulated list, which would necessarily constrain variation. Participants also responded to additional questions that are not included in the present report (see page 4 of the supplementary information for details).
At the end of each experience sampling day, participants received a modified day reconstruction survey93, in which they were presented with their self-generated brief description and social context from each prompt they completed during the day. For each prompt, participants were asked to provide additional details about the event, its social context, and accompanying affective features. Participants were requested to select three sampling events for which they provided a longer, more detailed description (>â200 words). Full details of the end-of-day survey are provided on pages 5â6 of the supplementary information. Data from these surveys are not reported here.
Physiological signal processing and feature extraction
Peripheral physiological signals were processed using an in-house pipeline coded in Python to accommodate the volume of physiological data collected (400â+âhours), as well as the variability in signal morphologies and artifacts produced during long-term ambulatory monitoring. Artifacts were identified using a series of quality checks, detailed for each signal below. Data affected by artifacts were excluded from analysis (i.e., artifacts were not corrected). This resulted in an average event-related data loss of 8.18% per day (SDâ=â3.35%) across ECG (5.79% average data loss; SDâ=â3.22%) and ICG (2.56% average daily data loss; SDâ=â2.03%) signals. The research team, including two trained psychophysiologists (J.B.W., K.S.Q.), met regularly throughout development of the processing pipeline to visually examine sample output and provide feedback on quality checks.
Electrocardiogram (ECG)
The ECG signal was processed following previous work94. First, the raw signal was passed through an elliptic bandpass filter to remove baseline and high frequency noise without affecting the waveform. Initial quality checks were then performed for each beat, checking for overall waveform shape, and acceptable minimum, maximum, and minimum-to-maximum values94. The default parameters, originally selected based on data collected in laboratory settings, were found to be too aggressive for this ambulatory data set, so we relaxed these thresholds based on manual review of the results by experts (J.B.W., K.S.Q.) for a random subset of the data. The parameters we implemented are reported in supplementary Table S2, along with the original recommendations94.
R-peak detection for ECG was performed using established methods95 and implemented using the BioSPPy package in Python96. Interbeat interval (IBI) was then derived as the average R-R interval. Additional quality checks (supplementary Table S2) were performed on each IBI series to ensure that values were within acceptable ranges (300â2000Â ms), and that expected beat-to-beat differences were consistent with normal beats and unlikely to be artifacts (following established benchmarks97). ECG data failing any quality check were excluded from analysis.
Respiratory sinus arrhythmia (RSA) was derived from the IBI series. RSA reflects high-frequency variability in IBI which occurs at the respiratory frequency and is often used as an estimate of parasympathetic influence on the heart57. RSA calculations were coded to mimic the processing steps of standard heart rate variability (HRV) analysis software (MindWare Technologies LTD, Westerville, Ohio), including: cubic interpolation of beat-to-beat IBI, detrending to minimize non-stationarity, tapering using a Hamming window, and lastly, fast Fourier transformation (FFT). RSA was calculated as the log of the area under the power spectrogram that lies between 0.12 and 0.4Â Hz.
Impedance cardiogram (ICG)
ICG feature detection was performed using an in-house implementation based on methods described in previous work94,98, and adjusted for ambulatory data using thresholds determined by comparing performance of the original algorithms with hand-scored impedance cardiographic data on a subset of participants. The first derivative of the basal impedance (Z0) signal, dZ/dt, was used as the basis for all analyses. The signal was segmented into time windows corresponding to 250Â ms before the ECG R-peak to 500Â ms after; eight such segments (i.e., 8 beats) were averaged together to form overlapping ensembles94. B points were detected in each ensemble by taking the first and second derivatives of the dZ/dt signal and comparing them with thresholds based on signal frequency (supplementary Table S2)98. Following standard procedures98, forward and reverse autoregressive modeling was then used to perform outlier detection and correction, such that all B points in a series were fit with the autoregressive model. X points were detected by examining the second derivative of the dZ/dt signal within each ensemble94. Segments of the ICG signal from which we could not detect B or X points and ICG data that corresponded to periods of unusable ECG data were excluded from analysis.
Two features we used were systolic time intervals, namely, pre-ejection period (PEP) and left ventricular ejection time (LVET). PEP represents the time interval in ms between the electrical event that signals the start of ventricular contraction to the opening of the aortic valve. PEP is often used as an inverse estimate of sympathetic control of the heart53,55,58. Here, PEP was calculated as the time in ms between the ECG R peak and the ICG B point (also referred to as PEPR99). LVET is the time interval in ms between the opening and closing of the aortic valve55, and was calculated as the time in ms between the ICG B point and X point. Quality checks (supplementary Table S2) were performed and we retained only values that occurred within an acceptable range (30â200Â ms for PEP, 100â500Â ms for LVET), and that did not result in changes in the gradient greater than 30Â ms from one ensemble to the next.
Two additional features we used were based on volumetric measures of cardiac function, namely, stroke volume (SV) and cardiac output (CO). SV represents the volume of blood ejected by the heart with each heartbeat (in mL), and here we calculated SV using Kubicekâs equation100. CO is the volume of blood circulated by the heart in L/min55.
Analysis
Within-person clustering analysis
We submitted data from each individual to a separate clustering analysis using Dirichlet Process-Gaussian Mixture Modeling (DP-GMM)49, 50. DP-GMM is a specialized variant of Gaussian Mixture Modeling (GMM)49. Data points (i.e., six-dimensional vectors of change scores for IBI, RSA, PEP, LVET, SV, and CO) were standardized prior to clustering by subtracting the mean and dividing by the standard deviation of each feature. Specific parameter values are reported in supplementary Table S3.
The first step in implementing standard GMM is to initialize a set of parameters for a predefined number of clusters represented in the data. Parameters (i.e., values) include each clusterâs location (i.e., mean), shape (i.e., covariance), and relative size (i.e., the mixture proportion or the prior probability of a point belonging to that cluster relative to others). Initially, these cluster parameters are chosen at random, with means and covariances selected from a Gaussian distribution that is proportional to the data, and with prior probabilities of equal value. GMM uses these initialized values to calculate the posterior probabilities of each point belonging to each cluster. Under mixture models, each point is given a posterior probability of belonging to every cluster; however, if patterns in the data are distinctive, posterior probabilities during later iterations will often be close to 1 for a single cluster and negligible for all others. Based on the posterior probabilities, the means and covariance matrices are updated iteratively using optimization techniques such as Expectation Maximization49 or Variational Inference49 until reaching convergence (i.e., until the values stop changing).
In the present analyses, we sought to avoid tuning and selecting the number of clusters for each participant. DP-GMM allowed us to accomplish this goal as, in contrast to standard GMM, DP-GMM discovers the number of clusters from input data through the use of a Dirichlet Process Prior101. The Dirichlet Process Prior enables a dataset to include an unknown (i.e., infinite) number of clusters, where the mixture proportion of each cluster is randomly assigned a value from a continuous probability distribution such that the mixture proportions must sum to 1. This random assignment process is repeated ad infinitum to create a vector of mixing proportions which then serves as the prior probabilities of a standard GMM. To place reasonable bounds on this process, the number of clusters made discoverable through Variational Inference was limited to the number of data points in the sample. Because DP-GMM iteratively optimizes parameters (means, covariances, and mixture proportions), clusters which initially acquire more points become more likely to continue acquiring points in the future. Consequently, as the model runs, some clusters grow while others shrink and are ultimately discarded.
As a fully Bayesian method, DP-GMM is more likely to avoid overfitting because it is self-regularized by the priors101. There are, however, model parameters that can constrain the number of clusters discovered. As detailed on page 28 of the supplementary information, we ran additional analyses to validate the values selected for these parameters. We also conducted parallel analyses on the example participants using a standard GMM approach, in which Bayesian information criterion (BIC) values were used to select the optimal number of clusters per participant. These analyses replicated our main findings with regard to the heterogeneity of clusters discovered and their relationship with mental features of affective experience (see pages 29â30 of the supplementary information for details).
Effect size estimates
We derived effect size estimates (similar to Cohenâs d56) for the mean change scores for IBI, RSA, PEP, LVET, SV, and CO for each cluster for each participant using the following formula:
where m is the weighted cluster mean for a given physiological feature, µ is the assumed population mean corresponding to no change (i.e., 0), and Ï is the standard deviation for the given physiological feature for that participant across all events included in the analysis. According to published recommendations for interpreting effect sizes54, we interpreted the magnitude of the weighted mean change score for each feature for each cluster as: |d|<â0.2 (negligible change); 0.2ââ¤|d|<â0.5 (small change); 0.5ââ¤|d|<â0.8 (medium change); |d|â¥â0.8 (large change).
Common patterns of change in physiological activity
We classified change score effect sizes as a negligible change (|d|<â0.2), decrease (dââ¤ââ0.2), or increase (dââ¥â0.2). Before classification, we rounded the effect size estimates to the nearest tenth (i.e., one decimal place) to facilitate comparison with published recommendations. Logical syntax was used to identify duplicate (i.e., repeated) strings of effect size classifications.
We performed further classifications of RSA and PEP change scores to assess overall trends in the change in PNS and SNS function, respectively. Larger change score effect sizes for RSA (i.e., greater variability in interbeat interval within the respiratory frequency band) indicate increased PNS function, such that: |d|<â0.2 (negligible PNS change); dââ¥â0.2 (PNS activation); and dââ¤ââ0.2 (PNS withdrawal). For PEP, larger change score effect sizes (i.e., longer pre-ejection periods) indicate SNS withdrawal, such that: |d|<â0.2 (negligible SNS change); dââ¥â0.2 (SNS withdrawal); and dââ¤ââ0.2 (SNS activation).
Change history
26 January 2021
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
References
Ekman, P. & Cordaro, D. T. What is meant by calling emotions basic. Emotion Rev. 3, 364â370 (2011).
Lench, H. C., Flores, S. A. & Bench, S. W. Discrete emotions predict changes in cognition, judgment, experience, behavior, and physiology: a meta-analysis of experimental emotion elicitations. Psychol. Bull. 137, 834â855 (2011).
Levenson, R. W. Basic emotion questions. Emotion Rev. 3, 379â386. https://doi.org/10.1177/1754073911410743 (2011).
Friedman, B. H. Feelings and the body: the Jamesian perspective on autonomic specificity of emotion. Biol. Psychol. 84, 383â393 (2010).
Levenson, R. W. The autonomic nervous system and emotion. Emotion Rev. 6, 100â112 (2014).
Roseman, I. J. A model of appraisal in the emotion system. In Appraisal processes in emotion: theory, methods, research (eds Scherer, K. R. et al.) 68â91 (Oxford University Press, Oxford, 2001).
Cowen, A. S. & Keltner, D. Self-report captures 27 distinct categories of emotion bridged by continuous gradients. Proc. Natl. Acad. Sci. 114, E7900âE7909 (2017).
Shaver, P., Schwartz, J., Kirson, D. & OâConnor, C. Emotion knowledge: Further exploration of a prototype approach. J. Pers. Soc. Psychol. 52, 1061â1086 (1987).
Scherer, K. R. On the nature and function of emotion: a component process approach. In Approaches to Emotion (eds Scherer, K. R. & Ekman, P.) 293â317 (Lawrence Erlbaum Associates, New York, 1984).
Siegel, E. H. et al. Emotion fingerprints or emotion populations? A meta-analytic investigation of autonomic features of emotion categories. Psychol. Bull. 144, 343â393. https://doi.org/10.1037/bul0000128 (2018).
Ekman, P., Levenson, R. W. & Friesen, W. V. Autonomic nervous system activity distinguishes among emotions. Science 221, 1208â1210 (1983).
Cacioppo, J. & Berntson, G. Relationship between attitudes and evaluative space: a critical review, with emphasis on the separability of positive and negative substrates. Psychol. Bull. 115, 401â423 (1994).
Russell, J. & Carroll, J. On the bipolarity of positive and negative affect. Psychol. Bull. 125, 3â30 (1999).
Adolphs, R. & Anderson, D. The neurobiology of emotion: a new synthesis (Princeton University Press, Princeton, 2018).
Barrett, L. F. Are emotions natural kinds?. Perspect. Psychol. Sci. 1, 28â58 (2006).
Barrett, L. F. Psychological construction: the Darwinian approach to the science of emotion. Emotion Rev. 5, 379â389 (2013).
Barrett, L. F. How emotions are made: the secret life the brain and what it means for your health, the law, and human nature (Houghton, Houghton Mifflin Harcourt, 2017).
Barrett, L. F. The theory of constructed emotion: an active inference account of interoception and categorization. Soc. Cognit. Affective Neurosci. https://doi.org/10.1093/scan/nsw154 (2017).
Ortony, A. & Clore, G. Can an appraisal model be compatible with psychological constructionism? In The psychological construction of emotion (eds Barrett, L. F. & Russell, J. A.) 305â333 (Guilford Press, Guilford, 2015).
Fridlund, A. J. The behavioral ecology view of facial displays, 25 years later. In The science of facial expression (eds Fernández-Dols, J.-M. & Russell, J. A.) 77â92 (Oxford University Press, Oxford, 2017).
Darwin, C. On the Origin of Species. (Harvard University Press, Harvard, 1859/2001).
Mayr, E. What makes biology unique? Considerations on the autonomy of a scientific discipline (Cambridge University Press, Cambridge, 2004).
Wilson-Mendenhall, C. D., Barrett, L. F. & Barsalou, L. W. Variety in emotional life: within-category typicality of emotional experiences is associated with neural activity in large-scale brain networks. Soc. Cognit. Affect. Neurosci. 10, 62â71. https://doi.org/10.1093/scan/nsu037 (2015).
Harmon-Jones, E., Harmon-Jones, C., Abramson, L. & Peterson, C. K. PANAS positive activation is associated with anger. Emotion 9, 183â196. https://doi.org/10.1037/a0014959 (2009).
Barrett, L. F. & Bliss-Moreau, E. Affect as a psychological primitive. Adv. Exp. Soc. Psychol. 41, 167â218. https://doi.org/10.1016/S0065-2601(08)00404-8 (2009).
Lindquist, K. A., Satpute, A. B., Wager, T. D., Weber, J. & Barrett, L. F. The brain basis of positive and negative affect: Evidence from a meta-analysis of the human neuroimaging literature. Cereb. Cortex 26, 1910â1922. https://doi.org/10.1093/cercor/bhv001 (2015).
Satpute, A., Kragel, P., Barrett, L., Wager, T. & Bianciardi, M. Deconstructing arousal into wakeful, autonomic and affective varieties. Neurosci. Lett. 693, 19â28 (2019).
Cacioppo, J. T., Berntson, G. G., Larsen, J. H., Poehlmann, K. M. & Ito, T. A. The psychophysiology of emotion. In The handbook of emotion (eds Lewis, R. & Haviland-Jones, J. M.) 173â191 (Guilford Press, Guilford, 2000).
Stemmler, G. Physiological processes during emotion. In The regulation of emotion (eds Phillipot, P. & Feldman, R. S.) 33â70 (Erlbaum, London, 2004).
Barrett, L. F., Adolphs, R., Marsella, S., Martinez, A. & Pollak, S. Emotional expressions reconsidered: challenges to inferring emotion in human facial movements. Psychol. Sci. Public Interest 20, 1â68 (2019).
Guillory, S. A. & Bujarski, K. A. Exploring emotions using invasive methods: review of 60 years of human intracranial electrophysiology. Soc. Cognit. Affect. Neurosci. 9, 1880â1889. https://doi.org/10.1093/scan/nsu002 (2014).
Lindquist, K. A., Wager, T. D., Kober, H., Bliss-Moreau, E. & Barrett, L. F. The brain basis of emotion: a meta-analytic review. Behav. Brain Sci. 35, 121â143. https://doi.org/10.1017/S0140525X11000446 (2012).
Satpute, A. B. & Lindquist, K. A. The Default Mode Networkâs role in discrete emotion. Trends Cognit. Sci. 23, 851â864 (2019).
Raz, G. et al. Functional connectivity dynamics during film viewing reveal common networks for different emotional experiences. Cognit. Affect. Behav. Neurosci. 16, 709â723 (2016).
Touroutoglou, A., Lindquist, K. A., Dickerson, B. C. & Barrett, L. F. Intrinsic connectivity in the human brain does not reveal networks for âbasicâ emotions. Soc. Cognit. Affect. Neurosci. 10, 1257â1265. https://doi.org/10.1093/scan/nsv013 (2015).
Clark-Polner, E., Johnson, T. D. & Barrett, L. F. Multivoxel pattern analysis does not provide evidence to support the existence of basic emotions. Cereb. Cortex 27, 1944â1948 (2017).
Durkee, P. K., Lukaszewski, A. W. & Buss, D. M. Pride and shame: Key components of a culturally universal status management system. Evol. Hum. Behav. 40, 470â478 (2019).
Keltner, D., Tracy, J. L., Sauter, D. & Cowen, A. What basic emotion theory really says for the twenty-first century study of emotion. J. Nonverbal Behav. 43, 195â201 (2019).
Wilhelm, F. H. & Grossman, P. Emotions beyond the laboratory: Theoretical fundaments, study design, and analytic strategies for advanced ambulatory assessment. Biol. Psychol. 84, 552â569 (2010).
Kuppens, P., Tuerlinckx, F., Russell, J. A. & Barrett, L. F. The relationship between valence and arousal in subjective experience. Psychol. Bull. 139, 917â940 (2013).
Conner, T. S. & Mehl, M. R. Ambulatory assessment: methods for studying everyday life. In Emerging trends in the social and behavioral sciences: an interdisciplinary, searchable, and linkable resource (eds Scott, R. et al.) 1â15 (Wiley, London, 2015).
Myrtek, M. Progress in ambulatory assessment: computer-assisted psychological and psychophysiological methods in monitoring and field studies (Hogrefe & Huber, Cambridge, 2001).
Myrtek, M., Aschenbrenner, E. & Brügner, G. Emotions in everyday life: an ambulatory monitoring study with female students. Biol. Psychol. 68, 237â255 (2005).
Salgado, S. & Kingo, O. S. How is physiological arousal related to self-reported measures of emotional intensity and valence of events and their autobiographical memories?. Conscious. Cogn. 75, 102811 (2019).
Blascovich, J. & Mendes, W. B. Challenge and threat appraisals. In Feeling and thinking: the role of affect in social contagion, Cambridge (ed. Forgas, J. P.) 59â82 (Cambridge University Press, Cambridge, 2001).
Seery, M. D. Challenge or threat? Cardiovascular indexes of resilience and vulnerability to potential stress in humans. Neurosci. Biobehav. Rev. 35, 1603â1610 (2011).
Tomaka, J., Blascovich, J., Kelsey, R. & Leitten, C. Subjective, physiological, and behavioral effects of threat and challenge appraisal. J. Pers. Soc. Psychol. 65, 248â260 (1993).
Wormwood, J. et al. Physiological indices of challenge and threat: A data-driven investigation of autonomic nervous system reactivity during an active coping stressor task. Psychophysiology 56, e13454 (2019).
Bishop, C. M. Pattern recognition and machine learning (Springer, Berlin, 2006).
Blei, D. M. & Jordan, M. I. Variational inference for Dirichlet process mixtures. Bayesian Anal. 1, 121â143 (2006).
Pedregosa, F. et al. Scikit-learn: machine learning in Python. J. Mach. Learn. Res. 12, 2825â2830 (2011).
Warriner, A. B., Kuperman, V. & Brysbaert, M. Norms of valence, arousal, and dominance for 13,915 English lemmas. Behav. Res. Methods 45, 1191â1207 (2013).
Cacioppo, J. T. et al. Autonomic cardiac control, II: Noninvasive indices and basal response as revealed by autonomic blockades. Psychophysiology 31, 586â598 (1994).
Lakens, D. Calculating and reporting effect sizes to facilitate cumulative science: A practical primer for t-tests and ANOVAs. Front. Psychol. 4, 863â863 (2013).
Sherwood, A. et al. Methodological guidelines for impedance cardiography. Psychophysiology 27, 1â23 (1990).
Cohen, J. The effect size index d: statistical power analysis for the behavioral sciences (Routledge Academic, Abingdon-on-Thames, 1988).
Berntson, G. G. et al. Heart rate variability: origins, methods and interpretive caveats. Psychophysiology 34, 623â647 (1997).
Berntson, G. G., Cacioppo, J. T. & Quigley, K. S. Autonomic cardiac control, I: Estimation and validation from pharmacological blockades. Psychophysiology 31, 572â585 (1994).
Berntson, G., Cacioppo, J. & Quigley, K. Autonomic determinism: The modes of autonomic control, the doctrine of autonomic space, and the laws of autonomic constraint. Psychol. Rev. 98, 459â487 (1991).
White, D. W. & Raven, P. B. Autonomic neural control of heart rate during dynamic exercise: revisited. J. Physiol. 592, 2491â2500 (2014).
Yamamoto, Y., Hughson, R. L. & Peterson, J. C. Autonomic control of heart rate during exercise studied by heart rate variability spectral analysis. J. Appl. Physiol. 71, 1136â1142 (1991).
Le Mau, T. et al. Uncovering emotional structure in 604 expert portrayals of experience. (under review).
Kleckner, I. R. & Quigley, K. S. Mapping the neurophysiological state of the body to affective experience. In The psychological construction of emotion (eds Barrett, L. F. & Russell, J. A.) 265â301 (Guilford Press, New York, 2015).
Sterling, P. Allostasis: a model of predictive regulation. Physiol. Behav. 106, 5â15 (2012).
Sterling, P. & Laughlin, S. Principles of neural design (MIT Press, London, 2015).
Barrett, L. F. & Simmons, W. K. Interoceptive predictions in the brain. Nat. Rev. Neurosci. 16, 419â429 (2015).
Hutchinson, J. B. & Barrett, L. F. The power of predictions: an emerging paradigm for psychological research. Curr. Directions Psychol. Sci. 28, 280â291 (2019).
Pezzulo, G., Rigoli, F. & Friston, K. Active Inference, homeostatic regulation and adaptive behavioural control. Prog. Neurobiol. 134, 17â35 (2015).
Seth, A. K. & Friston, K. J. Active interoceptive inference and the emotional brain. Philos. Trans. R. Soc. B Biol. Sci. 371, 20160007 (2016).
Wolpert, D. M. & Flanagan, J. R. Motor prediction. Curr. Biol. 11, R729âR732 (2001).
Shadmehr, R., Smith, M. A. & Krakauer, J. W. Error correction, sensory prediction, and adaptation in motor control. Annu. Rev. Neurosci. 33, 89â108 (2010).
Obrist, P. A. Cardiovascular psychophysiology: a perspective (Plenum, New York, 1981).
Obrist, P. A., Webb, R. A., Sutterer, J. R. & Howard, J. L. The cardiac-somatic relationship: some reformulations. Psychophysiology 6, 569â587 (1970).
Turner, J. R. & Carroll, D. Heart rate and oxygen consumption during mental arithmetic, a video game, and graded exercise: further evidence of metabolically-exaggerated cardiac adjustments?. Psychophysiology 22, 261â267 (1985).
Coan, J. A., Schaefer, H. S. & Davidson, R. J. Lending a hand: Social regulation of the neural response to threat. Psychol. Sci. 17, 1032â1039 (2006).
Framorando, D. & Gendolla, G. Itâs about effort: Impact of implicit affect on cardiovascular response is context dependent. Psychophysiology 56, e13436 (2019).
Harmon-Jones, E. & Peterson, C. K. Supine body position reduces neural response to anger evocation. Psychol. Sci. 20, 1209â1210 (2009).
MacCormack, J. & Lindquist, K. Feeling hangry? When hunger Is conceptualized as emotion. Emotion 19, 301â319 (2019).
Price, T. F. & Harmon-Jones, E. Approach motivational body postures lean toward left frontal brain activity. Psychophysiology 48, 718â722 (2011).
Mehl, M. R., Gosling, S. D. & Pennebaker, J. W. Personality in its natural habitat: Manifestations and implicit folk theories of personality in daily life. J. Pers. Soc. Psychol. 90, 862â877 (2006).
Kragel, P. A. & LaBar, K. S. Multivariate pattern classification reveals autonomic and experiential representations of discrete emotions. Emotion 13, 681â690. https://doi.org/10.1037/a0031820 (2013).
Stephens, C. L., Christie, I. C. & Friedman, B. H. Autonomic specificity of basic emotions: evidence from pattern classification and cluster analysis. Biol. Psychol. 84, 463â473 (2010).
83Picard, R. W., Vyzas, E. & Healey, J. Toward machine emotional intelligence: Analysis of affective physiological state. IEEE Transactions on Pattern Analysis & Machine Intelligence, 1175â1191 (2001).
Kragel, P. A. & LaBar, K. S. Multivariate neural biomarkers of emotional states are categorically distinct. Soc. Cognit. Affect. Neurosci. 10, 1437â1448 (2015).
Saarimäki, H. et al. Discrete neural signatures of basic emotions. Cereb. Cortex 26, 2563â2573 (2016).
Kragel, P. A., Knodt, A. R., Hariri, A. R. & LaBar, K. S. Decoding spontaneous emotional states in the human brain. PLoS Biol. 14, e2000106 (2016).
Khalaf, A. et al. Analysis of multimodal physiological signals within and between individuals to predict psychological challenge vs. threat. Expert Syst. App. 140, 112890 (2020).
Lindquist, K. A., Siegel, E. H., Quigley, K. S. & Barrett, L. F. The hundred-year emotion war: Are emotions natural kinds or psychological constructions? Comment on Lench, Flores, and Bench (2011). Psychol. Bull. 139, 255â263 (2013).
James, W. The physical basis of emotion. Psychol. Rev. 1, 516â529 (1894).
Qu, M., Zhang, Y., Webster, J. G. & Tompkins, W. J. Motion artifact from spot and band electrodes during impedance cardiography. IEEE Transactions on Biomedical Engineering, 1029â1036 (1986).
Tugade, M. M., Fredrickson, B. L. & Barrett, L. F. Psychological resilience and positive emotional granularity: examining the benefits of positive emotions on coping and health. J. Pers. 72, 1161â1190 (2004).
Widdershoven, R. L. A. et al. Effect of self-monitoring through experience sampling on emotion differentiation in depression. J. Affect. Disord. 244, 71â77 (2019).
Stone, A. A. et al. A population approach to the study of emotion: Diurnal rhythms of a working day examined with the Day Reconstruction Method. Emotion 6, 139â149 (2006).
Nabian, M. et al. An open-source feature extraction tool for the analysis of peripheral physiological data. IEEE J. Transl. Eng. Health Med. 6, 1â11 (2018).
95Hamilton, P. Open source ECG analysis. in Computers in Cardiology. 101â104 (IEEE).
96Carreiras, C. et al. BioSPPy: Biosignal Processing in Python. (2018).
Berntson, G. G., Quigley, K. S., Jang, J. F. & Boysen, S. T. An approach to artifact identification: Application to heart period data. Psychophysiology 27, 586â598 (1990).
Forouzanfar, M. et al. Toward a better noninvasive assessment of preejection period: a novel automatic algorithm for B-point detection and correction on thoracic impedance cardiogram. Psychophysiology 55, e13072 (2018).
Berntson, G. G., Lozano, D. L., Chen, Y. J. & Cacioppo, J. T. Where to Q in PEP. Psychophysiology 41, 333â337 (2004).
Kubicek, W., Karnegis, J., Patterson, R., Witsoe, D. & Mattson, R. Development and evaluation of an impedance cardiac output system. Aerosp. Med. 37, 1208â1212 (1966).
Rasmussen, C. E. The infinite Gaussian mixture model. in Advances in Neural Information Processing Systems. 554â560.
Author information
Authors and Affiliations
Contributions
J.B.W., K.S.Q., L.F.B., and J.D. designed the study. M.J.F., M.D., C.N., and K.H. collected the data. All authors developed the data pre-processing and analysis plan; Z.K., K.H., and M.J.F. analyzed the data. K.H., M.J.F., Z.K., and C.N. wrote the manuscript. All authors reviewed and revised the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hoemann, K., Khan, Z., Feldman, M.J. et al. Context-aware experience sampling reveals the scale of variation in affective experience. Sci Rep 10, 12459 (2020). https://doi.org/10.1038/s41598-020-69180-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-69180-y