设万维读者为首页 万维读者网 -- 全球华人的精神家园 广告服务 联系我们 关于万维
 
首  页 新  闻 视  频 博  客 论  坛 分类广告 购  物
搜索>> 发表日志 控制面板 个人相册 给我留言
帮助 退出
老钱的博客  
我不思,则我不在  
https://blog.creaders.net/u/6860/ > 复制 > 收藏本页
网络日志正文
老钱:关于“样本池”的普测方法 2020-05-27 09:03:18

老钱:关于“样本池”的普测方法

5/25/20


今天网上流转两篇文章,武汉千万人核酸检测,真实感染率曝光全民检测50万人已感染 武汉女斥政府脑残(视频)
 
武汉市的感染情况到底如何。这个问题不仅仅是武汉人关心,政府关心,也是全世界都很关心的问题。5月14日,武汉市掀起了“十日大会战任务”,要在10天之内,对全武汉的一千一百万人做一次核酸检测,彻底摸清武汉市的详细情况。


这么大的测试量,怎么能在十天完成。据说这是因为采取了混合样本池的方法。这个方法是英国科学家在《柳叶刀/Lancet》上发表的,为了提高大规模普查的效率和降低成本。
 
我思考了一下,想来做一个进一步理解和解释。
 
设想把一群人(比如说64人中的一半,32人)的每个人的采样样本混合在一起,成为一个样本池。因为32人的样品都搅合在一起了,只要其中有一个人是阳性,这个样本池整体就是阳性。如果检测结果是阴性,那么这个样本池里的任何一个人都必定是阴性。那么就不用对这32人一个一个地检测(但是仍然需要做了32个采样),而是仅一次测试就把这32个人都排除了。减少了31个测试。
 
对于剩下的另一半的32人呢?如法炮制。
 
再取一次16人的混合样本池,如果最好的结果是阴性。这16个人就一下子排除了。
第三次对剩下一半16人,做8人的混合样本池,最好的结果是阴性。这八个人就一次排除了。
第四次是对剩下8人,进行4人的混合样本池,最好的结果是阴性。这四个人就一次排除了。
第五次是对剩下4人,2人的混合样本池,如果是阴性。这两个人就排除了。
第六次,仅剩下最后两个人了。就要对每一个人都做测试。结果可能是3种。
两人都是阴性,那么这全部的64人就没有一人被感染!
如果一阴一阳,就是只有感染一人了;感染率就是1/64,即1.56%。
如果都是阳性,就是只有两人感染:感染率就是1/32,即3.12%。
 
这样的样本池法,最好的结果就是,对全部64人,采了64次样,但是仅仅做了7个测试,就完成了对全部64人的测试。每个人都必须被采样,这是必须的。但是,不是每一个人都需要做测试。最好的情况,就是只做了7次测试,就查清了64人。7/64,即是挨个检测的11%的检测量。
 
对于昂贵的耗料耗时耗人工的核酸检测来说,样本池法是非常有效,非常经济合算的。极大地降低了测试成本!
 
这样的样本池法,也叫做“二分法”。
 
上面是最幸运的情况。如果不幸呢?我们来看看“最不幸”的情况。
 
如果这个32人的“样本池”的测试结果是阳性,就必须对这个“样本池”里的人重新来一次采样和测试。可以一个人一个人的测试,找到具体谁是阳性,谁是阴性。不仅仅必须要对这32个人重新测试,还要对另一半的32人从头开始。
 
也可以,对这个阳性的“样本池”再来一次16人的混合样本池,最坏的结果是阳性。不仅仅必须对这16个人重新测试,还要对另一半的16人从头开始。
第三次是8人的混合样本池,最坏的结果是阳性。不仅仅必须对这8个人重新测试,还要对另一半的8人从头开始。
第四次是4人的混合样本池,最坏的结果是阳性。不仅仅必须对这4个人重新测试,还要对另一半的4人从头开始。
第五次是2人的混合样本池,如果结果是阳性。那么这最后四个人的每一个人都要采样测试。
 
这样的样本池法,最坏的结果就是,对全部64人,每一个人都被重复了4次采样和4次检测。一共采了256次样,做了256个测试。最终结果,最坏是64人都被感染了,百分之百的感染。最坏之下的最好的结果是,每个4人里只有一个人是阳性,另外3人都是阴性。就是四分之一,即25%的感染率。但是,仍然是要采256次样,做256个测试。
 
这样,样本池法反而增加了工作量,还不如64人,每一个人都测试一次了。
 
所以,样本池法的运用是有条件的,需要随机应变,灵活运用的。
 
而且,样本池太大了,可能降低了检测的灵敏度。我前面用的32人的样本池就太大了。这样“夸张”的样本池法,只是为了使得“减少测试量到11%”可以显易而见。如果32人中只有一个阳性的话,阳性的浓度立刻就降低到了1/32,很可能就检测不出来了。10人的样本池就会把一个阳性的浓度降低到1/10。多大的样本池是合适的呢?这就不是纸上谈兵,而应该从试验中总结出来。从微信上武汉大嫂的抱怨来看,在武汉的“大会战”中,是采用了10人的样本池。那个大嫂抱怨:“把我们十个人的样品不注名地放在一个容器里”。大嫂说,“那查个鬼啊!”。这是因为她不知道到样本池”的理论。
 
那么,我们就用10人“样本池”来分析。如果这个“样本池”是阴性,那么这10人就全部是阴性,就节约了9/10的测试量。如果是阳性呢,就要把这10人叫回来重新测试。
改为5人“样本池”,如果结果是阴性,这5个人都是阴性,不用再测试了。
对于其余5个人,一一测试,总测试量就是7次,7/10,也减少了3次。
如果再采取3人“样本池”,如果是阴性,这3人都是阴性,只要再分别测试剩下的两个人就完成了。这样的总测试量就是5次,减少了50%的工作量。
如果3人“样本池”是阳性,那么这5个人都要逐个测试。总工作量就是8次。8/10.
如果第一个5人“样本池”是阳性,就重复上面的分析,总测试量大约为50%到80%。因为他们有把握,剩下的5人应该全是阴性。因为,在决定采用“样本池”方法之前,他们已经做了前期试验,得到了平均感染率是5%。但是,为了避免漏侧,还是要对剩下的5人,做一个“样本池”测试。这样来,总测试量大约人人检测的为60%到90%。
 
到这里,我们要回过头来说,“武汉大嫂”的抱怨是完全正确的,即使采用样本池”方法,不可能把受测试的老百姓,反复的叫回来重新采样,更要被老百姓骂死了!一次采样的量必须够用测试3次。
 
所以,那个大嫂抱怨:“把我们十个人的样品不注名地放在一个容器里,那查个鬼啊!”。她是完全正确的。不管她懂不懂样本池”。
 
由此看来,样本池”的方法也只能用于第一步,顶多到第二步,到了4人“样本池”,就没有什么优越性了。
 
对于美国斯坦福大学研发出来的抗体快速检测来说,测试与采样基本上就是一件事儿。除了昂贵复杂的核酸检测,这样的混合样本池法,也就没有多少优越性了。
 
我不是生化医学领域里的人,只是一个工程师,只是有兴趣从简单算术来理解样本池”的理论。本文不谈如何政治问题,只是数学和逻辑。可能贻笑大方。欢迎专家批评指正。
 
 (全文到此结束)


浏览(755) (0) 评论(3)
发表评论
文章评论
作者:西岸 回复 gugeren 留言时间:2020-05-27 12:15:06

这是关于检测效率的问题,与检测精度没有关系,再精确的方式也存在检测效率问题。再说这是普测,目的不是针对治疗隔离,而是发现无症状阳性的比例。

从LA不久前对抗体的检测结果来看,美国实际被感染的人应该是目前被发现的数量的20-50倍,这显然是引起警报的,因为这意味很多人是无症状或者轻症,但都具备传染他人的能力。武汉的这种普测就能相对精确地知道被传染的比例。

回复 | 0
作者:西岸 留言时间:2020-05-27 12:03:03

这种做法并非什么新鲜的东西,刚来美国的时候给一个在海军实验室做地形模型的刚毕业的博士讲过适合他的项目的数据取样优化,设计模型就简化很多,就是这类方式,鬼子似乎不那么理解,挺奇怪的,这类方法中国大概工科生多数人都知道,哪怕不会应用也应该知道。

不过当时兄弟的英文也是不到能上课的水平,不过就是帮接待我们的美国朋友的忙,他是海军那厮的教授。

具体在关于这个病毒的检测使用这种方式,并非没有前提。第一就是假设呈阳性的人是少数,如果超过一半就没多大意义。

第二是假设病毒传染力很强,放在一起的样本会很容易发生交叉感染。

那么看看这次的病毒是否具有这两个特点就是了,明显是同时具备这两个特点的。

文科生啊。。。怎么夸呢?

回复 | 0
作者:gugeren 留言时间:2020-05-27 10:13:21

这种“二分法”在理论上是对的。

但是,它的基础是每次的测试都是正确的。就是说,既不会出现“假阳性”、也不会出现“假阴性”这两种情况;即使有,两者的比例也不能很高。

很可惜,中国测试盒的质量太烂了。中国出口或支援别国的测试盒被人家退货了。

回复 | 0
我的名片
老钱
注册日期: 2012-10-17
访问总量: 1,094,569 次
点击查看我的个人资料
Calendar
最新发布
· 老钱ZT:《历史的宫外孕》序言
· 老钱:关于Roof/屋顶
· 老钱:共产主义“打家劫舍”进行时
· 老钱:再谈“民粹主义”
· 老钱ZT:咋办?
· 老钱:这就是侵略战争!
· 老钱:再说两句朱令案
友好链接
· 润涛阎:润涛阎的博客
· 桑妮:桑妮
· 俺是吾丁:俺是吾丁的博客
· Jack之家:Jack之家
分类目录
【园艺】
· 老钱:关于Roof/屋顶
· 老钱:菊花2021
· 老钱的菊花
【大自然】
【世界】
· 老钱:袁红冰谈xi猪头
· 老钱:狼要吃羊,需要什么理由!
· 老钱:再谈是撤退还是出卖
· 老钱:南海问题(一)/从南京宝
· 老钱:再谈对新冠病毒的正确对策
· 老钱:祝福阿桑奇幸福健康
· 老钱:再谈佛系对策
· 老钱:有感于日本中小学校餐
· 老钱:我看方方日记(2)
· 老钱:病毒的佛系对策图解
【大陆】
· 老钱:这完全是人祸,这就是屠杀
· 老钱:竟无一人是男儿
· 老钱ZT:不能忘记的罪恶(1)
· 老钱:到底是谁杀害的中国人最多
· 老钱ZT: 谁需要敌人 谁需要战争
· 老钱:钟南山火线入帮? 闹剧,丑
【为美国忧】
· 老钱:再谈“民粹主义”
· 老钱:这就是侵略战争!
· 老钱:华人应该怎样自我定位?
· 老钱:什么人在反对这些提案?
· 老钱:谁被禁止购买美国房地产?
· 老钱:是“保护华人”,还是抗美援
· 老钱:这是美国的存亡之战!
· 老钱ZT方伟: 解密共运之路 从法
· 老钱ZT方伟:解密共运之路 从费
· 老钱ZT方伟: 揭开惊天秘密与美国
【永恒的六四】
· 老钱ZT: “狼心狗肺的东西”
· 老钱:林郑月娥让步了吗?
· 什么叫“境外势力”?
· 美丽的勇敢,勇敢的美丽 - 由6。
· “六.四”万岁!年轻人万岁!
【故里行】
· 《我们这一辈》
· “拆!拆!拆!”/故里行(七)
· 看财经频道/故里行(六)
· 大陆的两极/故里行(五)
· 火车上/故里行(四)
· 位虽卑,未敢忘国/故里行(三)
· 狗屎GDP/故里行(二)
· 南京一条小街/故里行(一)
【华人参政】
· 老钱:Michelle Au不是美国华人
· 老钱:极左派挑动制造族裔矛盾的
· 老钱:糊涂的华左朋友们应该醒醒
· 老钱:Ossoff竞选就是作弊!
· 老钱:答友人,这个疯狂的世界
· 老钱:与杨博士商榷
· 老钱:这就是民主党极左派所谓的
· 老钱:“全部选票”就是坚持“造假”
· 老钱:与杨贵葆先生商榷
· 老钱Z:一个清华同学的问题
【华人故事】
· 老钱:为什么有这么多糊涂的海外
· 老钱:疑点重重的一个自杀案
· 老钱: 谈谈我和华夏文摘
· 老钱:华人的困境和认知误区
· 老钱:她为苍生说过话
· 老钱:一个沥青厂,将矗立在我们
· 老钱:为傅雷收藏骨灰的江小燕的
· 老钱:纪念叶晓芸
· 老钱:我才遇到的一个网络诈骗
· 我们是史无前例的77级!
【在北美开车】
· 关于交通规则的差异
· 再说我的车库门
· 说说我的车库门
· 一次车祸(续)
· 一次车祸
【培养男子汉】
· 培养男子汉(四)天生我才必有用
· 培养男子汉(三)好学向上
· 培养男子汉(二)重视体育
· 培养男子汉(一)人格人品
【读书随笔】
· 老钱: 怎么离开这个世界?
· 读《松园旧事》- 中国二十世纪的
【电影评论】
· 老钱:我看《都挺好》
· 我看《芳华》
· “为‘影响’服务”-《我不是潘金莲
· 《武媚娘传奇》的观感
· 《山楂树之恋》
· 《铁娘子》获得了奥斯卡奖
【晚霞的凄凉】
· 晚霞的凄凉2,这是一个国家的,
· 晚霞的凄凉1,关于Alzheimer/老年
【谁的错】
· 老钱:为了父亲节
· 美国的问题(二)/ 关于平等的一
· 华裔不相信眼泪!
· 谁的错?2--浅谈税收与富人
· 谁的错?2--也说说华尔街与金融
· 谁的错?---浅谈全球化
【科学技术】
· 云计算的不久将来
· 一次云计算的讲座
· “立体打印机”,让我们睁大眼睛,
· 让我们睁大眼睛,高举双手。。。
· “好奇心”是推动历史的动力
【旅游观光】
· 老钱:莫斯科的地铁
· 老钱:列宾的画展-莫斯科行5
· 老钱:莫斯科的航空航天博物馆
· 老钱:莫斯科芭蕾舞大剧院
· 老钱:莫斯科的一个艺术博物馆
· 老钱:华灯初上莫斯科
· 日本印象(二)
· 日本印象(一)
· 在欧洲de走马观花及思索
· 在韩国坐马观花
【健康保健】
· 老钱:腮腺炎的故事
· 老钱:我们回到家了/自我隔离中
【动物世界】
· 濒临灭绝的鲨鱼
· 青蛙的故事(二)
· 青蛙的故事(一)
【环境保护】
· 老钱:沥青厂的再思考
· 大洋旋涡中心的垃圾场
【孩子教育】
· 浅议“虎妈”
【大选2012】
· 美国大选(7),选谁? ,关于欧
· 美国大选,选谁?as the same as
· 美国大选(5)选谁? 不能再选一只
· 美国大选(4)选谁? 不搞“阶级斗
· 美国大选(3)选谁?第一次总统竞
· 美国大选(2)选谁?关于平等的思
· 美国大选6选谁?有感于第二次总
· 美国大选1选谁?不要搞“社会主义
【布衣闲谈】
· 老钱ZT:《历史的宫外孕》序言
· 老钱:共产主义“打家劫舍”进行时
· 老钱ZT:咋办?
· 老钱:再说两句朱令案
· 老钱:我也说两句朱令案
· 老钱:不带他玩了,可以吧
· 老钱:美国乱象之根源
· 老钱:“不须放屁”
· 老钱:袁红冰谈xi猪头
· 老钱:李克强功过是非
存档目录
2024-04-14 - 2024-04-14
2024-03-02 - 2024-03-30
2024-02-03 - 2024-02-03
2024-01-05 - 2024-01-09
2023-12-07 - 2023-12-07
2023-11-12 - 2023-11-29
2023-10-12 - 2023-10-29
2023-09-05 - 2023-09-14
2023-08-06 - 2023-08-06
2023-07-16 - 2023-07-16
2023-06-01 - 2023-06-23
2023-01-17 - 2023-01-17
2022-12-31 - 2022-12-31
2022-11-05 - 2022-11-05
2022-10-29 - 2022-10-29
2022-09-06 - 2022-09-23
2022-08-10 - 2022-08-27
2022-07-25 - 2022-07-28
2022-06-07 - 2022-06-20
2022-05-10 - 2022-05-25
2022-04-08 - 2022-04-22
2022-03-02 - 2022-03-26
2022-01-15 - 2022-01-30
2021-12-05 - 2021-12-27
2021-11-12 - 2021-11-13
2021-09-03 - 2021-09-29
2021-08-19 - 2021-08-28
2021-07-10 - 2021-07-10
2021-06-09 - 2021-06-24
2021-04-09 - 2021-04-09
2021-03-13 - 2021-03-29
2021-02-13 - 2021-02-24
2021-01-03 - 2021-01-28
2020-12-08 - 2020-12-18
2020-11-02 - 2020-11-30
2020-10-07 - 2020-10-26
2020-09-01 - 2020-09-30
2020-08-05 - 2020-08-15
2020-07-14 - 2020-07-28
2020-06-12 - 2020-06-28
2020-05-05 - 2020-05-28
2020-04-16 - 2020-04-28
2020-03-02 - 2020-03-23
2020-02-06 - 2020-02-19
2020-01-12 - 2020-01-13
2019-12-03 - 2019-12-30
2019-11-03 - 2019-11-07
2019-10-10 - 2019-10-29
2019-09-01 - 2019-09-30
2019-08-02 - 2019-08-02
2019-07-24 - 2019-07-24
2019-06-17 - 2019-06-24
2019-05-04 - 2019-05-11
2019-04-08 - 2019-04-08
2019-03-11 - 2019-03-11
2019-02-20 - 2019-02-27
2019-01-07 - 2019-01-29
2018-12-05 - 2018-12-05
2018-11-22 - 2018-11-22
2018-10-01 - 2018-10-31
2018-09-19 - 2018-09-26
2018-07-21 - 2018-07-29
2018-06-13 - 2018-06-13
2018-05-01 - 2018-05-25
2018-03-26 - 2018-03-26
2018-02-09 - 2018-02-09
2018-01-02 - 2018-01-04
2017-12-16 - 2017-12-26
2017-11-07 - 2017-11-30
2017-10-06 - 2017-10-31
2017-09-11 - 2017-09-11
2017-08-15 - 2017-08-15
2017-06-23 - 2017-06-27
2017-05-29 - 2017-05-29
2017-04-04 - 2017-04-04
2017-02-08 - 2017-02-21
2017-01-04 - 2017-01-12
2016-12-03 - 2016-12-22
2016-11-01 - 2016-11-19
2016-10-04 - 2016-10-28
2016-09-23 - 2016-09-25
2016-08-07 - 2016-08-24
2016-07-08 - 2016-07-12
2016-04-07 - 2016-04-07
2016-03-01 - 2016-03-13
2016-02-17 - 2016-02-17
2015-12-02 - 2015-12-21
2015-11-03 - 2015-11-26
2015-10-03 - 2015-10-30
2015-09-30 - 2015-09-30
2015-08-09 - 2015-08-20
2015-06-17 - 2015-06-23
2015-04-02 - 2015-04-29
2015-03-05 - 2015-03-05
2014-10-07 - 2014-10-07
2014-09-18 - 2014-09-18
2014-08-01 - 2014-08-21
2014-07-29 - 2014-07-29
2014-06-08 - 2014-06-16
2014-05-06 - 2014-05-29
2014-04-08 - 2014-04-29
2014-03-11 - 2014-03-30
2014-02-04 - 2014-02-21
2014-01-06 - 2014-01-23
2013-12-10 - 2013-12-31
2013-11-04 - 2013-11-30
2013-10-05 - 2013-10-13
2013-09-02 - 2013-09-29
2013-07-06 - 2013-07-23
2013-05-26 - 2013-05-26
2012-12-01 - 2012-12-04
2012-11-02 - 2012-11-28
2012-10-18 - 2012-10-31
 
关于本站 | 广告服务 | 联系我们 | 招聘信息 | 网站导航 | 隐私保护
Copyright (C) 1998-2024. CyberMedia Network /Creaders.NET. All Rights Reserved.