Skip to content

模型Model

分层随机实验

Blocked randomized experiment · Stratified randomized experiment · Randomized block design

先按处理前资料固定区组,再在各块独立随机分配,以单位数权重合并块内效应,推导设计方差并检查不同处理比例与小块的估计边界。

把起点相似的单位放在一起,再在每组内部抽签,可以防止全部高起点单位碰巧进入同一臂。区组不是分析结束后随意贴上的标签:它规定哪些分配在实验中有机会发生,也决定均值差应该怎样加权。

形式陈述 ​

块内完全随机,块间独立 ​

在处理前将固定的 N 个单位划分为不相交区组 B1,…,BH,大小为 Nh≥2,总和为 N。每块固定处理人数 1≤n1h<Nh,令 n0h=Nh−n1h;块内从所有相应大小的子集等概率选择,各块抽签相互独立。

潜在结果表固定;每块各自使用等概率固定人数设计。记第 h 块平均效应为 τh,目标是所有单位的平均效应

τ=∑h=1HWhτh,Wh=Nh/N.

令 τ^h 为该块处理均值减对照均值,取

τ^B=∑hWhτ^h.

由块内Neyman 方差,有

EZτ^B=τ,VarZ(τ^B)=∑hWh2(S1h2n1h+S0h2n0h−Sτh2Nh).

三个有限总体方差都在块内计算,分母为 Nh−1。若每块每臂至少二人,便能使用

V^B=∑hWh2(s1h2n1h+s0h2n0h),

且其设计期望超出真实方差 ∑hWh2Sτh2/Nh≥0。小于二人的臂仍可定义均值差和真实方差,却没有对应的普通样本方差估计。

直觉

加权先由目标决定,再由设计计算波动 ​

每块的处理与对照都是同一块的互补样本,所以块内均值差无偏估计本块效应。期望线性性给出整体无偏性;块间独立使不同块估计量的协方差为零,再把每块方差乘上权重的平方即可。

Wh=Nh/N 不是精度调参。它表达每个单位在总体目标中占相同份量。若大块四人、小块两人,效应分别为三和零,单位平均效应是二;对两块简单平均会得到 1.5,那是另一个“块平均”目标。即使后者方差较小,也不能把估计目标悄悄替换。

这与分层调查抽样不同:调查从每层抽一份样本测同一个固定变量,本设计在每块把单位分到两个互补组,分别揭示两栏潜在结果。因此块内方差带有不可观测的效应异质项,不能直接套一栏抽样的分配公式。

例子与边界

两块八人的完整潜在结果 ​

下表只是用于复算的完整教学表;实际实验不会同时看到同一个人的两栏。

区组 单位 Y(0) Y(1)
A 1、2、3、4 0、1、2、3 1、1、4、4
B 5、6、7、8 10、11、12、13 11、13、12、14

每块选两人处理,WA=WB=1/2。两块平均效应均为一;A 块有 S1A2=3,S0A2=5/3,SτA2=2/3,方差为 13/6。B 块相应为 5/3,5/3,2/3,方差为 3/2。所以

Var(τ^B)=14(136+32)=1112.

总共有 (42)2=36 份等概率分配,逐份计算的平均估计为一,方差正是 11/12。若同八人改为完全随机选四人处理,共有七十份分配,方差为 211/14。在这个例子里,两个块的起点相差十,阻止块间人数失衡确实大幅减少了波动。

分层不是对所有固定表都降低方差 ​

另取四人的两栏都为 (0,1,0,1),按前两人和后两人分块,每块各一人处理。块内差分别独立取 ±1,两块平均的方差为 1/2。完全随机选两人处理时,有限总体方差 S2=1/3,均值差的方差为 S2/2+S2/2=1/3。

这份分块反而增大方差。结论取决于固定表、块的形成方式和比较的抽样框架,不能从“增加平衡约束”推出所有结果变量都更精确。

不同处理比例时,合并两臂均值会错位 ​

两块各四人,第一块的两栏结果都为零,第二块都为十,因此真实效应全为零。若第一块处理三人、第二块处理一人,合并处理均值为 2.5,合并对照均值为 7.5,差恒为 −5。块内差则都为零,τ^B=0。

原因是两臂在两个块的组成权重不同。只有所有块处理比例相同等适当条件下,合并两臂均值差才等于按 Nh/N 加权的块内差。区组标签不能在最后一步被丢掉。

推论与应用

检验同样只能在块内重抽签 ​

针对逐人无效应的尖锐零假设,调用Fisher 随机化检验:固定补全后的结果,在每块保留原处理人数,独立重抽,再计算预定的加权统计量。允许分配数为

M=∏h(Nhn1h).

八人例可完整枚举三十六次;把全部八人任意洗成四对四会使用七十份不同的分布,不再复现原实验。事后分层则是另一个分析程序,其分层人数可能随机变化,不能冒充事前分层设计。

实际点估计及可用的方差估计可在 O(N) 时间、O(H) 存储内按块累计。完整随机化检验的朴素时间为 O(MN);模拟时也必须分别在块内抽样。

当每块恰为一对、每对一人处理时,进入配对随机实验。它仍有精确设计方差,却要另用跨对差异构造保守方差估计。识别这个边界,比给未定义的块内样本方差补一个零更重要。

参考资料
关系图谱14 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

暂未标注直接上位概念。

下位 / 直接特例

类型化关系