A flexible Bayesian framework for unbiased estimation of timescales

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

[thumbnail of 2022.Zeraati.Bayesian_framework.pdf] 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 Administrator's edit/view item