如何在R中处理大数字？

Question

我想改变R语言计算的精度，例如我想用x = c(-2.5e+59, -5.6e+60)计算x^6。为了计算它，我需要更改R语言的精度，否则结果将为Inf，但我不知道如何操作。请注意保留原文中的HTML标签。

- user3430764

---

这更多关于R而非RStudio——如果问题在于超越双精度，那么就有其他问题需要担心（例如R链接的库的精度）。这也可能是那些指数幂运算（通常由语言后端实现并调用exp()）与重复乘法不同的浮点运算之一，在接近精度边缘时，这可能会产生重要后果。 - Livius

这个链接 http://stat.ethz.ch/R-manual/R-patched/library/base/html/zMachine.html 可以帮助你了解你的架构限制，这可以成为寻找能够克服任何内置限制的软件包的好出发点。 - hrbrmstr

感谢您的帮助和花费时间回答。 - user3430764

Answer 1

正如Livius在评论中指出的那样，这是R（事实上，大多数编程语言）中用二进制表示数字的问题。

为了处理非常大/小的浮点数，可以使用Rmpfr库：

install.packages("Rmpfr")
library("Rmpfr")
x <- c(-2.5e+59, -5.6e+60)
y <- mpfr(x, 6)  # the second number is how many precision **bits** you want - NB: not decimal places!
y^6
# 2 'mpfr' numbers of precision  6   bits 
# [1] 2.50e356 3.14e364

如果要处理比R能够处理的更大的数字（例如exp(1800)），您可以使用“Brobdingnag”包：

install.packages("Brobdingnag")
library("Brobdingnag")

## An example of a single number too large for R: 
10^1000.7
# [1] Inf

## Now using the Brobdingnag package:
10^as.brob(1000.7)
# [1] +exp(2304.2)

- Scott Ritchie

---

在Numerical Mathematics Task View中有一个专门介绍“多精度算术和符号数学”的部分。 - Joshua Ulrich

我仍然无法在R中显示像exp(1800)这样的数字，对吗？ - Jakob

@PeterPan 我已更新我的答案，包括有关数字大于 R 可处理的范围的信息。 - Scott Ritchie