首页 > 解决方案 > 如何使用 R 创建向量 (10^16 + 1, 10^16 + 2, ... , 10^16 + 1000)?

问题描述

我想创建以下向量 (10^16 + 1, 10^16 + 2, ... , 10^16 + 1000)。但如果我这样尝试:

daten <- rep(1e+16,1000)
daten + 1:1000

它给了我一个近似值:

[1] 10000000000000000 10000000000000002 10000000000000004 10000000000000004
[5] 10000000000000004 10000000000000006 10000000000000008 10000000000000008
[9] 10000000000000008 10000000000000010 10000000000000012 10000000000000012
[13] 10000000000000012 10000000000000014 10000000000000016 10000000000000016 
 ......

如何让 R 计算更精确?谢谢

标签: rprecisionmantissa

解决方案


Rmpfr软件包可让您根据需要获得尽可能高的精度。使用函数创建数字mpfr,指定值 + 您需要的精度位数。

library(Rmpfr)
daten <- mpfr(1:1000, 120)
daten <- daten + mpfr(1e16, 120)
print(daten[1:5])

这会产生以下结果:

5 'mpfr' numbers of precision  120   bits 
[1] 10000000000000001 10000000000000002 10000000000000003 10000000000000004 10000000000000005

推荐阅读