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

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

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

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

Определение

править

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

 .

Замечания

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

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

править

Для получения приближённых значений квантилей распределения хи-квадрат   существует целый ряд аппроксимаций. Их обзор и сравнение даны в статье: 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.