[논문 리뷰] Multiple Testing of Local Maxima for Detection of Peaks in Random Fields

Dan Cheng, Armin Schwartzman|arXiv (Cornell University)|2014. 05. 06.

Topological and Geometric Data Analysis참고 문헌 15인용 수 6

한 줄 요약

이 논문은 정적 마르코프 가우시안 잡음 하에서 다차원 랜덤 필드의 피크를 탐지하기 위해 국소 최대값을 테스트하는 위상적 다중 검정 프레임워크를 제안한다. Benjamini-Hochberg 절차를 사용하여 점차적 FDR 제어 및 검정력 일致성을 달성하며, 두 가지 p-값 계산 방법을 제공한다: 등방성 필드의 경우 정확한 높이 분포, 비등방성 필드의 경우 오버슈트 근사. 시뮬레이션과 나노현미경 데이터를 통해 검증된다.

ABSTRACT

A topological multiple testing scheme is presented for detecting peaks in images under stationary ergodic Gaussian noise, where tests are performed at local maxima of the smoothed observed signals. The procedure generalizes the one-dimensional scheme of Schwartzman et al. (2011) to Euclidean domains of arbitrary dimension. Two methods are developed according to two different ways of computing p-values: (i) using the exact distribution of the height of local maxima (Cheng and Schwartzman, 2014), available explicitly when the noise field is isotropic; (ii) using an approximation to the overshoot distribution of local maxima above a pre-threshold (Cheng and Schwartzman, 2014), applicable when the exact distribution is unknown, such as when the stationary noise field is non-isotropic. The algorithms, combined with the Benjamini-Hochberg procedure for thresholding p-values, provide asymptotic strong control of the False Discovery Rate (FDR) and power consistency, with specific rates, as the search space and signal strength get large. The optimal smoothing bandwidth and optimal pre-threshold are obtained to achieve maximum power. Simulations show that FDR levels are maintained in non-asymptotic conditions. The methods are illustrated in a nanoscopy image analysis problem of detecting fluorescent molecules against the image background.

연구 동기 및 목표

• 정적 마르코프 가우시안 잡음에 의해 오염된 영상에서 피크를 탐지하기 위한 다중 검정 절차를 개발한다.
• 일차원 피크 탐지 방법을 임의의 다차원 유클리드 영역으로 확장한다.
• 큰 영역과 신호 강도 조건 하에서 False Discovery Rate (FDR)의 점차적 제어 및 검정력 일치를 보장한다.
• 정확한 높이 분포(등방성 필드용)와 오버슈트 근사(비등방성 필드용)를 포함한 두 가지 p-값 계산 전략을 제공한다.
• 탐지 능력을 극대화하기 위해 스무딩 대역폭과 사전 임계값을 최적화한다.

제안 방법

• 피크 탐지 이전에 신호 대 잡음 비율을 향상시키기 위해 커널 스무딩을 적용한다.
• 사전 임계값 이상의 스무딩된 필드 내 국소 최대값을 식별하여 후보 피크를 정의한다.
• 전체 근본 가설(신호 없음) 하에서 각 후보 피크에서 p-값을 계산하며, 정확한 높이 분포 또는 오버슈트 근사를 사용한다.
• Benjamini-Hochberg 절차를 사용하여 p-값을 임계값 처리하고 FDR을 제어한다.
• 임의 행렬 이론의 결과를 활용하여 국소 최대값 높이의 정확한 및 근사 분포를 유도한다.
• 점차적 분석을 통해 스무딩 대역폭과 사전 임계값을 최적화하여 탐지 능력을 극대화한다.

실험 결과

연구 질문

• RQ1 국소 최대값 기반의 다중 검정 프레임워크는 다차원 랜덤 필드에서 점차적 FDR 제어를 달성할 수 있는가?
• RQ2 정확한 높이 분포가 알려지지 않은 상황에서 사전 임계값 처리는 FDR 제어와 검정력에 어떤 영향을 미치는가?
• RQ3 피크 탐지에서 탐지 능력을 극대화하기 위해 최적의 스무딩 대역폭과 사전 임계값은 무엇인가?
• RQ4 정확한 높이 분포와 오버슈트 근사는 등방성 및 비등방성 잡음 필드에서 성능가 비교해 볼 때 어떻게 다른가?
• RQ5 이 방법은 유한 표본, 비점차적 설정에서도 FDR 제어를 유지하는가?

주요 결과

• 제안된 방법은 큰 영역과 신호 강도 조건 하에서 점차적 FDR 제어를 달성한다.
• 신호 강도가 영역 크기의 로그 제곱근보다 더 빠르게 증가할 경우 검정력 일치가 입증된다.
• 등방성 필드의 경우, 국소 최대값의 정확한 높이 분포를 통해 사전 임계값 없이 최적의 p-값 계산이 가능하다.
• 비등방성 필드의 경우 정확한 분포가 알려져 있지 않기 때문에 오버슈트 근사를 사용하고 사전 임계값 처리와 조합함으로써 효과적인 p-값 계산이 가능하다.
• 시뮬레이션 결과는 비점차적, 유한 표본 조건에서도 FDR 수준이 유지됨을 확인한다.
• 이 방법은 나노현미경 데이터에서 형광 분자를 성공적으로 탐지하여 실제 영상 응용 분야에서의 실용성을 입증한다.