Калькулятор критерия Граббса: поиск грубых промахов
Проверьте серию измерений на грубые промахи по критерию Граббса (Смирнова–Граббса). Калькулятор найдёт подозрительное значение, сравнит статистику G с критическим значением, при необходимости повторит проверку и выдаст очищенные данные с новым средним и СКО.
Ход проверки
Что такое грубый промах
В серии повторных измерений иногда попадается значение, резко отличающееся от остальных: перепутали деление шкалы, неверно записали цифру, прибор дал сбой. Такой промах сильно искажает среднее и особенно стандартное отклонение, поэтому перед обработкой результатов его нужно найти и исключить. Но исключать «неудобные» значения на глаз нельзя — нужен объективный статистический критерий. Самый распространённый для нормально распределённых данных — критерий Граббса.
Как работает критерий
Находят значение, наиболее удалённое от среднего, и вычисляют статистику
G = max |xᵢ − x̄| / s
где x̄ — среднее, s — выборочное стандартное отклонение всей серии, включая подозрительное значение. Если G больше критического значения Gкр для данного объёма выборки n и уровня значимости α, значение признают выбросом. Критическое значение калькулятор вычисляет точно, через распределение Стьюдента:
Gкр = (n − 1)/√n · √( t² / (n − 2 + t²) ), t — квантиль Стьюдента уровня α/(2n) с n − 2 степенями свободы
После исключения выброса проверку можно повторить для оставшихся данных — так находят несколько промахов по очереди. Повторное применение снижает строгость процедуры, поэтому если «выбросов» набирается больше одного-двух на десяток измерений, стоит усомниться не в отдельных значениях, а в нормальности распределения или в самой методике.
Пример
Восемь измерений длины, мм: 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 |
|---|---|---|---|
| 3 | 1,153 | 1,154 | 1,155 |
| 4 | 1,462 | 1,481 | 1,496 |
| 5 | 1,671 | 1,715 | 1,764 |
| 6 | 1,822 | 1,887 | 1,973 |
| 7 | 1,938 | 2,020 | 2,139 |
| 8 | 2,032 | 2,127 | 2,274 |
| 9 | 2,110 | 2,215 | 2,387 |
| 10 | 2,176 | 2,290 | 2,482 |
| 11 | 2,234 | 2,355 | 2,564 |
| 12 | 2,285 | 2,412 | 2,636 |
| 13 | 2,331 | 2,462 | 2,699 |
| 14 | 2,372 | 2,507 | 2,755 |
| 15 | 2,409 | 2,548 | 2,806 |
| 16 | 2,443 | 2,586 | 2,852 |
| 17 | 2,475 | 2,620 | 2,894 |
| 18 | 2,504 | 2,652 | 2,932 |
| 19 | 2,531 | 2,681 | 2,968 |
| 20 | 2,557 | 2,708 | 3,001 |
| 21 | 2,580 | 2,734 | 3,031 |
| 22 | 2,603 | 2,758 | 3,060 |
| 23 | 2,624 | 2,780 | 3,087 |
| 24 | 2,644 | 2,802 | 3,112 |
| 25 | 2,663 | 2,822 | 3,135 |
| 26 | 2,681 | 2,841 | 3,158 |
| 27 | 2,698 | 2,859 | 3,179 |
| 28 | 2,714 | 2,876 | 3,199 |
| 29 | 2,730 | 2,893 | 3,218 |
| 30 | 2,745 | 2,908 | 3,236 |
| 35 | 2,812 | 2,978 | 3,316 |
| 40 | 2,868 | 3,036 | 3,381 |
| 50 | 2,957 | 3,128 | 3,482 |
| 60 | 3,027 | 3,200 | 3,560 |
| 80 | 3,132 | 3,306 | 3,673 |
| 100 | 3,210 | 3,384 | 3,754 |
Частые вопросы
Чем критерий Граббса лучше правила трёх сигм?
Правило «отбросить всё, что дальше 3s от среднего» не учитывает объём выборки. В серии из 5–10 измерений отклонение больше 3s математически почти невозможно даже для явного промаха — правило его не заметит. Критерий Граббса подбирает порог под конкретное n.
Можно ли исключить выброс, если он «правильный»?
Статистический критерий лишь указывает на подозрительное значение. Если есть физическая причина выброса (сбой, ошибка записи), его исключают с пометкой в отчёте. Если причины нет, а значение воспроизводится, возможно, это реальный эффект — тогда исключать его нельзя.
Какой уровень значимости выбрать?
Обычно α = 0,05. Уровень 0,01 строже: выбросом признаётся только очень далёкое значение, риск выбросить нормальное измерение меньше. Уровень 0,10 мягче и подходит для предварительного отсева.
Комментарии и вопросы