수치 적분 오류가 만든 이중극대 장기 모델에서의 비정상적 사후분포

초록

본 논문은 연속시간 미분방정식 기반 장기 모델에서 수치 적분(RK4 등) 오류가 파라미터 추정에 비정상적인 비식별성을 야기하여 사후분포가 두 개의 뚜렷한 모드(정답에 가까운 모드와 오류에 의한 모드)로 분리되는 현상을 이론과 시뮬레이션으로 입증한다. 또한 이러한 다중극성을 사전에 탐지하고, 적절한 MCMC 설정이나 적분 스킴 변경을 통해 회피할 수 있는 워크플로우를 제시한다.

상세 분석

이 연구는 선형 1차 ODE (dY/dt = \beta_0 - \beta_1 Y) 에 대해 RK4 적분을 적용했을 때, 동일한 관측 데이터에 대해 서로 다른 파라미터 조합 ((\beta_0,\beta_1))이 거의 동일한 수치 해를 생성한다는 사실을 밝혀냈다. 핵심은 ‘비식별성’이 전통적인 측정오차나 데이터 부족이 아니라, 수치 적분 자체가 모델의 확률적 부분이 되면서 발생한다는 점이다. 논문은 먼저 해석적으로 RK4 한 스텝을 전개하여
\