A/B实验

A/B实验本质上就是数理统计中假设检验的工程实践,所以我们从假设检验开始引入。

假设检验

假设检验基本概念

  • 原假设(null hypothesis)H0H_0:是实验者想要收集证据予以反对的假设。A/B实验中的原假设就是指“新策略没有效果”。
  • 备择假设(alternative hypothesis)H1H_1:是实验者想要收集证据予以支持的假设,与原假设互斥。A/B实验中的备择假设就是指“新策略有效果”。

利用反证法来检验假设,意味着我们要利用现有的数据,通过一系列方法证明原假设是错误的(伪),并借此证明备择假设是正确的(真)。这一套方法在统计学上被称作原假设显著性检验(NHST:null hypothesis significance testing)

第一类错误和显著性水平( α\alpha )

第一类错误指原假设实际上是正确的(真),但是我们假设检验的结论却显示”原假设错误,备择假设正确“。这一过程中我们拒绝了正确的原假设,所以第一类错误是“弃真”。

在实际操作中表现为:实验结论显示我的新策略有用,但实际上我的新策略没有用。

在统计学中,我们用显著性水平( α\alpha )来描述实验者犯第一类错误的概率。

统计显著性:统计显著性 = 1 - 第一类错误的概率( α\alpha ),也称“置信水平、置信度、置信系数”,它的存在是为了描述实验结果的可信度。当某个实验组的指标是显著的,说明这个实验结果大概率是可信的。这个概率一般设为95%,也就是说,系统有95%的信心确认这个实验结果是准确的。

p-value:P值就是在原假设为真的前提下随机抽取样本出现极端情况的概率。当p-value < 1 - 置信度水平,认为统计显著。

如果P值很小,说明这种情况发生的概率很小,但如果出现了,根据小概率原理,我们就有理由拒绝原假设。P值越小,说明实验发现的差异是因为抽样误差导致的概率越小,极大程度上还是由于本质上存在差异造成,我们拒绝原假设的理由越充分。

第二类错误( β\beta )和统计功效(statistics power)

第二类错误指原假设错误(伪),但是我们假设检验的结论却显示“原假设正确(真)、备择假设是错误的”,这一过程中我们接受了错误的原假设,所以第二类错误是“取伪”。

第二类错误在实际操作中表现为:我的新策略其实有效,但实验没能检测出来。

统计功效:在统计学中,统计功效 = 1 - 第二类错误的概率( β\beta ),统计功效在现实中表现为:我的新策略是有效的,我有多大概率在实验中检测出来。

TypeIAndIIError

在上面的图中,浅灰色区域就是我们犯第二类错误的域,其面积占右侧钟形图形面积的比率,就是我们犯第二类错误的概率。这一概率被记作 β\beta 。在图中的例子里, β\beta 的取值为47%(计算过程就不展示了)。也就是说,在上面的例子里,即便我的新策略B有效,我仍有47%的概率没法检验出它是有效的。而统计效力=1β=53%=1- \beta =53\%。这代表着我有53%的概率可以检测出新策略B是有效的。

接受原假设 拒绝原假设
原假设为真 正确:1α1-\alpha 第一类错误Type I(弃真):α\alpha
备择假设为真 第二类错误Type II(取伪):β\beta 正确:1β1-\beta

我们希望犯两类错误的概率都尽可能的小,但它们是矛盾的。假设检验的通常做法是,先保证第一类错误(“弃真”)的概率 α\alpha 有较小的值,在这个前提下,使犯第二类错误(”取伪“)的概率尽可能小。优先保证”弃真“的概率较小,容忍可能的”取伪“,在AB Test中可以理解为在没有充分的把握说明实验组的新策略一定好的情况下,相对保守一点使用对照组的原策略。

若要同时减少犯两类错误的概率,则必须增加样本量。

参数的假设检验

单个正态总体均值 μ\mu 的检验

X1,X2,,XnX_1,X_2,\cdots,X_n 是正态总体 N(μ,σ2)N(\mu,\sigma^2) 的样本,我们讨论关于均值 μ\mu 的假设检验问题。在应用上常见的是对如下形式的假设进行检验:

  1. H0:μ=μ0,H1:μμ0H_0:\mu=\mu_0,\quad H_1:\mu \neq \mu_0
  2. H0:μμ0,H1:μ>μ0H_0:\mu \leq \mu_0,\quad H_1:\mu > \mu_0
  3. H0:μμ0,H1:μ<μ0H_0:\mu \geq \mu_0,\quad H_1:\mu < \mu_0

其中 μ0\mu_0 是一个给定的数。

先考虑双边检验问题。以 Xˉ\bar{X} 估计 μ\mu ,考察 Xˉμ0\bar{X}-\mu_0 的情况。直观上 Xˉμ0|\bar{X}-\mu_0| 的值越小,与原假设越符合;反之, Xˉμ0|\bar{X}-\mu_0| 的值越大,与备择假设 H1H_1 越符合。因此,H0H_0 的拒绝域是 Xˉμ0|\bar{X}-\mu_0| 大于一个界限值的区域。

  1. 总体方差 σ2\sigma^2 已知时

    原假设 H0H_0 成立时,检验统计量 U=Xˉμ0σ/nN(0,1)U=\frac{\bar{X}-\mu_0}{\sigma/\sqrt{n}} \sim N(0,1) ,若给定显著性水平 α\alpha ,则有: P{Xˉμ0σ/n>uα2}=α P\Big\{ \frac{|\bar{X}-\mu_0|}{\sigma/\sqrt{n}} > u_\frac{\alpha}{2} \Big\}=\alpha 得到原假设 H0H_0 的拒绝域为 U>uα2|U|>u_{\frac{\alpha}{2}},或者: Xˉμ0>σnuα2 |\bar{X}-\mu_0|>\frac{\sigma}{\sqrt{n}}u_{\frac{\alpha}{2}} 在这个检验中,我们使用了统计量U=Xˉμ0σ/nU=\frac{\bar{X}-\mu_0}{\sigma/\sqrt{n}} ,其服从的分布为 N(0,1)N(0,1) 。若一检验方法使用了服从标准正态分布的统计量,则称为U检验或者称为Z检验

  2. 总体方差 σ2\sigma^2 未知时

    由于 σ2\sigma^2 未知,所以用样本标准差 SS 代替 σ\sigma 。在原假设 H0H_0 成立时检验统计量 T=Xˉμ0S/nt(n1)T=\frac{\bar{X}-\mu_0}{S/\sqrt{n}} \sim t(n-1) ,故有: P{Xˉμ0S/n>tα2(n1)}=α P\Big\{ \frac{|\bar{X}-\mu_0|}{S/\sqrt{n}} > t_\frac{\alpha}{2}(n-1) \Big\}=\alpha 得到原假设 H0H_0 的拒绝域为 T>tα2(n1)|T|>t_{\frac{\alpha}{2}}(n-1),或者: Xˉμ0>Sntα2(n1) |\bar{X}-\mu_0|>\frac{S}{\sqrt{n}}t_{\frac{\alpha}{2}}(n-1) 在上述检验中,我们使用了统计量T=Xˉμ0S/nT=\frac{\bar{X}-\mu_0}{S/\sqrt{n}} ,其服从的分布为 t(n1)t(n-1) 。若一检验方法使用了服从 tt 分布的统计量,则称为t检验

现在来考虑单边检验问题。

对于:H0:μμ0,H1:μ>μ0H_0:\mu \leq \mu_0,\quad H_1:\mu > \mu_0

仍以 Xˉ\bar{X} 估计 μ\mu ,考察 Xˉμ0\bar{X}-\mu_0 的值大于一个正界限值的区域。总体方差 σ2\sigma^2 已知时有: P{Xˉμ0σ/n>uα}=α P\Big\{ \frac{\bar{X}-\mu_0}{\sigma/\sqrt{n}} > u_\alpha \Big\}=\alpha 得到拒绝域为 U>uαU>u_\alpha 或者: Xˉμ0>σnuα \bar{X}-\mu_0>\frac{\sigma}{\sqrt{n}}u_\alpha 总体方差 σ2\sigma^2 未知时有: P{Xˉμ0S/n>tα(n1)}=α P\Big\{ \frac{\bar{X}-\mu_0}{S/\sqrt{n}} > t_\alpha(n-1) \Big\}=\alpha 得到拒绝域为 T>tα(n1)T>t_\alpha(n-1) 或: Xˉμ0>Sntα(n1) \bar{X}-\mu_0>\frac{S}{\sqrt{n}}t_\alpha(n-1) 对于:H0:μμ0,H1:μ<μ0H_0:\mu \geq \mu_0,\quad H_1:\mu < \mu_0

类似的,总体方差 σ2\sigma^2 已知时,拒绝域为: U<uαU<-u_\alpha 。总体方差 σ2\sigma^2 未知时,拒绝域为: T<tα(n1)T<-t_\alpha(n-1)

大样本非正态总体均值 μ\mu 的检验

对于大样本的非正态总体而言,根据中心极限定理,当样本量 nn 足够大时,其均值 μ\mu 仍然服从正态分布。如果总体方差 σ2\sigma^2 未知,可以用样本标准差 SS 代替 σ\sigma 。于是对于双边检验问题有: P{Xˉμ0S/n>uα2}=α P\Big\{ \frac{|\bar{X}-\mu_0|}{S/\sqrt{n}} > u_\frac{\alpha}{2} \Big\}=\alpha 得到原假设 H0H_0 的拒绝域为 U>uα2|U|>u_{\frac{\alpha}{2}},即: Xˉμ0>Snuα2 |\bar{X}-\mu_0|>\frac{S}{\sqrt{n}}u_{\frac{\alpha}{2}} 单边检验与之前类似,这里就不再赘述。

在实际AB实验的应用中,通常样本量都会是大于30的,符合大样本的条件。所以如果做均值的AB实验,就可以无脑选择U检验或者说Z检验,而无需考虑t检验,也不用在意总体的分布情况,只需要用样本标准差 SS 代替 σ\sigma 。t分布在 nn 很大的时候本身就会无限趋近于正态分布。

[!NOTE]

关于U检验和Z检验。国内教材基本都叫U检验,国外教材基本都叫Z检验。虽然U检验和Z检验是一回事儿,但是个人理解觉得,对于正态总体均值 μ\mu 的检验更适合叫做U检验,而对于大样本总体均值 μ\mu 的检验,更适合叫做Z检验。在实际AB实验中,样本量基本不会小于30,所以对于均值 μ\mu 检验的问题,无需考虑总体是什么分布,以及总体方差 σ2\sigma^2 是否已知的情况,通通采用Z检验即可。

在实际AB实验中,不仅会碰到均值 μ\mu 检验的问题,同时也经常碰到比率值或者说转化率类的检验问题。这类问题的总体实际上是服从二项分布的。在大样本的情况下同样可以转化为Z检验的问题。

设有 nn 重伯努利实验,其中事件 AA 发生的概率为 pp ,则事件 AA 发生次数的随机变量 XX 服从二项分布,即 XB(n,p)X \sim B(n,p) 。对于参数 pp 双边假设检验问题描述为: H0:p=p0,H1:pp0 H_0:p=p_0,\quad H_1:p \neq p_0 对于二项分布,其数学期望 E(X)=npE(X)=np ,方差 D(X)=np(1p)D(X)=np(1-p)

nn 很大时,根据中心极限定理,有: Z=XE(X)D(X)=np^np0np0(1p0)=p^p0p0(1p0)/nN(0,1) Z=\frac{X-E(X)}{D(X)}=\frac{n\hat{p}-np_0}{\sqrt{np_0(1-p_0)}}=\frac{\hat{p}-p_0}{\sqrt{p_0(1-p_0)/n}} \sim N(0,1) 若给定显著性水平为 α\alpha ,则有: p{p^p0p0(1p0)/n>zα2}=α p\Big\{ \frac{|\hat{p}-p_0|}{\sqrt{p_0(1-p_0)/n}}>z_{\frac{\alpha}{2}} \Big\} = \alpha 得到原假设 H0H_0 的拒绝域为 Z>zα2|Z|>z_{\frac{\alpha}{2}},或者: p^p0>p0(1p0)/nzα2 |\hat{p}-p_0|>\sqrt{p_0(1-p_0)/n} \cdot z_{\frac{\alpha}{2}} 单边检验与之前类似,这里就不再赘述。

results matching ""

    No results matching ""