Cargando…

Fitting a function to time-dependent ensemble averaged data

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...

Descripción completa

Detalles Bibliográficos
Autores principales: Fogelmark, Karl, Lomholt, Michael A., Irbäck, Anders, Ambjörnsson, Tobias
Formato: Online Artículo Texto
Lenguaje:English
Publicado: Nature Publishing Group UK 2018
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5934400/
https://www.ncbi.nlm.nih.gov/pubmed/29725108
http://dx.doi.org/10.1038/s41598-018-24983-y
_version_ 1783320105313107968
author Fogelmark, Karl
Lomholt, Michael A.
Irbäck, Anders
Ambjörnsson, Tobias
author_facet Fogelmark, Karl
Lomholt, Michael A.
Irbäck, Anders
Ambjörnsson, Tobias
author_sort Fogelmark, Karl
collection PubMed
description 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.
format Online
Article
Text
id pubmed-5934400
institution National Center for Biotechnology Information
language English
publishDate 2018
publisher Nature Publishing Group UK
record_format MEDLINE/PubMed
spelling pubmed-59344002018-05-10 Fitting a function to time-dependent ensemble averaged data Fogelmark, Karl Lomholt, Michael A. Irbäck, Anders Ambjörnsson, Tobias Sci Rep Article 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. Nature Publishing Group UK 2018-05-03 /pmc/articles/PMC5934400/ /pubmed/29725108 http://dx.doi.org/10.1038/s41598-018-24983-y Text en © The Author(s) 2018 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
spellingShingle Article
Fogelmark, Karl
Lomholt, Michael A.
Irbäck, Anders
Ambjörnsson, Tobias
Fitting a function to time-dependent ensemble averaged data
title Fitting a function to time-dependent ensemble averaged data
title_full Fitting a function to time-dependent ensemble averaged data
title_fullStr Fitting a function to time-dependent ensemble averaged data
title_full_unstemmed Fitting a function to time-dependent ensemble averaged data
title_short Fitting a function to time-dependent ensemble averaged data
title_sort fitting a function to time-dependent ensemble averaged data
topic Article
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5934400/
https://www.ncbi.nlm.nih.gov/pubmed/29725108
http://dx.doi.org/10.1038/s41598-018-24983-y
work_keys_str_mv AT fogelmarkkarl fittingafunctiontotimedependentensembleaverageddata
AT lomholtmichaela fittingafunctiontotimedependentensembleaverageddata
AT irbackanders fittingafunctiontotimedependentensembleaverageddata
AT ambjornssontobias fittingafunctiontotimedependentensembleaverageddata