[논문 리뷰] New procedures for discrete tests with proven false discovery rate control

Sebastian Döhler, Guillermo Durand외 1명|arXiv (Cornell University)|2017. 06. 18.

Statistical Methods in Clinical Trials인용 수 1

한 줄 요약

이 논문은 이산 검정 통계량을 가진 다중 가설 검정에 대한 새로운 FDR 제어 절차를 제안한다. 이론적 FDR 제어를 보장하기 위해 이산 분포의 구조와 신호 강도(π₀-적응)를 통합한 제안된 방법들은 고정된 가설 수에 대해 어떤 경우에도 이론적으로 FDR 제어를 보장하며, 시뮬레이션과 실제 데이터에서 기존 방법들에 비해 유의미하게 높은 검정력(power)을 보인다.

ABSTRACT

We consider the problem of testing simultaneously many null hypotheses when the test statistics have a discrete distribution. We present new modifications of the Benjamini-Hochberg procedure that incorporate the discrete structure of the data in an appropriate way. These new procedures are theoretically proved to control the false discovery rate (FDR) for any fixed number of null hypotheses. A strong point of our FDR controlling methodology is that it allows to incorporate at once the discreteness and the quantity of signal of the data (so called $\pi_0$-adaptation ). Finally, the power advantage of the new methods is demonstrated by using both numerical experiments and real data sets.

연구 동기 및 목표

• 검정 통계량이 이산 분포를 따를 경우 표준 FDR 절차의 성능 저하 문제를 해결하기 위해 다중 가설 검정의 과제를 다루는 것.
• 데이터 분포의 이산성 구조를 명시적으로 고려하는 FDR 제어 방법을 개발하는 것.
• 진정한 귀무가설의 비율(π₀)을 절차에 통합하여 검정력 향상을 위해 π₀-적응을 가능하게 하는 것.
• 임의의 고정된 수의 가설에 대해 이론적 FDR 제어를 보장하는 것 — 이는 渐近적 제어보다 더 강력한 보장이다.
• 수치 실험과 실제 데이터 응용을 통해 뛰어난 통계적 검정력을 입증하는 것.

제안 방법

• 이산 통계량의 특성을 반영하여 p-값 임계값을 조정하는 수정된 벤자민-호크베르그 절차를 제안한다.
• 이산 검정 통계량의 정확한 분포를 사용하여 FDR 제어를 유지하는 임계값을 계산한다.
• 진정한 귀무가설의 비율(π₀)의 추정치를 단계적 상향 절차에 통합하여 신호 강도에 적응하도록 한다.
• 각 이산 점에서의 확률 질량을 고려한 조건부 오류율 접근법을 적용하여 기각 임계값을 정밀하게 조정한다.
• 이산 누적분포함수를 사용하여 각 단계에서 임계값을 조정하는 수정된 단계적 상향 절차를 구현한다.
• 임의의 고정된 수의 가설에 대해 기대되는 잘못된 기각 비율이 α 이하로 유지됨을 증명함으로써 FDR 제어를 보장한다.

실험 결과

연구 질문

• RQ1 검정 통계량이 이산적일 경우 다중 검정에서 FDR 제어를 어떻게 유지할 수 있는가?
• RQ2 검정 통계량의 이산 분포를 고려함으로써 FDR 제어 절차를 개선할 수 있는가?
• RQ3 π₀-적응은 이산 다중 검정 상황에서 검정력에 얼마나 기여하는가?
• RQ4 이산 분포 하에서 임의의 고정된 수의 가설에 대해 이론적 FDR 제어를 보장할 수 있는가?
• RQ5 실제 데이터와 시뮬레이션된 이산 데이터에서 새로운 절차는 기존 방법들에 비해 어떤 정도의 검정력 향상을 보이는가?

주요 결과

• 제안된 절차는 임의의 고정된 수의 가설에 대해 이론적 FDR 제어를 달성하며, 이는 渐近적 제어보다 더 강력한 보장이다.
• 검정 통계량의 이산 분포를 통합함으로써 더 정확한 임계값과 함께 검정력 향상을 달성한다.
• 시뮬레이션 결과에서 표준 벤자민-호크베르그 절차 및 기타 이산 분포를 고려한 절차들에 비해 뚜렷한 검정력 우월성을 보였다.
• π₀-적응의 통합은 실제 신호 수준에 따라 절차가 동적으로 조정되어 탐지 능력을 향상시킨다.
• 실제 데이터 응용 결과에서 새로운 방법들이 더 많은 진정한 기각을 식별하면서도 FDR 제어를 유지함을 확인하였다.
• 다양한 이산 분포와 신호 강도에 걸쳐 뛰어난 안정성과 일관된 성능을 보이며 강건함을 입증하였다.