Как я могу обойти ошибку округления, которая вызывает бесконечный цикл в статистике Perl :: Descriptive?

Question

Я использую библиотеку Statistics::Descriptive в Perl для вычисления частотных распределений и решения проблемы с ошибкой округления с плавающей запятой.

Я передаю два значения, 0.205 и 0.205 (взятых из других чисел и 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 

} 

$ self-> 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

Я - статистика :: Описательный сопровождающий. Из-за его числового характера было сообщено много проблем округления. Я считаю, что этот конкретный был исправлен в более поздней версии той, которую вы использовали, которую я выпустил недавно, используя умножение для делений вместо + =.

Пожалуйста, используйте the most up-to-date version от CPAN, и это должно быть лучше.

Answer 2

Не совсем округления проблемы; Вы можете увидеть более точные значения с чем-то вроде

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

кажется мне, что есть ошибка в модуле, где она принимает диапазон выборки можно разделить на $ Перегородки штук; потому что плавающая точка не имеет бесконечной точности, это не всегда возможно. В вашем случае значения min и max являются точно смежными представляемыми значениями, поэтому не может быть более одного раздела. Я не знаю, для чего именно модуль использует разделы, поэтому я не уверен, каким может быть это воздействие. Еще одна проблема в модуле заключается в том, что он использует числа как хеш-ключи, которые неявно строит их, что немного округляет значение.

Вы можете иметь некоторый успех в отмывании данных через stringization перед подачей его к модулю:

$data = 0+"$data"; 

Это, по крайней мере убедиться, что два числа, что (с точностью печати по умолчанию) появляются равно фактически равны.

Answer 3

Это не должно вызывать бесконечный цикл. То, что вызовет бесконечный цикл, будет $self->{sample_range}/$partitions равно 0.