Статистический анализ данных (курс лекций, К.В.Воронцов)/2013

Материал из MachineLearning.

(Различия между версиями)
Перейти к: навигация, поиск
м (→Оценки)
м (→Оценки)
Строка 2: Строка 2:
= Оценки =
= Оценки =
{|class = "standard sortable"
{|class = "standard sortable"
-
! class="unsortable"| Студент !! #1 (1 балл) !! #2 (2 балла) !! Рецензирование #2 (1 балл) !! #3 (2 балла) !! Сумма
+
! class="unsortable"| Студент !! #1 (1 балл) !! #2 (2 балла) !! Рецензирование #2 (1 балл) !! #3 (2 балла) !! Рецензирование #3 (1 балл) !! Сумма
|-
|-
-
| align="center"| Березин Алексей || || || || ||
+
| align="center"| Березин Алексей || || || || || ||
|-
|-
-
| align="center"| Борисов Михаил || || || || ||
+
| align="center"| Борисов Михаил || || || || || ||
|-
|-
-
| align="center"| Гавриков Михаил || || || || ||
+
| align="center"| Гавриков Михаил || || || || || ||
|-
|-
-
| align="center"| Зак Евгений || || || || ||
+
| align="center"| Зак Евгений || || || || || ||
|-
|-
-
| align="center"| Исмагилов Тимур || || || || ||
+
| align="center"| Исмагилов Тимур || || || || || ||
|-
|-
-
| align="center"| Кондрашкин Дмитрий || || || || ||
+
| align="center"| Кондрашкин Дмитрий || || || || || ||
|-
|-
-
| align="center"| Куракин Александр || || || || ||
+
| align="center"| Куракин Александр || || || || || ||
|-
|-
-
| align="center"| Лобачева Екатерина || || || || ||
+
| align="center"| Лобачева Екатерина || || || || || ||
|-
|-
-
| align="center"| Любимцева Мария || || || || ||
+
| align="center"| Любимцева Мария || || || || || ||
|-
|-
-
| align="center"| Малышева Екатерина || || || || ||
+
| align="center"| Малышева Екатерина || || || || || ||
|-
|-
-
| align="center"| Меркулова Татьяна || || || || ||
+
| align="center"| Меркулова Татьяна || || || || || ||
|-
|-
-
| align="center"| Морозова Дарья || || || || ||
+
| align="center"| Морозова Дарья || || || || || ||
|-
|-
-
| align="center"| Нижибицкий Евгений || || || || ||
+
| align="center"| Нижибицкий Евгений || || || || || ||
|-
|-
-
| align="center"| Новиков Максим || || || || ||
+
| align="center"| Новиков Максим || || || || || ||
|-
|-
-
| align="center"| Огнева Дарья || || || || ||
+
| align="center"| Огнева Дарья || || || || || ||
|-
|-
-
| align="center"| Остапец Андрей || || || || ||
+
| align="center"| Остапец Андрей || || || || || ||
|-
|-
-
| align="center"| Потапенко Анна || || || || ||
+
| align="center"| Потапенко Анна || || || || || ||
|-
|-
-
| align="center"| Ромов Петр || || || || ||
+
| align="center"| Ромов Петр || || || || || ||
|-
|-
-
| align="center"| Фонарев Александр || || || || ||
+
| align="center"| Фонарев Александр || || || || || ||
|-
|-
-
| align="center"| Шаймарданов Ильдар || || || || ||
+
| align="center"| Шаймарданов Ильдар || || || || || ||
|}
|}
* Штраф за просрочку сдачи заданий начисляется из расчета 0.1 балла за сутки.
* Штраф за просрочку сдачи заданий начисляется из расчета 0.1 балла за сутки.

Версия 13:15, 4 октября 2013

Содержание

Оценки

Студент #1 (1 балл) #2 (2 балла) Рецензирование #2 (1 балл) #3 (2 балла) Рецензирование #3 (1 балл) Сумма
Березин Алексей
Борисов Михаил
Гавриков Михаил
Зак Евгений
Исмагилов Тимур
Кондрашкин Дмитрий
Куракин Александр
Лобачева Екатерина
Любимцева Мария
Малышева Екатерина
Меркулова Татьяна
Морозова Дарья
Нижибицкий Евгений
Новиков Максим
Огнева Дарья
Остапец Андрей
Потапенко Анна
Ромов Петр
Фонарев Александр
Шаймарданов Ильдар
  • Штраф за просрочку сдачи заданий начисляется из расчета 0.1 балла за сутки.
  • Задание считается сданным на момент получения проверяющим письма с отчётом (и кодом, если это указано в задании), при условии отсутствия необходимости внесения дополнений и исправлений.
  • Для допуска к экзамену необходимо сдать как минимум два задания, обязательно включая первое.

Задание 1. Исследование свойств одномерных статистических критериев на модельных данных

Необходимо провести исследование одного или нескольких классических критериев проверки статистических гипотез. Интерес представляет поведение достигаемого уровня значимости (p-value) как функции размера выборок и параметров распределения. В соответствии с индивидуальными параметрами задания необходимо указанным способом сгенерировать одну или несколько выборок из заданного распределения, выполнить проверку гипотезы при помощи соответствующего критерия, а затем многократно повторить эту процедуру для различных значений параметров. По результатам расчётов необходимо построить требуемые в задании графики, среди которых могут быть следующие:

  1. график зависимости достигаемого уровня значимости от значений параметров при однократном проведении эксперимента;
  2. график зависимости достигаемого уровня значимости одного или двух критериев от значений параметров, усреднённого по большому количеству повторений эксперимента (например, по 1000 повторений);
  3. график с эмпирическими оценками мощности одного или двух критериев для разных значений параметров.

В качестве оценки мощности принимается доля отвержений нулевой гипотезы среди всех проверок. То есть, если эксперимент повторялся k раз для каждого набора значений параметров, и в m из k случаев гипотеза была отвергнута на некотором фиксированном уровне значимости \alpha (примем \alpha=0.05), оценкой мощности будет отношение m/k.

Необходимо сдать: выполненный в Tex или Microsoft Word отчёт с описанием алгоритма, построенными графиками и выводами (объяснение полученных результатов моделирования, границы применимости критерия и т.д.), а также код на R, Матлабе или Питоне, при запуске которого на экран выводятся графики, соответствующие имеющимся в отчёте.

Задание принимается до 23:59 19.10.

Пример задания

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

X_1^n, \;\; X_{1i} \sim 0.9\cdot N(\mu_1,1)+ 0.1\cdot U\left[-5+\mu_1,5+\mu_1\right] — выборка длины n из смеси стандартного нормального N(\mu_1,1) и равномерного U\left[-5+\mu_1,5+\mu_1\right] распределений с весами 0.9 и 0.1 соответственно (при генерации выборки используется случайный датчик — если его значение не превосходит 0.9, то добавляем в выборку элемент, взятый из нормального распределения, иначе — элемент, взятый из равномерного).

X_2^n, \;\; X_{2i} \sim 0.9\cdot N(\mu_2,1)+ 0.1\cdot U\left[-5+\mu_2,5+\mu_2\right] — аналогичная выборка.

H_0\,:\; \mu_1=\mu_2, \;\; H_1\,:\; \mu_1\neq\mu_2.

\mu_1=0, \;\; \mu_2=-2\,:\,0.01\,:\,2, \;\; n=15\,:\,5\,:\,200.

При каждом значении \mu_2 выборки для разных значений n генерируются независимо.

Значения достигаемого уровня значимости при однократной генерации выборок.

Значения достигаемого уровня значимости, усрёднённые по 3000 экспериментам.

Значения эмпирических оценок мощности критерия при проведении 3000 экспериментов (\alpha=0.05).

Заметим, что однократная генерация выборок даёт достаточно нестабильные результаты, не позволяя точно оценить границы области, где нулевая гипотеза отклоняется, поэтому и необходимо усреднение по большому числу экспериментов.

Видно, что при достаточно большой разнице между средними и большом размере выборок наличие шума не мешает уверенно отклонять гипотезу однородности. Когда, наоборот, разница между средними невелика (меньше 0.2-0.5 в зависимости от размера выборок), мощность близка к нулю, а среднее значение достигаемого уровня значимости колеблется около 0.5, что логично, так как его распределение при справедливости нулевой гипотезы равномерно на [0,1].

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

Разность средних достигаемых уровней значимости на выборках без шума и с шумом.

Разность эмпирических оценок мощности на выборках без шума и с шумом.

Видно, что наличие шума всё меньше влияет на работу критерия с ростом объёма выборок и разницы между их средними. Тем не менее, в некоторых областях изменения параметров потеря мощности из-за 10% зашумления может составлять до 20%, а средний достигаемый уровень значимости может быть выше на 0.1.

Отметим, что приведённые количественные выводы справедливы только для шума рассматриваемой структуры.

Задания

Ссылки

Источник — «http://www.machinelearning.ru/wiki/index.php?title=%D0%A1%D1%82%D0%B0%D1%82%D0%B8%D1%81%D1%82%D0%B8%D1%87%D0%B5%D1%81%D0%BA%D0%B8%D0%B9_%D0%B0%D0%BD%D0%B0%D0%BB%D0%B8%D0%B7_%D0%B4%D0%B0%D0%BD%D0%BD%D1%8B%D1%85_%28%D0%BA%D1%83%D1%80%D1%81_%D0%BB%D0%B5%D0%BA%D1%86%D0%B8%D0%B9%2C_%D0%9A.%D0%92.%D0%92%D0%BE%D1%80%D0%BE%D0%BD%D1%86%D0%BE%D0%B2%29/2013»