[논문 리뷰] Rapid Acceleration of the Permutation Test via Slow Random Walks in the Permutation Group

Moo K. Chung, Yixian Wang외 2명|arXiv (Cornell University)|2018. 12. 17.

Bayesian Methods and Mixture Models인용 수 4

한 줄 요약

이 논문은 신경영상에서 정확한 순열 검정을 근사 없이 가속화하기 위해 순열군 위의 새로운 느린 랜덤 워크 전략을 제안한다. 대칭군의 대수적 구조를 활용함으로써, 노트북에서 1시간 이내에 전체 순열 검정을 수행할 수 있으며, 기존 방법보다 100배 이상 빠르다. 이는 HCP 및 ADNI 데이터셋에서 274명의 여성과 182명의 남성, 그리고 138쌍의 일란성 쌍둥이와 79쌍의 이란성 쌍둥이를 포함한 대규모 연구를 가능하게 한다.

ABSTRACT

The permutation test is an often used test procedure in brain imaging. Unfortunately, generating every possible permutation for large-scale brain image datasets such as HCP and ADNI with hundreds images is not practical. Many previous attempts at speeding up the permutation test rely on various approximation strategies such as estimating the tail distribution with known parametric distributions. In this study, we show how to rapidly accelerate the permutation test without any type of approximate strategies by exploiting the underlying algebraic structure of the permutation group. The method is applied to large number of MRIs in two applications: (1) localizing the male and female differences and (2) localizing the regions of high genetic heritability in the sulcal and gyral pattern of the human cortical brain.

연구 동기 및 목표

• 수백 명의 피험자가 포함된 대규모 뇌영상 연구에서 정확한 순열 검정의 계산적 병목 현상을 해결하기 위해.
• 파라메트릭 꼬리 추정이나 제한된 샘플링과 같은 근사 방법에 의존하지 않기 위해.
• 전체 순열 열거가 비현실적인 대규모 연구(예: HCP, ADNI)에서 정확한 추론을 가능하게 하기 위해.
• 군론적 랜덤 워크를 사용하여 이론적으로 탄탄하고 근사 없는 순열 검정 프레임워크를 제공하기 위해.
• 피질의 굴곡 및 고각 패턴에서의 성차이와 유전적 유사성의 유전적 유병률을 탐지하기 위해.

제안 방법

• 모든 $ (m+n)! $ 개의 순열을 열거할 필요 없이, 순열군 $\mathbb{S}_{m+n}$ 위에서 느린 랜덤 워크를 사용하여 순열을 체계적으로 탐색한다.
• 순열군의 대수적 구조를 활용하여, 전체 순열 공간을 훨씬 적은 샘플 수로 높은 확률로 커버할 수 있다.
• 전치와 순환의 복합을 통해 랜덤 워크를 구성함으로써, 에르고딕성과 균일 분포 수렴을 보장한다.
• 각 순열 샘플에서 검정 통계량(예: t통계량, 상관계수)을 계산하고, 더 극단적인 값을 유도하는 순열 비율을 통해 p-값을 추정한다.
• 열화도를 이용한 피질 표면에서의 디퓨전 매핑에 적용되며, 열전도는 굴곡 및 고각 패턴 간 국소 기하학적 유사도를 측정한다.
• 쌍둥이 상관계수 분석에 있어서, 근사 없이 10,000번의 랜덤 워크를 통해 평균 상관계수를 계산하여 헤리티빌리티 지수를 추정한다.

실험 결과

연구 질문

• RQ1 대규모 뇌영상 데이터셋에서 근사 없이 정확한 순열 검정을 크게 가속화할 수 있는가?
• RQ2 대칭군의 대수적 구조를 어떻게 활용하여 효율적으로 대표적인 순열을 생성할 수 있는가?
• RQ3 전체 순열이 계산적으로 비현실적인 대규모 연구(예: 274명의 여성, 182명의 남성)에서 정확한 추론을 가능하게 할 수 있는가?
• RQ4 이 비근사적 랜덤 워크 프레임워크를 통해 쌍둥이 상관계수와 헤리티빌리티를 얼마나 정확하게 추정할 수 있는가?
• RQ5 이 방법은 피질 표면에서의 디퓨전 매핑에서 유도된 복잡하고 비단조화적인 검정 통계량(예: 비단조화적 통계량)에 적용될 수 있는가?

주요 결과

• 이 방법은 노트북에서 40분 만에 50만 개의 순열을 계산하여, 기존 방법을 사용할 경우 18일 이상 소요되는 계산량을 처리하며, 100배의 속도 향상을 달성했다.
• 이 방법은 274명의 여성과 182명의 남성을 포함한 데이터셋에서 정확한 순열 검정를 수행한 최초의 사례로, 보상된 p < 0.05 기준으로 편도엽에서 유의미한 성차이를 국소화했다.
• 10,000번의 랜덤 워크를 통해 97초 만에 쌍둥이 상관계수 지ap을 계산하여, 소수점 3자리 이내로 수렴함을 확인했다.
• 헤리티빌리티 지수(HI) 지도는 고유전성 영역(HI > 0.7)의 70%가 중심 고각 및 굴곡 영역에 국소화되어 있음을 보여주었다.
• 성차이에 대한 최대 t통계량은 7.08이었고, 최소는 -6.44였으며, 임계값 ±4.27 및 4.48에서 보정된 p-값이 0.05 이하였다.
• 이 방법은 고차원 뇌영상 데이터에서 조차도 근사 없이 전체 통계적 추론을 달성할 수 있음을 입증했다.