r - 计算分数有时非常接近零的数值策略

Question

在 R 函数chisq.test()中有以下行：

PVAL <- (1 + sum(ss >= almost.1 * STATISTIC))/(B + 1)和

almost.1 <- 1 - 64 * .Machine$double.eps

这显然是一种计算调整，以避免为PVAL.

计算什么并不重要，但我们真正想要的是，一堆模拟的结果sum(ss >= STATISTIC)/ B在哪里，是一个要比较的固定值，是模拟的数量。我们正在计算大于ssSTATISTICBssSTATISTIC.

1分子和分母相加应该完成什么？

和

为什么我们需要乘以1 - 64 * .Machine$double.eps？