R 计算的百分比总和不等于一 - 精度问题答案

【问题标题】：R computed percentages do not sum up to one - precision issueR 计算的百分比总和不等于一 - 精度问题
【发布时间】：2015-12-08 08:51:24
【问题描述】：

我正在尝试处理 R 中的百分比，但我遇到了奇怪的问题。当我将向量的值转换为向量的sum 的百分比时，经常会发生它们加起来不等于一的情况。

小例子：

data <- rnorm(1000)*100
max <- 50
unlist(lapply(0:(1000/max-1), 
     function(i) 
        sum(
            data[(i*max+1):(i*(max+1))]
            /
            sum(data[(i*max+1):(i*(max+1))])
           )
        ))-1

它应该给出零向量，但是我得到了这个：

[1]  0.000000e+00  0.000000e+00 -1.110223e-16 -1.110223e-16  0.000000e+00 -1.110223e-16  0.000000e+00  0.000000e+00  0.000000e+00
[10]  0.000000e+00  0.000000e+00  2.220446e-16  0.000000e+00 -4.440892e-16  0.000000e+00  0.000000e+00  0.000000e+00  4.440892e-16
[19] -1.110223e-16  0.000000e+00

有什么补救办法吗？

【问题讨论】：

我的机器人朋友 Dason 用来说明 R 使用的浮点运算的例子是 sqrt(2)^2 == 2
@TylerRinker -- 或者，更清楚一点，print(sqrt(2)^2, digits=22)
另见 Patrick Burns 的优秀 R Inferno 的 Circle 1。

标签： r

【解决方案1】：

他们偏离了一个微不足道的数字。如果您想更改浮点运算中固有的这些微不足道的差异的显示方式，您可以使用 format 函数或其表亲之一，如 sprintf 或 formatC。这实际上是FAQ 7.31. 的一个实例如果您确实需要格式化方面的帮助，您应该描述一个特定的应用程序。如果你想强制查看零，你也可以使用 round()

round( unlist(lapply(0:(1000/max-1), 
 function(i) 
    sum(
        data[(i*max+1):(i*(max+1))]
        /
        sum(data[(i*max+1):(i*(max+1))])
       )
    ))-1  , digits=4)

【讨论】：

【解决方案2】：

一个更重要的问题是为什么你认为这些应该是0？

您正在使用浮点运算，并非所有数字都可以在您的计算机中精确表示。这涵盖（或相关）R FAQ 7.31，它解释了这种现象。

您可以忽略它（出于所有意图和目的，这些值为0）

> all.equal(tmp, rep(0, length(tmp))) ## tmp contain your numbers
[1] TRUE

或学习根据您的特定操作相应地处理它。一种方法是在某种程度上将它们四舍五入：

> round(tmp, 2)
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
> round(tmp, 3)
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
> round(tmp, 4)
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
> round(tmp, 5)
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

但这确实取决于你想用这些数字做什么。

【讨论】：