形式陈述
相邻时刻的极端观测若来自同一个冲击,还能当成多个独立机会吗?设 ( X t ) t ∈ Z 为严格平稳的随机过程 公理库 随机过程 Stochastic process · Random process 由同一随机实验产生、按时间或空间指标组织的一族随机变量。 :每个有限向量的联合分布在共同时间平移后不变。记共同分布函数 公理库 分布函数 Cumulative distribution function · CDF 实值分布在每个阈值左侧累积的概率,是非降右连续且端点极限为零和一的函数。 为 F ,M n = max 1 ≤ t ≤ n X t 。
若存在 θ ∈ [ 0 , 1 ] ,使对每个 τ > 0 及满足
n P ( X 0 > u n ) → τ 的门槛序列,都有
P ( M n ≤ u n ) ⟶ e − θ τ , 则称 θ 为该序列的极值指数。此处默认所需门槛序列存在;连续边缘分布通常可直接用高分位选取。
相同边缘的独立样本有极限 e − τ 。因此 θ = 1 对应独立情形的最大值强度,θ < 1 表示同样多的边缘超越被聚集成较少的有效极端机会。若边缘属于某个极值类型 公理库 Fisher–Tippett–Gnedenko 极值定理 Fisher-Tippett-Gnedenko theorem · Extreme types theorem 独立同分布最大值若存在线性归一化后的非退化极限,该极限必属于广义极值分布族。 G ,并且上述极限成立、θ > 0 ,则依赖序列按同一归一化的非退化最大值极限为 G θ 。当 θ = 0 时,定义只说明这些门槛下无超越概率趋于一;不能把形式上的 G 0 当成一条非退化实值极限分布。
严格平稳不同于仅有常均值、常方差与滞后协方差的弱平稳。二阶信息不足以定义本页的联合尾部极限。
直觉
想象一次大冲击连续留下两次很高的读数。单个时刻的超越概率没有告诉你,这两次读数其实只提供一次新的机会。最大值只关心有没有至少一簇出现,因此簇的到达强度决定“完全没有超越”的概率。
不过“极值指数等于平均簇大小的倒数”需要额外的簇极限与均值控制。定义本身只涉及最大值概率,并未自动给出一个唯一的有限样本分簇办法。
例子与边界
移动最大值精确给出二分之一
取 α > 0 ,令 Y t 独立同分布,P ( Y t ≤ y ) = exp ( − y − α ) 、y > 0 ,并定义
X t = max ( Y t , Y t − 1 ) . 它严格平稳,且相隔至少两个时间单位的变量独立 公理库 独立性 Statistical independence 从概率表理解独立性,区分两两、相互和条件独立,并用可计算反例澄清零协方差与条件均值的限度。 。边缘分布为 F ( x ) = exp ( − 2 x − α ) 。取
u n = ( 2 n / τ ) 1 / α , 便有 n P ( X 0 > u n ) = n ( 1 − e − τ / n ) → τ 。
但 n 个 X 的最大值实际只涉及 Y 0 , … , Y n ,所以
P ( M n ≤ u n ) = exp [ − ( n + 1 ) u n − α ] ⟶ e − τ / 2 . 故 θ = 1 / 2 。一个特别大的 Y t 同时使 X t 与 X t + 1 很大;相邻两个大冲击同时出现的概率在稀有事件极限下可以忽略,典型簇于是含两个超越。
图片加载失败 完全依赖可以使指数为零
令所有时刻都等于同一个随机变量:X t = Y 。它仍严格平稳。若 n P ( Y > u n ) → τ ,则
P ( M n ≤ u n ) = P ( Y ≤ u n ) → 1 = e − 0 τ . 所以 θ = 0 。增加观察时刻并未增加独立机会。这个例子没有短程分簇结构,不能把 1 / θ 草率解释为某个可直接估计的普通平均簇大小。
推论与应用
簇解释需要哪些额外控制
一种常见证明策略把 n 个时刻分成长块,中间留较短间隔。需要远隔长块的极端事件近似独立,且删掉的间隔总超越概率可忽略;还需要反聚簇控制,使一个极端冲击产生的簇在所选长块内结束,而不是无限延伸。
在这些条件下,若簇到达数趋向 Poisson ( θ τ ) ,每簇超越个数趋向正整数变量 C ,并有足够的一致可积性把期望传到极限,那么
τ = θ τ E C , θ = 1 / E C . 一致可积性不能省去:计数分布收敛不自动保证均值收敛,少量异常巨簇可能携带不消失的期望。极值指数存在也不自动替代全部混合、反聚簇与矩控制。
原始超越数可能不是 Poisson
对上面的双重移动最大值,由独立稀有 Bernoulli 和的Poisson 极限 公理库 Poisson 分布 Poisson distribution · 泊松分布 以强度参数描述固定窗口内稀有事件计数的离散分布。 ,底层超阈值创新数
K n = ∑ t = 1 n − 1 1 { Y t > u n } ⇒ Poisson ( τ / 2 ) . 两端创新的影响以概率趋零;两个相邻创新同时超越的概率由 n P ( Y 0 > u n ) 2 → 0 控制。因此原始超越数满足
N n := ∑ t = 1 n 1 { X t > u n } = 2 K n + o P ( 1 ) . 其极限是两倍 Poisson 变量,而不是均值 τ 的普通 Poisson 变量。极限均值为 τ 、方差为 2 τ ,且无超越概率为 e − τ / 2 。若忽略聚簇,只看边缘发生率,会把不确定性和最大值概率都算错。
指数不是相关系数,也不是独立性检验
θ 描述越来越高门槛下的极端行为。它不等于任何普通自相关系数;某些重尾甚至没有有限二阶矩。θ = 1 也只说明最大值极限没有因尾部聚簇而降低强度,不能反推出整个过程独立。
独立超越的 Poisson 点过程极限 公理库 极端超越的 Poisson 点过程极限 Poisson point-process limit for extremes · Extreme exceedance point process 正则变化独立样本的极端观测保留发生时间与归一化高度后,形成具有幂型高度强度的 Poisson 点过程极限。 提供没有聚簇时的基准。对依赖数据,先检查是否应当建模簇的到达,再处理每簇内部的多个读数,而不是仅把全部超越当成互相独立的点。
参考资料