Convergent cross mapping

Last updated

Convergent cross mapping (CCM) is a statistical test for a cause-and-effect relationship between two variables that, like the Granger causality test, seeks to resolve the problem that correlation does not imply causation. [1] While Granger causality is best suited for purely stochastic systems where the influences of the causal variables are separable (independent of each other), CCM is based on the theory of dynamical systems and can be applied to systems where causal variables have synergistic effects. As such, CCM is specifically aimed to identify linkage between variables that can appear uncorrelated with each other.

Contents

Theory

In the event one has access to system variables as time series observations, Takens' embedding theorem can be applied. Takens' theorem generically proves that the state space of a dynamical system can be reconstructed from a single observed time series of the system, . This reconstructed or shadow manifold is diffeomorphic to the true manifold, , preserving instrinsic state space properties of in .

Convergent Cross Mapping (CCM) leverages a corollary to the Generalized Takens Theorem [2] that it should be possible to cross predict or cross map between variables observed from the same system. Suppose that in some dynamical system involving variables and , causes . Since and belong to the same dynamical system, their reconstructions via embeddings and , also map to the same system.

The causal variable leaves a signature on the affected variable , and consequently, the reconstructed states based on can be used to cross predict values of . CCM leverages this property to infer causality by predicting using the library of points (or vice-versa for the other direction of causality), while assessing improvements in cross map predictability as larger and larger random samplings of are used. If the prediction skill of increases and saturates as the entire is used, this provides evidence that is causally influencing .

Cross mapping is generally asymmetric. If forces unidirectionally, variable will contain information about , but not vice versa. Consequently, the state of can be predicted from , but will not be predictable from .

Algorithm

The basic steps of convergent cross mapping for a variable of length against variable are:

  1. If needed, create the state space manifold from
  2. Define a sequence of library subset sizes ranging from a small fraction of to close to .
  3. Define a number of ensembles to evaluate at each library size.
  4. At each library subset size :
    1. For ensembles:
      1. Randomly select state space vectors from
      2. Estimate from the random subset of using the Simplex state space prediction
      3. Compute the correlation between and
    2. Compute the mean correlation over the ensembles at
  5. The spectrum of versus must exhibit convergence.
  6. Assess significance. One technique is to compare to computed from random realizations (surrogates) of .

Applications

CCM is used to detect if two variables belong to the same dynamical system, for example, can past ocean surface temperatures be estimated from the population data over time of sardines or if there is a causal relationship between cosmic rays and global temperatures. As for the latter it was hypothesised that cosmic rays may impact cloud formation, therefore cloudiness, therefore global temperatures. [3]

Extensions

Extensions to CCM include:

See also

Related Research Articles

In vector calculus, the Jacobian matrix of a vector-valued function of several variables is the matrix of all its first-order partial derivatives. When this matrix is square, that is, when the function takes the same number of variables as input as the number of vector components of its output, its determinant is referred to as the Jacobian determinant. Both the matrix and the determinant are often referred to simply as the Jacobian in literature.

In mathematics, conformal geometry is the study of the set of angle-preserving (conformal) transformations on a space.

In physics, Liouville's theorem, named after the French mathematician Joseph Liouville, is a key theorem in classical statistical and Hamiltonian mechanics. It asserts that the phase-space distribution function is constant along the trajectories of the system—that is that the density of system points in the vicinity of a given system point traveling through phase-space is constant with time. This time-independent density is in statistical mechanics known as the classical a priori probability.

The theory of functions of several complex variables is the branch of mathematics dealing with functions defined on the complex coordinate space , that is, n-tuples of complex numbers. The name of the field dealing with the properties of these functions is called several complex variables, which the Mathematics Subject Classification has as a top-level heading.

<span class="mw-page-title-main">Euler equations (fluid dynamics)</span> Set of quasilinear hyperbolic equations governing adiabatic and inviscid flow

In fluid dynamics, the Euler equations are a set of quasilinear partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

<span class="mw-page-title-main">Takens's theorem</span> Conditions under which a chaotic system can be reconstructed by observation

In the study of dynamical systems, a delay embedding theorem gives the conditions under which a chaotic dynamical system can be reconstructed from a sequence of observations of the state of that system. The reconstruction preserves the properties of the dynamical system that do not change under smooth coordinate changes, but it does not preserve the geometric shape of structures in phase space.

In mathematics, a Killing vector field, named after Wilhelm Killing, is a vector field on a Riemannian manifold that preserves the metric. Killing fields are the infinitesimal generators of isometries; that is, flows generated by Killing fields are continuous isometries of the manifold. More simply, the flow generates a symmetry, in the sense that moving each point of an object the same distance in the direction of the Killing vector will not distort distances on the object.

In mathematics, the jet is an operation that takes a differentiable function f and produces a polynomial, the truncated Taylor polynomial of f, at each point of its domain. Although this is the definition of a jet, the theory of jets regards these polynomials as being abstract polynomials rather than polynomial functions.

In general relativity, the metric tensor is the fundamental object of study. The metric captures all the geometric and causal structure of spacetime, being used to define notions such as time, distance, volume, curvature, angle, and separation of the future and the past.

Time-dependent density-functional theory (TDDFT) is a quantum mechanical theory used in physics and chemistry to investigate the properties and dynamics of many-body systems in the presence of time-dependent potentials, such as electric or magnetic fields. The effect of such fields on molecules and solids can be studied with TDDFT to extract features like excitation energies, frequency-dependent response properties, and photoabsorption spectra.

<span class="mw-page-title-main">Multiple integral</span> Generalization of definite integrals to functions of multiple variables

In mathematics (specifically multivariable calculus), a multiple integral is a definite integral of a function of several real variables, for instance, f(x, y) or f(x, y, z). Physical (natural philosophy) interpretation: S any surface, V any volume, etc.. Incl. variable to time, position, etc.

In mathematics, the (linear) Peetre theorem, named after Jaak Peetre, is a result of functional analysis that gives a characterisation of differential operators in terms of their effect on generalized function spaces, and without mentioning differentiation in explicit terms. The Peetre theorem is an example of a finite order theorem in which a function or a functor, defined in a very general way, can in fact be shown to be a polynomial because of some extraneous condition or symmetry imposed upon it.

In mathematics, a π-system on a set is a collection of certain subsets of such that

Holevo's theorem is an important limitative theorem in quantum computing, an interdisciplinary field of physics and computer science. It is sometimes called Holevo's bound, since it establishes an upper bound to the amount of information that can be known about a quantum state. It was published by Alexander Holevo in 1973.

<span class="mw-page-title-main">Causal sets</span> Approach to quantum gravity using discrete spacetime

The causal sets program is an approach to quantum gravity. Its founding principles are that spacetime is fundamentally discrete and that spacetime events are related by a partial order. This partial order has the physical meaning of the causality relations between spacetime events.

<span class="mw-page-title-main">Spacetime topology</span>

Spacetime topology is the topological structure of spacetime, a topic studied primarily in general relativity. This physical theory models gravitation as the curvature of a four dimensional Lorentzian manifold and the concepts of topology thus become important in analysing local as well as global aspects of spacetime. The study of spacetime topology is especially important in physical cosmology.

<span class="mw-page-title-main">Quantum non-equilibrium</span> Notions in the de Broglie-Bohm theory

Quantum non-equilibrium is a concept within stochastic formulations of the De Broglie–Bohm theory of quantum physics.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

<span class="mw-page-title-main">Causal fermion systems</span> Candidate unified theory of physics

The theory of causal fermion systems is an approach to describe fundamental physics. It provides a unification of the weak, the strong and the electromagnetic forces with gravity at the level of classical field theory. Moreover, it gives quantum mechanics as a limiting case and has revealed close connections to quantum field theory. Therefore, it is a candidate for a unified physical theory. Instead of introducing physical objects on a preexisting spacetime manifold, the general concept is to derive spacetime as well as all the objects therein as secondary objects from the structures of an underlying causal fermion system. This concept also makes it possible to generalize notions of differential geometry to the non-smooth setting. In particular, one can describe situations when spacetime no longer has a manifold structure on the microscopic scale. As a result, the theory of causal fermion systems is a proposal for quantum geometry and an approach to quantum gravity.

Empirical dynamic modeling (EDM) is a framework for analysis and prediction of nonlinear dynamical systems. Applications include population dynamics, ecosystem service, medicine, neuroscience, dynamical systems, geophysics, and human-computer interaction. EDM was originally developed by Robert May and George Sugihara. It can be considered a methodology for data modeling, predictive analytics, dynamical system analysis, machine learning and time series analysis.

References

  1. Sugihara, George; May, Robert; Ye, Hao; Hsieh, Chih-hao; Deyle, Ethan; Fogarty, Michael; Munch, Stephan (2012). "Detecting Causality in Complex Ecosystems". Science. 338 (6106): 496–500. Bibcode:2012Sci...338..496S. doi: 10.1126/science.1227079 . PMID   22997134. S2CID   19749064.
  2. Deyle, Ethan R.; Sugihara, George (2011). "Generalized Theorems for Nonlinear State Space Reconstruction". PLOS ONE. 6 (3): e18295. Bibcode:2011PLoSO...618295D. doi: 10.1371/journal.pone.0018295 . PMC   3069082 . PMID   21483839.
  3. Tsonis, Anastasios A.; Deyle, Ethan R.; Ye, Hao; Sugihara, George (2018), Tsonis, Anastasios A. (ed.), "Convergent Cross Mapping: Theory and an Example", Advances in Nonlinear Geosciences, Cham: Springer International Publishing, pp. 587–600, doi:10.1007/978-3-319-58895-7_27, ISBN   978-3-319-58895-7 , retrieved 2023-10-19
  4. Ye, Hao; Deyle, Ethan R.; Gilarranz, Luis J.; Sugihara, George (2015). "Distinguishing time-delayed causal interactions using convergent cross mapping". Scientific Reports. 5: 14750. Bibcode:2015NatSR...514750Y. doi:10.1038/srep14750. PMC   4592974 . PMID   26435402.
  5. Breston, Leo; Leonardis, Eric J.; Quinn, Laleh K.; Tolston, Michael; Wiles, Janet; Chiba, Andrea A. (2021). "Convergent cross sorting for estimating dynamic coupling". Scientific Reports. 11 (1): 20374. Bibcode:2021NatSR..1120374B. doi:10.1038/s41598-021-98864-2. PMC   8514556 . PMID   34645847. S2CID   238859361.

Further reading

Animations: