Skip to content

定理Theorem

Lee 选择界与分数截尾

Lee selection bounds · Lee trimming bounds · 单调选择下的锐界

在随机处理和单调选择下,对始终可见人群的平均效应给出可达截尾界,保留原子分数质量,并推导条件化后的正确人群权重。

形式陈述 ​

随机分配以后,留下来的人仍可能不同 ​

某项培训随机分配给求职者。我们希望比较工资,但只有找到工作的人才有工资记录。若培训提高了就业率,处理组中的有工资者就包含一批“只有接受培训才就业”的人,而对照组没有对应人群。随机化保证分配前可比,却不能保证按就业筛选后仍可比。

Lee界先明确一个可以在两个世界比较的目标:无论是否处理,结果都可见的人。它不要求认出每个具体的人,而是从两组可见比例推断这个人群在混合物中的份额。

采用潜在结果Y(0),Y(1)和潜在选择指标S(0),S(1)∈{0,1}。处理D∈{0,1}满足0<P(D=1)<1,观察

S=S(D),O=(D,S,SY(D)).

假定

(1)D⊥(Y(0),Y(1),S(0),S(1)),S(1)≥S(0)几乎处处.

第二项是选择单调性,不是结果单调性:允许个体的处理效应正负不一。令pd=P(S=1∣D=d),要求0<p0≤p1。定义始终可见事件A={S(0)=S(1)=1},目标为

(2)τA=E{Y(1)−Y(0)∣A}.

假定两组可见结果的绝对一阶矩有限。令F为Y∣D=1,S=1的条件分布,μ0=E(Y∣D=0,S=1),保留比例为

(3)a=p0p1∈(0,1].

从F中保留恰好a的概率质量,最低a份额的均值记为μ−(a),最高a份额的均值记为μ+(a)。允许在并列值处保留一部分质量。则τA的锐识别集合恰为

(4)[μ−(a)−μ0, μ+(a)−μ0].
直觉

三种潜在人群与一份混合分布 ​

式(1)只允许三种类型:始终可见(S(0),S(1))=(1,1)、仅处理后可见(0,1)、始终不可见(0,0)。它们的总体质量依次为

p0,p1−p0,1−p1.

所以对照组可见者全是A;处理组可见者中,A占p0/p1=a。若FA,FI分别为这两类人的处理结果分布,则

(5)F=aFA+(1−a)FI.

我们知道混合分布和份额,但不知道哪部分属于A。让最低的a份额属于A,就把A的处理均值压到最小;让最高份额属于A,则推到最大。这是求相容人群,不是说应把实际数据中某些人的结果永久删除。

有原子时怎样“恰好保留” ​

先看0<a<1。取下分位阈值c=inf{y:F(y)≥a},则

PF(Y<c)≤a≤PF(Y≤c).

全部保留Y<c,在Y=c处只保留缺少的质量a−PF(Y<c)。因此

(6)μ−(a)=EF[Y1{Y<c}]+c{a−PF(Y<c)}a.

上端取d=inf{y:F(y)≥1−a},全部保留Y>d,再从阈值原子中补足,得到

(7)μ+(a)=EF[Y1{Y>d}]+d{a−PF(Y>d)}a.

若阈值处没有原子,需要补的质量为零。这些公式允许负结果和无界支持,只要求绝对可积。a=1单列为μ−=μ+=EFY,不必把无穷分位端点代进公式。

为什么任何其他保留方式都不会越界 ​

任一相容的A分配可写为保留概率h(y)=P(A∣Y(1)=y,S(1)=1),其中0≤h≤1、EFh=a。它的处理均值为EF(Yh)/a。

令h−表示刚才的低端保留规则。对每个y都有

(y−c){h(y)−h−(y)}≥0.

因为c以下h−=1,两个因子都非正;c以上h−=0,两个因子都非负;在c处乘积为零。积分并使用EFh=EFh−=a,便有

EF(Yh)≥EF(Yh−).

高端规则同理给反向不等式。这个证明同时涵盖连续值、并列原子和分数保留,并没有假设数据无并列。

如何造出达到端点的完整世界 ​

给定任一合法h,在总体中先按p0,p1−p0,1−p1生成三种选择类型。在A中令Y(1)的分布为hdF/a,在仅处理可见类型中令它的分布为(1−h)dF/(1−a);后者只在a<1时需要。令A中的Y(0)具有观察到的对照可见分布,与其Y(1)如何耦合可任选,因为这里目标只用边缘均值。其余从未可见的潜在结果任选可积值。

最后独立产生D。这份完整律满足随机化与单调性,恢复全部观测律,并达到h指定的效应。取h−、h+就得到两端;对两个规则作凸混合仍满足EFh=a,给出每个内部值。于是式(4)确实是锐界。

例子与边界

截到一个原子的中间 ​

设p0=1/2,p1=3/4,故a=2/3。处理组可见结果分布为

PF(Y=0)=1/4,PF(Y=2)=1/2,PF(Y=6)=1/4,μ0=1.

低端先取全部1/4的零,再从结果二处取2/3−1/4=5/12的质量。高端先取全部1/4的六,再从二处取同样的5/12。因此

μ−=2(5/12)2/3=54,μ+=6(1/4)+2(5/12)2/3=72,(8)τA∈[1/4,5/2].

处理可见者均值为5/2,粗差为3/2,位于界内,却没有被点识别为始终可见人群的效应。若低端直接保留所有Y≤2,保留了3/4而非2/3,得到均值4/3;这是另一种份额,不能替代5/4。

图中每列高度是概率质量,取样界线穿过并列原子是数学上必要的分数分配;并不是把一个人的实际结果改小。

选择率相等,仍需单调性 ​

在单调性下,若p0=p1>0,则非负变量S(1)−S(0)的期望为零,因而它几乎处处为零。两臂可见者对应同一潜在人群,a=1,界缩成均值差。

去掉单调性,选择率相等却可能由反向转换抵消。构造四个各占1/4的类型:始终可见者有Y(0)=Y(1)=0;仅对照可见者有Y(0)=2;仅处理可见者有Y(1)=8;另有始终不可见者。随机化后p0=p1=1/2,两臂可见均值分别为一和四,粗差为三,但始终可见者真实效应为零。表面的p1≥p0只是一条必要的边缘检查,不证明逐人的单调性。

方向、空人群和目标不能混淆 ​

若合理假设是S(1)≤S(0),应截尾对照可见分布,保留p1/p0;处理可见均值已经属于A。若对照截尾均值为ν−,ν+,效应界是[μ1−ν+,μ1−ν−],相减会交换上下端。

若p0=0且假设S(1)≥S(0),则P(A)=0,式(2)的条件效应没有定义,不能报告一个“无限宽Lee区间”。若总体p1<p0,则它与所声明的向上单调模型矛盾;有限样本估计率逆序则可能只是抽样波动,需要另行推断,不能擅自交换处理标签。

本页目标不是整个总体ATE,也不是实际处理组全部可见者的效应。支持无界时仍可约束A的均值,恰是因为这个人群的两份潜在结果都在相应可见混合物中,而非因为未见者的结果被设成零。

推论与应用

先在层内截尾,再用目标人群权重平均 ​

设处理前协变量X取有限多个值,wx=P(X=x)>0,各层处理概率在零与一之间。将式(1)的独立性改为给定X的条件独立,保留全体的选择单调性。写pdx=P(S=1∣D=d,X=x),假定p0x>0。每层得到[Lx,Ux],但总体目标不是按wx简单平均。

由全期望公式和Bayes公式,始终可见者中的层权重为

vx=P(X=x∣A)=wxp0x∑zwzp0z,

所以准确界为

(9)[∑xvxLx, ∑xvxUx].

每层可独立选择达到端点的完整律,再以wx混合并按各层既定概率分配处理,证明两端同时可达。若某层p0x=0,它在目标A中权重为零,应排除这一未定义条件效应;只需总分母为正。

若处理还独立于X,不使用X的总体Lee界也有效。加入全部层内限制只会删去相容律,因此条件化界不会更宽。若分配率本身随X改变,直接混合两臂不能继续声称无条件随机化,比较粗界前须重新检查这一前提。

有限分布的可执行算法 ​

对支持点y1<⋯<yk及概率fj,从低端累加质量,最后一个点只取剩余所需质量;从高端反向做一次。排序需O(klog⁡k)比较,排好序后各做O(k)次运算,输出保留质量向量(rj),并检查

0≤rj≤fj,∑jrj=a,μ=∑jyjrja.

这三个等式和不等式构成可独立核验的截尾证书。原子概率为有理数时可全程精确计算。将经验频率代入可得到估计的界,但它本身不是置信区间;估计选择率、分位点及边界时的抽样误差需要独立处理,本页不把人口层面的锐性冒称有限样本覆盖。

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

拖动节点调整位置。

显示关系

显示:依赖

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