Skip to content

定理Theorem

分层抽样与 Neyman 分配

Stratified sampling · Neyman allocation

在独立层内抽样的精确方差上优化人数或费用分配,说明 Neyman 公式与实际整数约束。

形式陈述 ​

固定预算时,样本应该更多分给人数多的层,还是差异大的层?把总体划分为互不重叠且穷尽的 H 层,第 h 层大小为 Nh。在各层独立作简单随机不放回抽样,抽取 1≤nh≤Nh 人。先设各层 Nh≥2,定义总体层内方差 Sh2=∑i∈Uh(yi−Y¯h)2/(Nh−1),其中 Y¯h 是该层总体均值。总量估计及其设计方差为

t^=∑hNhy¯h,Vp(t^)=∑hNh2(1−nhNh)Sh2nh.

若每人的成本相同、总样本量固定为 n,暂时把 nh 当连续变量,各层 Sh>0,且最优解未触及各层人数上下限,则 Neyman 分配是

nh=nNhSh∑gNgSg.

若各层单位成本为 ch>0,在 ∑hchnh=C 下,最优比例改为 nh∝NhSh/ch。这两个约束不同,应先说明固定的是人数还是费用。

直觉

分层先确保每类总体都有计划地被代表;分配再决定每层测得多精细。层越大,它的均值误差对总体总量影响越大;层内波动越大,估计同样精度需要更多样本。

在方差公式中,有限总体修正对应的 −∑hNhSh2 不随分配变化。最小化只需考虑 ∑hNh2Sh2/nh。使用Lagrange 乘子法,加上约束的乘子项 λ(∑hnh−n),求导得到 −Nh2Sh2/nh2+λ=0,即 nh∝NhSh。成本约束同样给出分母的 ch。

乘子方程先给出候选点,全局最优性还可由Cauchy–Schwarz 不等式直接认证:

(∑hNh2Sh2nh)(∑hnh)≥(∑hNhSh)2.

等号恰在 nh∝NhSh 时成立。对成本约束,改用两组向量 NhSh/nh 与 chnh,同样得到等号条件 nh∝NhSh/ch。这证明的是所述连续内点问题;人数上下限生效时还须重新求解。

例子与边界

比例分配与最优分配的实算 ​

两层大小为 100,300,层内标准差为 10,20,总样本量 80。按人数比例分配得到 (20,60);Neyman 连续解为 (80/7,480/7)≈(11.43,68.57),可取邻近整数 (11,69) 后比较方差。

比例分配的总量方差是

1002(1−20/100)10220+3002(1−60/300)20260=520000.

整数分配 (11,69) 给出约 482648,更小。第二层的总体规模和内部波动都较大,因此即使已经人数较多,仍值得分配更高比例的测量资源。

若第二层每人的成本变成第一层的四倍,则成本最优的样本数比例为 100⋅10:300⋅20/2=1:3。最终人数还要用实际总预算缩放,不能沿用固定 80 人的约束。

必要的设计边界 ​

实际规划中的 Sh 通常来自历史资料或先导调查。若它估得不准,声称的最优性也随之改变。多目标调查还要照顾多个指标,而非只优化一个变量的方差。

若公式要求 nh>Nh,该层应全查,再对剩余预算重新分配。需要估计层内方差时,通常也须设置每层至少两个样本单位等下限。单人层直接全查,贡献零设计方差;零方差层按所需覆盖和最低样本量约束单独安排。整数、全查、最低样本量和访问成本都应进入最后的可执行方案。

推论与应用

这是设计阶段的优化:分层和分配在观察结果之前确定。事后按样本结果随意重组层,再套原设计方差,会改变随机性结构。若辅助变量已知总体总量,则分析阶段还可用校准进一步提高精度。

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

拖动节点调整位置。

显示关系

显示:依赖

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