From pre-processing to advanced dynamic modeling of pupil data.

The pupil of the eye provides a rich source of information for cognitive scientists, as it can index a variety of bodily states (e.g., arousal, fatigue) and cognitive processes (e.g., attention, decision-making). As pupillometry becomes a more accessible and popular methodology, researchers have proposed a variety of techniques for analyzing pupil data. Here, we focus on time series-based, signal-to-signal approaches that enable one to relate dynamic changes in pupil size over time with dynamic changes in a stimulus time series, continuous behavioral outcome measures, or other participants' pupil traces. We first introduce pupillometry, its neural underpinnings, and the relation between pupil measurements and other oculomotor behaviors (e.g., blinks, saccades), to stress the importance of understanding what is being measured and what can be inferred from changes in pupillary activity. Next, we discuss possible pre-processing steps, and the contexts in which they may be necessary. Finally, we turn to signal-to-signal analytic techniques, including regression-based approaches, dynamic time-warping, phase clustering, detrended fluctuation analysis, and recurrence quantification analysis. Assumptions of these techniques, and examples of the scientific questions each can address, are outlined, with references to key papers and software packages. Additionally, we provide a detailed code tutorial that steps through the key examples and figures in this paper. Ultimately, we contend that the insights gained from pupillometry are constrained by the analysis techniques used, and that signal-to-signal approaches offer a means to generate novel scientific insights by taking into account understudied spectro-temporal relationships between the pupil signal and other signals of interest.

Using measures of reading time regularity (RTR) to quantify eye movement dynamics, and how they are shaped by linguistic information.

In this article, we present the concept of reading time regularity (RTR) as a measure to capture reading process dynamics. The first study is concerned with examining one of the assumptions of RTR, namely, that process measures of reading, such as eye movement fluctuations and fixation durations, exhibit higher regularity when contingent on sequentially structured information, such as texts. To test this, eye movements of 26 German native speakers were recorded during reading-unrelated and reading-related tasks. To analyze the data, we used recurrence quantification analysis and sample entropy analysis to quantify the degree of temporal structure in time series of gaze steps and fixation durations. The results showed that eye movements become more regular in reading compared to nonreading conditions. These effects were most prominent when calculated on the basis of gaze step data. In a second study, eye movements of 27 native speakers of German were recorded for five conditions with increasing linguistic information. The results replicate the findings of the first study, verifying that these effects are not due to mere differences in task instructions between conditions. Implications for the concept of RTR and for future studies using these metrics in reading research are discussed.

Four Methods to Distinguish between Fractal Dimensions in Time Series through Recurrence Quantification Analysis.

Fractal properties in time series of human behavior and physiology are quite ubiquitous, and several methods to capture such properties have been proposed in the past decades. Fractal properties are marked by similarities in statistical characteristics over time and space, and it has been suggested that such properties can be well-captured through recurrence quantification analysis. However, no methods to capture fractal fluctuations by means of recurrence-based methods have been developed yet. The present paper takes this suggestion as a point of departure to propose and test several approaches to quantifying fractal fluctuations in synthetic and empirical time-series data using recurrence-based analysis. We show that such measures can be extracted based on recurrence plots, and contrast the different approaches in terms of their accuracy and range of applicability.

Changes in the Complexity of Limb Movements during the First Year of Life across Different Tasks.

Infants' limb movements evolve from disorganized to more selectively coordinated during the first year of life as they learn to navigate and interact with an ever-changing environment more efficiently. However, how these coordination patterns change during the first year of life and across different contexts is unknown. Here, we used wearable motion trackers to study the developmental changes in the complexity of limb movements (arms and legs) at 4, 6, 9 and 12 months of age in two different tasks: rhythmic rattle-shaking and free play. We applied Multidimensional Recurrence Quantification Analysis (MdRQA) to capture the nonlinear changes in infants' limb complexity. We show that the MdRQA parameters (entropy, recurrence rate and mean line) are task-dependent only at 9 and 12 months of age, with higher values in rattle-shaking than free play. Since rattle-shaking elicits more stable and repetitive limb movements than the free exploration of multiple objects, we interpret our data as reflecting an increase in infants' motor control that allows for stable body positioning and easier execution of limb movements. Infants' motor system becomes more stable and flexible with age, allowing for flexible adaptation of behaviors to task demands.

Multidimensional Cross-Recurrence Quantification Analysis (MdCRQA) - A Method for Quantifying Correlation between Multivariate Time-Series.

In this paper, Multidimensional Cross-Recurrence Quantification Analysis (MdCRQA) is introduced. It is an extension of Multidimensional Recurrence Quantification Analysis (MdRQA), which allows to quantify the (auto-)recurrence properties of a single multidimensional time-series. MdCRQA extends MdRQA to bi-variate cases to allow for the quantification of the co-evolution of two multidimensional time-series. Moreover, it is shown how a Diagonal Cross-Recurrence Profile (DCRP) can be computed from the MdCRQA output that allows to capture time-lagged coupling between two multidimensional time-series. The core concepts of these analyses are described, as well as practical aspects of their application. In the supplementary materials to this paper, implementations of MdCRQA and the DCRP as MatLab- and R-functions are provided.

Deriving inferential statistics from recurrence plots: A recurrence-based test of differences between sample distributions and its comparison to the two-sample Kolmogorov-Smirnov test.

Recurrence plots (RPs) have proved to be a very versatile tool to analyze temporal dynamics of time series data. However, it has also been conjectured that RPs can be used to model samples of random variables, that is, data that do not contain any temporal dependencies. In the current paper, we show that RPs can indeed be used to mimic nonparametric inferential statistics. Particularly, we use the case of the two-sample Kolmogorov-Smirnov test as a proof-of-concept, showing how such a test can be done based on RPs. Simulations on differences in mean, variance, and shape of two distributions show that the results of the classical two-sample Kolmogorov-Smirnov test and the recurrence-based test for differences in distributions of two independent samples scale well with each other. While the Kolmogorov-Smirnov test seems to be more sensitive in detecting differences in means, the recurrence based test seems to be more sensitive to detect heteroscedasticity and asymmetry. Potential improvements of our approach as well as extensions to tests with individual distributions are discussed.

Single-Word Recognition Need Not Depend on Single-Word Features: Narrative Coherence Counteracts Effects of Single-Word Features that Lexical Decision Emphasizes.

Research on reading comprehension of connected text emphasizes reliance on single-word features that organize a stable, mental lexicon of words and that speed or slow the recognition of each new word. However, the time needed to recognize a word might not actually be as fixed as previous research indicates, and the stability of the mental lexicon may change with task demands. The present study explores the effects of narrative coherence in self-paced story reading to single-word feature effects in lexical decision. We presented single strings of letters to 24 participants, in both lexical decision and self-paced story reading. Both tasks included the same words composing a set of adjective-noun pairs. Reading times revealed that the tasks, and the order of the presentation of the tasks, changed and/or eliminated familiar effects of single-word features. Specifically, experiencing the lexical-decision task first gradually emphasized the role of single-word features, and experiencing the self-paced story-reading task afterwards counteracted the effect of single-word features. We discuss the implications that task-dependence and narrative coherence might have for the organization of the mental lexicon. Future work will need to consider what architectures suit the apparent flexibility with which task can accentuate or diminish effects of single-word features.

Synchronized arousal between performers and related spectators in a fire-walking ritual.

Collective rituals are present in all known societies, but their function is a matter of long-standing debates. Field observations suggest that they may enhance social cohesion and that their effects are not limited to those actively performing but affect the audience as well. Here we show physiological effects of synchronized arousal in a Spanish fire-walking ritual, between active participants and related spectators, but not participants and other members of the audience. We assessed arousal by heart rate dynamics and applied nonlinear mathematical analysis to heart rate data obtained from 38 participants. We compared synchronized arousal between fire-walkers and spectators. For this comparison, we used recurrence quantification analysis on individual data and cross-recurrence quantification analysis on pairs of participants' data. These methods identified fine-grained commonalities of arousal during the 30-min ritual between fire-walkers and related spectators but not unrelated spectators. This indicates that the mediating mechanism may be informational, because participants and related observers had very different bodily behavior. This study demonstrates that a collective ritual may evoke synchronized arousal over time between active participants and bystanders. It links field observations to a physiological basis and offers a unique approach for the quantification of social effects on human physiology during real-world interactions.

Effectiveness of an online recovery training for employees exposed to blurred boundaries between work and non-work: Bayesian analysis of a randomised controlled trial.

BACKGROUND: Blurred work-non-work boundaries can have negative effects on mental health, including sleep. OBJECTIVES: In a randomised control trial, we aimed to assess the effectiveness of an online recovery training programme designed to improve symptoms of insomnia in a working population exposed to blurred boundaries. METHODS: 128 participants with severe insomnia symptoms (Insomnia Severity Index ≥15) and working under blurred work and non-work conditions (segmentation supplies <2.25) were randomly assigned to either the recovery intervention or a waitlist control group (WLC). The primary outcome was insomnia severity, assessed at baseline, after 2 months (T2) and 6 months (T3). FINDINGS: A greater reduction in insomnia was observed in the intervention compared with the WLC group at both T2 (d=1.51; 95% CI=1.12 o 1.91) and T3 (d=1.63; 95% CI=1.23 to 2.03]. This was shown by Bayesian analysis of covariance (ANCOVA), whereby the ANCOVA model yielded the highest Bayes factor (BF 10=3.23×e60] and a 99.99% probability. Likewise, frequentist analysis revealed significantly reduced insomnia at both T2 and T3. Beneficial effects were found for secondary outcomes including depression, work-related rumination, and mental detachment from work. Study attrition was 16% at T2 and 44% at T3. CONCLUSIONS: The recovery training was effective in reducing insomnia symptoms, work related and general indicators of mental health in employees exposed to blurred boundaries, both at T2 and T3. CLINICAL IMPLICATIONS: In addition to demonstrating the intervention's effectiveness, this study exemplifies the utilisation of the Bayesian approach in a clinical context and shows its potential to empower recipients of interventional research by offering insights into result probabilities, enabling them to draw informed conclusions. TRIAL REGISTRATION NUMBER: German Clinical Trial Registration (DRKS): DRKS00006223, https://drks.de/search/de/trial/DRKS00006223.

Mother-infant social gaze dynamics relate to infant brain activity and word segmentation.

The 'social brain', consisting of areas sensitive to social information, supposedly gates the mechanisms involved in human language learning. Early preverbal interactions are guided by ostensive signals, such as gaze patterns, which are coordinated across body, brain, and environment. However, little is known about how the infant brain processes social gaze in naturalistic interactions and how this relates to infant language development. During free-play of 9-month-olds with their mothers, we recorded hemodynamic cortical activity of ´social brain` areas (prefrontal cortex, temporo-parietal junctions) via fNIRS, and micro-coded mother's and infant's social gaze. Infants' speech processing was assessed with a word segmentation task. Using joint recurrence quantification analysis, we examined the connection between infants' ´social brain` activity and the temporal dynamics of social gaze at intrapersonal (i.e., infant's coordination, maternal coordination) and interpersonal (i.e., dyadic coupling) levels. Regression modeling revealed that intrapersonal dynamics in maternal social gaze (but not infant's coordination or dyadic coupling) coordinated significantly with infant's cortical activity. Moreover, recurrence quantification analysis revealed that intrapersonal maternal social gaze dynamics (in terms of entropy) were the best predictor of infants' word segmentation. The findings support the importance of social interaction in language development, particularly highlighting maternal social gaze dynamics.

Distinguishing between struggling and skilled readers based on their prosodic speech patterns in oral reading: An exploratory study in grades 2 and 4.

The purpose of this study was to examine if prosodic patterns in oral reading derived from Recurrence Quantification Analysis (RQA) could distinguish between struggling and skilled German readers in Grades 2 (n = 67) and 4 (n = 69). Furthermore, we investigated whether models estimated with RQA measures outperformed models estimated with prosodic features derived from prosodic transcription. According to the findings, struggling second graders appear to have a slower reading rate, longer intervals between pauses, and more repetitions of recurrent amplitudes and pauses, whereas struggling fourth graders appear to have less stable pause patterns over time, more pitch repetitions, more similar amplitude patterns over time, and more repetitions of pauses. Additionally, the models with prosodic patterns outperformed models with prosodic features. These findings suggest that the RQA approach provides additional information about prosody that complements an established approach.

A Multivariate Method for Dynamic System Analysis: Multivariate Detrended Fluctuation Analysis Using Generalized Variance.

Fractal fluctuations are a core concept for inquiries into human behavior and cognition from a dynamic systems perspective. Here, we present a generalized variance method for multivariate detrended fluctuation analysis (mvDFA). The advantage of this extension is that it can be applied to multivariate time series and considers intercorrelation between these time series when estimating fractal properties. First, we briefly describe how fractal fluctuations have advanced a dynamic system understanding of cognition. Then, we describe mvDFA in detail and highlight some of the advantages of the approach for simulated data. Furthermore, we show how mvDFA can be used to investigate empirical multivariate data using electroencephalographic recordings during a time-estimation task. We discuss this methodological development within the framework of interaction-dominant dynamics. Moreover, we outline how the availability of multivariate analyses can inform theoretical developments in the area of dynamic systems in human behavior.

Maximal strength measurement: A critical evaluation of common methods-a narrative review.

Measuring maximal strength (MSt) is a very common performance diagnoses, especially in elite and competitive sports. The most popular procedure in test batteries is to test the one repetition maximum (1RM). Since testing maximum dynamic strength is very time consuming, it often suggested to use isometric testing conditions instead. This suggestion is based on the assumption that the high Pearson correlation coefficients of r ≥ 0.7 between isometric and dynamic conditions indicate that both tests would provide similar measures of MSt. However, calculating r provides information about the relationship between two parameters, but does not provide any statement about the agreement or concordance of two testing procedures. Hence, to assess replaceability, the concordance correlation coefficient (ρ c) and the Bland-Altman analysis including the mean absolute error (MAE) and the mean absolute percentage error (MAPE) seem to be more appropriate. Therefore, an exemplary model based on r = 0.55 showed ρ c = 0.53, A MAE of 413.58 N and a MAPE = 23.6% with a range of -1,000-800 N within 95% Confidence interval (95%CI), while r = 0.7 and 0.92 showed ρ c = 0.68 with a MAE = 304.51N/MAPE = 17.4% with a range of -750 N-600 N within a 95% CI and ρ c = 0.9 with a MAE = 139.99/MAPE = 7.1% with a range of -200-450 N within a 95% CI, respectively. This model illustrates the limited validity of correlation coefficients to evaluate the replaceability of two testing procedures. Interpretation and classification of ρ c, MAE and MAPE seem to depend on expected changes of the measured parameter. A MAPE of about 17% between two testing procedures can be assumed to be intolerably high.

The interplay between posture control and memory for spatial locations.

Three experiments examined interactions between posture control in upright stance and a concurrent location memory task. Healthy young participants stood upright and memorized the locations of dots presented on a computer screen. In the retrieval phase, they indicated whether arrows presented on the screen would pass through any of the memorized locations. Postural sway variability was measured either during the retention period or during retrieval. Relative to not performing the memory task, postural sway variability increased in the retention period when the eyes were closed, but remained unaffected when the eyes were open. During retrieval, postural sway variability was reduced relative to the no-memory-task condition. Results were interpreted in terms of dual-task costs associated with maintaining multiple frames of reference.

Modeling items for text comprehension assessment using confirmatory factor analysis.

Reading is a complex cognitive task with the ultimate goal of comprehending the written input. For longer, connected text, readers generate a mental representation that serves as its basis. Due to limited cognitive resources, common models of discourse representation assume distinct processing levels, each relying on different processing mechanisms. However, only little research addresses distinct representational levels when text comprehension is assessed, analyzed or modelled. Moreover, current studies that tried to relate process measures of reading (e.g., reading times, eye movements) to comprehension did not consider comprehension as a multi-faceted, but rather a uni-dimensional construct, usually assessed with one-shot items. Thus, the first aim of this paper is to use confirmatory factor analysis (CFA) to test whether comprehension can be modelled as a uni-or multi-dimensional concept. The second aim is to investigate how well widely used one-shot items can be used to capture comprehension. 400 participants read one of three short stories of comparable length, linguistic characteristics, and complexity. Based on the evaluation of three independent raters per story, 16 wh-questions and 60 yes/no-statements were compiled in order to retrieve information at micro and inference level, and 16 main contents were extracted to capture information at the macro level in participants' summaries. Still, only a fraction of these items showed satisfactory psychometric properties and factor loadings - a blatant result considering the common practice for item selection. For CFA, two models were set up that address text comprehension as either a one-dimensional construct (a uni-factor model with a single comprehension factor), or a three-dimensional construct reflecting the three distinct representational levels (three correlated first-order factors). Across stories and item types, model fit was consistently better for the three-factor model providing evidence for a multi-dimensional construct of text comprehension. Our results provide concrete guidance for the preparation of comprehension measurements in studies investigating the reading process.

Interpersonal Physiological Synchrony Predicts Group Cohesion.

A key emergent property of group social dynamic is synchrony-the coordination of actions, emotions, or physiological processes between group members. Despite this fact and the inherent nested structure of groups, little research has assessed physiological synchronization between group members from a multi-level perspective, thus limiting a full understanding of the dynamics between members. To address this gap of knowledge we re-analyzed a large dataset (N = 261) comprising physiological and psychological data that were collected in two laboratory studies that involved two different social group tasks. In both studies, following the group task, members reported their experience of group cohesion via questionnaires. We utilized a non-linear analysis method-multidimensional recurrence quantification analysis that allowed us to represent physiological synchronization in cardiological interbeat intervals between group members at the individual-level and at the group-level. We found that across studies and their conditions, the change in physiological synchrony from baseline to group interaction predicted a psychological sense of group cohesion. This result was evident both at the individual and the group levels and was not modified by the context of the interaction. The individual- and group-level effects were highly correlated. These results indicate that the relationship between synchrony and cohesion is a multilayered construct. We re-affirm the role of physiological synchrony for cohesion in groups. Future studies are needed to crystallize our understanding of the differences and similarities between synchrony at the individual-level and synchrony at the group level to illuminate under which conditions one of these levels has primacy, or how they interact.

Dataset of physiological, behavioral, and self-report measures from a group decision-making lab study.

This paper presents data from a study conducted in twenty groups of three participants each. Data were collected from sixty participants during a lab visit which was video recorded from several angles. Upon arrival to the lab and following informed consent, participants were told that they would be a part of a group decision-making task and were given instructions for a procedure titled "the desert survival task" Lafferty and Pond (1974). Participants were then connected to several electrodes on their upper body and palm for the collection of their electrocardiogram, respiration and electrodermal activity throughout the group task. Participants then performed the task together. The collection of physiological data from all group members was conducted simultaneously and in synchrony with the video recording. The video recordings of the group interactions were later coded by trained psychology students for positive affective behaviors made by participants (smiling and laughing) throughout the group task. Self-report measures (trait anxiety and social phobia) were collected prior to the group task from all participants. This multimodal dataset thus integrates behavioral, self-report, and physiological measures from group members, which are important for understanding group dynamics. These data will allow verification, replications, and additional analyses of the data from new perspectives.

Group-level physiological synchrony and individual-level anxiety predict positive affective behaviors during a group decision-making task.

Joint performance can lead to the synchronization of physiological processes among group members during a shared task. Recently, it has been shown that synchronization is indicative of subjective ratings of group processes and task performance. However, different methods have been used to quantify synchronization, and little is known about the effects of the choice of method and level of analysis (individuals, dyads, or triads) on the results. In this study, participants performed a decision-making task in groups of three while physiological signals (heart rate and electrodermal activity), positive affective behavior, and personality traits were measured. First, we investigated the effects of different levels of analysis of physiological synchrony on affective behavior. We computed synchrony measures as (a) individual contributions to group synchrony, (b) the average dyadic synchrony within a group, and (c) group-level synchrony. Second, we assessed the association between physiological synchrony and positive affective behavior. Third, we investigated the moderating effects of trait anxiety and social phobia on behavior. We discovered that the effects of physiological synchrony on positive affective behavior were particularly strong at the group level but nonsignificant at the individual and dyadic levels. Moreover, we found that heart rate and electrodermal synchronization showed opposite effects on group members' display of affective behavior. Finally, trait anxiety moderated the relationship between physiological synchrony and affective behavior, perhaps due to social uncertainty, while social phobia did not have a moderating effect. We discuss these results regarding the role of different physiological signals and task demands during joint action.

Bed-Sharing in Couples Is Associated With Increased and Stabilized REM Sleep and Sleep-Stage Synchronization.

BACKGROUND/OBJECTIVES: Sharing the bed with a partner is common among adults and impacts sleep quality with potential implications for mental health. However, hitherto findings are contradictory and particularly polysomnographic data on co-sleeping couples are extremely rare. The present study aimed to investigate the effects of a bed partner's presence on individual and dyadic sleep neurophysiology. METHODS: Young healthy heterosexual couples underwent sleep-lab-based polysomnography of two sleeping arrangements: individual sleep and co-sleep. Individual and dyadic sleep parameters (i.e., synchronization of sleep stages) were collected. The latter were assessed using cross-recurrence quantification analysis. Additionally, subjective sleep quality, relationship characteristics, and chronotype were monitored. Data were analyzed comparing co-sleep vs. individual sleep. Interaction effects of the sleeping arrangement with gender, chronotype, or relationship characteristics were moreover tested. RESULTS: As compared to sleeping individually, co-sleeping was associated with about 10% more REM sleep, less fragmented REM sleep (p = 0.008), longer undisturbed REM fragments (p = 0.0006), and more limb movements (p = 0.007). None of the other sleep stages was significantly altered. Social support interacted with sleeping arrangement in a way that individuals with suboptimal social support showed the biggest impact of the sleeping arrangement on REM sleep. Sleep architectures were more synchronized between partners during co-sleep (p = 0.005) even if wake phases were excluded (p = 0.022). Moreover, sleep architectures are significantly coupled across a lag of ± 5min. Depth of relationship represented an additional significant main effect regarding synchronization, reflecting a positive association between the two. Neither REM sleep nor synchronization was influenced by gender, chronotype, or other relationship characteristics. CONCLUSION: Depending on the sleeping arrangement, couple's sleep architecture and synchronization show alterations that are modified by relationship characteristics. We discuss that these alterations could be part of a self-enhancing feedback loop of REM sleep and sociality and a mechanism through which sociality prevents mental illness.

Switching between reading tasks leads to phase-transitions in reading times in L1 and L2 readers.

Reading research uses different tasks to investigate different levels of the reading process, such as word recognition, syntactic parsing, or semantic integration. It seems to be tacitly assumed that the underlying cognitive process that constitute reading are stable across those tasks. However, nothing is known about what happens when readers switch from one reading task to another. The stability assumptions of the reading process suggest that the cognitive system resolves this switching between two tasks quickly. Here, we present an alternative language-game hypothesis (LGH) of reading that begins by treating reading as a softly-assembled process and that assumes, instead of stability, context-sensitive flexibility of the reading process. LGH predicts that switching between two reading tasks leads to longer lasting phase-transition like patterns in the reading process. Using the nonlinear-dynamical tool of recurrence quantification analysis, we test these predictions by examining series of individual word reading times in self-paced reading tasks where native (L1) and second language readers (L2) transition between random word and ordered text reading tasks. We find consistent evidence for phase-transitions in the reading times when readers switch from ordered text to random-word reading, but we find mixed evidence when readers transition from random-word to ordered-text reading. In the latter case, L2 readers show moderately stronger signs for phase-transitions compared to L1 readers, suggesting that familiarity with a language influences whether and how such transitions occur. The results provide evidence for LGH and suggest that the cognitive processes underlying reading are not fully stable across tasks but exhibit soft-assembly in the interaction between task and reader characteristics.