r - 控制R中打印输出的小数位数

Question

R 中有一个选项可以控制数字显示。例如:

options(digits=10)

应该给出 10 位数字的计算结果，直到 R session 结束。在R的帮助文件中，digits参数的定义如下:

digits: controls the number of digits to print when printing numeric values. It is a suggestion only. Valid values are 1...22 with default 7

所以，它说这只是一个建议。如果我希望始终显示 10 位数字而不是更多或更少怎么办？

我的第二个问题是，如果我想显示超过 22 位数字，即进行更精确的计算，例如 100 位数字，该怎么办？是否可以使用基础 R，或者我是否需要额外的包/函数？

编辑:感谢jmoy的建议，我尝试了sprintf("%.100f",pi)并且它给出了

[1] "3.1415926535897931159979634685441851615905761718750000000000000000000000000000000000000000000000000000"

有 48 位小数。这是 R 可以处理的最大限制吗？

Answer 1

这只是一个建议，因为您可以很容易地编写一个忽略选项值的打印函数。内置打印和格式化功能确实使用 options 值作为默认值。

关于第二个问题，由于 R 使用有限精度算术，因此您的答案在小数点后 15 或 16 位之后就不准确，因此通常不需要更多。 gmp和 rcdd包处理多精度算术(通过与 gmp 库的接口(interface))，但这主要与大整数有关，而不是 double 的更多小数位。

Mathematica或Maple将允许您根据需要提供任意数量的小数位。

编辑:
考虑小数位和有效数字之间的差异可能会很有用。如果您正在进行的统计测试依赖于超过 15 位有效数字的差异，那么您的分析几乎肯定是垃圾。

另一方面，如果您只处理非常小的数字，那就不是什么问题，因为 R 可以处理小到 .Machine$double.xmin 的数字(通常为 2e- 308)。

比较这两种分析。

x1 <- rnorm(50, 1, 1e-15)
y1 <- rnorm(50, 1 + 1e-15, 1e-15)
t.test(x1, y1)  #Should throw an error

x2 <- rnorm(50, 0, 1e-15)
y2 <- rnorm(50, 1e-15, 1e-15)
t.test(x2, y2)  #ok

在第一种情况下，数字之间的差异仅出现在许多有效数字之后，因此数据“几乎恒定”。在第二种情况下，虽然数字之间的差异大小相同，但与数字本身的大小相比，它们是很大的。

<小时/>

正如 e3bo 所提到的，您可以使用 Rmpfr 包来使用多精度 float 。

mpfr("3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825")

与常规( double )数字向量相比，这些向量使用起来更慢且占用更多内存，但如果您遇到条件较差的问题或不稳定的算法，它们可能会很有用。