Zeraati, R, Engel, TA, Levina, A (March 2022) A flexible Bayesian framework for unbiased estimation of timescales. Nature Computational Science, 2 (3). pp. 193-204. ISSN 2662-8457
PDF
2022.Zeraati.Bayesian_framework.pdf Available under License Creative Commons Attribution. Download (4MB) |
Abstract
Timescales characterize the pace of change for many dynamic processes in nature. They are usually estimated by fitting the exponential decay of data autocorrelation in the time or frequency domain. Here we show that this standard procedure often fails to recover the correct timescales due to a statistical bias arising from the finite sample size. We develop an alternative approach to estimate timescales by fitting the sample autocorrelation or power spectrum with a generative model based on a mixture of Ornstein–Uhlenbeck processes using adaptive approximate Bayesian computations. Our method accounts for finite sample size and noise in data and returns a posterior distribution of timescales that quantifies the estimation uncertainty and can be used for model selection. We demonstrate the accuracy of our method on synthetic data and illustrate its application to recordings from the primate cortex. We provide a customizable Python package that implements our framework via different generative models suitable for diverse applications.
Item Type: | Paper |
---|---|
Subjects: | bioinformatics > computational biology neurobiology > neuroscience organism description > animal > mammal > primates |
CSHL Authors: | |
Communities: | CSHL labs > Engel lab |
SWORD Depositor: | CSHL Elements |
Depositing User: | CSHL Elements |
Date: | 1 March 2022 |
Date Deposited: | 13 Apr 2022 22:05 |
Last Modified: | 15 Nov 2023 19:25 |
PMCID: | PMC9835171 |
URI: | https://repository.cshl.edu/id/eprint/40582 |
Actions (login required)
Administrator's edit/view item |