Python Pandas - 与 Scipy 的优势比(P 值 = 0？)

Question

我正在尝试计算 2 x 2 网格的优势比(如下)

    Lower       Higher
A   11772336    18837138
B   3624890     4263509

我已经完成了这段代码来获取优势比:

import scipy.stats as stats
table = df.values
reversedTable = table[::-1]

oddsratio, pvalue = stats.fisher_exact(reversedTable)
print("OddsR: ", oddsratio, "p-Value:", pvalue)

输出为:

OddsR:  1.3604416468907712 p-Value: 0.0

如果我使用在线计算器，我得到的 P 值为 <0.0001

有没有办法让我的代码中的 P 值不被舍入？还是我做错了什么？

非常感谢

Answer 1

在您的例子中，对数优势比的 z 分数为 382.246，相关的双尾 p 值为 3.169 × 10^-31731。这小于可以用 64 位 float 表示的最小正值 (2⁻¹⁰⁷⁴)，这是 fisher_exact 返回的 p 值。