跳到正文
格致开物MATHWIKI

辅助统计量

辅助统计量(ancillary statistic)是在指定模型中分布不随未知参数改变的统计量。它仍是随机的,也可能有助于检查模型或进行条件推断;“辅助”不表示其观测值没有用途。必须先说明哪些量未知,因为同一个统计量对一个参数可以辅助,对扩大的参数向量却未必辅助。

设两个独立读数 X1,X2∼N(μ,4),未知量只有共同位置 μ。两数之差 D=X1−X2∼N(0,8) 不依赖 μ,所以是这个位置模型的辅助统计量。共同位置被相减消去,读数间的随机差异则保留下来。

平移后的同一差值

取教学读数 (5,9),均值 M=7、差值 D=−4。两读数都增加十,变成 (15,19),均值变成十七,而差仍为负四。

两组读数共同平移十,均值随之平移而差值始终为负四
箭头连接的是同一对读数的平移。示意说明位置被相减消去;辅助性还需用整个抽样分布确认。

把模型写成 Xi=μ+σZi,其中 Zi 独立标准正态、已知 σ=2,则 D=σ(Z1−Z2)。独立正态之差均值为零、方差为二,所以得到前面的 N(0,8)。这证明的是对每个 μ∈ℝ 的分布相同,不只是两组观察值碰巧相减相等。

一般定义是:对模型 {Pθ},统计量 A(X) 辅助,当且仅当其诱导分布 PθA 对所有参数都相同。它是一个边缘分布条件,不自动断言它与样本均值或其他统计量独立。[1]

正态模型中均值与差的独立证明

对一般已知 σ>0,令 M=(X1+X2)/2、D=X1−X2。反变换为 X1=M+D/2,X2=M−D/2, Jacobian行列式的绝对值为一。平方项展开为 (x1−μ)2+(x2−μ)2=2(m−μ)2+d22. 因此联合密度恰好分解为 fM,D(m,d)=1πσe−(m−μ)2/σ2×12πσe−d2/(4σ2). 两因子分别是 N(μ,σ2/2) 与 N(0,2σ2) 的密度,所以 M,D 独立。这是对正态模型的直接计算;并不是从“差是辅助统计量”这一句话单独推出独立性。巴苏定理给出另一种在完备充分条件下证明独立性的办法。

主例中,M∼N(μ,2),而差的抽样分布为 N(0,8)。例如在该模型下,差的绝对值至少为四的概率为 2Φ(−48)≈0.1573, 其中 Φ 是标准正态分布函数。这个尾概率不随共同位置改变,不能拿它来区分 μ=0 和 μ=10。

参数扩大时的变化

如果把方差也改为未知,D 的方差是 2σ2,其分布已依赖参数向量 (μ,σ2)。此时它不是完整两参数模型的辅助统计量,尽管它的分布仍不依赖位置坐标。

变式。 保持方差已知,但改为 σ=3,则 D∼N(0,18),仍对未知位置辅助;不能继续套用方差八。辅助性是在一个固定模型内考察参数无关性,不要求更换已知实验精度后也有同一分布。

辅助性不保证任意独立性

取两次独立均匀读数 Xi∼U(μ,μ+1)。极差 R=|X1−X2| 仍因平移消去而具有参数无关的分布。可是,给定 R=r∈(0,1) 后,均值只能落在 μ+r2≤M≤μ+1−r2 之内,而且条件密度在该区间均匀。这个区间随 r 改变;例如极差接近一时,均值必靠近区间中点。因此 M 与 R 不独立。

该条件分布可用 (M,X1−X2) 的面积变换得到:联合密度在两个对应符号的可行区域内为常数,固定绝对差后,均值的允许线段长度恰为 1−r。这个反例把辅助分布与联合独立性明确分开,也说明条件推断为什么可能利用辅助量观察到的大小,而不把它解释为关于位置的直接证据。

参考文献

  1. ↑ Will Fithian,Completeness, Ancillarity, and Basu's Theorem,“Ancillarity”及“Basu's theorem”,关于参数无关分布与独立性条件的区别。