Physics > Data Analysis, Statistics and Probability
[Submitted on 9 Jul 2021 (v1), last revised 22 Jul 2022 (this version, v3)]
Title:Efficient Bayesian estimation of the generalized Langevin equation from data
View PDFAbstract:Modeling non-Markovian time series is a recent topic of research in many fields such as climate modeling, biophysics, molecular dynamics, or finance. The generalized Langevin equation (GLE), given naturally by the Mori-Zwanzig projection formalism, is a frequently used model including memory effects. In applications, a specific form of the GLE is most often obtained on a data-driven basis. Here, Bayesian estimation has the advantage of providing both suitable model parameters and their credibility in a straightforward way. It can be implemented in the approximating case of white noise, which, far from thermodynamic equilibrium, is consistent with the fluctuation-dissipation theorem. However, the exploration of the posterior, which is done via Markov chain Monte Carlo sampling, is numerically expensive, which makes the analysis of large data sets unfeasible. In this work, we discuss an efficient implementation of Bayesian estimation of the GLE based on a piecewise constant approximation of the drift and diffusion functions of the model. In this case, the characteristics of the data are represented by only a few coefficients, so that the numerical cost of the procedure is significantly reduced and independent of the length of the data set. Further, we propose a modification of the memory term of the GLE, leading to an equivalent model with an emphasis on the impact of trends, which ensures that an estimate of the standard Langevin equation provides an effective initial guess for the GLE. We illustrate the capabilities of both the method and the model by an example from turbulence.
Submission history
From: Clemens Willers [view email][v1] Fri, 9 Jul 2021 17:23:47 UTC (859 KB)
[v2] Thu, 2 Jun 2022 16:31:27 UTC (1,050 KB)
[v3] Fri, 22 Jul 2022 14:44:56 UTC (1,050 KB)
Current browse context:
physics.data-an
Change to browse by:
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.