在R中得到0而不是一个精确的结果

Question

我怎样才能得到实际的精确结果而不是圆形的结果呢？

result = ((0.61)**(10435)*(0.39)**(6565))/((0.63)**(5023)*(0.60)**(5412)*(0.37)**(2977)*(0.40)**(3588))

退出：

NaN

因为分母是0

Answer 1

我认为对数是处理大幂指数的有力工具(参见https://mathworld.wolfram.com/Logarithm.html中的属性)。

您可以尝试首先在数学表达式上使用log，然后依次应用exp，即，

result <- exp((10435*log(0.61)+6565*log(0.39)) - (5023*log(0.63)+5412*log(0.60)+ 2977*log(0.37)+3588*log(0.40)))

这给

> result
[1] 0.001219116

Answer 2

R不能处理如此大的指数，因为它将收敛到超过其精度的0。精确并不是无限的。为了满足您的需要，您需要一个任意的精确包，比如Rmpfr。

library(Rmpfr)

precision <- 120

result <- (mpfr(0.61, precision)**10435 * mpfr(0.39, precision)**6565) /
  (mpfr(0.63, precision)**5023 * mpfr(0.60, precision)**5412 * mpfr(0.37, precision)**2977 * mpfr(0.40, precision)**3588)

print(result)

输出：

1 'mpfr' number of precision  120   bits 
[1] 0.0012191160601483692718001967190171336975