V N Smelyanskiy et al 2009 New J. Phys. 11 053012 doi:10.1088/1367-2630/11/5/053012
Recovering 'lost' information in the presence of noise: application to rodent–predator dynamics
V N Smelyanskiy1,5, D G Luchinsky2,3, M M Millonas4 and P V E McClintock3
Show affiliationsA Hamiltonian approach is introduced for the reconstruction of trajectories and models of complex stochastic dynamics from noisy measurements. The method converges even when entire trajectory components are unobservable and the parameters are unknown. It is applied to reconstruct nonlinear models of rodent–predator oscillations in Finnish Lapland and high-Arctic tundra. The projected character of noisy incomplete measurements is revealed and shown to result in a degeneracy of the likelihood function within certain null-spaces. The performance of the method is compared with that of the conventional Markov chain Monte Carlo (MCMC) technique.
-
GENERAL SCIENTIFIC SUMMARY
Introduction and background. Data measured from complex (multidimensional, nonlinear) stochastic dynamical systems are often very hard to interpret because seemingly essential information is missing. For example the data set may be incomplete and/or very little may be known about the underlying dynamical system. For historic time series where measurements cannot be repeated, missing information can be considered 'lost'. Examples arise in climate evolution, predator–prey communities, and coupled matter–radiation systems. The problem is usually compounded by random fluctuations (noise).Main results. The paper describes a solution of this enduring problem based on a novel Hamiltonian formalism derived from stochastic dynamics using maximum likelihood estimation. The key result is our discovery of the Hamiltonian equations underlying the problem through analogy with the Wentzel–Friedlin theory of large fluctuations. We apply the formalism to analyse rodent–predator oscillations for which the predator population could not be observed. Contrary to earlier belief, noise-corrupted measurements of the prey dynamics alone are sufficient for computing both the likelihood of the predator population and the ecological model parameters. The method has advantages over the conventional Markov-chain Monte-Carlo technique.
Wider implications. The dynamics of our exemplar predator–prey community has much in common with e.g. models of economic cycles, the susceptible-exposed-infectious-removed-susceptible (SEIRS) model in epidemiology, and competing modes in lasers. So the method will be of immediate relevance to a broad interdisciplinary community. There is also the promise of reconstructing long-past historical events, e.g. the 1918 influenza epidemic, or even the Black Death.
- PACS
- Subjects
- Dates
-
Issue 5 (May 2009)
Received 17 February 2009
Published 27 May 2009
-
Recovering 'lost' information in the presence of noise: application to rodent–predator dynamics
V N Smelyanskiy et al 2009 New J. Phys. 11 053012
-
The covariant entropy bound in gravitational collapse
Sijie Gao and José P. S. Lemos JHEP04(2004)017
-
A preparative comparison using the 'harmonization' method on standards of synthetic natural gas (CCQM-P87)
G Vargha et al 2008 Metrologia 45 08004
-
On the covariant quantization of type-II superstrings
Sebastian Guttenberg et al JHEP06(2004)030
-
Final Report of CCQM P49: Natural gas types IV and V
Adriaan M H van der Veen et al 2005 Metrologia 42 08004
-
Factors influencing graphene growth on metal surfaces
E Loginova et al 2009 New J. Phys. 11 063046
-
Generation and detection of a spin entanglement in nonequilibrium quantum dots
Stefan Legel et al 2008 New J. Phys. 10 045016
-
Out-of-equilibrium electromagnetic radiation
Julien Serreau JHEP05(2004)078
-
The loop group of E8 and K-theory from 11d
Allan Adams and Jarah Evslin JHEP02(2003)029
-
Final report on CCQM-P48: Uranium isotope ratio measurements in simulated biological/environmental materials
C R Quétel et al 2007 Metrologia 44 08010
Users also read
What's this?- Preparing US community greenhouse gas inventories for climate action plans
- Regional climate consequences of large-scale cool roof and photovoltaic array deployment
- Areal changes of land ecosystems in the Alaskan Yukon River Basin from 1984 to 2008
Related review articles
What's this?- The effective temperature
- Integral characteristics: a key to understanding structure formation in stochastic dynamic systems
- Universal randomness