[论文解读] Exponential-type Inequalities Involving Ratios of the Modified Bessel Function of the First Kind and their Applications

Prakash Balachandran, Weston D. Viles+1 more|arXiv (Cornell University)|Nov 6, 2013

Random Matrices and Applications参考文献 22被引用 8

一句话总结

本文建立了关于和 $ H(\nu,x) = \sum_{k=1}^{\infty} \frac{I_{\nu+k}(x)}{I_{\nu}(x)} $ 的紧致上界与下界，利用对比例 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 的指数型界，改进了修正贝塞尔函数的渐近逼近，并推导出斯凯勒姆分布的更紧致的集中不等式。关键贡献在于 $ \nu = 0 $ 附近 $ H(\nu,x) $ 边界中体现的 $ \sqrt{x} $ 缩放行为，从而实现了改进的非渐近误差控制。

ABSTRACT

The modified Bessel function of the first kind, $I_ν(x)$, arises in numerous areas of study, such as physics, signal processing, probability, statistics, etc. As such, there has been much interest in recent years in deducing properties of functionals involving $I_ν(x)$, in particular, of the ratio ${I_{ν+1}(x)}/{I_ν(x)}$, when $ν,x\geq 0$. In this paper we establish sharp upper and lower bounds on $H(ν,x)=\sum_{k=1}^{\infty} {I_{ν+k}(x)}/{I_ν(x)}$ for $ν,x\geq 0$ that appears as the complementary cumulative hazard function for a Skellam$(λ,λ)$ probability distribution in the statistical analysis of networks. Our technique relies on bounding existing estimates of ${I_{ν+1}(x)}/{I_ν(x)}$ from above and below by quantities with nicer algebraic properties, namely exponentials, to better evaluate the sum, while optimizing their rates in the regime when $ν+1\leq x$ in order to maintain their precision. We demonstrate the relevance of our results through applications, providing an improvement for the well-known asymptotic $\exp(-x)I_ν(x)\sim {1}/{\sqrt{2πx}}$ as $x ightarrow \infty$, upper and lower bounding $\mathbb{P}\left[W=ν ight]$ for $W\sim Skellam(λ_1,λ_2)$, and deriving a novel concentration inequality on the $Skellam(λ,λ)$ probability distribution from above and below.

研究动机与目标

• 推导 $ H(\nu,x) = \sum_{k=1}^{\infty} \frac{I_{\nu+k}(x)}{I_{\nu}(x)} $ 的更紧致的非渐近界，其中 $ \nu, x \geq 0 $，这对统计网络分析至关重要。
• 解决现有基于几何级数的 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 界的局限性，这些界无法捕捉 $ \nu = 0 $ 附近的 $ \sqrt{x} $ 缩放行为。
• 改进当 $ x \to \infty $ 时的渐近逼近 $ \exp(-x)I_\nu(x) \sim \frac{1}{\sqrt{2\pi x}} $，尤其针对小的 $ \nu $。
• 推导斯凯勒姆分布 $ \text{Skellam}(\lambda,\lambda) $ 的新非渐近集中不等式，对 $ \mathbb{P}[|W| \geq \nu] $ 提供显式的上下界。
• 为所有 $ \nu \in \mathbb{N} $ 和 $ \lambda_1, \lambda_2 \geq 0 $ 的 $ W \sim \text{Skellam}(\lambda_1, \lambda_2) $ 提供 $ \mathbb{P}[W = \nu] $ 的显式上下界。

提出的方法

• 该方法通过使用优化速率的指数函数，对比例 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 进行上下界控制，尤其在 $ \nu + 1 \leq x $ 时有效。
• 这些指数界随后被用于控制无穷级数 $ H(\nu,x) $，取代在 $ \nu = 0 $ 附近无法保持 $ \sqrt{x} $ 缩放行为的几何级数近似。
• 在 $ \nu + 1 \leq x $ 的区域内，构造了 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 的边界 $ \mathcal{L}(\nu,x) $ 和 $ \mathcal{U}(\nu,x) $，以确保在渐近极限下保持紧致性。
• 通过恒等式 $ \mathbb{P}[|W| > \nu] = 2H(\nu, 2\lambda) \exp(-2\lambda) I_\nu(2\lambda) $ 将结果应用于斯凯勒姆分布，从而推导出集中不等式。
• 该方法利用修正贝塞尔函数 $ I_\nu(x) $ 的已知性质，包括其级数表示和渐近行为，以导出代数上可处理的界。
• 理论界通过数值验证，显示其在 $ x < 10 $ 时显著优于标准渐近界 $ \frac{1}{\sqrt{2\pi x}} $。

实验结果

研究问题

• RQ1 能否推导出对 $ H(\nu,x) = \sum_{k=1}^{\infty} \frac{I_{\nu+k}(x)}{I_{\nu}(x)} $ 的更紧致界，使其在 $ \nu \approx 0 $ 附近保留 $ \sqrt{x} $ 缩放行为，而现有基于几何级数的方法无法做到？
• RQ2 如何优化 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 的指数型界，以提高 $ \exp(-x)I_\nu(x) $ 的渐近逼近精度？
• RQ3 能否为 $ \text{Skellam}(\lambda,\lambda) $ 分布推导出具有显式上下界控制的非渐近集中不等式，以描述 $ \mathbb{P}[|W| \geq \nu] $？
• RQ4 对于所有 $ \nu \in \mathbb{N} $ 和 $ \lambda_1, \lambda_2 \geq 0 $，$ W \sim \text{Skellam}(\lambda_1, \lambda_2) $ 的 $ \mathbb{P}[W = \nu] $ 的精确上下界是什么？
• RQ5 序列 $ \{ \exp(-x)I_\nu(x) \}_{\nu=-\infty}^{\infty} $ 的归一化条件是否可推广至广义格点 $ \mathbb{N} + \alpha $？若可推广，归一化常数是多少？

主要发现

• 本文建立了 $ \frac{I_{\nu+1}(x)}{I_\nu(x)} $ 的紧致上下界 $ \mathcal{L}(\nu,x) $ 和 $ \mathcal{U}(\nu,x) $，在 $ \nu + 1 \leq x $ 时采用优化的指数速率，确保在渐近区域的精度。
• 这些界导出的 $ H(\nu,x) $ 估计在 $ \nu \approx 0 $ 时随 $ x $ 呈 $ \sqrt{x} $ 缩放，纠正了先前基于几何级数方法在 $ \sqrt{x} $ 缩放上的不足。
• 改进后的界导致对 $ \exp(-x)I_\nu(x) $ 的逼近更紧致，优于标准渐近界 $ \frac{1}{\sqrt{2\pi x}} $，尤其在 $ x < 10 $ 时表现更优。
• 推导出斯凯勒姆分布 $ \text{Skellam}(\lambda,\lambda) $ 的新集中不等式，对所有 $ \nu \in \mathbb{N} $ 和 $ \lambda \geq 0 $，提供了 $ \mathbb{P}[|W| \geq \nu] $ 的显式上下界。
• 基于所推导的 $ H(\nu,x) $ 界及斯凯勒姆分布的质量归一化性质，提供了 $ W \sim \text{Skellam}(\lambda_1, \lambda_2) $ 的 $ \mathbb{P}[W = \nu] $ 的显式上下界。
• 数值验证表明，新界在小 $ x $ 时对 $ \exp(-x)I_\nu(x) $ 的逼近显著优于渐近界 $ \frac{1}{\sqrt{2\pi x}} $，同时保持渐近一致性。