r - 将浮点精度添加到qnorm/pnorm吗？

Question

我将有兴趣从当前级别计算qnorm/pnorm时增加浮点数限制，例如:

x <- pnorm(10) # 1
qnorm(x) # Inf
qnorm(.9999999999999999444) # The highst limit I've found that still return a <<Inf number

在合理的时间内可以做到吗？如果是这样，怎么办？

Answer 1

如果参数位于右上角，则应该可以通过计算1-p来获得更好的精度。像这样:

> x = pnorm(10, lower.tail=F)
> qnorm(x, lower.tail=F)
10

我希望(尽管我不确定)，pnorm()函数引用的是C或Fortran例程，该例程停留在硬件支持的任何浮点大小上。重新安排您的问题可能更好，因此不需要精度。

然后，如果要处理非常大的z值，则可以使用log.p = T:

> qnorm(pnorm(100, low=F, log=T), low=F, log=T)
100

抱歉，这不是您要找的东西。但我认为它将具有更高的可扩展性-pnorm在高z值(毕竟为e ^(-x ^ 2))下如此快地命中1，即使添加更多位，它们也会很快用完。