공통 대각선 효과 검정을 위한 마코프 기저 구축

본 논문은 교차표 분석에서 대각선 셀들의 빈도가 비대각선에 비해 현저히 크거나 작을 때 사용되는 quasi‑independence 모델에 대한 간소화된 버전인 공통 대각선 효과 모델(Common Diagonal Effect Model, CDEM)을 조건부 검정하는 방법을 제시한다. CDEM은 모든 대각선 셀에 동일한 추가 파라미터 γ를 부여함으로써, 각 대각선 셀마다 별도의 파라미터를 두는 quasi‑independence 모델보다 파라미터 수를 크게 줄인다. 이 모델의 충분통계는 행합, 열합, 그리고 대각선 셀들의 총합 x_S 로 구성된다. 이러한 충분통계가 고정된 상황에서, 표들의 전체 공간(섬유, fiber)을 연결하는 마코프 연쇄를 설계하려면 적절한 마코프 기저가 필요하다. 기존 연구에서는 행·열 합만 고정된 경우, “square‑free” 2차 이동(형태 +1,−1,−1,+1)만으로도 섬유를 연결할 수 있음을 보였다. 그러나 대각선 합까지 고정하면, 이러한 2차 이동만으로는 섬유가 분리되는 현상이 발생한다. 이는 “subtable sum problem”의 일종으로, 특히 대각선 셀 집합 S가 구조적 영(zero) 셀인 경우와 유사하지만, 일반적인 양·음 값이 허용된다는 점에서 차이가 있다. 따라서 저자는 새로운 이동들을 도입해 이 문제를 해결한다. 새로운 이동은 총 여섯 종류로 구분된다. 1. **Type I**: 기존의 2차 square‑free 이동으로, 네 개의 서로 다른 행·열 인덱스를 사용한다. 2. **Type II I**: 3차 이동으로, 세 개의 서로 다른 대각선 인덱스(i, i′, i″)를 선택해 두 개의 양·음 쌍을 교환한다. 이 이동은 “dispensable”이며, 최소 기저를 구성할 때 각 (i,i′,i″) 삼중에 대해 두 개만 선택하면 된다. 3. **Type III**: 3차 이동으로, 두 대각선 인덱스와 하나의 비대각선 인덱스를 결합한다. 4. **Type IV**: 3차 이동으로, 두 행 인덱스와 두 열 인덱스를 모두 서로 다르게 선택한다. 이 이동은 “indispensable”이며, 모든 섬유 연결에 필수적이다. 5. **Type V**: 4차 이동으로, 두 대각선 셀을 포함하고 나머지 두 셀을 비대각선에 배치한다(비‑square‑free). 6. **Type VI**: 4차 square‑free 이동으로, 네 개의 서로 다른 행·열 인덱스를 사용한다. **Theorem 1**에서는 min(R,C)≥3, max(R,C)≥4인 경우 위의 Type I–VI 이동이 대각선 합 제약을 만족하는 모든 표들을 연결한다는 것을 증명한다. 증명은 두 표 X와 Y의 차이 행렬 Z=X−Y에 대해 L1 거리(|Z|₁) 감소를 보이는 이동을 단계별로 찾아가는 방식이다. Z의 부호 패턴을 전부 경우별로 분석하고, 각 경우에 적절한 Type I–VI 이동을 적용하면 |Z|₁이 감소한다. 결국 Z가 영 행렬이 되므로 X와 Y가 연결된다. 이 과정에서 Type II I와 Type IV, V, VI 이동이 “indispensable”함을 확인하고, 최소 기저를 구성하기 위해서는 Type II I 이동을 각 (i,i′,i″)에 대해 두 개만 선택하면 충분함을 도출한다. 따라서 **Corollary 1**에서는 최소 마코프 기저가 Type I, II, IV, V, VI 및 각 삼중에 대한 두 개의 Type II I 이동으로 이루어진다고 명시한다. 알고리즘적 구현에서는 위에서 정의한 마코프 기저를 이용해 MCMC 샘플링을 수행한다. 기대 빈도는 IPF(iterative proportional fitting) 알고리즘을 통해 quasi‑independence 모델과 CDEM 각각에 대해 계산한다. 검정 통계량은 로그우도비 G²=2∑_{i,j} x_{ij} log( m̂_QI_{ij} / m̂_S_{ij} ) 로 정의하고, 10,000개의 표를 8,000번의 버닝인 후 수집한다. 실제 데이터(예: 병리학자 진단 교차표)에서는 관측된 G² 값이 13.5505이며, χ² 근사 p‑값은 0.003585이다. MCMC 기반 p‑값은 0.00379로 거의 일치한다. 따라서 CDEM은 5% 유의수준에서 기각된다. 두 번째 예시에서도 동일한 절차를 적용해 CDEM의 적합도를 평가한다. 논문의 주요 기여는 다음과 같다. - **마코프 기저의 명시적 구성**: 대각선 합 제약을 포함한 로그선형 모델에 대해, 기존 2차 이동만으로는 불가능했던 연결성을 확보하기 위해 3차·4차 이동을 포함한 여섯 종류의 이동을 제시하고, 이들이 최소 마코프 기저를 이룸을 증명하였다. - **최소성 증명**: 불필요한 이동을 배제하고, 각 이동의 “indispensable”·“dispensable” 특성을 구분함으로써 최소 기저를 명확히 정의하였다. 이는 실제 MCMC 구현 시 연산 효율성을 크게 향상시킨다. - **조건부 검정 프레임워크 제공**: 알제브라적 접근과 MCMC 샘플링을 결합해, 대각선 효과가 있는지 여부를 전통적인 χ² 근사와 비교해 보다 정확한 p‑값을 추정할 수 있는 실용적인 절차를 제시하였다. - **일반화 가능성**: 대각선 합 제약은 구조적 영 셀 문제의 일반화 형태이므로, 본 연구에서 제시한 기저 설계 방법은 다른 서브테이블 합 고정 모델이나 구조적 영 셀을 포함하는 다양한 로그선형 모델에도 적용 가능하다. 결론적으로, 본 논문은 대각선 효과를 단일 파라미터로 요약하는 CDEM에 대한 조건부 적합도 검정을 알제브라적 통계와 MCMC 방법을 통해 체계화했으며, 이를 통해 복잡한 로그선형 모델의 검정 문제를 실용적으로 해결할 수 있음을 보여준다.