Question

有人可以向我解释一下all.equal的容差参数吗？

手册说（?all.equal）：

tolerance：数字≥0。小于容差的差异不是   考虑。

scale = NULL（默认值）的数值比较首先完成   计算两个数值向量的平均绝对差值。   如果这小于公差或不是有限的绝对差异   使用，否则相对差异按平均绝对值缩放   差。

示例：

all.equal(0.3, 0.26, tolerance=0.1)

返回Mean relative difference: 0.1333333

为什么这里会返回平均相对差异？两个数值向量的平均绝对差值不小于容差吗？

0.3 - 0.26 = 0.04 < 0.1

谢谢！

Answer 1

如果target大于tolerance，则似乎会检查relative error <= tolerance。也就是说，abs(current-target)/target <= tolerance in：

all.equal(target, current, tolerance)

例如：

all.equal(3, 6, tolerance = 1)
# TRUE --> abs(6-3)/3 <= 1

相反，如果target小于tolerance，则all.equal会使用mean absolute difference。

all.equal(0.01, 4, tolerance = 0.01)
# [1] "Mean absolute difference: 3.99"

all.equal(0.01, 4, tolerance = 0.00999)
# [1] "Mean relative difference: 399"

all.equal(4, 0.01, tolerance = 0.01)
# [1] "Mean relative difference: 0.9975"

但是，这是不文档所述的内容。为了进一步了解为什么会发生这种情况，让我们看看all.equal.numeric中的相关代码段：

# take the example: all.equal(target=0.01, current=4, tolerance=0.01)
cplx <- is.complex(target) # FALSE
out <- is.na(target) # FALSE
out <- out | target == current # FALSE

target <- target[!out] # = target (0.01)
current <- current[!out] # = current (4)

xy <- mean((if(cplx) Mod else abs)(target - current)) # else part is run = 3.99

# scale is by default NULL
what <- if (is.null(scale)) {
    xn <- mean(abs(target)) # 0.01
    if (is.finite(xn) && xn > tolerance) { # No, xn = tolerance
        xy <- xy/xn
        "relative"
    }
    else "absolute" # this is computed for this example
}
else {
    xy <- xy/scale
    "scaled"
}

上面代码中检查的所有内容（仅显示OP中示例的必要部分）是：删除任何NA和相等的值（target和{{ 1}}）来自current和target。然后将current计算为xy和target的平均绝对差值。但是，确定它是current还是relative取决于absolute部分。此处what未检查任何条件。它取决于xy上{em>仅 xn。

因此，总而言之，OP粘贴的部分（为方便起见粘贴在这里）：

如果此（意思是，平均绝对差）小于容差或不是有限的，则使用绝对差异，否则相对差异按平均绝对差异缩放。 / p>

似乎错误/误导。

all.equal（）的容差参数如何工作？

1 个答案: