Abstract

Socioemotional processes engaged in daily life may afford and/or constrain individuals’ emotion regulation in ways that affect psychological health. Recent findings from experience sampling studies suggest that persistence of negative emotions (emotion inertia), the strength of relations among an individual’s negative emotions (density of the emotion network), and cycles of negative/aggressive interpersonal transactions are related to psychological health. Using multiple bursts of intensive experience sampling data obtained from 150 persons over one year, person-specific analysis, and impulse response analysis, this study quantifies the complex and interconnected socioemotional processes that surround individuals’ daily social interactions and on-going regulation of negative emotion in terms of recovery time. We also examine how this measure of regulatory inefficiency is related to interindividual differences and intraindividual change in level of depressive symptoms. Individuals with longer recovery times had higher overall level of depressive symptoms. Also, during periods where recovery time of sadness was longer than usual, individuals’ depressive symptoms were also higher than usual, particularly among individuals who experienced higher overall level of stressful life events. The findings and analysis highlight the utility of a person-specific network approach to study emotion regulation, how regulatory processes change over time, and potentially how planned changes in the configuration of individuals’ systems may contribute to psychological health.

1. Introduction

Lifespan developmental theories view persons as complex dynamic systems, with feelings, thoughts, and actions that are interconnected and that change over time. Individual development is the product of numerous dynamic processes that span multiple levels of analysis, multiple domains of functioning, and multiple time-scales [13]. At their core, most developmental theories conceive individuals’ development as the output of a complex dynamic system that reorganizes and changes as individuals transition between life phases and are affected by life events.

In line with basic tenets of developmental theory, conceptual work and empirical studies have elaborated the view that individuals’ on-going socioemotional processes are a dynamic system wherein emotions and social behaviors interact to produce or influence psychological health and development [46]. Emotions facilitate social behaviors [7], and social interactions regulate emotions [8], and these system dynamics are associated with general psychological health and with depression [911]. More specifically, the cyclic cause-effect structures nested within interconnected networks of emotions and social behaviors may afford and/or constrain individuals’ emotion regulation in ways that affect psychological health. This paper presents new empirical work examining how the dynamics of individuals’ daily emotional and social experiences are related to both interindividual differences and intraindividual changes in depressive symptoms. By leveraging intensive experience sampling data and new time-series based network methods, we examine how inefficiency of socioemotional regulation processes is linked to individuals’ experience of depressive symptoms.

1.1. Emotion Regulation, Depression, and Feedback Loops

Substantial evidence suggests that inefficient regulation of negative emotions puts individuals at risk for and is a marker of depression [12]. Consistent with identification of the persistence of negative mood as a core symptom of depression [13], multiple experience sampling studies illustrate how greater emotional inertia, the extent of carryover from moment to moment, of negative emotion (i.e., within-person autocorrelation) is associated with depression and psychological maladjustment [14, 15]. Expanding on these studies, researchers have begun examining the temporal dynamics of multiple emotions simultaneously. Taking advantage of the conceptual and methodological advantages of time-series based network methods [16], Pe and colleagues [17], for example, examined how the structure of relations among 11 negative and positive emotions was related to psychological health. They found that individuals with major depressive disorder had stronger temporal relations among emotions, greater density of the emotion network, especially greater density of the network of negative emotions, compared to controls (see also [18]). Interpretation focuses on how network density indicates the emotion system’s resistance to change and presence of spirals of mutually reinforcing negative emotions. The underlying idea is that dense networks are likely to contain feedback loops in which any disruption in negative emotions will reverberate and persist, a form of inefficient emotion regulation [19]. In contrast, sparse networks are unlikely to contain feedback loops, so that any newly introduced negative emotion dissipates before influencing other parts of the network, a form of efficient emotion regulation.

Building on this work, we further expand the relevant network of variables to include aspects of both emotional and social experience. Theoretical models and empirical findings already highlight the variety of interpersonal mechanisms that may be contributing to maladaptive emotion regulation and risk for depression [20, 21]. For example, excessive seeking of reassurance [22], negative evaluation from friends [23], and excessive self-disclosure regardless of social context [6] have all been linked to depression. Such behaviors often solicit negative reactions from social partners, including domineering behavior, dismissiveness, or rejection [22, 2427] that, when an individual has strong emotional reactivity to close relationships, can cycle back to produce more negative emotions [28]. Indeed, experience sampling studies show that social interactions characterized by less warm and more submissive behaviors are accompanied by more negative emotions [29, 30]. Particularly problematic configurations of emotional and social experiences involve feedback loops that perpetuate experience of negative emotions, for example, when negative emotions lead to maladaptive social behaviors and interpersonal interactions lead to negative emotions (e.g., a cycle of negative/aggressive interpersonal transactions [31, 32]). In sum, studies of daily life suggest links between interpersonal and emotional experiences, with the possibility that specific types of dynamics, namely, feedback loops, are associated with individuals’ experience of depressive symptoms.

To illustrate more directly how feedback loops may contribute to regulation, network-based models of two hypothetical individuals’ socioemotional dynamics are shown in Figure 1. In Individual A’s network (Figure 1(a)), higher sadness leads to lower happiness (a temporal relation where sadness influences subsequent happiness at -0.6), lower happiness leads to lower social engagement (a temporal relation where happiness influences subsequent communal behavior at +0.6), and lower social engagement leads to higher sadness (a temporal relation where communal behavior influences subsequent sadness at -0.6). The overall effect of this cyclic structure is positive, as indicated by multiplication of the three temporal relations (-0.6 0.6 (-0.6) = 0.216). Hence the three temporal relations together form a positive feedback loop, a structure that sustains changes in sadness. In contrast, in Individual B’s network (Figure 1(b)), there is a negative feedback loop because social engagement is positively associated with sadness. The overall effect of this cyclic structure is negative (-0.6 0.6 0.6 = -0.216); hence the three temporal relations together form a negative feedback loop. Here, the structure of the relations facilitates regulation of sadness. The impact of the feedback loops can be examined mathematically using impulse response analysis ([33]; details introduced in the Method section). As shown in the accompanying temporal profile of sadness, Individual A recovers from an increase in sadness (sadness = 1.0 at t = 1) by t = 18 (Figure 1(c)). In contrast, Individual B recovers from an increase in sadness (sadness = 1.0 at t = 1) by t = 15 (Figure 1(d)). This comparison illustrates that the positive feedback loop extends the recovery time of sadness by approximately 20%, an inefficient emotion regulation process.

1.2. Experience Sampling, Person-Specific Networks, and Recovery Time

Identifying feedback loops is relatively straightforward in low-dimensional systems or through experimental manipulation of controllable systems [34]. The potential complexity, however, compounds as the systems become larger (more variables) and more complicated (more realistic). Study of high-dimensional, multivariate human systems can be facilitated by (a) collection of intensive experience sampling (time-series) data and a combination of methods that (b) identify the network structure from those empirical time-series, (c) characterize performance of the network (e.g., efficiency of emotion regulation) by impulse response analysis, and (d) describe interindividual and intraindividual differences in systems through regression-based modeling.

1.2.1. Intensive Experience Sampling

Technological advances in mobile computing provide an infrastructure that allows for unprecedented opportunity to obtain the temporally dense and comprehensive experience sampling needed for studying individuals as high-dimensional, multivariate dynamic systems [35]. Studies wherein individuals provide many reports on their emotions and interpersonal behaviors as they go about their daily lives are beginning to obtain the types of multivariate time-series data needed to identify and model the complex feedback loops involved in emotion regulation. Multiple time-scale or “measurement burst” study designs, wherein data are collected at both micro- and macrotimescales (hours and months; [36, 37]) provide new opportunities to observe how the moment-to-moment processes governing short-term behavior (e.g., socioemotional dynamics) evolve over the long-term.

1.2.2. Modeling Socioemotional Process as a Person-Specific Network

The time series data collected in experience sampling studies facilitates use of network methods for studying within-person processes [16, 19]. In these models, the temporal relations among variables in the time series data are estimated by person-specific multivariate time-series methods and depicted as a network of paths between nodes [38]. For example, the intensive experience sampling data shown in Figure 2(a) was used to derive the network of relations shown in Figure 2(b). Each emotion and social behavior variable is represented as a node (labeled circles) in the network, and the statistical relations between each pair of variables is represented as an edge (arrows). Influences of variables on other variables are represented as directional paths that indicate how changes in one variable influence other variables subsequently. The directionality of edge indicates causal relations (e.g., the edge pointing from happiness to anger indicates that changes in happiness are likely to lead to changes in anger). The sign, strength, and temporal lag of the relations are indicated by color (green = positive, red = negative), line width (wider = stronger), and line-type (dashed = lag-1, solid = contemporaneous), respectively. Altogether the 13 nodes and connecting edges in the network shown in Figure 2(b) provide a model for how this specific individual’s socioemotional system functioned during the 21-day period during which they provided the data.

1.2.3. Recovery Time as a Description of Individuals’ Emotion Regulation

The person-specific network depicted in Figure 2(b) is notably larger and more complex than the networks depicted in Figure 1. While the network does provide better coverage of the socioemotional space (13 versus 3 variables), identification and interpretation of the embedded feedback loops are substantially more difficult. Often, the structure of larger networks is quantified using summary measures such as network-density and node-centrality (see [39]). However, because these metrics usually involve summing the absolute value of edges, they do not differentiate between positive (excitatory) or negative (inhibitory) feedback loops, i.e., loops that have opposite regulatory function (note that the networks in Figure 1 both have density = 2.4). More direct quantification of the emotion regulatory implications of the network structure can be obtained using impulse response analysis [33]. An “impulse” is given to a specific node and the behavior of the system observed through simulation over many time steps [4042]. For example, in Figure 2(c), we see the behavior of the network shown in Figure 2(b) after an impulse is delivered to the sadness node. The impulse filters through connected nodes (due to the temporal relations) before returning to equilibrium. The recovery time depends on the existence and configuration of positive and/or negative loops within the network (see, e.g., [4345]), with positive feedback loops extending recovery time and negative feedback loops shortening recovery time. This simulation approach provides new opportunities to measure the efficiency of individuals’ emotion regulation.

1.3. The Present Study: Associations between Recovery Time and Depressive Symptoms

In this study we use multiple “bursts” of intensive experience sampling data obtained from 150 persons over one year, person-specific network analysis, and impulse response analysis methods to derive a new network metric, recovery time, to quantify emotion regulation efficiency and emphasize consideration of the complex and interconnected socioemotional processes surrounding day-to-day social interactions and regulation of emotions. We then examine how recovery time for sadness is related to interindividual differences and intraindividual change in the experience of depressive symptoms. We hypothesize that the network’s efficiency to regulate sadness is related to level of depressive symptoms because persistent sadness is a prominent feature of depression [13]. Specifically, we expect that individuals with longer recovery time will experience more depressive symptoms (between-person association). During bursts where recovery times are longer than usual, the prototypical person will also experience more depressive symptoms than usual (within-person association).

Acknowledging that emotion regulation is related to and influenced by interpersonal behaviors (see [6] for review), we include as many variables as possible in the person-specific networks. Available data obtained in the context of the normal social interactions in daily lives include some negative emotions (sadness, anger, and ashamed), some positive emotions (happy, proud), interpersonal behaviors and perceptions (communion, agency), perceived benefits for self and other, perceived control, and self-esteem. Rather than specifying how each variable contributes to regulation of sadness, we explicitly take a holistic view and use the temporal relations among all these variables to derive how an individual’s regulation of sadness is related to the experience of depressive symptoms. Given that change in depressive symptoms may also be related to individuals’ immediate life context [4648], we control for differences in stressful life events.

2. Method

Our analysis makes use of data from the Intraindividual Study of Affect, Health, and Interpersonal Behavior (iSAHIB), a multiple time-scale experience sampling study designed for articulation and study of process-oriented theory and methods [49]. An analysis tutorial is available in the supplementary material (available here), as well as at https://quantdev.ssri.psu.edu/tutorials.

2.1. Participants

The iSAHIB sample consists of 150 adults (50% women), recruited from The Pennsylvania State University and surrounding community and stratified by gender and age to cover the full adult life span. Participants ranged in age from 18 to 89 years (MAge = 47.10, SDAge = 18.76) and had obtained between 2 and 24 years of formal education (MEduc = 16.36, SDEduc = 3.90), with 91% self-identifying as Caucasian (4% African American, 1% Asian American, and 4% Mixed or Other ethnicity). Most individuals identified as heterosexual (93%) with 6% identifying as bisexual/gay/lesbian. After participants were recruited, informed of the intensive nature of the assessments, and self-selected into the study, they began the assessment protocol. Over the course of a year, they provided extensive reports about their lives through a combination of web-based (completed during visits to the laboratory) and smartphone-based (completed multiple times per day during regular daily life) questionnaires.

Although the participants are drawn from the general community, they exhibit a range of maladjustment. Borderline personality disorder symptoms were measured using the Personality Assessment Inventory-Borderline Features Scale [50], a 24-item Likert-scale questionnaire, with items measured on a 1 to 4 Likert scale. Participants were prompted to “give your own opinion of yourself” on 4 dimensions, including affective instability, identity problems, negative relations, and self-harm. Item responses were summed to obtain a composite measure of disorder severity, which ranged from 27 to 72 (MPAIBFS = 44.6, SDPAIBFS = 10.0). Child abuse and trauma experiences were captured with the Child Abuse and Trauma Scale [51], a 14-item questionnaire, with items measured on a 1 to 5 Likert scale. Participants were promoted to “respond to the question in terms of the person or persons who had the primary responsibility for your upbringing as a child” on 4 dimensions, including physical, verbal, sexual, and emotional abuse. Item responses were summed to obtain a composite measure of the severity of child abuse and trauma, which ranged from 14 to 59 (MCATS = 23.3, SDCATS = 9.3). Interpersonal problems were measured using the Interpersonal Problems Circumplex [52], a 32-item questionnaire, with items measured on a 0 to 4 Likert scale. Participants were prompted to rate themselves on 8 dimensions, including domineering, vindictive, cold, socially avoidant, nonassertive, exploitable, overly nurturing, and intrusive. Item responses were summed to obtain a composite measure of the severity of interpersonal problems, which ranged from 6 to 73 (MIPCS = 33.6, SDIPCS = 15.7).

2.2. Experience Sampling Procedure

Participants completed three 21-day “measurement-bursts” spaced at approximately even intervals over one year. During each 21-day burst, individuals used a study-provided smartphone with a customized application to provide event-contingent reports of face-to-face social interactions that lasted longer than five minutes. Each interaction-level report included descriptions of when and where the interaction occurred, whom the interaction was with, how they and their partners behaved, and how they felt afterward. To facilitate compliance, the smartphones were programmed to chime a prompt if the participant did not provide a report for any two-hour span between 8 am and 8 pm. Data flow was monitored in real-time, a process that enabled the research staff to make periodic “check-in” calls that supported, motivated, and helped participants (e.g., solving technical problems) provide high-quality data. Over the entire study period, participants provided multivariate time-series with, on average, 427.4 observations (SDT = 145.7, Range = 88 to 869) during 422 measurement burst periods (of a possible 450, because of some sample attrition). At the beginning and end of each burst, individuals visited the laboratory, received training or debriefing, picked-up or dropped-off smartphones, and completed demographic, health, personality, and other questionnaires. Participants were compensated $500 for completing the entire protocol.

2.3. Measures

The present analysis uses all available emotion and interpersonal variables in the experience sampling protocol to derive network representations of individuals as high-dimensional, multivariate dynamic systems and relates specific aspects of network function (recovery time following an impulse of sadness) to burst-level data on individuals’ depressive symptoms and experience of life events.

2.3.1. Emotion and Interpersonal Behavior

After every social interaction (event-contingent sampling), participants were prompted with a series of items that they responded to using a ‘touch-point continuum’ (e.g., slider-type interface that allowed for continuous range between left anchor “Not at All” and right anchor “Very Much”) that was digitally coded on a 0 to 100 scale (numbers not visible to participants). After each social interaction, individuals reported on five emotions. Individuals’ shame was measured using the item, “How ASHAMED do you feel right now?” (“Not at all … Very much”). Across all social interactions, individuals rated very low levels of shame (M = 6.45, SD = 9.42). Individuals’ anger was measured using one item, “How ANGRY do you feel right now?” (M = 8.19, SD = 12.85). Individuals’ sadness was measured using the item, “How SAD do you feel right now?” (M = 10.63, SD = 15.51). Individuals’ happiness was measured using the item, “How HAPPY do you feel right now?” (M = 63.82, SD = 22.00). Individuals’ pride was measured using the item, “How PROUD do you feel right now?” (M = 54.14, SD = 27.23).

Individuals also reported on a variety of interpersonal and intrapersonal factors. Social partners’ interpersonal communion and agency [53] was measured by asking participants to “Rate how THE OTHER PERSON acted”, and two sliders with end-point anchors of “Distant … Friendly” and “Submissive … Dominant” for communion-other (M = 80.44, SD = 15.96) and agency-other (M = 56.79, SD = 18.20), respectively. Parallel measurement of the participant’s own interpersonal behavior, communion-self (M = 82.09, SD = 14.78) and agency-self (M = 55.08, SD = 17.23), was prompted by the item “Rate how YOU acted”, and the same two sliders. Individuals’ current impression of benefit for self (benefit-self) was measured by the item, “How useful was this interaction for YOU?” (“Very costly Very Beneficial”; M = 64.11, SD = 20.46). In parallel, individuals’ current impression of benefit for the person they just interacted with (benefit-other) was measured by the item, “How useful was this interaction for THE OTHER PERSON?” (M = 65.51, SD = 19.13). Individuals’ perceived control was measured using the item, “I have control over the things happening to me right now." (“Not at all... Very much”; M = 70.81, SD = 21.67), and self-esteem with the item, “I have high self-esteem right now." (M = 67.09, SD = 23.77).

For illustration, one individual’s multivariate, 13-dimensional time-series data from one burst are shown in Figure 2(a). As can be seen, the value of each variable fluctuated from interaction to interaction across the course of study, some more than others, with some moving in synchrony (e.g., anger and shame, cross-correlation = 0.34), some moving in opposite directions (e.g., anger and happiness, cross-correlation = -0.22) and some with minimal interrelations (e.g., agency-other and agency-self, cross-correlation = -0.001). Notable are the length of this time-series (T > 150), the extent of intraindividual variability, the “stationarity” of the series (i.e., fluctuating rather than drifting up or down over time), and the level of synchrony among pairs of variables.

2.3.2. Depressive Symptoms and Life Events

Prior to each of the three bursts, individuals’ recent experience of depressive symptoms was measured using the 20-item Center of Epidemiologic Studies Depression Scale (CESD; [48]). Participants were prompted with the stem, “How often have you felt this way during the past week?” followed by a list of symptoms (e.g., loss of appetite, restless sleep, feeling lonely, being happy). Each of these items required indication in 1 of 4 checkboxes labeled “Rarely or none of the time (less than 1 day), Some or little of the time (1-2 days), Occasionally or moderate amount of time (3-4 days), Most or all of the time (5-7 days)”. Item responses on a 0 to 3 scale were summed to obtain a composite measure of the severity of depressive symptoms for each burst (M = 10.00, SD = 8.25). A CESD score of 16 is a recommended cut-off to screen for clinical depression [54]. In the current sample, the number of participants scoring ≥16 was 67 participants at burst 1, 57 participants at burst 2, and 47 participants at burst 3. In the context of a general community (versus a clinical) sample, we chose to operationalize depression dimensionally (level of depressive symptoms) rather than categorically (depressive disorder present/absent). Clinical psychology and psychiatry are moving toward dimensional models of psychopathology due to the limitations of categorical models of mental disorder [55]. It is common to employ dimensional measures of depressive symptoms in clinical research, including studies using the CESD (e.g., [56]) and studies examining associations among emotion networks and depressive symptoms (e.g., [18]). Dimensional models of depression are better predictors of functioning than categorical ones [57] and are part of a broader dimension of internalizing symptoms [58].

Given that burst-to-burst changes in depressive symptoms may be related to individuals’ immediate life context [4648], we controlled for differences in recent life events. Individuals’ recent life experiences were probed using 12 items adapted from life events scales [59, 60]. Participants were prompted with the stem, “Since the last time we saw you, [or at the first visit, “In the last 6 months,”] did you experience a change in the following and, if so, how much did it affect you?” followed by a list that included significant life events (e.g., change in relationship status, loss of loved one, hospitalization). Each of these items required participants to check 1 of 5 checkboxes labeled “Did not experience, Not affected, A little bit, Somewhat, A lot”. Responses coded on a 0 to 4 scale were summed to obtain a composite measure of the impact of recent life events for everyone at each burst (M = 6.18, SD = 6.08).

2.4. Data Preparation and Analysis

There were three stages in the data analysis. In the first stage, uSEM [61] was used to construct person-specific networks that describe the configuration and temporal relations underlying each person’s 13-dimensional multivariate time-series data at each burst (422 networks). In the second stage, the behavior of these networks was quantified using impulse response analysis [33] to obtain network-specific recovery times for sadness, a measure of emotion regulation. In the third stage, we used multilevel models to examine how recovery time was related to interindividual differences and intraindividual changes in depressive symptoms, controlling for differences and changes in life events.

2.4.1. Data Preparation

Before analysis, the data were examined for suitability of application. The general guiding principle for application of uSEM is that the multivariate time-series data can be treated as weakly stationary (with means and variance-covariance structure that is constant over time; see details on preprocessing in [7]). Visual inspection of each individual’s data and testing of polynomial trends (see [62]) suggested that the data were reasonably stationary (an exemplar participant’s time-series is shown in Figure 2(a)). To focus analysis on intraindividual regulation dynamics, the 13-variable time-series for each burst for each person were standardized into a z-metric (M = 0, SD = 1), thereby effectively removing burst-to-burst and person-level differences in level and variance (see, e.g., [63]).

2.4.2. Construction of Networks

Individual data from each burst were then modeled as a multinode dynamic network using a unified Structural Equation Model (uSEM, [61]). In brief, the multivariate observed time-series is modeled as the output of a latent variable time series ,where is a factor loading matrix and is a time-series of residuals with variance-covariance structure given by a matrix , that is assumed diagonal. The temporal relations among the set of latent constructs in (the circles in Figure 2(b)) are then modeled aswhere is a vector of the lag-1 version of the multivariate latent time-series; is a matrix of regression parameters that describe the contemporaneous relations among the latent variables (solid arrows in Figure 2(b)), is a matrix of regression parameters that describe the lag-1 relations (auto- and cross-regressions) among the latent variables (dashed arrows in Figure 2(b)), and is a multivariate “shock” or input time series. Together, the contemporaneous relations in and auto- and cross-regressive relations in indicate the causal influences among variables through which exogenous input is processed and diffused (i.e., dynamic regulation). At the practical level, the uSEM model is estimated using an iterative search process wherein a series of models are constructed and tested for improvements in fit. At each step, Lagrange Multiplier tests (modification indices; Sörbom, 1989) are used to select the path that facilitates maximum improvement in fit. This element is then freed, the model reestimated, and a new set of modification indices calculated, iteratively adding paths until further addition does not significantly improve model fit. The model expansion was constrained so that only and blocks of the model parameter matrix were freed, thus keeping the time-series structure of the model intact. Bidirectional paths in the contemporaneous relations are avoided by including all potential autoregression relations in the initial model and by deeming the opposite path unavailable when any given element in was freed. In experience sampling study designs being used to collect multivariate, intensive longitudinal data, the item pool has often been optimized to minimize participant burden. In cases where each latent construct has been measured by only one item, the factor loading matrix is configured as an identity matrix , and all elements of , the variance-covariance matrix of , and are fixed = 0.

Person-specific models for each burst were estimated by adapting R code from the Group Iterative Multiple Model Estimation package (GIMME; [64]). The adaptation was mainly to ensure the model fitting procedure will produce an interpretable result, including allowing no more than one direction of contemporaneous relation being fitted between two variables, and setting the autoregression to be freed in the initial iteration in the fitting procedure. Once the person-specific models for each burst were obtained, and matrices were extracted and drawn as network graphs using the qgraph package [65]. Conceptually, the resulting network describes how behavior is organized and proceeds at the micro time-scale. A sample network is shown in Figure 2(b), where, for example, sadness was predictive of lower other’s communal behavior at the next observation, and other’s communal behavior was also predicting of lower sadness at the next observation (red dashed line from other’s communion to sadness and the opposite direction in Figure 2(b)). These two edges together form a positive feedback loop between sadness and other’s communion.

2.4.3. Impulse Response Analysis and Recovery Time

Each of the 422 networks (150 persons x 3 bursts, minus some attrition) was then summarized with respect to the recovery time of sadness, quantified as the number of time steps until the level of sadness returns to near equilibrium (e.g., within 0.01 of the asymptote) after a hypothetical perturbation. Formally, the impulse response simulation model is constructed by converting the uSEM into a vector autoregression model and doing one step ahead forecasting (see (3); [33, 61, 66]).

In our case, the system is set in motion by sending an initial impulse to the sadness node (sadness = 1.0 at t = 1) and computing how the system evolves over 150 time steps (to guarantee a sufficient length for all nodes to return to equilibrium). The time profile obtained from the impulse response analysis of the network in Figure 2(b) is shown in Figure 2(c). Recovery time, defined as time to return within ±0.01 of equilibrium, was then derived through a backward search to accommodate oscillation in the time profiles. Specifically, we searched backward from the end of the time profile, to identify the time-step, denoted as k, where the level of a specific variable was first outside the ±0.01 boundary. Recovery time was then quantified for the sadness node as k, the number of time steps from perturbation to equilibrium. Within-person differences across bursts are illustrated in Figure 3. The recovery time of sadness for this individual changed across bursts, starting at k1 = 3 in the first burst, increasing to a k2 = 8 in the second burst, before returning to k3 = 2. Because the distribution of recovery times is skewed (cannot go below zero), scores were log-transformed before being used in the third stage of analysis.

2.4.4. Associations between Recovery Time and Depressive Symptoms

Because the equilibrium represents the average value of sadness and the normal sadness level is rather low (M = 10.63, SD = 15.51), disruption of sadness from the equilibrium is most likely to result in an increased level of sadness, which is undesirable subjectively. Thus, positive feedback loops around sadness which sustain this disruption are also undesirable. This informed our hypothesis that emotion regulation inefficiency of sadness embedded in individuals’ socioemotional networks is related to higher (individual differences in) and increased (intraindividual changes in) depressive symptoms (ICC = 0.65, skewness = 1.25), controlling for life events. Making use of and accommodating the nested nature of the multiple-burst longitudinal data (422 bursts nested within 150 persons), hypotheses were examined within a multilevel modeling framework [67]. Following usual practice, the predictor variables were split into time-invariant (person-level means; OverallLifeEventsi, OverallRecoveryTimei) and time-varying (burst-to-burst deviations, BurstLifeEventsib, and BurstRecoveryTimeib) components. Relations among the extended set of variables were then examined using 2-level models of the form where the repeated measures of depressive symptoms for individual i at burst b, , are modeled as a function of person-specific intercepts, , that indicate baseline level of depressive symptoms; person-specific coefficients, and , that indicate the extent of within-person associations between burst-specific log recovery time or life events, respectively, and depressive symptoms; and coefficient, , that capture how life events moderate the within-person association between burst-specific log recovery time and depressive symptoms. Person-specific coefficients were simultaneously modeled as a function of person-level predictorswhere to are sample-level parameters and to are residual unexplained between-person differences that are assumed multivariate normal with variances , , and and covariances , , and . The model was fit to the data using the nlme package in R [68], with incomplete data (0.2%) treated as missing at random. Person-level predictors were sample-centered to facilitate interpretation of model parameters as representing effects for the prototypical person (as described by the average demographics above). All possible interactions were tested, but, to maintain parsimony in the presentation of the final models, iteratively trimmed to remove those that were nonsignificant (α = 0.05) and not directly related to the hypotheses (always retaining the component main effects and lower-order interactions). Also, a variety of random effects structures were tested, with the different configurations having little influence on the fixed effects parameters or interpretations (i.e., no changes in significant effects). For parsimony, we present and interpret only the final model, which included random effects u0i and u2i.

3. Results

3.1. Socioemotional Networks and Recovery Time

Person-specific socioemotional networks were derived from the 13-variate time-series obtained from each person during each 21-day burst of measurement using uSEM. Of the total 422 network models, 411 fit the data well, as indicated by at least three of the following fit criteria: RMSEAs 0.08, SRMRs 0.08, CFIs 0.95, NNFI 0.95 (see Beltz et al., 2013). Models from 11 bursts with relatively short time-series (T = 31 to 62, compared to Taverage = 145.7) did not fit well and were set aside. The lagged and contemporaneous relations among variables were extracted and used in an impulse response analysis to calculate person- and burst-specific recovery times for sadness. Recovery time, interpreted as a measure of (in)efficiency of emotion regulation, ranged from 1 to 57 (Msad.RT = 4.63, SDsad.RT = 5.32).

3.2. Associations between Recovery Time and Depressive Symptoms

Results from the multilevel model examining between-person and within-person associations between log recovery time of sadness and level of depressive symptoms are shown in Table 1. Level of depressive symptoms for a prototypical individual in an average burst was 9.84 ( = 9.84, p < 0.001) on a 0 to 60 scale. As expected, between-person differences in log recovery time of sadness were associated with differences in level of depressive symptoms, with, as shown in Figure 4(a), longer recovery times linked to higher level of depressive symptoms (γ01 = 2.29, p = 0.031), even after controlling for the significant effect of stressful life events (γ02 = 0.70, p < 0.001). The within-person association was also significant. During a burst where log recovery time of sadness was longer than usual, the prototypical person had a higher level of depressive symptoms (γ10 = 1.26, p = 0.013). However, this association was moderated by individuals’ overall exposure to life events (γ12 = 0.32, p = 0.005). As shown in Figure 4(b), for an individual with high (+1 SD; blue line) exposure to life events, within-person changes in recovery time for sadness were strongly linked to depressive symptoms (implied within-person association β1i = 2.80, 95% CI = [0.22, 5.01]), while for an individual with low (-1SD; red line) exposure to life events, within-person changes in recovery time for sadness were not linked to depressive symptoms (implied within-person association β1i = -0.28, 95% CI = [-2.86, 2.05]). Post hoc probing using the Johnson-Neyman method ([69]; implemented using the probemod package in R; [70]) indicated a significant within-person link between recovery time and depressive symptoms when level of (sample-centered) overall life-events was greater than 2.49.

Post hoc analysis controlling for age, gender, and education (in years) found the same pattern of results. The within-person association between log recovery time and depressive symptoms was significant (γ10 = 1.35, p = 0.009), while the between-person association became marginally significant (γ01 = 1.96, p = 0.065). Additionally, a separate post hoc analysis controlled for scores on the Borderline Features, Child Abuse and Trauma, and Interpersonal Problem Circumplex scales. The within-person association between log recovery time and depressive symptoms was significant (γ10 = 1.16, p = 0.02), while the between-person association became nonsignificant, but the direction of the association held (γ01 = 0.93, p = 0.30).

The autoregression of sadness of course plays a substantial role in recovery time of sadness. Thus, to check whether the results were only driven by a single variable rather than the feedback loops embedded in the larger network, we reran the impulse response analysis with the autoregression of sadness set to zero and calculated the recovery times again. Between-person differences in this log recovery time were no longer associated with depressive symptoms (γ01 = 1.14, p = 0.31), but the within-person association remained robust. In bursts where this log recovery time was longer, depressive symptoms were higher (γ10 = 1.22, p = 0.01), highlighting the effect of temporal relations aside from the autoregression (emotional inertia) of sadness.

4. Discussion

This paper examined emotion regulation from a within-person, process-oriented, and network perspective. Data obtained across multiple “bursts” of intensive experience sampling data over one year were used to construct person-specific networks that described the complex and interconnected socioemotional processes that surround individuals’ day-to-day social interactions and on-going regulation of negative emotion. Impulse response analysis was used to describe and quantify the efficiency of individuals’ regulation in terms of recovery time, and multilevel models were used to examine how recovery time was associated with between-person differences and within-person change in individuals’ experience of depressive symptoms.

Working from basic principles, individuals were viewed as holistic interactive dynamic systems with a broad range of emotions and interpersonal behaviors that influence how they regulate negative emotions. We found that the behavior of the person-specific networks, in particular the recovery time of the illustrative negative emotion (sadness), was related to both between-person differences and within-person changes in depressive symptoms. In line with hypotheses, individuals with longer recovery times had higher overall level of depressive symptoms, even after controlling for recent life events. Also, during periods where recovery time of sadness was longer than usual, depressive symptoms were also higher than usual, particularly among individuals who experienced a higher overall level of life events. These results indicated both between-person and within-person links between regulatory inefficiency and depressive symptoms.

4.1. Person-Specific Networks, Emotion Regulation, and Depressive Symptoms

In line with prior work [17, 18, 71], we found evidence that emotional experiences in daily life are temporally related to each other (e.g., sadness, happiness) and to interpersonal behaviors (communion, agency). Distinct from prior studies using a multilevel modeling approach, the person-specific analysis approach allowed idiosyncratic temporal relations between any of the two nodes without constraining it to a sample-level mean. In the same way that multilevel models allow for interindividual differences in the relations among 2 or 3 variables, the person-specific network approach allows for heterogeneity in the structure of relations among many variables. Indeed, of the 411 networks we obtained, none were identical. Each individual and each burst provided a unique configuration of temporal relations, and this provided for examination of both interindividual differences and intraindividual changes in the socioemotional (network) dynamics.

Viewed as holistic representations of individual functioning, all 411 socioemotional network configurations were evaluated with respect to how well that configuration facilitated regulation of low arousal negative emotion, specifically, recovery time of sadness. Generally, regulatory processes, including those involved in emotion regulation, are engaged to bring a system back to equilibrium [15]. In prior work, efficiency to return to equilibrium has been examined through analysis of univariate time-series data. Specifically, the extent of autocorrelation or autoregression, emotion inertia, is quantified using experience sampling of single emotions and interpreted as an indicator of emotion dysregulation [14, 15]. The present study, along with other recent work [17], expands the investigation of emotion inertia and regulation by explicitly acknowledging that other aspects of daily life (e.g., interpersonal relations, control, self-esteem) can afford or constrain emotion regulation.

The between-person findings that individuals with longer recovery times for sadness had higher overall level of depressive symptoms confirms the interpretation of the recovery time metric as a measure of regulatory (in)efficiency and aligns with prior work demonstrating the relation between depression (or other psychological maladjustment) and inefficient regulation of negative emotions, and emotion inertia [12, 14, 15]. The interpretation is further bolstered by the within-person findings. Indeed, during periods where recovery time of sadness was longer than usual, depressive symptoms were also higher than usual. This is a new finding, facilitated by a multiple burst study design that provides for a more direct test of the within-person links between emotion (dys)regulation and depressive symptoms that has not been available in prior (single-burst or cross-sectional study) studies. Interestingly, after controlling for between-person differences in demographic variables as well as broad maladjustment scales, the within-person association between recovery time of sadness and depressive symptoms remains robust, indicating the within-person association deserves more consideration for future examination.

Importantly, this within-person association was moderated by overall level of life events. The fact that the link between regulatory inefficiency (as indicated by longer recovery time) and depressive symptoms was especially strong among individuals who experienced multiple, impactful life events suggests that the within-person links may be easier to observe when individuals are under duress, for example, in a “testing-the-limits” context (see [72]). To illustrate, consider the diagnosis of heart disease. Individuals are typically subjected to an exercise “stress test” wherein their cardiovascular reactivity and regulation is observed as their bodies are pushed towards their physiological limits (e.g., running on a treadmill). The general idea is to produce a situation where dysfunctionality can be more easily observed. Differences in cardiovascular function are not so apparent when individuals are resting or going about their daily lives. Under “stress” conditions, however, differences in functionality become clear and can be diagnosed and subsequently treated. Our moderation results suggest that “stress test” paradigms may also be useful in the study of emotion regulation in daily life. While it may require monitoring individuals for long durations, identification of natural experiments wherein individuals’ adaptive capacities are being pushed to their limits should provide further opportunity to observe differences in emotion regulation.

4.2. Impulse Response Analysis and Person-Specific Intervention

The moderation result also highlights opportunity to use impulse response analysis in studies of intervention. The results of this study suggest that the combination of intensive experience sampling data, network methods, and impulse response analysis could inform the design and deployment of person-specific prevention or intervention [16, 73]. First, based on the information presented in the network, clinicians can discover maladaptive feedback loops and design a targeted treatment plan. For example, consider a case where, after a few weeks of monitoring, the clinician finds that their patient’s socioemotional network contains a feedback loop wherein sadness leads to distant, unfriendly behavior (low communion), which in turn leads to more sadness. This configuration implies that persistence of negative mood might be alleviated through some social skills training, thus breaking the link between sadness and withdrawn, quarrelsome social behavior. Clinicians might also test or demonstrate the probable effect of a particular treatment plan by using impulse response analysis to simulate how different network constructions (current versus ideal) lead to different recovery times. For patient networks that contain multiple maladaptive feedback loops, impulse response analysis could be used to compare potential efficiency of different treatment plans. Of course, experience sampling data collected prior to, during, and after treatment would be especially useful for evaluating, if indeed, the network configuration changed in ways that facilitate functionality and health.

4.3. Limitations and Outlook

The results of this study must be interpreted with respect to some limitations in design and implementation. First, although stratified by age and gender, the sample of persons who provided data were somewhat homogenous. This study was based on a convenience sample that lived in or near a university town and consisted of individuals who were willing to participate in an intensive experience sampling protocol. Before generalizing to the larger population, it will be particularly useful to engage with other populations. Clinical populations, very old persons, or individuals who recently experienced particular types of life events (trauma), for example, may allow for more detailed study of network structure and regulation processes that are under duress. Since our sample had a homogeneous demographic profile (e.g., greater than 90% white, cisgender, and heterosexual), and given differences in depression rates among marginalized populations, it is important to verify these results in more diverse and/or special populations.

The study design used here followed persons intensively during three 21-day periods over one year, with each burst of measurement producing multivariate time-series of about 150 observations, a length that facilitated construction of 13-variable networks. Even so, the number of measurements available did not allow us to study intraindividual change in person-specific network structure within-burst. More dense sampling and/or longer time-series would provide possibilities to segment each time-series into multiple windows and either study how the network changes over shorter time-scales or obtain estimates of the reliability of the recovery time estimates (see, e.g., [74]). This will also be an opportunity to probe deeper into the relation of specific pairs of nodes (e.g., synchronization) and characteristics of a specific node (e.g., recovery time). Our hope is that as new technologies (e.g., wearable sensors) begin delivering more intensive and longer data streams, the opportunities to model more nuanced intraindividual change in network behavior will expand.

The data used here come from a multidisciplinary study that was designed around a select set of substantive domains and sampling procedures. Acknowledging that many actions, thoughts, and feelings are involved in the ongoing emotion regulation processes individuals engage in throughout daily life, we constructed the person-specific networks using 13 continuous-value variables assessed in an experience sampling questionnaire. While this set of variables matched our intent to model individuals as high-dimensional, multivariate dynamic systems, future studies might be either more selective or more inclusive in determining which variables to assess repeatedly and include as part of the dynamic system. More focused networks might be constructed using only emotions. Broader networks might also look to include a broader range of self-perceptions (e.g., self-worth) and cognitions (e.g., perseveration).

As in other work where network methods are being applied to experience sampling data [71, 75, 76], our analysis made use of a three-stage procedure. There are, of course, some risks in using the output from one analysis as input for the next analysis because the uncertainties present in earlier analyses are ignored in subsequent analysis. Generally, it would be better to estimate all the models simultaneously in a single model. To our knowledge, however, this is not yet possible. Therefore, the results should still be interpreted cautiously and conservatively.

We chose to perturb sadness and characterized the impulse response curve by recovery time because persistent sad mood is associated with depression. There could be other ways to perturb the system [77] and alternative characterizations of the impulse response analysis (area under the impulse response analysis curve; [78]). Future research could further examine various aspects of system behaviors when different nodes are perturbed (e.g., a node of social behavior).

4.4. Conclusion

Building upon previous work examining emotion regulation process with network approaches, this paper merged intensive experience sampling data and time-series based network methods to construct person-specific socioemotional networks. The evidence of interconnected networks showed that emotions and social behaviors are indeed working together interaction by interaction. Using recovery time to quantify regulatory efficiency of the socioemotional network, we provide further empirical evidence that the regulatory efficiency of the socioemotional dynamics is associated with depressive symptoms. The evidence of this association showed that the interconnected newtork of emotions and social behaviors are indeed contributing to emotion regulation.

Data Availability

The R code to conduct analytical steps (e.g., uSEM model fitting, impulse response analysis, network visualization) and the visualization of all person-specific networks used to support the findings of this study are included within the supplementary information file(s).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Thanks are due to to the study participants for providing a detailed glimpse of their daily lives for such an extended period of time and to the many research assistants who helped obtain such rich data. This work was supported by the National Institute on Health (RC1 AG035645, R01 HD076994, R24 HD041025, UL1 TR002014, and T32 DA017629), National Science Foundation I/UCRC Center for Healthcare Organization Transformation (CHOT, NSF I/UCRC award #1624727), National Institute on Aging Grant (T32 AG049676), an ISSBD-JJF Mentored Fellowship for Early Career Scholars, and the Penn State Social Science Research Institute.

Supplementary Materials

There are two files in the Supplementary material. The “tutorial” documents the R code to conduct analytical steps (e.g., uSEM model fitting, impulse response analysis, network visualization) and the “411 networks merged file” consists of visualization of all person-specific networks. (Supplementary Materials)