形式陈述
设 r n → ∞ 且
r n ( T n − θ ) → d Z , 其中箭头表示依分布收敛 理路 依分布收敛 Convergence in distribution · Weak convergence 分布函数在极限分布连续点处收敛的随机变量收敛概念。 ,T n , θ ∈ R k 。若映射 g : U → R m 在 θ 的某个开邻域 U 上 Borel 可测,并在 θ 处 Fréchet 可微,即存在一个线性映射作为局部导数 理路 多元函数导数 Derivative in several variables · Jacobian derivative Fréchet 导数在有限维 Euclidean 坐标中的 Jacobian、梯度与方向导数表示。 ,使线性近似的误差相对于输入偏移趋于零,则
r n { g ( T n ) − g ( θ ) } → d D g ( θ ) Z . 为使 g ( T n ) 对每个有限 n 都有定义,约定在 U 外将 g 延拓为一个预先固定的有限向量,例如零向量。此延拓仍 Borel 可测。由给定的缩放极限和 r n → ∞ 可得 T n → θ 依概率,故 Pr ( T n ∉ U ) → 0 ;这些例外事件上的取值不影响上述极限。仅在一点可微不自动保证整个邻域上的可测性。
这里 D g ( θ ) 是导数的 m × k 矩阵(Jacobian),把输入的微小误差转成输出误差。速率 r n 是已知极限的一部分,常见取值为 n 。
特别地,若 T n 具有渐近正态性 理路 估计量的渐近正态性 Asymptotic normality of estimators 估计误差在明确中心、速率与协方差下趋于正态律,以及线性化、标准化和边界失效机制。 ,即
n ( T n − θ ) → d N k ( 0 , V ) , 则
n { g ( T n ) − g ( θ ) } → d N m ( 0 , D g ( θ ) V D g ( θ ) T ) . 协方差公式使用多元正态分布的线性变换规律 理路 多元正态分布 Multivariate normal distribution · Multivariate Gaussian distribution · Jointly Gaussian vector 以所有线性组合都正态刻画联合高斯向量,并由特征函数连接线性构造、退化支撑、全维密度与条件分布。 。结论是对每个固定参数点的渐近陈述。若 Jacobian 秩不足,极限协方差退化;不能把零矩阵误报为“估计没有误差”。
直觉
一阶推导
由 Fréchet 可微性的定义,
g ( T n ) − g ( θ ) = D g ( θ ) ( T n − θ ) + r ( T n − θ ) , 其中 ‖ r ( h ) ‖ / ‖ h ‖ → 0 。符号 O P ( 1 ) 表示随机量以高概率保持有界,o P ( 1 ) 表示依概率趋零。分布收敛蕴含 r n ( T n − θ ) = O P ( 1 ) ,因而 T n → θ 依概率,并有
r n r ( T n − θ ) = o P ( 1 ) . 首项由线性连续映射传播 Z ,余项再由Slutsky 定理 理路 Slutsky 定理 Slutsky's theorem 依分布收敛随机量与依概率收敛常量组合时,和、积与合法商保持相应分布极限。 消失。这说明可微性需要统一余项控制,仅有各方向导数不足。
例子与边界
独立同分布 Bernoulli 样本的比例 p ^ = n − 1 ∑ i X i 满足
n ( p ^ − p ) → d N ( 0 , p ( 1 − p ) ) , 0 < p < 1. 取 logit g ( p ) = log p / ( 1 − p ) ,导数为 g ′ ( p ) = 1 / [ p ( 1 − p ) ] ,故
n { g ( p ^ ) − g ( p ) } → d N ( 0 , 1 p ( 1 − p ) ) . 若 n = 500 , p ^ = 0.4 ,logit 估计约 − 0.405 ,估计标准误为
1 500 ( 0.4 ) ( 0.6 ) ≈ 0.0913 . 以标准正态分位数 1.96 作近似,logit 区间为 − 0.4055 ± 1.96 ( 0.0913 ) ,约 [ − 0.5844 , − 0.2265 ] 。对两个端点应用反变换 h ( t ) = e t / ( 1 + e t ) ,得到概率区间约 [ 0.358 , 0.444 ] 。这说明方差传播先在 logit 坐标完成,再把端点送回概率空间。
有限样本可能有 p ^ = 0 或 1 ,此时 logit 无定义。对固定内部点 p ,这些事件的概率 ( 1 − p ) n + p n 趋零,可在其上另作有限赋值以陈述极限定理;实际遇到全成功或全失败数据时,仍需另选区间程序。仅有“端点在 ( 0 , 1 ) ”并不保证覆盖良好。
导数退化时
若 g ′ ( θ ) = 0 ,一阶 Delta 方法只给退化零极限。若标量 g 二阶可微、g ′ ( θ ) = 0 且 r n ( T n − θ ) ⇒ Z ,则二阶展开给
r n 2 { g ( T n ) − g ( θ ) } → d 1 2 g ″ ( θ ) Z 2 . 例如 T n = X ¯ 、E [ X ] = 0 、n X ¯ ⇒ N ( 0 , σ 2 ) ,取 g ( t ) = t 2 ,则
n X ¯ 2 → d σ 2 χ 1 2 , 这里用到标准正态平方服从自由度一的卡方分布 理路 卡方分布 Chi-square distribution · Chi-squared distribution · χ² distribution 若干独立标准正态变量平方和的分布,以自由度记录独立平方方向的数量。 ;所得并不是 n 尺度的正态极限。
边界与失败情形
映射在参数点不可微时,普通 Delta 方法失效。g ( t ) = | t | 在零点把正态极限变成折叠正态,需要方向可微 Delta 方法。
参数在边界时,T n 的原始极限本身可能非正态;不能先假定正态再传播。logit 导数在 p → 0 , 1 发散,也提示内部点近似不统一。
plug-in Jacobian D g ( T n ) 替换 D g ( θ ) 需要导数连续与 T n 一致。数值上导数巨大时,渐近方差会放大,有限样本近似可能很差。
比值给出必须保留协方差的多元例子:g ( a , b ) = a / b 在 b ≠ 0 时梯度为 ( 1 / b , − a / b 2 ) 。若分子、分母由同一样本估计,渐近方差含交叉协方差;把两个标准误分别相除会漏掉依赖,并可能得到完全错误的尺度。
若映射 g n 随样本量变化,固定函数版本不能直接套用。增长维数、变化截断阈值或正则化强度都要求统一 Taylor 余项;形式上对每个 n 可微,不保证 r n 倍余项仍依概率消失。
推论与应用
Delta 方法传播比值、对数、风险差和预测函数的标准误。多元情形保留协方差交叉项,不能只逐分量变换方差。
它提供一阶局部近似,不给有限样本偏差;曲率较强时可用二阶修正、Bootstrap 理路 Bootstrap Bootstrap resampling · Nonparametric bootstrap 用条件重采样律近似抽样误差:完整计算均值的有限分布,证明有限方差下一致性,并由最大值的持久原子判定失败。 或直接枢轴方法核验。
Delta 方差的 plug-in 估计还需要原协方差估计一致。若 D g ( θ ) 有非满秩,输出某些方向没有一阶波动,报告矩阵逆或标准化统计量前应先检查秩;二阶项可能在这些方向上成为主导。
对置信区间,先在线性化坐标构造再反变换,与直接在原坐标使用对称区间通常只首阶等价。有限样本端点、覆盖和不对称性可能不同,应根据参数空间和数值稳定性选择并说明。
比值的固定分母条件不能悄悄随样本消失
若 b n = c / n 、a n = r ∗ b n ,其中 c ≠ 0 ,且分子分母的根号样本量误差趋于非退化联合正态 ( Z A , Z B ) ,则观测比值趋于 ( r ∗ c + Z A ) / ( c + Z B ) (样本分母为零时预先赋固定有限值;由缩放分母的连续极限,这个事件概率趋零)。分母的相对误差不趋零,不能把随机分母统一替换成 b n ;对每个固定非零分母成立的一阶结论没有自动覆盖这条变化参数序列。
Fieller 方法 理路 Fieller 比值置信集合 Fieller confidence set · Fieller ratio interval · 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。