Опубликован: 16.12.2009 | Уровень: для всех | Доступ: платный
Лекция 4:

Статистический анализ числовых величин (непараметрическая статистика)

Неустойчивость параметрических методов отбраковки резко выделяющихся результатов наблюдений

При обработки реальных экономических данных, полученных в процессе наблюдений, измерений, расчетов, иногда один или несколько результатов наблюдений резко выделяются, т.е. далеко отстоят от основной массы данных. Такие резко выделяющиеся результаты наблюдений часто считают содержащими грубые погрешности, соответственно называют промахами или выбросами. В рассматриваемых случаях возникает естественная мысль о том, что подобные наблюдения не относятся к изучаемой совокупности, поскольку содержат грубую погрешность, а получены в результате ошибки, промаха. В метрологии об этом явлении говорят так: "Грубые погрешности и промахи возникают из-за ошибок или неправильных действий оператора (его психо-физиологического состояния, неверного отсчета, ошибок в записях или вычислениях, неправильного включения приборов и т.п.), а также при кратковременных резких изменений проведения измерений (вибрации, поступления холодного воздуха, толчка прибора оператором и т.п.). Если грубые погрешности и промахи обнаруживают в процессе измерений, то результаты, содержащие их, отбрасывают. Однако чаще всего их выявляют только при окончательной обработке результатов измерений с помощью специальных критериев оценки грубых погрешностей" [7, с.46-47].

Есть два подхода к обработке данных, которые могут быть искажены грубыми погрешностями и промахами:

  1. отбраковка резко выделяющихся результатов наблюдений, т.е. обнаружение наблюдений, искаженных грубыми погрешностями и промахами, и исключение их из дальнейшей статистической обработки;
  2. применение устойчивых (робастных) методов обработки данных, На результаты работы которых мало влияет наличие небольшого числа грубо искаженных наблюдений

В настоящем пункте обсуждаются методы отбраковки.

Наиболее изучена ситуация, когда результаты наблюдений - числа x_1., x_2.,\dots, x_n., резко выделяется один результат наблюдения, для определенности, максимальный xmax .

Простейшая вероятностно-статистическая модель такова [8]. При нулевой гипотезе H_0 результаты наблюдения x_1., x_2.,\dots, x_n рассматриваются как реализация независимых одинаково распределенных случайных величин числа X_1., X_2.,\dots, X_n. с функцией распределения F(x). При альтернативной гипотезе H_1 случайные величины X_1., X_2., \dots, X_n. также независимы, X_1., X_2.,\dots, X_{n-1} имеют распределение F(x) , а X_n - распределение G(x), оно "существенно сдвинуто вправо" относительно F(x), например, G(x)=F(x - A) , где A достаточно велико. Если альтернативная гипотеза справедлива, то при A \to \infty вероятность равенства

X_n=\max(X_1, X_2, \dots, X_n)

стремится к 1, поэтому естественно применять решающее правило следующего вида:

если x_{max}.> d, то принять H_1,

если x_{max}.\le d, то принять H_0, (1)

где d - параметр решающего правила, который следует определять из вероятностно-статистических соображений.

При справедливости нулевой гипотезы

P\{max_1 \le i \n X_i >d\}=\{F(d)\}^n

Статистический критерий проверки гипотезы H_0, основанный на решающем правиле вида (1), имеет уровень значимости \alpha, если

P\{\max_{1 \le i \le n}X_i >d\}=1-\{F(d)\}^n= \alpha

т.е.

F(d)=\sqrt[n]{1-\alpha} ( 2)

Из соотношения (2) определяют граничное значение d=d(\alpha, n) в решающем правиле (1).

При больших n и малых \alpha

F(d)=\sqrt[n]{1- \alpha}=1-\frac{\alpha}{n}+O \left( \frac{\alpha^2}{n^2} \right) ( 3)

поэтому в качестве хорошего приближения к d(\alpha, n) рассматривают (1-\alpha/n) - квантиль распределения F(x) .

Пусть правило отбраковки задано в соответствии с выражениями (1) и (2) с некоторой функцией распределения F, однако выборка берется из функции распределения G, мало отличающейся от F в смысле расстояния Колмогорова

\rho (F,G)=\sup_x|F(x)-G(x)| \le \delta ( 4)

С помощью соотношения (3) получаем, что величина \gamma= G(d) для d из уравнения (2) находится между \gamma_1=\max(0, 1-\frac{\alpha}{n}- \delta) и \gamma_2=\min(1-\frac{\alpha}{n}+ \delta, 1). Уровень значимости критерия, построенного для F, при применении к наблюдениям из G есть 1- \gamma^n и может принимать любые значения в отрезке [1- \gamma_2^n; 1-\gamma_1^n] . В частности, при \delta = 0,01, \alpha =0,05, n = 5 возможные значения уровня значимости заполняют отрезок [0; 0,1], т.е. уровень значимости может быть в 2 раза выше номинального, а если n возрастает до 30, то максимальный уровень значимости есть 0,297, т.е. почти в 6 раз выше номинального. При дальнейшем росте n верхняя граница для уровня значимости, как нетрудно видеть, приближается к 1.

Рассмотрим и другой вопрос - насколько правило отбраковки с уровнем значимости \alpha для G может отличаться от такового для F при справедливости неравенства (4). С использованием соотношения (3) заключаем, что из

G(d)=1- \frac{\alpha}{n} ( 5)

следует, что \gamma_1 \le F(d) \le \gamma_2 где \gamma_1 и \gamma_2 выписаны выше. Решение уравнения (5) может принимать любое значение в отрезке [F^{-1}(\gamma_1); F^{-1}(\gamma_2)]. В частности, при \alpha=0.05 и n = 5 для стандартного нормального распределения F имеем d(\alpha, n) = 2,319, при \delta =0,01 решение уравнения (5) может принимать любое значение в отрезке [2,054; + \infty], при \delta =0,005 - любое значение в [2,170; 2,576].

При использовании любого другого расстояния между функциями распределения выводы о неустойчивости правил отбраковки также справедливы. Отметим, что проведенные рассмотрения выполнены в рамках "общей схемы устойчивости" (см. ниже лекцию об устойчивости статистических процедур).

Рассмотренные примеры показывают, что при конкретном значении \delta = 0,01 в неравенстве (4) весьма неустойчивы как уровни значимости при фиксированном правиле отбраковки, так и параметр d правила отбраковки при фиксированном уровне значимости. Обсудим, насколько реалистично определение функции распределения с точностью \delta \le 0.01

Есть два подхода к определению функции распределения результатов наблюдений: эвристический подбор с последующей проверкой с помощью критериев согласия и вывод из некоторой вероятностной модели.

Пусть с помощью критерия согласия Колмогорова проверяется гипотеза о том, что выборка взята из распределения F. Пусть функции распределения F и G удовлетворяют соотношению (4). Пусть на самом деле выборка взята из распределения G, а не F. При каких \delta не удастся различить F и G? Для определенности, при каких \delta гипотеза согласия с F будет приниматься не менее чем в 50% случаев?

Критерий согласия Колмогорова основан на статистике

\gamma_n=\sqrt{n} \rho (F_n, H) ( 6)

где расстояние \rho между функциями распределения определено выше в формуле (4); H - та функция распределения, согласие с которой проверяется, а F_n - эмпирическая функция распределения (т.е. F_n(х) равно доле наблюдений, меньших х, в выборке объема n ). Как показал А.Н. Колмогоров в 1933 г., функция распределения случайной величины \gamma_n при росте объема выборки n сходится к некоторой функции распределения К(х), которую ныне называют функцией Колмогорова. При этом К(1,36)= 0,95 и К(0,83)=0,50.

Поскольку выборка взята из распределения G, то с вероятностью 0,50

\rho (F_n, G) < 0.83/\sqrt n ( 7)

(при больших n ). Тогда для рассматриваемой выборки с учетом неравенства (4) и неравенства треугольника для расстояния Колмогорова и симметричности этого расстояния имеем

\rho (F_n, F) \le \rho (F_n, G)+ \rho (G, F)= \rho (F_n, G)+ \rho (F, G) < 0.83/\sqrt n + \delta

Если

0.83/ \sqrt n + \delta \le 1.36/ \sqrt n

т.е.

\delta \sqrt n \le 0.53 ( 8)

то, согласно формуле (6), гипотеза согласия принимается по крайней мере с той же вероятностью, с которой выполнено неравенств (7), т.е. с вероятностью не менее 0,50. Для \sigma=0.01 это условие выполняется при n \le 2809. Таким образом, для определения функции распределения с точностью \sigma \le 0.01 с помощью критерия согласия Колмогорова необходимо несколько тысяч наблюдений, что для большинства эконометрических задач нереально.

При втором из названных выше подходов к определению функции распределения ее конкретный вид выводится из некоторой системы аксиом, в частности, из некоторой модели порождения соответствующей случайной величины. Например, из модели суммирования вытекает нормальное распределение, а из мультипликативной модели перемножения - логарифмически нормальное распределение. Как правило, при выводе используется предельный переход. Так, из Центральной Предельной Теоремы теории вероятностей вытекает, что сумма независимых случайных величин может быть приближена нормальным распределением. Однако более детальный анализ, в частности, с помощью неравенства Берри-Эссеена (см. предыдущий пункт) показывает, что для гарантированного достижения точности \sigma \le 0.01 необходимо более полутора тысяч слагаемых. Такого количества слагаемых реально, конечно, указать почти никогда нельзя. Это означает, что при решении практических эконометрических задач теория дает возможность лишь сформулировать гипотезу о виде функции распределения, а проверять ее надо с помощью анализа реальной выборки объема, как показано выше, не менее нескольких тысяч.

Таким образом, в большинстве реальных ситуаций определить функцию распределения с точностью \sigma \le 0.01 невозможно.

Итак, показано, что правила отбраковки, основанные на использовании конкретной функции распределения, являются крайне неустойчивыми к отклонениям от нее распределения элементов выборки, а гарантировать отсутствие подобных отклонений невозможно. Поэтому отбраковка по классическим правилам математической статистики не является научно обоснованной, особенно при больших объемах выборок. Указанные правила целесообразно применять лишь для выявления "подозрительных" наблюдений, вопрос об отброаковке которых должен решаться из соображений соответствующей предметной области, а не из формально-математических соображений.

Выше для простоты изложения рассмотрен лишь случай полностью известного распределения F, для которого изучено правило отбраковки, заданное формулами (1) и (2). Аналогичные выводы о крайней неустойчивости правил отбраковки справедливы, если "истинное распределение" принадлежит какому-либо параметрическому семейству, например, нормальному, Вейбулла-Гнеденко, гамма.

Параметрическим методам отбраковки, основанным на моделях тех или иных параметрических семейств распределений, посвящены тысячи книг и статей. Приходится признать, что они имеют в основном внутриматематический интерес. При обработке реальных данных следует применять устойчивые методы, в частности, непараметрические.

Михаил Агапитов
Михаил Агапитов

Не могу найти  требования по оформлению выпускной контрольной работы по курсу профессиональной переподготовки "Менеджмент предприятия"

Подобед Александр
Подобед Александр

Я нажал кнопку "начать курс" и почти его уже закончил, но для получения диплома на бумаге, нужно его же оплатить? Как оплатить? 

Вячеслав Гримальский
Вячеслав Гримальский
Россия
Михаил Байков
Михаил Байков
Россия, Москва, Московский Авиационный Институт, 2009