PerlのStatistics :: Descriptiveで無限ループを引き起こす丸め誤差を回避するにはどうしたらいいですか？

Question

私は頻度分布を計算し、浮動小数点の丸め誤差の問題に対処するためにPerlでStatistics::Descriptiveライブラリを使用しています。

私は統計モジュールに0.205と0.205の2つの値（他の数値とsprintf'dから取得）を渡し、頻度分布を計算するよう依頼しますが、無限ループに陥っています。私はそれがやっていることがわかりますデバッガでステップスルー

：私は期待通り（範囲は最大 - 最小である）

my $interval = $self->{sample_range}/$partitions; 

my $iter = $self->{min}; 

while (($iter += $interval) < $self->{max}) { 

    $bins{$iter} = 0; 

    push @k, $iter; ##Keep the "keys" unstringified 

} 

$自己> sample_rangeは2.77555756156289e-17ではなく0を返します。これは、ループ（（min + = range）< max）が（すべての目的と目的のために）無限ループに入ることを意味します。

DB < 8> print $ self - > {max};
0.205
DB < 9> print $ self - > {min};
0.205
DB < 10> print $ self - > {max} - $ self - > {min};
2.77555756156289e-17

これは丸めの問題のようです。私はこれを私の側で修正する方法を考えることはできません、そして、私はライブラリの編集が良いアイデアであるとは確信していません。私は回避策または代替の提案を探しています。

乾杯、 ニール

Answer 1

私はStatistics :: Descriptive maintainerです。数値的性質のため、多くの丸め問題が報告されています。私は、この特定のものが、最近リリースされたものを使用しているものに、+ =の代わりに除算を使って、後のバージョンで修正されたと考えています。

CPANのthe most up-to-date versionを使用してください。

Answer 2

ない正確に丸め問題。あなたは何か欠陥は、それがサンプルの範囲は$パーティション片に分割することができますを前提とモジュールでありますように私には見える

printf("%.18g %.18g", $self->{max}, $self->{min}); 

などで、より正確な値を見ることができます。浮動小数点は無限の精度を持たないため、これは必ずしも可能ではありません。あなたの場合、最小値と最大値は正確に隣接した表現可能な値なので、複数のパーティションは存在できません。私はモジュールがパーティションを使用していることを正確にはわからないので、この影響がどのようなものかはわかりません。 モジュールのもう1つの問題は、数値をハッシュキーとして使用していることです。 は暗黙のうちに値をわずかに丸めます。

あなたはモジュールに を供給する前にstringizationを通して、あなたのデータを洗濯中にいくつかの成功を持っていることがあります。これは、少なくとも（デフォルトの印刷精度と）同じ表示されていることを二つの数字が実際にあることを保証します

$data = 0+"$data"; 

を等しい。

Answer 3

無限ループを引き起こしてはいけません。そのループが無限になる原因は、$self->{sample_range}/$partitionsが0の場合です。