pith. sign in

arxiv: 1706.04961 · v3 · pith:ELGI2KHKnew · submitted 2017-06-15 · ❄️ cond-mat.stat-mech · cond-mat.soft· physics.data-an

Fast Bayesian inference of the multivariate Ornstein-Uhlenbeck process

classification ❄️ cond-mat.stat-mech cond-mat.softphysics.data-an
keywords bayesianprocessornstein-uhlenbeckestimateestimatesmatricesmethodmultivariate
0
0 comments X
read the original abstract

The multivariate Ornstein-Uhlenbeck process is used in many branches of science and engineering to describe the regression of a system to its stationary mean. Here we present an $O(N)$ Bayesian method to estimate the drift and diffusion matrices of the process from $N$ discrete observations of a sample path. We use exact likelihoods, expressed in terms of four sufficient statistic matrices, to derive explicit maximum a posteriori parameter estimates and their standard errors. We apply the method to the Brownian harmonic oscillator, a bivariate Ornstein-Uhlenbeck process, to jointly estimate its mass, damping, and stiffness and to provide Bayesian estimates of the correlation functions and power spectral densities. We present a Bayesian model comparison procedure, embodying Ockham's razor, to guide a data-driven choice between the Kramers and Smoluchowski limits of the oscillator. These provide novel methods of analyzing the inertial motion of colloidal particles in optical traps.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Stochastic problems in pulsar timing

    astro-ph.HE 2026-04 unverdicted novelty 5.0

    Analytical solutions to Langevin equations for red noise and GWB in pulsars show that an Ornstein-Uhlenbeck spin frequency model is inconsistent with stationary signals, while an overdamped oscillator model and a two-...