Following the pivotal work of Sevastyanov, who considered branching processes with homogeneous Poisson immigration, much has been done to understand the behaviour of such processes under different types of branching and immigration mechanisms. Recently, the case where the times of immigration are generated by a non-homogeneous Poisson process was considered in depth. In this work, we try to demonstrate how one can use the framework of point processes in order to go beyond the Poisson process. As an illustration, we show how to transfer techniques from the case of Poisson immigration to the case where it is spanned by a determinantal point process.
論文ID : 2411.12474タイトル : Multi-type branching processes with immigration generated by point processes著者 : Martin Minchev, Maroussia Slavtchova-Bojkova分類 : math.PR(確率論)発表時期 : 2024年11月(arXiv プレプリント)論文リンク : https://arxiv.org/abs/2411.12474 本論文は、Sevastyanovによる同次ポアソン移民分岐過程に関する先駆的研究に基づき、ポアソン過程を超えた一般的なポイントプロセスの枠組みが分岐過程の移民モデリングに適用される方法を検討している。著者は、ポアソン移民の場合の技術的方法を行列式ポイントプロセス(Determinantal Point Process, DPP)で駆動される移民の場合に一般化する方法を示し、分岐過程理論に新しい理論的枠組みと分析ツールを提供している。
中核的問題 : 従来の分岐過程移民モデルは主にポアソン過程に基づいているが、ポアソン過程は等分散性仮定と時間斉次性の制限を持ち、実際の応用では限界がある重要性 : 分岐過程は生物学、疫学、金融学など広範な分野に応用されており、より柔軟な移民メカニズムは現実現象をより良く記述できる既存の限界 :
ポアソン過程の等分散性仮定(分散が平均に等しい)は実際にはしばしば成立しない 時間斉次性の制限はモデルの適用可能性を制限する 一般的なポイントプロセス移民を扱うための統一的な理論的枠組みが欠けている 一般的なポイントプロセスに基づく分岐過程移民の理論的枠組みを確立する ポアソン過程の制限を超える数学的ツールを提供する 実際の応用に対してより柔軟なモデリング選択肢を提供する 理論的枠組み : ポイントプロセスのLaplace汎関数に基づく統一的な理論的枠組みを確立し、一般的なポイントプロセスで駆動される分岐過程移民を分析する生成関数の導出 : DPP、Cox過程、および分数ポアソン過程移民の場合の確率生成関数の明示的表現を提供する漸近的挙動の分析 : DPP移民多型分岐過程が亜臨界、臨界、および超臨界の場合における極限定理を証明するモーメント計算 : DPP移民の場合の過程の平均と共分散の正確な表現を導出する誤り修正 : 文献17, Theorem 8 および36, Theorem 3 の技術的誤りを修正する一般的なポイントプロセスΦ \Phi Φ で駆動される移民時刻を持つ多型連続時間Markov分岐過程Z ( t ) = ( Z 1 ( t ) , … , Z d ( t ) ) Z(t) = (Z_1(t), \ldots, Z_d(t)) Z ( t ) = ( Z 1 ( t ) , … , Z d ( t )) の漸近的性質を研究する。ここで:
Z i ( t ) Z_i(t) Z i ( t ) は時刻t t t における第i i i 型粒子の数を表す移民はポイントプロセスΦ \Phi Φ のジャンプ時刻T 1 < T 2 < ⋯ T_1 < T_2 < \cdots T 1 < T 2 < ⋯ で発生する 各移民は独立同分布の粒子群I ( k ) ∼ I I^{(k)} \sim I I ( k ) ∼ I をもたらす ポイントプロセスΦ = ∑ i ≤ κ δ X i \Phi = \sum_{i \leq \kappa} \delta_{X_i} Φ = ∑ i ≤ κ δ X i を定義し、その主要な特性は以下を含む:
強度測度 : Λ ( B ) = E [ Φ ( B ) ] \Lambda(B) = E[\Phi(B)] Λ ( B ) = E [ Φ ( B )] Laplace汎関数 : L Φ ( f ) = E [ e − ∫ f d Φ ] = E [ e − ∑ i ≤ κ f ( X i ) ] L_\Phi(f) = E[e^{-\int f d\Phi}] = E[e^{-\sum_{i \leq \kappa} f(X_i)}] L Φ ( f ) = E [ e − ∫ fd Φ ] = E [ e − ∑ i ≤ κ f ( X i ) ] 結合強度 : ρ k ( x 1 , … , x k ) \rho_k(x_1, \ldots, x_k) ρ k ( x 1 , … , x k ) 移民を伴う分岐過程の生成関数は以下のように表現できる:
G Z ( t ) ( s ) = L Φ ( f t ) , f t ( x ) = − ln ( G Z × ( t − x ) ( s ) ) 1 { x ≤ t } G_{Z(t)}(s) = L_\Phi(f_t), \quad f_t(x) = -\ln(G_{Z_\times(t-x)}(s))1_{\{x \leq t\}} G Z ( t ) ( s ) = L Φ ( f t ) , f t ( x ) = − ln ( G Z × ( t − x ) ( s )) 1 { x ≤ t }
ここでZ × Z_\times Z × は移民のない基礎分岐過程である。
( Λ , K ) (\Lambda, K) ( Λ , K ) -DPPの場合、結合強度は:
ρ n ( x 1 , … , x n ) = det ( K ( x i , x j ) ) 1 ≤ i , j ≤ n \rho_n(x_1, \ldots, x_n) = \det(K(x_i, x_j))_{1 \leq i,j \leq n} ρ n ( x 1 , … , x n ) = det ( K ( x i , x j ) ) 1 ≤ i , j ≤ n
Laplace汎関数は:
L Φ ( f ) = 1 + ∑ n ≥ 1 ( − 1 ) n n ! ∫ R + n D ( x ) ∏ i = 1 n ( 1 − e − f ( x i ) ) Λ ( d x ) L_\Phi(f) = 1 + \sum_{n \geq 1} \frac{(-1)^n}{n!} \int_{\mathbb{R}_+^n} D(x) \prod_{i=1}^n (1-e^{-f(x_i)}) \Lambda(dx) L Φ ( f ) = 1 + ∑ n ≥ 1 n ! ( − 1 ) n ∫ R + n D ( x ) ∏ i = 1 n ( 1 − e − f ( x i ) ) Λ ( d x )
導向測度がη \eta η のCox過程の場合:
L Φ ( f ) = E [ exp ( − ∫ R + ( 1 − e − f ( x ) ) η ( d x ) ) ] L_\Phi(f) = E\left[\exp\left(-\int_{\mathbb{R}_+} (1-e^{-f(x)}) \eta(dx)\right)\right] L Φ ( f ) = E [ exp ( − ∫ R + ( 1 − e − f ( x ) ) η ( d x ) ) ]
パラメータ( β , λ ) (\beta, \lambda) ( β , λ ) の分数ポアソン過程の場合:
L Φ ( f ) = E [ exp ( − λ ∫ R + ( 1 − e − f ( t ) ) d Y β ( t ) ) ] L_\Phi(f) = E\left[\exp\left(-\lambda \int_{\mathbb{R}_+} (1-e^{-f(t)}) dY_\beta(t)\right)\right] L Φ ( f ) = E [ exp ( − λ ∫ R + ( 1 − e − f ( t ) ) d Y β ( t ) ) ]
ここでY β Y_\beta Y β はβ \beta β -安定従属過程の逆過程である。
( Λ , K ) (\Lambda, K) ( Λ , K ) -DPP移民の場合:
平均 :
E [ Z i ( t ) ] = ∫ ( 0 , t ] K ( x , x ) E [ Z × , i ( t − x ) ] Λ ( d x ) E[Z_i(t)] = \int_{(0,t]} K(x,x) E[Z_{\times,i}(t-x)] \Lambda(dx) E [ Z i ( t )] = ∫ ( 0 , t ] K ( x , x ) E [ Z × , i ( t − x )] Λ ( d x ) 共分散 :
Cov ( Z i ( t ) , Z j ( t ) ) = ∫ ( 0 , t ] K ( x , x ) E [ Z × , i ( t − x ) Z × , j ( t − x ) ] Λ ( d x ) \text{Cov}(Z_i(t), Z_j(t)) = \int_{(0,t]} K(x,x) E[Z_{\times,i}(t-x)Z_{\times,j}(t-x)] \Lambda(dx) Cov ( Z i ( t ) , Z j ( t )) = ∫ ( 0 , t ] K ( x , x ) E [ Z × , i ( t − x ) Z × , j ( t − x )] Λ ( d x ) − ∫ ( 0 , t ] 2 K 2 ( x , y ) E [ Z × , i ( t − x ) ] E [ Z × , j ( t − y ) ] Λ ( d x ) Λ ( d y ) - \int_{(0,t]^2} K^2(x,y) E[Z_{\times,i}(t-x)]E[Z_{\times,j}(t-y)] \Lambda(dx)\Lambda(dy) − ∫ ( 0 , t ] 2 K 2 ( x , y ) E [ Z × , i ( t − x )] E [ Z × , j ( t − y )] Λ ( d x ) Λ ( d y ) ρ \rho ρ を生成行列A A A のPerron-Frobenius根とし、∫ ( 0 , ∞ ) e − ρ x K ( x , x ) Λ ( d x ) < ∞ \int_{(0,\infty)} e^{-\rho x} K(x,x) \Lambda(dx) < \infty ∫ ( 0 , ∞ ) e − ρ x K ( x , x ) Λ ( d x ) < ∞ ならば、R d \mathbb{R}^d R d 値確率変数W W W が存在して:
Z ( t ) e ρ t → d W ( t → ∞ ) \frac{Z(t)}{e^{\rho t}} \xrightarrow{d} W \quad (t \to \infty) e ρt Z ( t ) d W ( t → ∞ )
定常DPPでλ ( x ) ∼ λ ∞ e δ x \lambda(x) \sim \lambda_\infty e^{\delta x} λ ( x ) ∼ λ ∞ e δ x の場合:
超臨界の場合 (δ > max { ρ , 0 } \delta > \max\{\rho, 0\} δ > max { ρ , 0 } ):
Z i ( t ) e δ t → L 2 A i = K ∗ λ ∞ ∫ 0 ∞ e − δ x E [ Z × , i ( x ) ] d x \frac{Z_i(t)}{e^{\delta t}} \xrightarrow{L^2} A_i = K^* \lambda_\infty \int_0^\infty e^{-\delta x} E[Z_{\times,i}(x)]dx e δ t Z i ( t ) L 2 A i = K ∗ λ ∞ ∫ 0 ∞ e − δ x E [ Z × , i ( x )] d x 臨界成長 (ρ > 0 , δ = ρ \rho > 0, \delta = \rho ρ > 0 , δ = ρ ):
Z i ( t ) t e δ t → L 2 A i ′ = K ∗ λ ∞ ⟨ u , E [ I ] ⟩ v i \frac{Z_i(t)}{te^{\delta t}} \xrightarrow{L^2} A_i' = K^* \lambda_\infty \langle u, E[I] \rangle v_i t e δ t Z i ( t ) L 2 A i ′ = K ∗ λ ∞ ⟨ u , E [ I ]⟩ v i 臨界の場合 (ρ = 0 , δ = 0 \rho = 0, \delta = 0 ρ = 0 , δ = 0 ):
Z ( t ) t → d Y v \frac{Z(t)}{t} \xrightarrow{d} Y v t Z ( t ) d Y v
ここでY ∼ Γ ( K ( 0 , 0 ) λ ∞ β , 1 / Q ) Y \sim \Gamma(K(0,0)\lambda_\infty \beta, 1/Q) Y ∼ Γ ( K ( 0 , 0 ) λ ∞ β , 1/ Q ) 亜臨界の場合 (ρ < 0 , δ = 0 \rho < 0, \delta = 0 ρ < 0 , δ = 0 ):
Z ( t ) → d X Z(t) \xrightarrow{d} X Z ( t ) d X
ここでX X X のLaplace変換は式(6.17)で与えられる。Laplace汎関数を通じて一般的なポイントプロセス移民を扱うための統一的方法を確立する 従来の方法で区間内の点数分布を具体的に計算する必要性を回避する 支配収束定理とマルチンゲール理論を巧妙に組み合わせ、DPPの複雑な相関構造を処理する Hadamard不等式を利用して行列式の増長を制御し、積分の可制御性を確保する 文献における固有ベクトル正規化に関する技術的誤りを特定し修正する 異なる文献における固有ベクトル定義の相違を明確にする 本論文は以下の重要な研究に基づいている:
Sevastyanov (1957) : 同次ポアソン移民分岐過程の基礎的研究Mitov et al. (2018) : 非同次ポアソン移民の深入的研究Slavtchova-Bojkova et al. (2022-2023) : 多型の場合への拡張Rabehasaina & Woo (2021) : 特性関数法の応用既存研究と比較して、本論文の主な利点は:
ポアソン過程を超える一般的な理論的枠組みを提供する より簡潔で統一的な導出方法を提供する 文献の技術的誤りを修正する 一般的なポイントプロセス移民に基づく分岐過程の理論的枠組みを確立した DPP移民の場合における完全な漸近理論を証明した 実用的な計算ツールと方法を提供した 技術的条件 : いくつかの結果は較強いモーメント条件と可積性仮定を必要とする計算の複雑性 : DPPの行列式構造は高次元の場合に計算が複雑である実際の応用 : 具体的な応用事例と数値実験が欠けているより一般的な可分解(reducible)の場合を研究する 数値計算方法とアルゴリズムを開発する 生物学、疫学における具体的な応用を探索する 理論的深さ : 厳密で完全な数学理論を提供し、この分野の重要な空白を埋める方法的革新 : Laplace汎関数法の応用は従来の複雑な導出過程を簡素化する技術的厳密性 : 証明過程は詳細で、技術処理は適切である実用的価値 : 実際のモデリングに対してより柔軟なツールを提供する応用指向 : 理論導出に偏り、実際の応用の展示が欠けている数値的側面 : 計算方法と数値実験が提供されていない可読性 : 非専門家にとって技術的敷居が高い学術的貢献 : 分岐過程理論に重要な理論的拡張を提供する方法論的価値 : Laplace汎関数法は他の関連研究にインスピレーションを与える可能性がある応用の可能性 : 実際の問題のモデリングに新しい理論的基礎を提供する移民事象間の相関性を考慮する必要がある人口動力学モデリング 排斥性または集約性の特性を持つ疫病伝播モデル 金融市場のジャンプ拡散過程モデリング 論文は50篇の重要な文献を引用しており、主に以下を含む:
Athreya & Ney (1972): 分岐過程の古典的教科書 Sevastyanov (1957, 1971): 移民分岐過程の基礎的研究 Baccelli et al. (2024): ポイントプロセス理論の参考文献 Mitov et al. (2018): 最近の関連重要研究 総合評価 : これは分岐過程理論において重要な貢献をなした高品質な理論数学論文である。理論に偏り応用展示が欠けているが、確立された理論的枠組みは重要な学術的価値と応用の可能性を持っている。