4.1. Расчет оптимального интервала на гистограмме
Расчет оптимального количества интервалов и диапазона значений, попадающих в интервал гистограммы распределения, является важной составляющей базовой статистики и может служить основой для выделения групп равномерно распределенных данных и выработки классификации на их основе. В данных методических рекомендациях рассматривается несколько простых подходов расчета оптимального количества интервалов и диапазона изменения значений. Первый подход основан на эвристической формуле Стерджесса (14)Herbert A. S. The choice of a class interval // Journal of the American statistical association. 1926. Vol. 21, №153. Р. 65–66.:
n = log2 N + 1 = 3.33 ∙ log10 N + 1, (14)
где N – объем анализируемой выборки (количество значений); n – количество интервалов разбиения данных.
В работеШторм Р. Теория вероятностей. Математическая статистика. Статистический контроль качества. М.: Мир, 1970. 368 с., описывающей применение статистического анализа при контроле качества продукции, применяют уравнение Брукса–Каррузера (15):
n = 5 ∙ log10 N, (15)
Наиболее простым методом оценки количества интервалов является (16)Heinhold I., Gaede K. W. Ingeniur statistic. München; Wien, Springler Verlag, 1964. 352 p.:
, (16)
При больших объемах данных рекомендуется использовать уравнение (17)Ченцов Н. Н. Статистические решающие правила и оптимальные выводы. М.: Наука, 1972. 520 с.:
, (17)
После расчета количества интервалов разбиения гистограммы распределения проводят вычисления ширины интервала по уравнению (18):
, (18)
где max(X) – максимальное значение исследуемой переменной X; min(X) – минимальное значение исследуемой переменной X; n – количество интервалов на гистограмме.
Соответственно, с использованием уравнения (18) вычисляются первый и последующие интервалы по уравнению (19):
Xi+1 = Xi + ∆Xhist, (19)
где Xi+1 – правая граница i-го интервала, не превышающая max(X); Xi – левая граница i-го интервала, начало отсчета которой равно min(X).
Высота столбца диаграммы определяется по уравнению (20):
, (20)
где N – общее количество исследований; Ni – количество исследований в i-м интервале; Hi – высота интервала.
Более точная (и она же более сложная) оценка размера интервала может быть получена при оценках статистической мощностиЛемешко Б. Ю., Чимитова Е. В. О выборе числа интервалов в критериях согласия типа c2// Заводская лаборатория. Диагностика материалов. 2003. Т. 69, №1. С. 61–67., данный подход не рассматривается в настоящих методических рекомендациях.