Skip to content

定理Theorem

极端超越的 Poisson 点过程极限

Poisson point-process limit for extremes · Extreme exceedance point process

正则变化独立样本的极端观测保留发生时间与归一化高度后,形成具有幂型高度强度的 Poisson 点过程极限。

形式陈述 ​

只知道最大值会丢掉其他大观测的时间与高度。怎样同时保留全部极端点?设 X1,X2,… 为来自同一总体的IID 正变量,尾函数 F―∈RV−α,α>0,并选择 an>0 使 nF―(an)→1。

在状态空间 E=[0,1]×(0,∞] 上定义随机点测度

Nn=∑i=1nδ(i/n,Xi/an).

取高度坐标 1/y 后,E 同胚于 [0,1]×[0,∞),是第二可数局部紧 Hausdorff 空间。这里使用允许重数、紧集上点数有限的Radon 点测度。它与经验概率测度都把观测变成点质量,但每点质量为一,没有除以 n,并且同时记录时间与高度。

由正则变化,Nn 依分布趋向 Poisson 随机测度 N,强度为

μ(dt,dy)=dtαy−α−1dy,μ([0,1]×{∞})=0.

极限随机测度的定义是:对互不相交、强度有限的 Borel 区域 A1,…,Am,N(Aj) 相互独立,且各自服从均值 μ(Aj) 的Poisson 分布。若要从 vague 收敛进一步得到 Nn(Aj) 的联合计数收敛,可取这些区域相对紧且边界强度为零;在这里,相对紧意味着高度与零保持正距离。

空间使用局部有限点测度的 vague 拓扑:通过连续紧支撑测试函数的积分比较测度,测试支撑必须离高度零有正距离。高度端点 ∞ 可用坐标 1/y 紧化。随机点测度的分布则按这个取值空间上的概率律弱收敛理解。

直觉

把 n 个观测压进单位时间,再把高度除以典型最大值尺度 an。普通观测都挤向高度零,只有少量极端点仍留在可见区域。每个时刻贡献可见点的概率约为 1/n,独立的小概率机会累加后便出现 Poisson 规律。

高度越低,可见点越多,因此强度不是有限总质量的概率密度。这个极限专门观察远离零的稀有区域,而不是要求把全部 n 个点同时变成有限个点。

例子与边界

一个矩形的强度算到底 ​

取精确 Pareto 尾 P(X>x)=x−2、x≥1,则 an=n。考虑

A=(0,1/2]×(1,2].

这是前半段时间内,归一化高度超过一但不超过二的点。强度为

μ(A)=12∫122y−3dy=12(1−14)=38.

有限样本中,相关计数是 Binomial(⌊n/2⌋,3/(4n)),恰好趋向 Poisson(3/8)。因此该区域没有极端点的概率趋于 e−3/8≈0.6873。

若另取 B=(1/2,1]×(2,∞],强度为 1/8,两区计数极限独立。即便两个区域共享时间区间,只要高度部分不相交,联合稀有类别的极限计数仍独立;有限样本中同一个观测只能落进一类造成的依赖会在极限消失。

为什么必须避开高度零 ​

对每个 ε>0,

μ([0,1]×[ε,∞])=ε−α<∞.

但让 ε↓0,强度发散。不能把包括高度零的整幅图当成一个有限均值 Poisson 计数区域,也不能将强度 μ 当作总质量为一的概率分布。

如果观测有高值聚簇,独立 Poisson 点极限也可能失效。比如移动最大值让一个大创新产生相邻两个超越;归一化时间后它们可挤到同一极限时刻,形成簇点过程,见极值指数。

推论与应用

单区域的均值先由尾比值得到 ​

对矩形 A=(s,t]×(x,y],其中 0<x<y≤∞,约定 ∞−α=0,正则变化给

nP(X1/an∈(x,y])=n[F―(anx)−F―(any)]⟶x−α−y−α.

时间区间约有 n(t−s) 个观测,因此 ENn(A)→(t−s)(x−α−y−α)=μ(A)。有限个端点取整的影响趋零。

联合计数为何变独立 ​

取若干不交矩形,记 qn,i,j 为第 i 个点落入第 j 个区域的概率。单个点至多进入一个区域,不同 i 则独立。对任意固定实系数向量 θ,把计数的线性组合视为实随机变量,其指数矩可写成

Eexp⁡(∑jθjNn(Aj))=∏i=1n[1+∑jqn,i,j(eθj−1)].

所有区域离零有共同正距离,所以 maxi,jqn,i,j=O(1/n),而 ∑iqn,i,j→μ(Aj)。取对数,用 log⁡(1+z)=z+O(z2);二次余量的总和是 O(1/n)。极限为

exp⁡[∑jμ(Aj)(eθj−1)],

恰好分解成相互独立的 Poisson 指数矩。对任意固定向量 v 令 θ=tv,便在 t=0 的邻域得到线性组合的矩母函数收敛;矩母函数连续性定理及 Cramér–Wold 判据因此给出整个有限计数向量的联合收敛。

要从有限区域计数走到点测度收敛,还需控制整族点测度的概率律。令 Kj=[0,1]×[1/j,∞],每个窗口上的期望计数满足 Cj:=supnENn(Kj)<∞。给定 δ>0,为各窗口选整数 mj,使 Cj/mj<δ2−j;Markov 不等式和可数并集界保证,概率至少 1−δ 时所有窗口都同时满足 Nn(Kj)≤mj。这样的点测度集合在 vague 拓扑下相对紧:每个紧窗口内只有有界个点,可逐窗抽取,并允许相邻点合并为重数。

这个局部有限点测度空间在 vague 拓扑下是 Polish 空间,所以上述紧控制可由Prokhorov 定理转为概率律的弱收敛子列。连续性矩形上的联合计数唯一确定每个子列极限,故整个序列趋向所述 Poisson 随机测度。

最大值与第 k 大值都是空点数问题 ​

最大值不超过 anx,等价于区域 [0,1]×(x,∞] 内没有点。因此

P(Mn/an≤x)→e−x−α,x>0.

这重新给出 Fréchet 极限,但还保留了更多信息。第 k 大观测不超过 anx,等价于该区域至多有 k−1 个点,所以

P(X(n−k+1)/an≤x)⟶e−x−α∑j=0k−1x−αjj!.

同一强度测度因而同时控制最高点、多个极端次序统计量及指定时间窗口中的超越次数。它给出已知尾部归一化下的分布近似;若 α 或 an 需要估计,参数不确定性不会自动包含在 Poisson 随机波动里。

参考资料
  • Richard L. Smith and Ishay Weissman, 《Extreme Values》,§1.2.4 的时间—高度点过程与 Poisson 强度;§3.5 的二维超越过程。
关系图谱27 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系