r - 2 つの非常に類似した値の差の最初の有効数字を見つける

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