Калькулятор критерия Граббса: поиск грубых промахов

Проверьте серию измерений на грубые промахи по критерию Граббса (Смирнова–Граббса). Калькулятор найдёт подозрительное значение, сравнит статистику G с критическим значением, при необходимости повторит проверку и выдаст очищенные данные с новым средним и СКО.

Через пробел, точку с запятой или с новой строки; можно вставить столбец из Excel. От 3 до 10 000 значений.

Найдено выбросов—
Среднее без выбросов—
СКО без выбросов—
Исходные данные—

Ход проверки

Данные без выбросов

Рекламный блок «после калькулятора» (ads/after-calc.html)

Что такое грубый промах

В серии повторных измерений иногда попадается значение, резко отличающееся от остальных: перепутали деление шкалы, неверно записали цифру, прибор дал сбой. Такой промах сильно искажает среднее и особенно стандартное отклонение, поэтому перед обработкой результатов его нужно найти и исключить. Но исключать «неудобные» значения на глаз нельзя — нужен объективный статистический критерий. Самый распространённый для нормально распределённых данных — критерий Граббса.

Как работает критерий

Находят значение, наиболее удалённое от среднего, и вычисляют статистику

G = max |xᵢ − x̄| / s

где x̄ — среднее, s — выборочное стандартное отклонение всей серии, включая подозрительное значение. Если G больше критического значения Gкр для данного объёма выборки n и уровня значимости α, значение признают выбросом. Критическое значение калькулятор вычисляет точно, через распределение Стьюдента:

Gкр = (n − 1)/√n · √( t² / (n − 2 + t²) ), t — квантиль Стьюдента уровня α/(2n) с n − 2 степенями свободы

После исключения выброса проверку можно повторить для оставшихся данных — так находят несколько промахов по очереди. Повторное применение снижает строгость процедуры, поэтому если «выбросов» набирается больше одного-двух на десяток измерений, стоит усомниться не в отдельных значениях, а в нормальности распределения или в самой методике.

Рекламный блок «внутри статьи» (ads/in-article.html)

Пример

Восемь измерений длины, мм: 10,2; 10,4; 10,3; 10,5; 10,3; 10,4; 11,6; 10,2. Среднее 10,4875, стандартное отклонение 0,4612. Дальше всех от среднего значение 11,6: G = (11,6 − 10,4875) / 0,4612 ≈ 2,412. Для n = 8 и α = 0,05 критическое значение 2,127 — значит, 11,6 является промахом. Без него среднее 10,329, s = 0,111. Повторная проверка: максимальное отклонение даёт G ≈ 1,541, а критическое значение для n = 7 равно 2,020 — больше выбросов нет. Обратите внимание, как изменилось СКО: с 0,46 до 0,11 — один промах завышал разброс в четыре раза.

Таблица критических значений

Двусторонний критерий: проверяется отклонение в любую сторону от среднего.

nα = 0,10α = 0,05α = 0,01
31,1531,1541,155
41,4621,4811,496
51,6711,7151,764
61,8221,8871,973
71,9382,0202,139
82,0322,1272,274
92,1102,2152,387
102,1762,2902,482
112,2342,3552,564
122,2852,4122,636
132,3312,4622,699
142,3722,5072,755
152,4092,5482,806
162,4432,5862,852
172,4752,6202,894
182,5042,6522,932
192,5312,6812,968
202,5572,7083,001
212,5802,7343,031
222,6032,7583,060
232,6242,7803,087
242,6442,8023,112
252,6632,8223,135
262,6812,8413,158
272,6982,8593,179
282,7142,8763,199
292,7302,8933,218
302,7452,9083,236
352,8122,9783,316
402,8683,0363,381
502,9573,1283,482
603,0273,2003,560
803,1323,3063,673
1003,2103,3843,754

Частые вопросы

Чем критерий Граббса лучше правила трёх сигм?

Правило «отбросить всё, что дальше 3s от среднего» не учитывает объём выборки. В серии из 5–10 измерений отклонение больше 3s математически почти невозможно даже для явного промаха — правило его не заметит. Критерий Граббса подбирает порог под конкретное n.

Можно ли исключить выброс, если он «правильный»?

Статистический критерий лишь указывает на подозрительное значение. Если есть физическая причина выброса (сбой, ошибка записи), его исключают с пометкой в отчёте. Если причины нет, а значение воспроизводится, возможно, это реальный эффект — тогда исключать его нельзя.

Какой уровень значимости выбрать?

Обычно α = 0,05. Уровень 0,01 строже: выбросом признаётся только очень далёкое значение, риск выбросить нормальное измерение меньше. Уровень 0,10 мягче и подходит для предварительного отсева.

Рекламный блок «внизу страницы» (ads/bottom.html)