Skip to content

定理Theorem

Delta 方法

施工中 · 内容逐步修订

Delta method

用可微映射的一阶 Taylor 展开传播估计量的渐近分布与协方差。

形式陈述 ​

设 rn→∞ 且

rn(Tn−θ)→dZ,

其中箭头表示依分布收敛,Tn,θ∈Rk。若映射 g:U→Rm 在 θ 的某个开邻域 U 上 Borel 可测,并在 θ 处 Fréchet 可微,即存在一个线性映射作为局部导数,使线性近似的误差相对于输入偏移趋于零,则

rn{g(Tn)−g(θ)}→dDg(θ)Z.

为使 g(Tn) 对每个有限 n 都有定义,约定在 U 外将 g 延拓为一个预先固定的有限向量,例如零向量。此延拓仍 Borel 可测。由给定的缩放极限和 rn→∞ 可得 Tn→θ 依概率,故 Pr(Tn∉U)→0;这些例外事件上的取值不影响上述极限。仅在一点可微不自动保证整个邻域上的可测性。

这里 Dg(θ) 是导数的 m×k 矩阵(Jacobian),把输入的微小误差转成输出误差。速率 rn 是已知极限的一部分,常见取值为 n。

特别地,若 Tn 具有渐近正态性,即

n(Tn−θ)→dNk(0,V),

则

n{g(Tn)−g(θ)}→dNm(0,Dg(θ)VDg(θ)T).

协方差公式使用多元正态分布的线性变换规律。结论是对每个固定参数点的渐近陈述。若 Jacobian 秩不足,极限协方差退化;不能把零矩阵误报为“估计没有误差”。

直觉

一阶推导 ​

由 Fréchet 可微性的定义,

g(Tn)−g(θ)=Dg(θ)(Tn−θ)+r(Tn−θ),

其中 ‖r(h)‖/‖h‖→0。符号 OP(1) 表示随机量以高概率保持有界,oP(1) 表示依概率趋零。分布收敛蕴含 rn(Tn−θ)=OP(1),因而 Tn→θ 依概率,并有

rnr(Tn−θ)=oP(1).

首项由线性连续映射传播 Z,余项再由Slutsky 定理消失。这说明可微性需要统一余项控制,仅有各方向导数不足。

例子与边界

独立同分布 Bernoulli 样本的比例 p^=n−1∑iXi 满足

n(p^−p)→dN(0,p(1−p)),0<p<1.

取 logit g(p)=logp/(1−p),导数为 g′(p)=1/[p(1−p)],故

n{g(p^)−g(p)}→dN(0,1p(1−p)).

若 n=500,p^=0.4,logit 估计约 −0.405,估计标准误为

1500(0.4)(0.6)≈0.0913.

以标准正态分位数 1.96 作近似,logit 区间为 −0.4055±1.96(0.0913),约 [−0.5844,−0.2265]。对两个端点应用反变换 h(t)=et/(1+et),得到概率区间约 [0.358,0.444]。这说明方差传播先在 logit 坐标完成,再把端点送回概率空间。

有限样本可能有 p^=0 或 1,此时 logit 无定义。对固定内部点 p,这些事件的概率 (1−p)n+pn 趋零,可在其上另作有限赋值以陈述极限定理;实际遇到全成功或全失败数据时,仍需另选区间程序。仅有“端点在 (0,1)”并不保证覆盖良好。

导数退化时 ​

若 g′(θ)=0,一阶 Delta 方法只给退化零极限。若标量 g 二阶可微、g′(θ)=0 且 rn(Tn−θ)⇒Z,则二阶展开给

rn2{g(Tn)−g(θ)}→d12g″(θ)Z2.

例如 Tn=X¯、E[X]=0、nX¯⇒N(0,σ2),取 g(t)=t2,则

nX¯2→dσ2χ12,

这里用到标准正态平方服从自由度一的卡方分布;所得并不是 n 尺度的正态极限。

边界与失败情形 ​

映射在参数点不可微时,普通 Delta 方法失效。g(t)=|t| 在零点把正态极限变成折叠正态,需要方向可微 Delta 方法。

参数在边界时,Tn 的原始极限本身可能非正态;不能先假定正态再传播。logit 导数在 p→0,1 发散,也提示内部点近似不统一。

plug-in Jacobian Dg(Tn) 替换 Dg(θ) 需要导数连续与 Tn 一致。数值上导数巨大时,渐近方差会放大,有限样本近似可能很差。

比值给出必须保留协方差的多元例子:g(a,b)=a/b 在 b≠0 时梯度为 (1/b,−a/b2)。若分子、分母由同一样本估计,渐近方差含交叉协方差;把两个标准误分别相除会漏掉依赖,并可能得到完全错误的尺度。

若映射 gn 随样本量变化,固定函数版本不能直接套用。增长维数、变化截断阈值或正则化强度都要求统一 Taylor 余项;形式上对每个 n 可微,不保证 rn 倍余项仍依概率消失。

推论与应用

Delta 方法传播比值、对数、风险差和预测函数的标准误。多元情形保留协方差交叉项,不能只逐分量变换方差。

它提供一阶局部近似,不给有限样本偏差;曲率较强时可用二阶修正、Bootstrap或直接枢轴方法核验。

Delta 方差的 plug-in 估计还需要原协方差估计一致。若 Dg(θ) 有非满秩,输出某些方向没有一阶波动,报告矩阵逆或标准化统计量前应先检查秩;二阶项可能在这些方向上成为主导。

对置信区间,先在线性化坐标构造再反变换,与直接在原坐标使用对称区间通常只首阶等价。有限样本端点、覆盖和不对称性可能不同,应根据参数空间和数值稳定性选择并说明。

比值的固定分母条件不能悄悄随样本消失 ​

若 bn=c/n、an=r∗bn,其中 c≠0,且分子分母的根号样本量误差趋于非退化联合正态 (ZA,ZB),则观测比值趋于 (r∗c+ZA)/(c+ZB)(样本分母为零时预先赋固定有限值;由缩放分母的连续极限,这个事件概率趋零)。分母的相对误差不趋零,不能把随机分母统一替换成 bn;对每个固定非零分母成立的一阶结论没有自动覆盖这条变化参数序列。

Fieller 方法改为反演候选比值的线性对比检验,并保留有限区间、两条射线、半直线或全实轴。其精确版本明确要求联合高斯且估计量协方差已知;若换成估计协方差或一般渐近实验,应重新证明检验校准,不能把“没有先相除”本身当成覆盖证明。

参考资料
  • John Duchi(授课)、Sky Cao(记录),Stats 300b: Theory of Statistics, Lecture 2,Stanford University,讲义首页标注 Winter 2018,第 2–3 页 Theorem 3–4(一阶及二阶 Delta 方法)。
  • A. W. van der Vaart, Asymptotic Statistics, Cambridge University Press, 1998,§3.1。
  • R. J. Serfling, Approximation Theorems of Mathematical Statistics, Wiley, 1980,§1.7。
  • Thomas S. Ferguson, A Course in Large Sample Theory, Chapman & Hall, 1996,Ch. 2。
关系图谱16 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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