English

Fast Bayesian inference of the multivariate Ornstein-Uhlenbeck process

Statistical Mechanics 2018-08-01 v3 Soft Condensed Matter Data Analysis, Statistics and Probability

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)O(N) Bayesian method to estimate the drift and diffusion matrices of the process from NN 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.

Keywords

Cite

@article{arxiv.1706.04961,
  title  = {Fast Bayesian inference of the multivariate Ornstein-Uhlenbeck process},
  author = {Rajesh Singh and Dipanjan Ghosh and R. Adhikari},
  journal= {arXiv preprint arXiv:1706.04961},
  year   = {2018}
}

Comments

add published version

R2 v1 2026-06-22T20:19:59.943Z