[논문 리뷰] A test for dependence between two point processes on the real line

Patrick Rubin‐Delanchy, Nicholas A. Heard|arXiv (Cornell University)|2014. 08. 17.

Stochastic processes and statistical mechanics참고 문헌 29인용 수 4

한 줄 요약

이 논문은 실수선 상의 두 1차원 점과정 간의 의존성에 대한 정확한 일반화된 우도비 검정을 제안한다. 이 검정은 A 과정의 사건 이후 B 과정의 강도 변화를 통해 자극 효과를 모델링한다. 검정은 균일한 포아송 과정에서의 변화점 검출과의 이중성에 기반하여 유한표본 정확한 p-값을 도출하며, 점점이 근사적으로 가중 칼모고로프-스미르노프 검정과 동치가 된다. 구현은 R 패키지 'mppa'에 제공된다.

ABSTRACT

Many scientific questions rely on determining whether two sequences of event times are associated. This article introduces a likelihood ratio test which can be parameterised in several ways to detect different forms of dependence. A common finite-sample distribution is derived, and shown to be asymptotically related to a weighted Kolmogorov-Smirnov test. Analysis leading to these results also motivates a more general tool for diagnosing dependence. The methodology is demonstrated on data generated on an email network, showing evidence of information flow using only timing information. Implementation code is available in the R package `mppa'.

연구 동기 및 목표

• 실수선 상의 두 점과정 간의 자극 또는 상관관계를 탐지하기 위한 통계적으로 엄밀하고 정확한 검정을 개발하는 것, 특히 신경 신호 전파나 네트워크 통신 데이터와 같은 맥락에서 적용 가능한 것.
• 비균일한 강도 조건 하에서도 유효하며, p-값 산출에 점점이론 근사에 의존하지 않는 방법을 제공하는 것.
• 검정 문제와 균일한 포아송 과정에서의 변화점 검출 간의 연결 고리를 확립하여 수학적 분석을 단순화하는 것.
• 파rameterization 을 통해 자극, 상관관계 또는 억제와 같은 다양한 형태의 의존성을 탐지할 수 있는 융통성 있는 프레임워크를 제공하는 것.
• 실제 데이터(예: 이메일 네트워크나 신경 기록)에 적용 가능한 실용적이고 구현 가능한 검정을 제공하고, R 코드를 통해 구현 가능하게 하는 것.

제안 방법

• 검정은 B 과정의 강도를 조각별로 일정한 값으로 모델링하며, A 과정의 각 사건 이후 τ 시간 이내에는 더 높은 강도 λ₁, 이후에는 더 낮은 강도 λ₂를 가진다. 귀무가설 하에서는 λ₁ = λ₂이다.
• 대립가설 하에서 λ₁ > λ₂일 때와 귀무가설 하에서 λ₁ = λ₂일 때의 우도를 비교하기 위해 일반화된 우도비 통계량을 사용하며, 검정 통계량은 A 과정의 사건으로부터 τ 이내에 발생한 B 과정의 사건 수로부터 유도된다.
• 문제를 균일한 포아송 과정에서의 변화점 검출 문제로 변환하는 이중성(보조정리 1)을 활용하여, 정확한 유한표본 분포 결과를 도출할 수 있다.
• 검정 통계량은 점점이 근사적으로 가중 칼모고로프-스미르노프 검정과 동치임을 보이며, 잘 알려진 비모수적 적합도 검정과의 연결 고리를 제공한다.
• 최적의 τ와 우도비를 추정하는 알고리즘을 통해 절차를 구현하며, 대립가설 하에서 λ₁과 λ₂에 대한 일致한 추정량을 제공한다.
• B 과정의 기저 강도 r(t)가 알려지지 않은 경우에도 검정이 유효하며, r(t)가 유계이면서 적분 가능하다면 검정 결과는 그대로 유지된다.

실험 결과

연구 질문

• RQ1 유한 표본에서도 실수선 상의 두 점과정 간의 자극 행동을 정확한 우도비 검정이 탐지할 수 있는가?
• RQ2 한 과정의 기저 강도가 알려지지 않거나 비균일한 경우, 두 점과정 간의 의존성을 어떻게 진단할 수 있는가?
• RQ3 두 점과정 간 의존성에 대한 우도비 검정 통계량의 유한표본 분포는 무엇인가?
• RQ4 제안된 검정은 가중 칼모고로프-스미르노프 검정과 같은 기존 비모수적 검정과 어떻게 관련이 있는가?
• RQ5 이 방법은 단순한 자극 외에도 억제나 상관관계와 같은 다른 형태의 의존성을 탐지할 수 있는가?

주요 결과

• 검정은 유한표본에서 정확한 p-값을 제공하므로, 작고 제한된 데이터셋에서도 신뢰할 수 있으며 점점이론 근사에 의존하지 않는다.
• 우도비 검정 통계량은 점점이 근사적으로 가중 칼모고로프-스미르노프 검정과 동치이며, 잘 알려진 비모수적 절차와의 연결 고리를 제공한다.
• 정규성 조건 하에서 제시된 모델 하에서 자극 효과를 탐지하기 위해 균일하게 가장 강력한 검정이다.
• 기저 강도 함수가 비균일하더라도, 그 함수가 유계이면서 가측 가능하다면 검정은 여전히 유효하다.
• 최적의 τ 및 강도 파rameter의 일致한 추정은 표본 크기가 증가함에 따라 검정의 검정력을 유지함을 보장한다.
• 이 방법은 R 패키지 'mppa'에 구현되어 있어, 이메일 네트워크의 이벤트 시점이나 신경 기록과 같은 실제 데이터에 적용 가능하다.