Skip to content

定理Theorem

样本分位数的渐近线性化

Sample quantile asymptotics · Bahadur representation for a sample quantile

在目标分位附近密度连续且为正时,经验分位的随机逆可线性化,余项由局部经验分布振荡控制。

形式陈述 ​

经验分布是一条阶梯函数,不能像普通光滑函数那样直接求逆导数。样本分位数的正态近似究竟怎样得到?设 X1,…,Xn 是来自 F 的独立同分布样本,固定 p∈(0,1),目标分位 q=F←(p) 唯一,且 F 在 q 的某个邻域有连续密度 f,f(q)>0。

令经验分布函数为 Fn,采用左广义逆约定

q^n=Fn←(p)=X(⌈np⌉),

其中括号秩是次序统计量。则有渐近线性展开

q^n−q=p−Fn(q)f(q)+oP(n−1/2)=1n∑i=1np−1{Xi≤q}f(q)+oP(n−1/2).

因此

n(q^n−q)⇒N(0,p(1−p)f(q)2).

余项结论比“左右两边有相同极限分布”更强:它说明二者在同一份样本上的差,比通常的 n−1/2 误差还小。

直觉

若经验 CDF 在真实分位 q 处偏高,说明落在 q 左边的样本过多,达到目标概率 p 的经验位置就向左移动。于是分位误差带有 p−Fn(q) 的负反馈符号。

CDF 的纵向误差要除以局部斜率 f(q) 才变成横向位置误差。密度小,附近样本稀疏,同样的概率误差会引起更大的位置移动。但要把随机的 q^n 代回经验过程,还必须证明这个小邻域内没有额外的一阶随机振荡。

例子与边界

四分之三分位的指数例子 ​

率一指数分布有 F(x)=1−e−x。取 p=3/4,则 q=log⁡4、f(q)=1/4,所以

n(q^n−log⁡4)⇒N(0,3).

线性化项为 4[3/4−Fn(log⁡4)]。如果 n=400 的样本在真实分位左侧有 310 个观测,则 Fn(q)=0.775,一阶预测的分位误差为 4(0.75−0.775)=−0.10。这只是一阶近似,但方向与比例都由公式清楚给出。

密度为零时,速度真的改变 ​

在 [−1,1] 上取分布函数 F(x)=(1+x3)/2,区间外分别补为零和一。中位数是零,密度 f(x)=3x2/2 在零处消失。

此时对固定 t,

F(tn−1/6)=12+t32n.

用“经验中位数不超过某高度”等价于该高度以下的二项计数至少达到中位秩,再对这个二项计数作中心极限定理,可得

P(n1/6q^n≤t)⟶Φ(t3).

极限不是正态,速度是 n−1/6 而非 n−1/2。把 f(0)=0 塞进分母,只会写出无意义的表达式,不能替代这项新的尺度分析。

推论与应用

第一步:把经验分位限制在根号 n 邻域 ​

记 q±=q±M/n。密度连续且为正,故对大 n,F(q+)−p≥f(q)M/(2n),p−F(q−)≥f(q)M/(2n)。

事件 q^n>q+ 蕴含 Fn(q+)<p。而 Fn(q+) 的方差至多 1/(4n),Chebyshev 界便给 P(q^n>q+)≤C/M2。左侧同理。先让 n 大,再让 M 大,得到

n(q^n−q)=OP(1).

第二步:证明局部经验振荡消失 ​

对固定 M,定义

An(h)=n[(Fn−F)(q+h/n)−(Fn−F)(q)],|h|≤M.

要证的是 sup|h|≤M|An(h)|→0 依概率。先固定一个 h:增量是长度 O(n−1/2) 小区间内的中心化样本比例,其方差满足

Var(An(h))≤P{X 位于 q 与 q+h/n 之间}≤CMn−1/2→0.

再把 [−M,M] 分成 J 个小格,格点为 hj。对固定有限 J,把各格点的 Chebyshev 界用并集界相加,所有格点的 An(hj) 同时依概率趋零。若 hj≤h≤hj+1,由 Fn 单调且邻域内 f≤C,

An(hj)−2CMJ≤An(h)≤An(hj+1)+2CMJ.

因此先让 n→∞,再让 J→∞,就控制了整个局部上确界。这个有限网格加单调夹逼的步骤,补上了不能直接对随机位置套固定点极限的缺口。

结合第一步的根号 n 紧性,可把 h=n(q^n−q) 代入,得到

n[(Fn−F)(q^n)−(Fn−F)(q)]=oP(1).

第三步:反演阶梯误差并应用中心极限定理 ​

以趋于一的概率,经验分位位于有连续密度的邻域,样本在那里没有并列值,故

0≤Fn(q^n)−p<1/n.

同时 F(q^n)−F(q)=f(q)(q^n−q)+oP(n−1/2)。把局部振荡式整理为

Fn(q^n)−p=Fn(q)−p+f(q)(q^n−q)+oP(n−1/2),

即可解出形式陈述的线性展开。最后,nFn(q) 是参数 p 的二项计数,中心极限定理给其中心化极限;Slutsky 定理消去余项,得到分位正态极限。

展开中的单观测项正是分位泛函的影响函数。在这里,影响函数能用于抽样方差的理由,是随机逆与局部振荡的余项已经得到证明,而不只是算出了污染方向的导数。

参考资料
关系图谱17 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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