[논문 리뷰] Intertwining of birth-and-death processes

Jan M. Swart|arXiv (Cornell University)|2010. 04. 30.

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

한 줄 요약

이 논문은 디아코니스와 미클로의 생애·사멸 과정에 대한 확률적 커플링을 확장하여, 원래 과정을 아래와 위에서 확률적으로 둘러싸는 제3의 과정을 도입함으로써, 세 과정이 모두 수준 N에 동시에 도달하도록 보장한다. 주요 기여는 생성자 고유값을 이용한 대칭적 통합 구조를 제안하는 것으로, 중간 과정이 고유값 순서로 정렬된 생애율을 가진 두 개의 순수 생애 과정 사이에 확률적으로 끼워지게 되며, 이는 첫 번째 도착 시간이 서로 독립적인 지수분포의 합으로 분포함을 증명한다.

ABSTRACT

It has been known for a long time that for birth-and-death processes started in zero the first passage time of a given level is distributed as a sum of independent exponentially distributed random variables, the parameters of which are the negatives of the eigenvalues of the stopped process. Recently, Diaconis and Miclo have given a probabilistic proof of this fact by constructing a coupling between a general birth-and-death process and a process whose birth rates are the negatives of the eigenvalues, ordered from high to low, and whose death rates are zero, in such a way that the latter process is always ahead of the former, and both arrive at the same time at the given level. In this note, we extend their methods by constructing a third process, whose birth rates are the negatives of the eigenvalues ordered from low to high and whose death rates are zero, which always lags behind the original process and also arrives at the same time.

연구 동기 및 목표

• 원래의 생애·사멸 과정에 대해 디아코니스와 미클로의 커플링 방법을 확장하여 하한을 제공하는 과정을 포함시키는 것.
• 원래 과정을 고유값 순서로 정렬된 생애율을 가진 두 개의 순수 생애 과정 사이에 대칭적인 확률적 샌드위치로 끼우는 것.
• 첫 번째 도착 시간이 독립적인 지수분포의 합으로 분포함을 증명하는 것. 이때 대수적 라플라스 변환 방법이 아닌 확률적 커플링 기법을 사용한다.
• 마르코프 과정의 통합 관계를 통해 고전적인 첫 번째 도착 시간 분포 결과에 대한 구성적 확률적 증명을 제공하는 것.

제안 방법

• 정지된 생성자의 비영 고유값의 음수인 $\lambda_i$를 기반으로, 생애율 $b^+_i = \lambda_{N+1-i}$를 가지는 빠른 순수 생애 과정 $X^+$를 구성한다.
• 생애율을 증가순서로 정렬한 $b^-_i = \lambda_i$를 가지는 느린 순수 생애 과정 $X^-$를 구성한다.
• 확률 커널 $K^-$와 $K^+$를 정의하여, $K^-(x, \{0,\ldots,x\}) = 1$, $K^+(x, \{0,\ldots,x\}) = 1$, 그리고 $K^-(N,N) = K^+(N,N) = 1$를 만족시키며, $N$에서의 흠직을 보장한다.
• 통합 관계 수립: $GK^+ = K^+G^+$ 및 $GK^- = K^-G^-$, 여기서 $G$, $G^+$, $G^-$는 각각 $X$, $X^+$, $X^-$의 생성자이다.
• 통합 성질을 활용하여 $X$, $X^+$, $X^-$를 커플링함으로써, 모든 $t \geq 0$에 대해 $X^-_t \leq X_{t \wedge \tau_N} \leq X^+_t$를 만족시키며, 세 과정이 거의 확실히 동시에 $N$에 도달하도록 보장한다.
• 페르로-프로베니우스 정리와 생성자의 스펙트럼 성질을 활용하여, 구성된 과정에서 생애율이 양수이자 엄격히 증가함을 보장한다.

실험 결과

연구 질문

• RQ1 원래의 생애·사멸 과정이 고유값 순서로 정렬된 생애율을 가진 두 개의 순수 생애 과정 사이에 확률적으로 끼워지게 하여, 둘 다 수준 $N$에 동시에 도달하도록 하는 확률적 커플링을 구성할 수 있는가?
• RQ2 라플라스 변환 방법이 아닌 대칭적 커플링 프레임워크를 통해 생애·사멸 과정의 첫 번째 도착 시간 분포를 유도할 수 있는가?
• RQ3 정지된 생성자의 고유값은 첫 번째 도착 시간 분포에 어떤 역할을 하는가? 그리고 이 고유값들은 보조 순수 생애 과정의 생애율로 어떻게 실현될 수 있는가?
• RQ4 통합 관계를 상하 둘 다를 끼치는 순수 생애 과정을 포함하도록 확장할 수 있는가? 이때 동일한 흡수 시간을 유지할 수 있는가?
• RQ5 커널 $K^-$와 $K^+$는 생성자의 스펙트럼 분해와 어떻게 관련되어 있으며, 그들의 확률적 양성과 $N$에서의 흡수와의 일관성을 보장하는 요소는 무엇인가?

주요 결과

• 수준 $N$으로의 첫 번째 도착 시간 $\tau_N$은 서로 독립적인 지수분포의 합으로 분포하며, 그 비율은 $\lambda_1 < \cdots < \lambda_N$이다. 여기서 $\lambda_i$는 정지된 생성자의 비영 고유값의 음수이다.
• 모든 $t \geq 0$에 대해 $X^-_t \leq X_{t \wedge \tau_N} \leq X^+_t$를 만족시키며, 세 과정이 거의 확실히 동시에 $N$에 도달하는 커플링이 존재한다.
• 원래 생애·사멸 과정의 생성자 $G$는 $N+1$개의 서로 다른 고유값 $0 = -\lambda_0 > -\lambda_1 > \cdots > -\lambda_N$을 가지며, 고유값 $\lambda_i$는 보조 순수 생애 과정에서 생애율로 실현된다.
• 커널 $K^+$는 원래 과정 $X$를 생애율 $\lambda_N, \ldots, \lambda_1$을 가진 빠른 순수 생애 과정 $X^+$로 매핑하며, 통합 관계 $GK^+ = K^+G^+$를 만족한다.
• 커널 $K^-$는 원래 과정 $X$를 생애율 $\lambda_1, \ldots, \lambda_N$을 가진 느린 순수 생애 과정 $X^-$로 매핑하며, 통합 관계 $GK^- = K^-G^-$를 만족한다.
• 이 구성은 고유값 $\lambda_i$가 엄격히 증가하고 양수임을 보장하며, $\lambda_i = b^-_i$ 및 $\lambda_{N+1-i} = b^+_i$임을 확인함으로써, 커플링 내에서 스펙트럼 순서가 유지됨을 확인한다.