跳到正文
格致开物MATHWIKI

最小充分统计量

最小充分统计量(minimal sufficient statistic)是能够从任意其他充分统计量恢复的充分摘要。“最小”比较的是信息包含关系,不是程序占用的字节数,也不是表达式里写了几个坐标。两个摘要若能可测地互相恢复,表示的是同一层信息。

以独立正态样本 X1,…,Xn∼N(μ,σ2) 为例,取 n≥2,未知参数为 μ∈ℝ、σ2>0。均值告诉我们样本中心,却没有说明点离中心多远。把中心相同而分散程度不同的样本合并,会损失关于未知方差的信息。

相同均值保留了哪些差别

取两组教学样本 x=(0,2,4),y=(1,1,4). 两者均值都是二,但 ∑ixi2=20,∑i(xi−2)2=8,∑iyi2=18,∑i(yi−2)2=6. 在完整的两参数正态模型中,它们的密度之比为 exp⁡(−1/σ2)。这个比值随方差变化,所以这两个样本不能在不损失参数信息的前提下合并成同一个“均值二”。

两个均值相同的样本在数轴上的分散程度不同,平方离差和分别为八与六
第二组的数值一出现两次,图中用两个标记表示。共同均值不能区分图中的两种分散程度。

恢复关系与正式定义

固定模型与样本空间。统计量 T 最小充分,是指它本身充分,且对每个充分统计量 S,都存在不依赖未知参数的可测函数 a,使 T(X)=a(S(X))Pθ-几乎处处,对每个 θ. 这不是要求能从 T 恢复原样本。恰恰相反,它要求每种充分摘要都不能进一步丢掉 T 所表示的信息。通常的有限或欧氏模型中,似然比判据可以识别这种摘要。[1]

对具有共同正密度的常规模型,一个常用判据是:若两个样本点给出相同的 T 值,恰好等价于其密度之比不依赖参数,则 T 是最小充分统计量。这里使用这一判据;以下完整证明它在两参数正态模型中的等价条件,而不把单个正态例子当作所有模型的测度论证明。

正态模型的两个必要摘要

记 S(x)=∑ixi,Q(x)=∑ixi2. 联合密度可以展开为 pμ,σ2(x)=(2πσ2)−n/2×exp⁡[−Q(x)−2μS(x)+nμ22σ2]. 它只通过 (S,Q) 依赖样本,由因子分解定理可知该二元摘要充分。

比较任意两个样本 x,y,令 ΔS=S(x)−S(y)、ΔQ=Q(x)−Q(y),则 log⁡pμ,σ2(x)pμ,σ2(y)=μΔSσ2−ΔQ2σ2. 如果比值与全部 μ,σ2 无关,先固定任意正方差并改变 μ,就必须有 ΔS=0;再改变方差,必须有 ΔQ=0。反之,两差都为零时比值恒为一。因此比值与参数无关,当且仅当两个样本的 (S,Q) 相同,判据给出最小充分性。

样本均值与平方离差和 X¯=Sn,V=∑i(Xi−X¯)2=Q−S2n 可以从 (S,Q) 恢复;反过来 S=nX¯、Q=V+nX¯2。所以 (X¯,V) 是等价的最小充分摘要。若采用样本方差 V/(n−1),固定样本量下也没有改变这层信息。

从主例核对判据

第一组样本的 (S,Q)=(6,20),第二组为 (6,18)。它们只在第二个坐标不同;代入比值公式,得到前面的 e−1/σ2。若把第一组重排为 z=(4,0,2),则两个坐标都不变,密度之比恒为一。被舍弃的是独立同分布模型中无关的排列次序,而不是所有样本差异。

变式。 假如方差已知为四,只有均值未知,比值中的 −ΔQ/8 已经是一个与未知均值无关的常数。此时只需 ΔS=0,所以 S 或 X¯ 就是最小充分统计量。最小性随参数空间改变,不能脱离“哪些参数未知”来谈。

充分性与最小性也不等于完备性。前两者讨论摘要保留了哪些信息,完备性则对该摘要所有可积函数的期望施加条件,作用不同。存在参数相关支撑或密度为零的模型时,不能把上面的正密度比值直接写成 0/0 再消去;应先检查支撑和适用的似然等价判据。

参考文献

  1. ↑ Will Fithian,Sufficiency,“Minimal sufficiency”,关于可恢复关系及似然比等价的判据。连续情形的零测集按支配模型的通常版本处理。