3

我正在尝试重现导致数据集的计算data.ref。我想通过将参考数据与我的计算结果进行比较来测试我当前的实现效果如何,data.my. 由于数据的每一列都应该在列内具有可比较的大小,但不一定在列之间,我一直在看

(data.ref - data.my) / data.ref

将错误放在可比较的范围内。然而,由于数据最终将被四舍五入,我真正想做的只是快速而肮脏地检查数据有多少有效数字值得一致。也就是说,由于我期望data.ref并且data.my彼此非常接近,我想回答这个问题:每对对应条目不同的第一个有效数字是多少?

是否有执行此操作的 R 函数?

4

1 回答 1

1

ceiling(log10(abs(data.ref, data.my)))似乎可以解决问题。

例子:

> data.my <- c(20, 30, 32, 32.01, 32.012)
> data.ref <- rep(32, length(data.my))
> ceiling(log10(abs(data.my - data.ref)))
[1]    2    1 -Inf   -2   -1
于 2015-10-15T18:15:43.113 回答