- •Работа № 1. Выборки и их представление Основные понятия
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •Работа №2. Предельные теоремы.
- •Теорема Бернулли
- •1. Выполнение в пакете Statistica
- •2) Выполнение в пакете spss
- •Сжатие распределения с ростом числа слагаемых.
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •Усиленный закон больших чисел.
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss.
- •Теорема Гливенко основная теорема статистики
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss.
- •Центральная предельная теорема
- •Одинаково распределенные слагаемые.
- •1) Выполнение в пакете statistica
- •2. Выполнение в пакете spss.
- •Различно распределенные слагаемые
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •Работа № 3. Оценки
- •Постановка конкретной задачи.
- •Теоретическое сравнение оценок
- •Статистическое сравнение оценок
- •Задание для самостоятельной работы
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •Работа n4. Доверительные границы и интервалы
- •Определения и построение интервалов
- •Уровень доверия
- •Интервалы для параметров нормального распределения
- •Задание на самостоятельную работу
- •1. Выполнение в пакете statistica
- •2. Выполнение в пакете spss
- •Работа № 5. Критерий "хи-квадрат" проверки гипотез
- •Проверка простой гипотезы о вероятностях
- •Проверка сложной гипотезы о вероятностях
- •Проверка гипотезы о типе распределения
- •1. Выполнение в пакете statistica
- •Примеры проверки простой гипотезы о распределении
- •1. Выполнение в пакете statistica
- •Проверка гипотезы о независимости признаков (таблица сопряженности признаков)
- •1. Выполнение в пакете statistica
- •Проверка гипотезы об однородности выборок
- •1. Выполнение в пакете statistica
- •Задание
- •Работа № 6. Различение двух простых гипотез. Различение при фиксированном объеме наблюдений
- •Последовательное различение двух простых гипотез (последовательный анализ Вальда).
- •Задание
- •1. Выполнение в пакете staTiStica
- •Оценка вкладов. Если гипотеза ha отклоняется, следует оценить вклады aj уровней фактора; оценка
- •1. Выполнение в пакете statistica
- •Работа № 8. Линейный регрессионный анализ
- •1. 1. Простая линейная регрессия
- •2.2. Множественная регрессия
- •3. 3. Нелинейная зависимость
- •4.4. Нелинейная зависимость (обобщение)
- •Литература
Проверка сложной гипотезы о вероятностях
Пусть A1, ...,Am - m исходов некоторого опыта, n - число независимых повторений опыта, 1,...,m - числа появлений исходов. Проверяемая гипотеза Н предполагает, что вероятности исходов P(Ai) являются известными функциями pi(a) k-мерного параметра a = (a1,...,ak), т.е.
Н: Р(Аi) = pi(a), i = 1, ..., m,
но значение а неизвестно.
Для проверки гипотезы Н определим статистику
(6)
По теореме Фишера, если Н верна, то при n распределение статистики Х2 асимптотически подчиняется распределению хи-квадрат с числом степеней свободы f = m -1- k, и потому отклоняем Н, если
h,
(7)
где h = Q(1-, f) - квантиль уровня 1- распределения хи-квадрат с числом степеней свободы f; такой порог обеспечивает выбранный уровень вероятности P(отклонить Н / Н) ошибки 1-го рода. Если (7) не выполняется, делаем вывод, что наблюдения не противоречат гипотезе. Распределению хи-квадрат с f = m -1- k степенями свободы асимптотически подчиняется также статистика
,
(8)
где
- оценка максимального правдоподобия
для а,
и потому в (7) может быть использована
статистика (8) вместо (6). Процедура (7)
может быть записана иначе: если
P{f2 X2} (9)
то гипотеза Н отклоняется.
Проверка гипотезы о типе распределения
Пусть требуется проверить гипотезу о том, что выборка x1, ..., xn извлечена из совокупности, распределенной по некоторому закону, известному с точностью до k-мерного параметра а=(а1,...,аk). Оказываются теоретически обоснованными следующие действия: разобьем весь диапазон наблюдений на m интервалов, определим значения i -число наблюдений в i-м интервале, получим значение оценки минимизацией (6) или методом максимального правдоподобия, определим вероятности pi( ) попадания в i-й интервал, вычислим (6) или (8) и примем решение по (7).
Пример 1. Проверка нормальности. Проверим гипотезу о нормальном законе распределения размеров головок заклепок, сделанных на одном станке, по выборке объема n = 200; измерения приведены в таблице 1 [1, с.15]. Оценками для а (среднего) и (стандартного отклонения) являются
и
.
Таблица 1.
Диаметры 200 головок заклепок, мм |
|||||||||
13.39 |
13.33 |
13.56 |
13.38 |
13.43 |
13.37 |
13.53 |
13.40 |
13.25 |
13.37 |
13.28 |
13.34 |
13.50 |
13.38 |
13.38 |
13.45 |
13.47 |
13.62 |
13.45 |
13.39 |
13.53 |
13.58 |
13.32 |
13.27 |
13.42 |
13.40 |
13.57 |
13.46 |
13.33 |
13.40 |
13.57 |
13.36 |
13.43 |
13.38 |
13.26 |
13.52 |
13.35 |
13.29 |
13.48 |
13.43 |
13.40 |
13.39 |
13.50 |
13.52 |
13.39 |
13.39 |
13.46 |
13.29 |
13.55 |
13.31 |
13.29 |
13.33 |
13.38 |
13.61 |
13.55 |
13.40 |
13.20 |
13.31 |
13.46 |
13.13 |
13.43 |
13.51 |
13.50 |
13.38 |
13.44 |
13.62 |
13.42 |
13.54 |
13.31 |
13.58 |
13.41 |
13.49 |
13.42 |
13.45 |
13.34 |
13.47 |
13.48 |
13.59 |
13.20 |
14.56 |
13.55 |
13.44 |
13.50 |
13.40 |
13.48 |
13.29 |
13.31 |
13.42 |
13.32 |
13.48 |
13.43 |
13.26 |
13.58 |
13.38 |
13.48 |
13.45 |
13.29 |
13.32 |
13.24 |
13.38 |
13.34 |
13.14 |
13.31 |
13.51 |
13.59 |
13.32 |
13.52 |
13.57 |
13.62 |
13.29 |
13.23 |
13.37 |
13.64 |
13.30 |
13.40 |
13.58 |
13.24 |
13.32 |
13.52 |
13.50 |
13.43 |
13.58 |
13.63 |
13.48 |
13.34 |
13.37 |
13.18 |
13.50 |
13.45 |
13.60 |
13.38 |
13.33 |
13.57 |
13.28 |
13.32 |
13.40 |
13.40 |
13.33 |
13.20 |
13.44 |
13.34 |
13.54 |
13.40 |
13.47 |
13.28 |
13.41 |
13.39 |
13.48 |
13.42 |
13.46 |
13.28 |
13.46 |
13.37 |
13.53 |
13.43 |
13.30 |
13.45 |
13.40 |
13.45 |
13.40 |
13.33 |
13.39 |
13.56 |
13.46 |
13.26 |
13.35 |
13.42 |
13.36 |
13.44 |
13.41 |
13.43 |
13.51 |
13.51 |
13.24 |
13.34 |
13.28 |
13.37 |
13.54 |
13.43 |
13.35 |
13.52 |
13.23 |
13.48 |
13.48 |
13.54 |
13.41 |
13.51 |
13.44 |
13.36 |
13.36 |
13.53 |
13.44 |
13.69 |
13.66 |
13.32 |
13.26 |
13.51 |
13.38 |
13.46 |
13.34 |
