找到兩個非常相似的值之間的差異的第一個有效數字

Question

我想重現導致數據集data.ref的計算。我想通過將參考數據與我的計算結果data.my進行比較來測試我目前的實施做得如何。由於數據的每一列應有柱內可比的大小，而不是列之間不一定，我一直在尋找

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

把誤差在可比的規模。但是，由於數據最終會四捨五入，我真正想做的只是快速骯髒地檢查數據有多少有價值的數據。也就是說，因爲我期望data.ref和data.my彼此非常接近，所以我想回答這個問題：每對相應條目的不同之處，第一個有效數字是多少？

有沒有一個R函數來做到這一點？

Answer 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