[논문 리뷰] Markov chain Monte Carlo test of toric homogeneous Markov chains

Akimichi Takemura, Hisayuki Hara|arXiv (Cornell University)|2010. 04. 20.

Markov Chains and Monte Carlo Methods참고 문헌 13인용 수 7

한 줄 요약

이 논문은 토릭 동차 마르코프 체인(THMC) 모델, 즉 표준 동차 마르코프 체인의 캡슐화가 되는 곡선형 지수가족 모델에 대한 조건부 적합도 검정을 위한 마르코프 기반을 개발한다. 임의의 길이의 이 상태 체인과 유한한 상태 공간을 가진 세 시점 체인에 대해 완전한 마르코프 기반을 제공하며, MCMC를 통한 정확한 조건부 추론이 가능하게 하며, 수치 예제를 통해 약한 점근적 근사로 인해 동질성 가설이 기각됨을 보여준다.

ABSTRACT

Markov chain models are used in various fields, such behavioral sciences or econometrics. Although the goodness of fit of the model is usually assessed by large sample approximation, it is desirable to use conditional tests if the sample size is not large. We study Markov bases for performing conditional tests of the toric homogeneous Markov chain model, which is the envelope exponential family for the usual homogeneous Markov chain model. We give a complete description of a Markov basis for the following cases: i) two-state, arbitrary length, ii) arbitrary finite state space and length of three. The general case remains to be a conjecture. We also present a numerical example of conditional tests based on our Markov basis.

연구 동기 및 목표

• 표준 동차 마르코프 체인의 캡슐화 지수가족인 토릭 동차 마르코프 체인(THMC) 모델에 대한 정확한 조건부 검정을 개발하는 것.
• 작은 표본 설정에서 마르코프 체인 모델 검정에 대한 대규모 표본 카이제곱 근사의 한계를 해결하는 것.
• 일반적으로 동차 마르코프 체인 모델이 곡선형 지수가족이 되는 경우, THMC 모델에 대한 완전한 마르코프 기반을 도출하는 것.
• 유도된 마르코프 기반을 사용하여 마르코프 체인 몬테카를로(MCMC)를 통한 정확한 추론을 가능하게 하는 것.
• 일반적인 S와 T에 대한 마르코프 기반의 구조적 성질을 조사하고, 간선-기호 패턴 그래프에서의 확장된 단순 루프를 제안하는 것.

제안 방법

• 자유 매개변수 γi와 βij를 가진 캡슐화 지수가족으로서 THMC 모델을 구성하여 총 확률에 대해서만 정규화가 가능하도록 한다.
• 마르코프 기반 방법론(Diaconis와 Sturmfels)을 사용하여 귀무가설인 동질성 하에서 충분통계량을 유지하는 움직임을 도출한다.
• Takemura와 Aoki의 거리 감소 논증을 적용하여 S=2 및 T≥3인 경우 마르코프 기반의 완전성을 증명한다.
• 소규모 사례(예: S=3, T=5)에서 계산적 검증을 위해 컴퓨터 대수 시스템(4ti2)을 사용한다.
• 간선-기호 패턴 그래프에서의 '확장된 단순 루프'를 정의하고 분석하여 마르코프 기반의 일반적 구조를 추측한다.
• 유도된 마르코프 기반을 사용하여 MCMC 샘플링을 수행하고, 피어슨 카이제곱과 같은 검정 통계량의 정확한 p-값을 계산한다.

실험 결과

연구 질문

• RQ1 상태 공간이 두 개이고 길이가 임의일 때, 토릭 동차 마르코프 체인 모델에 대한 완전한 마르코프 기반은 무엇인가요?
• RQ2 체인이 세 시점이고 유한한 상태 공간 S를 가질 때, THMC 모델에 대한 완전한 마르코프 기반은 무엇인가요?
• RQ3 점근적 근사가 신뢰할 수 없을 경우, MCMC를 사용하여 THMC 모델에 대한 정확한 조건부 검정을 어떻게 수행할 수 있나요?
• RQ4 일반적인 S와 T에 대한 마르코프 기반의 구조적 성질은 무엇이며, 간선-기호 패턴 그래프에서의 확장된 단순 루프로 기술할 수 있나요?
• RQ5 주어진 S에 대해, 마르코프 기반의 모든 움직임을 생성하는 데 필요한 확장된 단순 루프의 수에 대한 균일한 상한이 존재하는가요?

주요 결과

• T≥3인 두 상태 토릭 동차 마르코프 체인에 대해 완전한 마르코프 기반을 도출하였으며, 이는 단순 루프에 기반한 움직임으로 구성된다.
• 유한한 상태 공간 S를 가진 세 시점 체인에 대해 완전한 마르코프 기반을 도출하였으며, 이는 최대 세 개의 확장된 단순 루프를 포함하는 움직임을 사용한다.
• 마리화나 사용 데이터에 대한 MCMC 기반 정확한 검정에서 정확한 p-값은 0.0184로, 점근적 p-값 0.0212보다 略적으로 낮아, THMC 모델이 기각됨을 시사한다.
• 세포 빈도가 낮아서 점근적 카이제곱 근사가 데이터에 대해 부적절함을 확인하였으며, 이는 정확한 방법의 필요성을 강조한다.
• 저자들은 각 S에 대해, 마르코프 기반의 모든 움직임이 최대 kS개의 확장된 단순 루프의 합으로 구성될 수 있는 상한 kS가 존재한다고 추측한다. 이때 kS = S는 강력한 후보이다.
• S=2인 경우 마르코프 기반은 최대 두 개의 확장된 단순 루프의 합으로 이루어지고, S=3, T=4인 경우 움직임은 최대 세 개의 루프를 포함한다.