7. Проверка, что выборка осуществлялась по
случайному закону
Будем использовать критерий согласия Пирсона
. (1)
Если, как это обычно имеет место,
,
то указанную формулу можно преобразовать к виду
более удобному для вычислений
. (2)
Проверяемая гипотеза состоит в том, что выбор 150 изделий из нормальной генеральной совокупности с μ = 60 и σ = 10 был произведён по случайному закону. Так как известны параметры нормальной генеральной совокупности, то критерий χ2 может только проверить случайный характер выборки.
В качестве уровня значимости выберем α = 0,05.
Для того, чтобы получить визуальное представление о степени соответствия нашей выборки нормальной кривой, воспользуемся гистограммой частот 150 изделий (см. рис. Б) с наложенной на неё нормальной кривой с параметрами μ = 60 и σ = 10. Вычисленные значения ординат нормальной кривой приведены в таблице 7.
Таблица 7
Ординаты нормальной кривой: μ=60; σ=10; n=150; ∆Х=8
|
|
Х
|
Z (x/σ) |
T (z2/2) |
|
Yс
|
|
28 |
-32 |
-3,2 |
5,12 |
0,000598 |
00,29 |
|
36 |
-24 |
-2,4 |
2,88 |
0,056140 |
02,69 |
|
44 |
-16 |
-1,6 |
1,28 |
0,278040 |
13,31 |
|
52 |
0-8 |
-0,8 |
0,32 |
0,726150 |
34,76 |
|
60 |
000 |
00,0 |
0 |
1,000000 |
47,87 |
|
68 |
008 |
00,8 |
0,32 |
0,726150 |
34,76 |
|
76 |
016 |
01,6 |
1,28 |
0,278040 |
13,31 |
|
84 |
024 |
02,4 |
2,88 |
0,056140 |
02,69 |
|
92 |
032 |
03,2 |
5,12 |
0,000598 |
00,29 |
- середина
теоретического интервала
.
Умножение на
,
а не на
вызвано
тем, что площадь гистограммы должна быть равна 100, а не 1 и являться суммой
площадей прямоугольников с основанием ΔХ.
Увеличение числа интервалов и изменение их границ вызвано симметричностью нормальной кривой относительно μ.
При рассмотрении рис. Г может создаться впечатление, что между нормальной кривой и гистограммой наблюдается несоответствие. Однако при этом следует помнить, что выборка содержит только 150 изделий, и наличие даже существенного расхождения не следует считать слишком неожиданным.
В таблице 8 приведены данные наблюдений с
правилами распределения площади под кривой нормального распределения.
Фактический процент наблюдения находится путём простого подсчёта частот mi
статистического ряда (таблица 3), попадающих в заданные интервалы
±
zs и соотнесённых с объёмом выборки n = 150. В нашем случае наблюдается хорошее
согласие процентов наблюдений.
Таблица 8
Сопоставление опытных данных с правилами распределения площади под кривой нормального распределения
|
z |
|
Теоретический процент наблюдений |
Фактический процент наблюдений |
|
1 |
от 49,83 до 71,61 |
68,27 |
70,00 |
|
2 |
от 38,94 до 82,50 |
95,45 |
95,33 |
|
3 |
от 28,05 до 93,99 |
99,73 |
100,00 |
Рис. Г. Нормальная кривая и гистограмма 150 изделий, взятых по случайному закону из приближённо нормальной генеральной совокупности (μ = 60, σ = 10)
Теоретические частоты
для
рассматриваемого распределения 150 изделий, взятых из генеральной совокупности
с μ
= 60 и σ
= 10, приведены в таблице 9. Порядок вычисления объясняют заголовки каждого из
столбцов этой таблицы.
Очевидно, что значений
1-го
столбца следует выбирать нечётное количество с учётом величины μ
так,
чтобы серединное значение
равнялось значению
μ.
Таблица 9
Вычисление теоретических частот для 150 изделий, взятых по случайному закону из приближённо нормальной генеральной совокупности (μ=60; σ=10; n=150;ΔX=8)
|
|
|
X
( |
|
Интегральные относительные частоты |
Первая разность Δ |
m* (n·Δ) |
|
… |
-∞ |
-∞ |
-∞ |
-0,50000 |
0,00226 |
0,34 |
|
|
32 |
-28 |
-2,8 |
-0,49744 |
|
|
|
36 |
|
|
|
|
0,02019 |
3,03 |
|
|
40 |
-20 |
-2,0 |
-0,47725 |
|
|
|
44 |
|
|
|
|
0,09232 |
13,85 |
|
|
48 |
-12 |
-1,2 |
-0,38493 |
|
|
|
52 |
|
|
|
|
0,22951 |
34,43 |
|
|
56 |
-4 |
-0,4 |
-0,15542 |
|
|
|
60 |
|
|
|
|
0,31084 |
46,63 |
|
|
64 |
4 |
0,4 |
0,15542 |
|
|
|
68 |
|
|
|
|
0,22951 |
34,43 |
|
|
72 |
12 |
1,2 |
0,38493 |
|
|
|
76 |
|
|
|
|
0,09232 |
13,85 |
|
|
80 |
20 |
2,0 |
0,47725 |
|
|
|
84 |
|
|
|
|
0,02019 |
3,03 |
|
|
88 |
28 |
2,8 |
0,49744 |
|
|
|
… |
|
|
|
|
0,00226 |
0,34 |
|
|
∞ |
∞ |
∞ |
0,50000 |
|
|
|
Σ |
… |
… |
… |
… |
0,99940 |
149,93 |
Значения в первой и последней строчках 2, 3 и
4-го столбцов равны соответственно -∞ и +∞, так как нормальное
распределение теоретически простирается от -∞ до +∞. Числа,
указанные в 5-м столбце, взяты из таблицы 2.3, приведённой в Приложении 2.
Наличие знака "минус" у первых цифр этого столбца, соответствующих
значениям
, меньшим нуля,
обусловлено тем, что при Хр < μ
< 0.
В 6-м столбце приведены первые разности.
Так, например,
|
Хс |
Первая разность |
Относительная частота |
|
52 |
-0,15542 - (-0,38493) |
= 0,22951 |
|
60 |
0,15542 - (-0,15542) |
= 0,31084 |
|
68 |
0,38493 - 0,15542 |
= 0,22951 |
Вычисленное по формуле (1) значение χ2
приведено
в таблице 10. Как следует из этой таблицы, значение χ2
= 1,2433. В этой таблице приведено семь групп
(интервалов), то есть на две группы меньше, чем в таблице 9. Сокращение числа
групп осуществлено за счёт объединения первых двух, а также двух последних
групп. Это сделано потому, что значения m* в первых двух и последних двух группах
весьма малы. Очень часто применяется правило, указывающее, что критерий χ2
может
применяться в тех случаях, когда каждая теоретическая группа содержит по
крайней мере пять наблюденных значений, а общее количество наблюденных значений
составляет по крайней мере 50. Введение этих ограничений имеет вполне
определённую цель: гарантировать, что распределение наблюденных значений mi
относительно значений теоретических частот
будет
настолько близко к нормальному, что применение при оценке вероятностей таблиц
для χ2
будет
вполне обоснованным.
Хотя таблица 10 содержит семь групп, имеется
только шесть степеней свободы (ν = S - 1 = 7 - 1 = 6), так
как на теоретические частоты накладывается одно ограничение:
.
В таблице для χ2 из
Приложения 3 для ν = 6 находим
При
этом область принятия будет определяться соотношением χ2
< 12,592, а область отклонения - соотношением χ2
≥
12,592.
Таблица 10
Вычисление χ2 для распределения частот при ширине группового интервала, равной 8, и выборке из 150 изделий, взятых по случайному закону из нормальной генеральной совокупности:μ=60; σ =10
|
Границы теоретических интервалов |
m |
m* |
m - m* |
(m - m*)2 |
|
|
-∞ ÷ 40 |
4 |
3,37 |
0,63 |
0,3969 |
0,1178 |
|
40 ÷ 48 |
13 |
13,85 |
-0,85 |
0,7225 |
0,0522 |
|
48 ÷ 56 |
32 |
34,43 |
-2,43 |
5,9049 |
0,1715 |
|
56 ÷ 64 |
46 |
46,63 |
-0,63 |
0,3969 |
0,0085 |
|
64 ÷ 72 |
35 |
34,43 |
0,57 |
0,3249 |
0,0094 |
|
72 ÷ 80 |
15 |
13,85 |
1,15 |
1,3225 |
0,0955 |
|
80 ÷ ∞ |
5 |
3,37 |
1,63 |
2,6569 |
0,7884 |
|
Всего |
150 |
149,93 |
0,07 |
… |
1,2433 |
Так как вычисленное значение χ2
составляет
1,2433, то оно попадает в область принятия, в связи с чем нет оснований для
отклонения гипотезы о том, что выборка осуществлялась по случайному закону.
8. Проверка гипотезы о нормальном распределении
генеральной совокупности
Если параметры распределения неизвестны, но можно считать, что выборка взята по случайному закону, то может возникнуть желание проверить гипотезу о том, что генеральная совокупность является нормальной. При этом снова примем α = 0,05.
Так как параметры распределения неизвестны,
используем совместные оценки максимума правдоподобия, исчисленные на основании
группировки (таблица 4). Такими оценками являются
для
μ
и
s2 для σ2.
По формулам
и
,
находим:
В таблице 11 показано вычисление теоретических
частот, а из таблицы 12 по формуле (2) находим χ2 = 0,59.
Таблица 11
Вычисление нормальных частот для 150 изделий,
взятых по случайному закону из генеральной совокупности с неизвестными
параметрами:
=60,72; s=10,89
|
|
|
X
( |
х /s |
Интегральные относительные частоты |
Δ |
m* (n·Δ) |
|
… |
-∞ |
-∞ |
-∞ |
-0,50000 |
0,02872 |
004,31 |
|
|
40 |
-20,72 |
-1,903 |
-0,47128 |
|
|
|
44 |
|
|
|
|
0,09228 |
013,84 |
|
|
48 |
-12,72 |
-1,168 |
-0,37900 |
|
|
|
52 |
|
|
|
|
0,21260 |
031,89 |
|
|
56 |
-4,72 |
-0,433 |
-0,16640 |
|
|
|
60 |
|
|
|
|
0,28431 |
042,66 |
|
|
64 |
3,28 |
0,301 |
0,11791 |
|
|
|
68 |
|
|
|
|
0,23292 |
034,94 |
|
|
72 |
11,28 |
1,036 |
0,35083 |
|
|
|
76 |
|
|
|
|
0,11081 |
016,62 |
|
|
80 |
19,28 |
1,770 |
0,46164 |
|
|
|
… |
|
|
|
|
0,03836 |
005,74 |
|
|
∞ |
∞ |
∞ |
0,50000 |
|
|
|
Σ |
… |
… |
… |
… |
1,00000 |
150,00 |