Adaptive, locally-linear models of complex dynamics

The dynamics of complex systems generally include high-dimensional, non-stationary and non-linear behavior, all of which pose fundamental challenges to quantitative understanding. To address these difficulties we detail a new approach based on local linear models within windows determined adaptively from the data. While the dynamics within each window are simple, consisting of exponential decay, growth and oscillations, the collection of local parameters across all windows provides a principled characterization of the full time series. To explore the resulting model space, we develop a novel likelihood-based hierarchical clustering and we examine the eigenvalues of the linear dynamics. We demonstrate our analysis with the Lorenz system undergoing stable spiral dynamics and in the standard chaotic regime. Applied to the posture dynamics of the nematode $C. elegans$ our approach identifies fine-grained behavioral states and model dynamics which fluctuate close to an instability boundary, and we detail a bifurcation in a transition from forward to backward crawling. Finally, we analyze whole-brain imaging in $C. elegans$ and show that the stability of global brain states changes with oxygen concentration.Comment: 25 pages, 16 figure

Capturing the nonlinear dynamics of animal behavior with Applications to the Nematode C. elegans

From microorganisms to humans, animals behave by making complex changes in their shape and posture over time with remarkable flexibility. To deal with the complexity of animal behavior existing analysis methods view it as a discrete time process, which is composed of transitions between a finite number of stereotyped motifs, such as walking or reaching. This viewpoint, however, ignores the fact that most behavior is not stereotyped. There is, therefore, a need for a perspective that captures the continuous complexity of animal behavior and offers detailed insights into general principles underlying its generation and control. In my Ph.D. thesis, I propose a new approach of analyzing animal behavior, based on the idea that it is fundamentally a continuous time spatiotemporal dynamical system. I develop methods to transform behavioral recordings into a geometrical object called the "behavioral state space". As an organism moves, the corresponding behavioral state traces out a continuous trajectory in the state space, such that the geometry and topology of the trajectories encode quantitative and qualitative properties of behavior. Finally, I characterize an organism\u27s behavioral dynamics in terms of the topological invariants estimated from the local Jacobians of the state space trajectories. The invariants capture essential aspects of a dynamical system, such as the number of degrees of freedom, symmetries in the governing equations of motion, and measures of predictability and variability. I use the tools and concepts developed the above to perform a detailed characterization of continuous dynamics of freely behaving \textit{C. elegans} worms.Okinawa Institute of Science and Technology Graduate Universit

A Markovian dynamics for C.elegansC. elegans behavior across scales

How do we capture the breadth of behavior in animal movement, from rapid body twitches to aging? Using high-resolution videos of the nematode worm $C. elegans$, we show that a single dynamics connects posture-scale fluctuations with trajectory diffusion, and longer-lived behavioral states. We take short posture sequences as an instantaneous behavioral measure, fixing the sequence length for maximal prediction. Within the space of posture sequences we construct a fine-scale, maximum entropy partition so that transitions among microstates define a high-fidelity Markov model, which we also use as a means of principled coarse-graining. We translate these dynamics into movement using resistive force theory, capturing the statistical properties of foraging trajectories. Predictive across scales, we leverage the longest-lived eigenvectors of the inferred Markov chain to perform a top-down subdivision of the worm's foraging behavior, revealing both ``runs-and-pirouettes'' as well as previously uncharacterized finer-scale behaviors. We use our model to investigate the relevance of these fine-scale behaviors for foraging success, recovering a trade-off between local and global search strategies.Comment: 28 pages, 14 figure

Capturing the Continuous Complexity of Behavior in C. elegans

Animal behavior is often quantified through subjective, incomplete variables that may mask essential dynamics. Here, we develop a behavioral state space in which the full instantaneous state is smoothly unfolded as a combination of short-time posture dynamics. Our technique is tailored to multivariate observations and extends previous reconstructions through the use of maximal prediction. Applied to high-resolution video recordings of the roundworm \textit{C. elegans}, we discover a low-dimensional state space dominated by three sets of cyclic trajectories corresponding to the worm's basic stereotyped motifs: forward, backward, and turning locomotion. In contrast to this broad stereotypy, we find variability in the presence of locally-unstable dynamics, and this unpredictability shows signatures of deterministic chaos: a collection of unstable periodic orbits together with a positive maximal Lyapunov exponent. The full Lyapunov spectrum is symmetric with positive, chaotic exponents driving variability balanced by negative, dissipative exponents driving stereotypy. The symmetry is indicative of damped, driven Hamiltonian dynamics underlying the worm's movement control.Comment: 26 pages, 14 figure

25th annual computational neuroscience meeting: CNS-2016

The same neuron may play different functional roles in the neural circuits to which it belongs. For example, neurons in the Tritonia pedal ganglia may participate in variable phases of the swim motor rhythms [1]. While such neuronal functional variability is likely to play a major role the delivery of the functionality of neural systems, it is difficult to study it in most nervous systems. We work on the pyloric rhythm network of the crustacean stomatogastric ganglion (STG) [2]. Typically network models of the STG treat neurons of the same functional type as a single model neuron (e.g. PD neurons), assuming the same conductance parameters for these neurons and implying their synchronous firing [3, 4]. However, simultaneous recording of PD neurons shows differences between the timings of spikes of these neurons. This may indicate functional variability of these neurons. Here we modelled separately the two PD neurons of the STG in a multi-neuron model of the pyloric network. Our neuron models comply with known correlations between conductance parameters of ionic currents. Our results reproduce the experimental finding of increasing spike time distance between spikes originating from the two model PD neurons during their synchronised burst phase. The PD neuron with the larger calcium conductance generates its spikes before the other PD neuron. Larger potassium conductance values in the follower neuron imply longer delays between spikes, see Fig. 17.Neuromodulators change the conductance parameters of neurons and maintain the ratios of these parameters [5]. Our results show that such changes may shift the individual contribution of two PD neurons to the PD-phase of the pyloric rhythm altering their functionality within this rhythm. Our work paves the way towards an accessible experimental and computational framework for the analysis of the mechanisms and impact of functional variability of neurons within the neural circuits to which they belong

Capturing the Continuous Complexity of Behavior in C. elegans

Animal behavior is often quantified through subjective, incomplete variables that may mask essential dynamics. Here, we develop a behavioral state space in which the full instantaneous state is smoothly unfolded as a combination of short-time posture dynamics. Our technique is tailored to multivariate observations and extends previous reconstructions through the use of maximal prediction. Applied to high-resolution video recordings of the roundworm C. elegans , we discover a low-dimensional state space dominated by three sets of cyclic trajectories corresponding to the worm’s basic stereotyped motifs: forward, backward, and turning locomotion. In contrast to this broad stereotypy, we find variability in the presence of locally-unstable dynamics, and this unpredictability shows signatures of deterministic chaos: a collection of unstable periodic orbits together with a positive maximal Lyapunov exponent. The full Lyapunov spectrum is symmetric with positive, chaotic exponents driving variability balanced by negative, dissipative exponents driving stereotypy. The symmetry is indicative of damped, driven Hamiltonian dynamics underlying the worm’s movement control

Capturing the continuous complexity of behaviour in Caenorhabditis elegans

Animal behaviour is often quantified through subjective, incomplete variables that mask essential dynamics. Here, we develop a maximally predictive behavioural-state space from multivariate measurements, in which the full instantaneous state is smoothly unfolded as a combination of short-time posture sequences. In the off-food behaviour of the roundworm Caenorhabditis elegans, we discover a low-dimensional state space dominated by three sets of cyclic trajectories corresponding to the worm’s basic stereotyped motifs: forward, backward and turning locomotion. We find similar results in the on-food behaviour of foraging worms and npr-1 mutants. In contrast to this broad stereotypy, we find variability in the presence of locally unstable dynamics with signatures of deterministic chaos: a collection of unstable periodic orbits together with a positive maximal Lyapunov exponent. The full Lyapunov spectrum is symmetric with positive, chaotic exponents driving variability balanced by negative, dissipative exponents driving stereotypy. The symmetry is indicative of damped–driven Hamiltonian dynamics underlying the worm’s movement control

Maximally predictive states: From partial observations to long timescales

Isolating slower dynamics from fast fluctuations has proven remarkably powerful, but how do we proceed from partial observations of dynamical systems for which we lack underlying equations? Here, we construct maximally predictive states by concatenating measurements in time, partitioning the resulting sequences using maximum entropy, and choosing the sequence length to maximize short-time predictive information. Transitions between these states yield a simple approximation of the transfer operator, which we use to reveal timescale separation and long-lived collective modes through the operator spectrum. Applicable to both deterministic and stochastic processes, we illustrate our approach through partial observations of the Lorenz system and the stochastic dynamics of a particle in a double-well potential. We use our transfer operator approach to provide a new estimator of the Kolmogorov-Sinai entropy, which we demonstrate in discrete and continuous-time systems, as well as the movement behavior of the nematode worm C. elegans