[論文レビュー] Inter-arrival time distribution for the non-homogeneous Poisson process

G. Yakovlev, John B. Rundle+2 more|arXiv (Cornell University)|Jul 27, 2005

Advanced Frequency and Time Standards被引用数 11

ひとこと要約

本稿では、有限または無限時間区間における非定常ポアソン過程（NHPP）の到着間隔分布について、正確な解析的表現を導出する。平均値関数 Λ(t) を用い、到着時刻密度を統合することで、観測窓効果を考慮した閉形式の確率密度関数（pdf）が得られ、プロセスの有限時間記憶性に起因する指数分布からの逸脱が示される。

ABSTRACT

We derive an analytical expression of the inter-arrival time distribution for a non-homogeneous Poisson process (NHPP). This expression is exact and is applicable to any time interval, finite or infinite. As an illustration, we present simulation results for three different intensity functions.

研究の動機と目的

• 有限または無限時間区間における非定常ポアソン過程（NHPP）の到着間隔分布に対する正確な解析的表現を導出すること。
• 特に有限観測窓下における到着間隔分布の正確な解が不足しているという問題を解決すること。
• 有限観測区間が、定常ポアソン過程で見られる指数分布から逸脱する原因となる記憶効果を定量化すること。
• 定数、減衰、増加の3種類の強度関数に対して、数値シミュレーションを用いて導出した式の妥当性を検証すること。

提案手法

• 不交差区間における増分の独立性を用い、イベント数 n および到着時刻 t_k のすべての可能性について和をとることで、生存関数 G(x) = P(S > x) を導出する。
• 数値実験における到着間隔のシミュレーションに逆変換法を適用し、NHPPの強度関数 λ(t) と整合性を保つ。
• 時間から累積イベントレートへのマッピングを可能にするために、平均値関数 Λ(t) = ∫₀ᵗ λ(u)du を適用し、時間区間内でのイベント確率の正確な計算を可能にする。
• 生存関数 G(x) を微分することで到着間隔の確率密度関数（pdf）を導出し、λ(t) と Λ(t) を含む二重積分表現が得られる。
• 有限区間の境界効果を考慮するため、最初の到着間隔に対して別個の項を追加し、分布の正規化および完全性を保証する。
• 10⁸回の実現について、λ(t) = 1/c、c/(1+t)、および ct の3つの強度関数に対して、理論的結果とシミュレーション結果を比較することで、解析的解の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 有限時間区間 T で観測された非定常ポアソン過程の正確な到着間隔分布は何か？
• RQ2 有限観測窓 T が、無限時間の場合と比較して到着間隔分布にどのように影響を与えるか？
• RQ3 有限区間が引き起こす記憶効果は、定常ポアソン過程で見られる指数分布からの逸脱をどの程度引き起こすか？
• RQ4 時間依存の強度関数と観測窓の境界制約の両方を考慮した解析的表現を導出できるか？

主な発見

• 有限区間 T における正確な到着間隔確率密度関数（pdf）は、pdf(x) = (1/Λ(T)) ∫₀ᵀ⁻ˣ λ(y)λ(y+x)exp{−[Λ(y+x)−Λ(y)]}dy + [λ(x)exp{−Λ(x)}]/Λ(T) で与えられ、境界効果を考慮している。
• 定常ケース λ(t) = 1/c では、導出された pdf は pdf(x) = (T + c − x)/(cT) exp(−x/c) に簡略化され、有限 T に起因する指数分布からの逸脱が示される。
• λ(t) = 1/c、c/(1+t)、および ct のすべての強度関数について、シミュレーションが理論的 pdf と非常に良好に一致しており、理論と数値結果の一致が確認された。
• 有限観測窓は、時刻 y 以降に T−y より長い到着間隔は観測できないため、プロセスに記憶を導入する。これは、定数 λ(t) に対しても指数分布から逸脱する分布を生じさせる。
• 導出した式は適切に正規化されており、G(0) = 1 かつ G(T) = 0 であるため、生存関数の定義と整合的であることが確認された。
• 本手法により、定常ポアソン過程の指数的到着間隔分布が、任意の決定的強度関数を持つ非定常ケースへと一般化された。