Квантили распределения хи-квадрат

Кванти́ли распределе́ния хи-квадра́т — числовые характеристики, широко используемые в задачах математической статистики таких как построение доверительных интервалов, проверка статистических гипотез и непараметрическое оценивание.

Квантиль хи-квадрат — это число (величина хи-квадрат), при котором функция распределения хи-квадрат равна заданной (затребованной) вероятности а.

Равенство функции распределения хи-квадрат вероятности а означает, что с вероятностью а будут наблюдаться значения хи-квадрат, не большие, чем найденный (определенный согласно функции распределения) квантиль хи-квадрат. Таким образом, найти квантиль означает разграничить распределения хи-квадрат согласно заданной вероятности а.

Определение

Пусть  — функция распределения хи-квадрат с степенями свободы, и . Тогда -квантилем этого распределения называется число такое, что

.

Замечания

  • Прямо из определения следует, что случайная величина, имеющая распределение хи-квадрат с степенями свободы, не превышает значение с вероятностью и превышает его с вероятностью .
  • Функция строго возрастает для любого . Следовательно, определена её обратная функция , и
.
  • Функция не имеет простого представления. Однако, возможно вычислить её значения численно.

Аппроксимация квантилей

Для получения приближённых значений квантилей распределения хи-квадрат существует целый ряд аппроксимаций. Их обзор и сравнение даны в статье: Zar (1978)[1]. Две из них приведены ниже.

  • Аппроксимация Корниша — Фишера[2]:
,

где:

,

Здесь — обратная функция Лапласа, называемая также нормальной квантильной функцией или пробит-функцией, а — обратная функция ошибок. Для этих функций, в свою очередь, имеются аппроксимации. Например, использовалась следующая аппроксимация:

при
при .

Более точное приближение даёт использование аппроксимации Виницкого[3] для обратной функции ошибок:

где — подгоночный параметр. Относительная погрешность формулы Виницкого с параметром не превышает 0.002[3] для всех ненулевых значений

  • Аппроксимация Голдштейна[4]:
,

где d определяется аналогично, а коэффициенты a, b,c приведены в таблице

a b c
1.0000886 -0.2237368 -0.01513904
0.4713941 0.02607083 -0.008986007
0.0001348028 0.01128186 0.02277679
-0.008553069 -0.01153761 -0.01323293
0.00312558 0.005169654 -0.006950356
-0.0008426812 0.00253001 0.001060438
0.00009780499 -0.001450117 0.001565326

Таблица квантилей

Нижеприведённая таблица получена с помощью функции chi2inv Архивная копия от 4 декабря 2009 на Wayback Machine пакета MATLAB.

Также квантили можно получить с помощью других программных средств:

Чтобы получить значение , необходимо найти строку, соответствующую нужному , и колонку, соответствующую нужному . Искомое число находится в таблице на их пересечении.
Например:

;
.

См. также

Примечания

  1. Zar J.H. Approximations for the percentage points of the chi-squared distribution (англ.) // Journal of the Royal Statistical Society. Series C (Applied Statistics). — Wiley, 1978. — Vol. 27, iss. 3. — P. 280-290. Архивировано 12 января 2024 года.
  2. Golberg H., Levine H. Approximate formulas for the percentage points and normalization of t and // AMS. 1945. V.17. P. 216—225.
  3. 1 2 Winitzki S. A handy approximation for the error function and its inverse (англ.). — 2008.
  4. Goldstein R.B. Chi-square quantiles, Algorithm 451 // Commun. Assoc. Comp. 1973. V. 16. P. 483—485.