다변량 오른슈타인‑우렌벡 과정의 초고속 베이지안 추정

초록

다변량 오른슈타인‑우렌벡(OU) 과정을 N개의 이산 관측값만으로 O(N) 시간에 드리프트와 확산 행렬을 정확히 추정하는 베이지안 방법을 제시한다. 네 개의 충분통계 행렬을 이용해 MAP 추정식과 표준 오차를 도출하고, 브라운 운동 조화진동자를 사례로 질량·감쇠·강성 추정을 수행한다. 또한 모델 선택을 위한 베이지안 증거와 오컴 팩터를 이용해 Kramers와 Smoluchowski 한계 모델을 데이터 기반으로 구분한다.

상세 분석

본 논문은 연속적인 가우시안·마코프 과정인 다변량 OU 프로세스의 파라미터 추정을, 관측 데이터의 길이 N에 대해 선형 시간 복잡도 O(N)으로 수행할 수 있는 베이지안 프레임워크를 구축한다. 핵심은 전이 확률 밀도와 정규분포 형태의 정상 상태 분포가 정확히 계산될 수 있다는 점이다. 저자들은 이 확률 모델을 이용해 전체 로그 사후 확률을 네 개의 충분통계 행렬 T₁~T₄(각각 한 시점 간 자기·상호 공분산 및 초기 공분산)로 표현한다. 이러한 행렬은 데이터 스트림을 한 번 순회하면서 O(N) 연산으로 누적 가능하므로, 메모리와 계산량이 크게 절감된다.

사후 확률을 미분해 MAP 해를 구하면, 회귀 행렬의 지수 형태 Λ* = T₂ T₃⁻¹, 그리고 공분산 Σ* = (1/N)