Skip to content

定理Theorem

Neyman 随机化方差

Neyman randomization variance · Finite-population treatment-effect variance

在固定潜在结果和完全随机分配下推导均值差的精确方差,分离不可观测的效应异质项,并证明常用方差估计只在设计期望上保守。

同一批人重新抽签,处理组与对照组的均值差会改变。Neyman 方差回答这种改变有多大:无需先假定每个人的效应相同,却必须把两栏潜在结果和互补的分组同时放进计算。

形式陈述 ​

固定表与完全随机分配 ​

取 N≥2 个单位,固定有限实数潜在结果 Yi(1),Yi(0)。预先固定 1≤n1≤N−1,令 n0=N−n1;所有大小为 n1 的处理组等概率。这是不放回设计,两组的单位索引互为补集。

记 Y¯(z)=N−1∑iYi(z)、τi=Yi(1)−Yi(0)、τ=Y¯(1)−Y¯(0),并定义有限总体量

Sz2=1N−1∑i{Yi(z)−Y¯(z)}2,Sτ2=1N−1∑i(τi−τ)2.

观测均值差为

τ^=1n1∑iZiYi(1)−1n0∑i(1−Zi)Yi(0).

对分配取期望和方差,有精确恒等式

EZτ^=τ,VarZ(τ^)=S12n1+S02n0−Sτ2N.

当两臂人数都至少为二,令 sz2 为该臂已观测结果、分母为 nz−1 的样本方差。通常使用

V^=s12n1+s02n0,EZV^−VarZ(τ^)=Sτ2N≥0.

这里“保守”说的是重复分配的平均值。它不是每份数据上 V^ 都大于真实方差,也不单独保证正态置信区间具有有限样本覆盖。

直觉

两组互补产生哪一项协方差 ​

由设计对称性,EZi=n1/N。令 f=n1/N,对 i≠j 有

Var(Zi)=f(1−f),Cov(Zi,Zj)=−f(1−f)N−1.

负协方差来自固定总人数:选中一个单位,会稍微减少另一个被选中的机会。对任何固定数列 ai,用协方差双线性展开可得

Var(∑iZiai)=n1n0NSa2,Sa2=1N−1∑i(ai−a¯)2.

把估计量写成

τ^=∑iZi{Yi(1)n1+Yi(0)n0}−1n0∑iYi(0).

后项固定。取 ai=Yi(1)/n1+Yi(0)/n0,得到

Var(τ^)=n0Nn1S12+n1Nn0S02+2NS10,

其中 S10=(N−1)−1∑i{Yi(1)−Y¯(1)}{Yi(0)−Y¯(0)}。再用 Sτ2=S12+S02−2S10,便得到形式陈述中的公式。不能把两臂当作独立抽样,直接遗漏最后的异质项。

为什么样本方差能估计各栏,却补不回异质项 ​

每个处理组是 Y(1) 这一固定栏的简单随机样本,对照组同理。用成对差恒等式

sz2=1nz(nz−1)∑i<j:Zi=Zj=z{Yi(z)−Yj(z)}2

以及同臂成对纳入概率 nz(nz−1)/[N(N−1)],取期望得到 Esz2=Sz2。两项相加即证明保守期望公式。

Sτ2 需要每个人的两栏结果;一个实验不能同时观察它们。因此在不增加效应结构时,不能照公式把该项直接代入估计。恒定个体效应使 Sτ2=0,此时 V^ 无偏;异质性通常只支持保守性,而非精确恢复方差。

例子与边界

四个单位的全部六次抽签 ​

令 Y(0)=(1,2,3,4)、Y(1)=(2,4,6,8),两人处理、两人对照。目标 τ=5/2,总体方差为 S02=5/3、S12=20/3、Sτ2=5/3,所以真实分配方差为 15/4。

处理组 τ^ V^
1、2 −1/2 5/4
1、3 1 5
1、4 5/2 37/4
2、3 5/2 13/4
2、4 4 5
3、4 11/2 5/4

六个估计量平均为 5/2,对该均值的平方偏差平均为 15/4。六个方差估计平均为 25/6,超出真实方差 5/12=Sτ2/N。但第一行的 V^=5/4 明显小于 15/4,已经否定逐次保守的解释。

公式适用的设计与目标 ​

若处理人数随机变化、各单位分配机会不等,或只接受某些平衡分配,上述 Zi 协方差就改变了。不能只将实际人数填进完全随机公式。分层设计应使用块内方差加权;配对设计每臂每块只有一人,更不能使用未定义的块内样本方差。

即使误差分布看起来接近正态,τ^±1.96V^ 也只是需要另外条件的近似程序。有限总体中心极限定理通常要求没有少数单位主导波动;这里的精确均值、方差恒等式本身没有给出这样的尾部结论。

推论与应用

一次计算和一份完整解释 ​

实际数据只需分别累计各臂人数、结果和及平方偏差,计算 τ^,V^。使用稳定的在线中心化递推可避免“大平方和减去大均值平方”的相消;时间 O(N)、额外空间 O(1)。完整潜在结果表只存在于教学例或模拟中,此时才能枚举 (Nn1) 次分配核对真实方差。

报告应写明目标是当前 N 人的平均效应,随机性来自完全随机分配,以及方差估计的保守含义。若还想推断另一个总体,需要抽样或外推假设;随机抽签没有自动把这批单位变成总体的代表样本。

迁移核对。 将例中每个人的处理结果改为 Yi(0)+2,异质项消失,S12=S02=5/3,真实方差变为 5/3,而非零。每人效应相同不等于每次组均值差都相同:两组原本的水平仍可能不平衡。

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

拖动节点调整位置。

显示关系

显示:依赖

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

使用的工具

被这些条目使用