동적 마코프 기저를 통한 다차원 교차표 샘플링

본 논문은 다차원 교차표의 정확한 조건부 검정을 위해 필요한 표본 생성 문제를 다룬다. 전통적으로는 고정된 주변합(마진)만을 만족하는 표의 집합 T에 대해, 모든 가능한 차이(이동)들을 모아 만든 마코프 기저를 사전에 계산한 뒤, 이를 이용해 마코프 체인을 수행한다. 그러나 마코프 기저는 일반적으로 매우 큰 크기를 가지며, 특히 구조적 영(0)이나 셀별 상하한이 존재하는 경우 그 계산이 조합론적으로 불가능해진다. 이러한 한계를 극복하고자 저자들은 ‘동적 마코프 기저’를 제안한다. 동적 기저의 핵심은 현재 테이블 n*에 대해 필요한 이동을 즉시 생성한다는 점이다. 이를 위해 먼저 셀 인덱스를 임의의 순열 δ에 따라 재배열하고, 순차적으로 각 셀에 할당할 수 있는 정수값 구간을 선형계획법(LP)으로 계산한다. 구간 계산은 실제 정수 해가 아닌 실수 해를 구한 뒤, 바닥값과 천장값을 취해 정수 구간으로 근사한다. 이렇게 하면 각 단계마다 가능한 값들의 집합이 명확히 정의되고, 그 안에서 사전 정의된 이산 분포(f)로 값을 샘플링한다. 알고리즘 1은 이러한 절차를 구체화한 것으로, 셀을 하나씩 채워 나가면서 현재까지 선택된 값들이 전체 제약을 만족하는지 실시간으로 검증한다. 만약 어느 단계에서 허용 구간이 비어 있으면(즉, 현재까지 선택된 값들로는 더 이상 유효한 표를 만들 수 없을 때) 알고리즘은 즉시 종료하고, 다른 순열 δ 혹은 다른 초기값을 시도한다. 이 과정에서 생성된 표 nδ는 Tδ(δ에 의해 재배열된 T)의 한 원소이며, nδ와 동일한 순열에 의해 재배열된 다른 표들 사이의 차이는 ‘국소 이동(local move)’이라 부른다. 이러한 국소 이동들의 집합은 대칭성을 만족하고, 전체 T를 연결하는 마코프 기저 역할을 한다. 즉, 전체 마코프 기저를 미리 구할 필요 없이, 매 반복마다 현재 상태에 맞는 최소한의 이동 집합만을 동적으로 구성한다. 논문은 이론적 정당성을 두 단계로 제시한다. 첫째, 순열 δ에 대해 정의된 이웃 집합 nbd_Tδ(n*δ) = Tδ는 T와 일대일 대응함을 보이며, 따라서 이 집합이 생성하는 이동은 T 내의 모든 가능한 변환을 포함한다. 둘째, 모든 가능한 순열 δ에 대해 얻어지는 이동들의 합집합이 전통적인 마코프 기저와 동등한 연결성을 제공함을 증명한다. 이로써 동적 마코프 기저가 마코프 체인의 에르고딕성(ergodicity)을 보장한다는 것이 입증된다. 실험 부분에서는 두 가지 실제 데이터를 사용한다. 첫 번째는 구조적 영을 포함한 3차원 표이며, 기존 소프트웨어(예: 4ti2, Macaulay2)로는 마코프 기저를 구할 수 없을 정도로 복잡한 제약을 가진다. 동적 기저를 적용한 결과, 수천 번의 체인 반복만으로 정확한 p‑값을 추정할 수 있었으며, 표본의 독립성 및 수렴성을 다양한 진단 지표(예: Gelman‑Rubin, ESS)로 확인했다. 두 번째는 8차원, 2×2×…×2 형태의 매우 희소한 표로, 셀 수가 256개에 불과하지만 제약이 복잡해 전통적인 방법이 전혀 적용되지 못한다. 동적 기저는 이 경우에도 평균 0.2초 내에 하나의 유효한 이웃을 생성했으며, 10⁶번 이상의 체인 반복을 수행해도 메모리 사용량이 일정하게 유지되는 장점을 보였다. 또한 저자들은 R 패키지 형태로 구현된 코드를 부록에 제공한다. 주요 함수는 (1) 순열 δ 생성, (2) 셀별 실수 구간 계산을 위한 LP 솔버 호출, (3) 이산 분포에 따라 셀값 샘플링, (4) 마코프 체인 진행을 위한 메트로폴리스‑헤스팅스 수용률 계산이다. 사용자는 데이터와 제약만 입력하면 자동으로 동적 마코프 기저 기반 MCMC를 실행할 수 있다. 결론적으로, 이 논문은 마코프 기저의 계산적 병목을 완전히 회피하면서도, 고차원·제약이 복잡한 교차표에 대한 정확한 조건부 검정을 가능하게 하는 새로운 프레임워크를 제시한다. 통계적 검정, 베이지안 사후 샘플링, 그리고 데이터 보안(통계적 비공개) 분야에서 구조적 영과 셀 상하한이 흔히 나타나는 상황에 특히 유용할 것으로 기대된다.