Your browser has javascript turned off or blocked. This will lead to some parts of our website to not work properly or at all. Turn on javascript for best performance.

The browser you are using is not supported by this website. All versions of Internet Explorer are no longer supported, either by us or Microsoft (read more here: https://www.microsoft.com/en-us/microsoft-365/windows/end-of-ie-support).

Please use a modern browser to fully experience our website, such as the newest versions of Edge, Chrome, Firefox or Safari etc.

Anders Irbäck. Photo.

Anders Irbäck

Professor

Anders Irbäck. Photo.

Fitting a function to time-dependent ensemble averaged data

Author

  • Karl Fogelmark
  • Michael A. Lomholt
  • Anders Irbäck
  • Tobias Ambjörnsson

Summary, in English

Time-dependent ensemble averages, i.e., trajectory-based averages of some observable, are of importance in many fields of science. A crucial objective when interpreting such data is to fit these averages (for instance, squared displacements) with a function and extract parameters (such as diffusion constants). A commonly overlooked challenge in such function fitting procedures is that fluctuations around mean values, by construction, exhibit temporal correlations. We show that the only available general purpose function fitting methods, correlated chi-square method and the weighted least squares method (which neglects correlation), fail at either robust parameter estimation or accurate error estimation. We remedy this by deriving a new closed-form error estimation formula for weighted least square fitting. The new formula uses the full covariance matrix, i.e., rigorously includes temporal correlations, but is free of the robustness issues, inherent to the correlated chi-square method. We demonstrate its accuracy in four examples of importance in many fields: Brownian motion, damped harmonic oscillation, fractional Brownian motion and continuous time random walks. We also successfully apply our method, weighted least squares including correlation in error estimation (WLS-ICE), to particle tracking data. The WLS-ICE method is applicable to arbitrary fit functions, and we provide a publically available WLS-ICE software.

Department/s

  • Computational Biology and Biological Physics
  • eSSENCE: The e-Science Collaboration

Publishing year

2018-12-01

Language

English

Publication/Series

Scientific Reports

Volume

8

Issue

1

Document type

Journal article

Publisher

Nature Publishing Group

Topic

  • Probability Theory and Statistics

Status

Published

ISBN/ISSN/Other

  • ISSN: 2045-2322