【万能求助】核酸检测误诊可能,概率论以及贝叶斯公式(原创)
所在版块:社会百科 发贴时间:2021-11-11 12:23  评分:

用户信息
复制本帖HTML代码
高亮: 今天贴 X 昨天贴 X 前天贴 X 
原创文章,理性讨论,不抬杠不撕逼。 

太长不看版:
1.新加坡眼app提到4例PCR疑似假阳,都来自来福士诊所。
2.贝叶斯公式证明如果新加坡普查, 一个人PCR阳性但真正得病的概率为0.44
3.如果怀疑自己PCR假阳,应该怎么办?


1. 4例PCR疑似假阳,都来自来福士诊所,血清全部呈现阴性。

有网友收到PCR阳性,当天马上去再测一次变阴。 阴变阳有可能是后来得的, 阳变阴逻辑上证明误诊可能,尤其是在这么短的时间内再测,有人还不止测了一次。4例都是来福士,应当引起注意了。


2. 贝叶斯公式以及误诊概率

新加坡总人口545万(2021年6月公开数据),新冠总确诊数227681(2021/11/11新浪数据),所以估计目前患病比例为0.04.  以事件A表示PCR呈现阳性,以A’表示PCR呈现阴性(事件A的补集).  PCR检测准确率数据不好找,假设一个人确实患病(设为事件C, P(C)=0.04),他做PCR呈现阳性的概率为P(A|C)=0.95.   假设一个人没有得病(C’,事件C的补集),他做PCR呈现阴性的概率为P(A’|C’)=0.95  (根据诊所报告写的confidence level 0.95).

其中P(A|C)=0.95, P(C)=0.04, P(A|C’)=1-P(A|C)=0.05,  P(C’)=1-P(C)=0.96

所以如果新加坡在所有人群中普查,某人得到诊所诊断PCR呈现阳性,他真正得病的概率(P(C|A)由贝叶斯公式

P(C|A)=P(A|C)*P(C)/[P(A|C)*P(C)+P(A|C’)*P(C’)]=0.44

在回国人群中,因为都是是无症状的,如果PCR呈现阳性,真正得病得概率仅为0.44(理论上还要小,因为这个包含了有症状的那些人).  这就是为什么中国发现PCR阳性要反复确认。 
新加坡这里只要PCR阳性,马上发隔离令,这是不是也是一种简单粗暴和医疗资源浪费呢?
(这个模型比较粗糙,值得商榷,欢迎提供资料继续讨论)


3. 如果怀疑自己PCR是假阳,应该怎么办? 

其实我是来问这个问题的,医学手段上有什么方法可以确定的证明自己是PCR假阳,因为现在一旦PCR阳性,大使馆就要求康复三个月才能回去中国,如果是假阳,太耽误事情了。


参考资料:
1.   新加坡眼app文章. 狮城茶馆->《疑似新冠误诊,近期回国注意了》 11.8
2. 《概率论与数理统计》浙江大学第三版 



该帖荣获当日十大第8,奖励楼主4分以及6华新币,时间:2021-11-11 22:00:22。
该帖荣获当日十大第8,奖励楼主4分以及6华新币,时间:2021-11-12 22:00:19。
.
欢迎来到华新中文网,踊跃发帖是支持我们的最好方法!

 相关帖子 我要回复↙ ↗回到正文
【万能求助】核酸检测误诊可能,概率论以及贝叶斯公式(原创) 漫漫长路   (2350 bytes , 1733reads )
楼主,我用你的公式算出了一个标准答案:95%. Blzrd   (215 bytes , 13reads )
你是比较靠谱的 漫漫长路   (175 bytes , 22reads )
你要算的是测过pcr的并且阳性的这些人确实患病的概率, Blzrd   (129 bytes , 11reads )
不是, 看原帖第二条, 漫漫长路   (122 bytes , 15reads )
楼主的概率是瞎算 目测是样本空间概率模糊所致 功夫熊猫   (8 bytes , 14reads )
这里出阳性报告不想你说的那么简单和粗暴 yayayayayaya   (214 bytes , 22reads )
moh_qsm@moh.gov.sg 联系了新加坡卫生部,说把详细和报告情况发到这里,他们会去调查 tiantang   (0 bytes , 21reads )
如果新加坡在所有人群中普查,某人得到诊所诊断PCR呈现阳性,他真正得病的概率0.44 Coye1214   (90 bytes , 21reads )
这就是概率论的美妙之处,反直觉。 漫漫长路   (0 bytes , 14reads )
楼主,你预设的几个概率没一个对的。 TPMT   (583 bytes , 28reads )
P(C)不是患病的概率,而是随机在新加坡抽选一人,他是冠病患者的概率。 漫漫长路   (79 bytes , 15reads )
PCR准确率是什么意思? Confidence level 0.95又指的是什么? ESPRITEE   (0 bytes , 11reads )
置信水平吧, 和confidence level一样 漫漫长路   (194 bytes , 14reads )
好好养病吧,别想太多 dancingdoll   (125 bytes , 11reads )
我就想问下 NoGrabTaxi   (35 bytes , 15reads )
条件概率的性质 漫漫长路   (0 bytes , 18reads )
还是 NoGrabTaxi   (16 bytes , 13reads )
我有朋友也是这种情况,但是大使馆不承认假阳性 tiantang   (251 bytes , 15reads )
多次pcr或十四天后pcr,抽血,血清也 Currant   (72 bytes , 12reads )
最有力的证明就是时间,目前的设置是3个月, Blzrd   (57 bytes , 13reads )
大使馆应该不会承担这个风险 tiantang   (48 bytes , 19reads )
这些人想回国和大使馆穷对付,是痴心妄想 熊走瞧   (0 bytes , 19reads )
呵呵 dancingdoll   (18 bytes , 14reads )
楼主是数学达人生物不通,pcr 阳性就是阳性, 熊走瞧   (159 bytes , 20reads )
概率 Currant   (0 bytes , 9reads )
概念和你说那个是两回事吧。 Currant   (0 bytes , 15reads )
楼主虽然知道贝叶斯公式但是似乎不懂生物和医学原理 dancingdoll   (502 bytes , 19reads )
第三点描述毫无因果啊 原因应该是art specificity很高所以假阳可能低所以不用做pcr了 功夫熊猫   (102 bytes , 20reads )
你没理解我的意思 dancingdoll   (126 bytes , 14reads )
PCR测的是抗原,也就是病毒的核酸序列。不是抗体。 TPMT   (53 bytes , 26reads )
第二条是针对楼主列出来的四个病例里面IgG, IgM测试阴性的 dancingdoll   (18 bytes , 13reads )
Sorry,可能我断句有问题。误读了楼主意思。 TPMT   (10 bytes , 26reads )
看到了 dancingdoll   (150 bytes , 14reads )
感谢科普 漫漫长路   (38 bytes , 10reads )
检测报告上应该有阳性信号是在第几个cycle出现的 dancingdoll   (165 bytes , 12reads )
结论 sakuragi   (31 bytes , 18reads )
忘记是哪里看的 sixsixsix   (475 bytes , 95reads )
请问本地人去什么医院 dancingdoll   (12 bytes , 17reads )