破产问题与简单序贯检测

本文是对Simple Sequential A/B Testing的解读。
该方法归属序贯检测类,可以用于伯努利分布场景,随着抽样持续进行,判断接受零假设或备择假设(关于序贯检测)。

理论依据——赌徒破产问题

一个赌徒携带d元进入赌场,每次下注一元,胜利额外赢得一元,失败则输掉一元。假如赌博是公平的,赌徒在第n轮剩余0元的概率(记作R_{n,d})是多少?

概率计算

使用R_{n,d}来表示赌徒起始d元在第n轮输光的概率。
总局:n,输:(n + d) / 2,胜:(n - d) / 2。

  1. 赌场可以借钱给赌徒的情况,此时为简单的排列组合关系。

R_{n,d} = \binom{n}{(n + d)/ 2} * 2^{-n}

  1. 赌场不肯借钱给赌徒情况。

R_{n,d} = \frac{d}{n}\binom{n}{(n + d)/ 2} * 2^{-n}

这是简单序贯检测的基础,将在下面推导。

随机游走和选举定理

在此对第二种情况进行计算。

  1. 取纵坐标为钱,横坐标为赌博轮次,则过程为一维随机游走。举例:


    赌徒场景2
  2. 此时从A出发到K,满足条件路径数(中途不与0轴接触),等于从K出发到A的路径数。场景转换为从0开始,到d结束。


    转换
  3. 在满足中途不触碰0轴条件,第一步必须为正。则从J到A路径数等于从A到K路径数。


    镜像

    如图所示,每一条从J到第一个接触横轴的点的路径,总有一条从K出发到相同点的映射路径,后续都走相同路径。因此J到A经过横轴路径数等于K到A的路径数。

  4. 设从0出发,有p步向上,q步向下。根据3,满足不触碰横轴概率为:
    (\binom{p + q - 1} {p - 1} - \binom{ p + q - 1}{p}) / \binom{p + q}{p} = \frac{p - q}{p + q}

此公式被称为选举定理

p + q = n, p - q = d带入,原问题最终结果为\frac{d}{n}\binom{n}{(n + d)/ 2} * 2^ {-n}

简单序贯检测

简单序贯检测是基于上面的结论产生的。

与赌徒破产问题的关系

均等流量的转化率型ab测试场景中,实验组、对照组随着时间推移,都会产生转化。在零假设下,下一次转化发生在实验组或对照组的概率是相等的。因此可以转换为赌徒破产问题。

每轮次:下一次转化发生;
每轮获胜者:下一次转化所在的组;
破产的轮次n:实验组、对照组转换数之和;
赌徒初始筹码d:因为转换数积累是从零开始的,不能直接套用。可以认为赌徒初始资本为0,输到-d时破产。

假设检验设计

以下仅介绍单尾情况,双尾的扩展请参考原文。

原假设

零假设(H0):实验组转化率等于对照组转化率(对应未破产情况);
备择假设(H1):实验组转化率小于对照组转换率(对应破产情况)。

假阳性控制(alpha)

在零假设假设下,第n轮赌徒破产概率为: R_{n,d} = \frac{d}{n}\binom{n}{(n + d)/ 2} * 2^{-n}
在N轮及N轮之前,赌徒破产概率为:\sum_{n = 1} ^{N}R_{n,d}
利用此公式控制假阳性,通过控N和d的选取,可以控制假阳性水平。

假阴性控制(beta)

控制假阴性,需要预设期望的最小观测效果(mde,此处选相对效果)。当实际提升效果等于MDE时,假阴性概率等于预设值。
在备择假设下,当实验组提升了mde时,下一次转化发生在实验概率为P_t = \frac {1 + mde} {2 + mde},对照组概率为P_c = \frac{1}{2 + mde}
第n轮破产概率:R_{n,d} = \frac{d}{n}\binom{n}{(n + d)/ 2} * P_c ^ {(n - d) / 2}* P_t^{(n + d) / 2}
在N轮及N轮之前,赌徒破产概率为:\sum_{n = 1} ^{N}R_{n,d}
取上面的概率为功效(power),通过次N和d的选取,可以控制假阴性水平。

选择结束条件

联立上述不等式,使N和d满足:

\sum_{n = 1}^N \frac{d}{n}\binom{n}{(n + d)/ 2} * 2^{-n} < \alpha

\sum_{n = 1}^N \frac{d}{n}\binom{n}{(n + d)/ 2} * P_c^{(n - d) / 2}* P_t^{(n + d) / 2} > 1 - \beta

不等式很难直接求解,可通过计算机遍历可能的N和d,找到合适的值。

具体流程
  1. 实验开始,实验组、对照组从0开始计数;
  2. 每有一个转化,对应的组计数+1,并进行判断;
  3. 对照转化数 - 实验组转化数 >= d,接受备择假设;
  4. 对照转化数 + 实验组转化数 = N,接受原假设。

优缺点

优点

  1. 基于弱假设,易于理解和证明;
  2. 过程易于操作;
  3. 在低转换率时,需要样本量小于固定水平检验需要的样本量;
  4. 不存在“偷看”问题。

缺点

  1. 只适用于近似伯努利分布场景;
  2. 结束条件难以直接计算,需要通过计算机遍历查找;
  3. 在高转化率时,需求样本量大于固定水平检验需要的样本量;
  4. 无法直接给出置信区间和P值。
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,588评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,456评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,146评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,387评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,481评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,510评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,522评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,296评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,745评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,039评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,202评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,901评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,538评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,165评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,415评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,081评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,085评论 2 352

推荐阅读更多精彩内容