Математическая статистика. Проверка гипотезы о виде закона распределения. Практикум
.pdf
2) принять в качестве оценки параметра показательного распределения величину, обратную выборочной средней:
* 1 ; xв
3) найти вероятность попадания X в частичные интервалы по формуле
Pi P xi X xi 1 e xi e xi 1 ; 4) вычислить теоретические частоты
ni nPi ,
где n l ni – объём выборки;
(1.7)
[xi ; xi 1 )
(1.8)
i 1
5)сравнить эмпирические и теоретические частоты с помощью критерия Пирсона.
При этом число степеней свободы k l 1 r l 2 , так как пока-
зательное распределение определяется одним параметром и этот параметр оценивается по выборке, т.е. r 1; l – число интервалов.
1.3.3. Проверка гипотезы о равномерном распределении генеральной совокупности
Для этого надо:
1) оценить параметры a и b – концы интервала, в котором наблюдались возможные значения X, по формулам
a* x S |
3 , b* x S 3 , |
(1.9) |
где a* и b* – оценки параметров;
2) найти дифференциальную функцию (плотность вероятности) предполагаемого распределения
f x |
1 |
; |
(1.10) |
|
b* a* |
||||
|
|
|
3) найти теоретические частоты:
11
|
n1 nP1 |
n f |
x x a* n |
|
1 |
|
x1 a* ; |
|||||||||
|
* |
a |
* |
|||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
b |
|
|
||
|
|
|
|
n |
|
|
|
xi xi 1 , |
i 2, 3, 4, ..., s 1 ; (1.11) |
|||||||
|
|
|
|
|
|
|||||||||||
n2 |
n3 |
... nl 1 |
b* a* |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
nl n |
|
|
1 |
|
|
b* xl 1 ; |
|
|
||||
|
|
|
|
|
b |
* |
|
* |
|
|
||||||
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|||
4) сравнить эмпирические и теоретические частоты с помощью критерия Пирсона, приняв число степеней свободы k l 3 , где l – число интервалов, на которые разбита выборка.
12
2. ПРИМЕРЫ ВЫПОЛНЕНИЯ ИНДИВИДУАЛЬНЫХ ЗАДАНИЙ
Пример 1
Имеютсярезультаты50 измеренийнекоторойслучайнойвеличиныХ: 36, 30, 34, 32, 40, 31, 32, 33, 30, 32, 35, 37, 35, 30, 32, 30, 31, 38, 31, 30, 32, 42, 30, 35, 37, 32, 31, 45, 30, 33, 31, 35, 37, 30, 33, 32, 41, 38, 39,
36, 35, 33, 30, 43, 34, 45, 32, 39, 31, 34.
Принять l = 5.
1. Составим вариационный ряд.
Пусть значение случайной величины xj встречается nj раз.
j |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
xj |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
45 |
nj |
9 |
6 |
8 |
4 |
3 |
5 |
2 |
3 |
2 |
2 |
1 |
1 |
1 |
1 |
2 |
m
nj 50 n.
j 1
Составим интервальный вариационный ряд.
Наименьшеезначение варианты xmin x1 30 , наибольшее значение
x |
x |
45 . Длина частичного интервала |
h |
xm x1 |
45 30 3. |
|||||||
|
||||||||||||
max |
m |
|
|
|
|
|
|
|
|
l |
5 |
|
Расчет приведен в табл. 2.1. |
|
|
|
|
||||||||
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Таблица 2.1 |
|
|
|
|
|
|
|
|
||||||
|
|
|
Частичный |
Середина частичного |
Сумма частот вариант |
|
||||||
Номер |
|
интервала |
|
интервала |
Плотность |
|||||||
|
интервал |
|
|
|
n j |
|
ni |
|||||
|
|
xi xi 1 |
|
|
|
|
частоты |
|||||
интервала i |
|
[xi ; xi+1) |
xi* |
|
|
ni |
|
h |
||||
|
|
|
|
|||||||||
|
|
|
|
2 |
|
|
|
[xi ;xi 1) |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
[30; 33) |
31,5 |
|
|
23 |
7,67 |
|
|||
|
2 |
|
[33; 36) |
34,5 |
|
|
12 |
4 |
|
|||
|
3 |
|
[36; 39) |
37,5 |
|
|
|
7 |
|
2,33 |
|
|
|
4 |
|
[39; 42) |
40,5 |
|
|
|
4 |
|
1,33 |
|
|
|
5 |
|
[42; 45) |
43,5 |
|
|
|
4* |
1,33 |
|
||
____________
* Интервалы [42; 45) и [45; ) объединены.
13
2. Построим эмпирическую функцию распределения.
Для того чтобы построить эмпирическую функцию распределения
F xi nni ,
можно считать, что все значения непрерывной случайной величины, попавшие в интервал [xi; xi+1), совпадают со значением xi :
F x0 P X 30 0 ;
F x |
P X 33 23 0,46 ; |
|
|||
1 |
|
50 |
|
|
|
|
|
|
|
||
F x |
P X 36 23 12 35 0,70 ; |
||||
2 |
|
|
50 |
50 |
|
|
|
|
|
||
F x |
P X 39 35 7 |
42 0,84 ; |
|||
3 |
|
|
50 |
50 |
|
|
|
|
|
||
F x |
P X 42 42 4 |
46 0,92 ; |
|||
4 |
|
|
50 |
50 |
|
|
|
|
|
||
F x |
P X 45 n nx 45 50 2 0,96 ; |
||||
5 |
|
|
n |
|
50 |
|
|
50 |
|
||
F x P X 45 |
1. |
|
|
||
Таким образом: |
50 |
|
|
|
|
|
|
|
|
||
|
|
|
0 |
при |
x 30; |
|
|
|
|
при |
30 x 33; |
|
|
0,46 |
|||
|
|
|
|
при |
33 x 36; |
|
F x |
0,70 |
|||
|
0,84 |
при |
36 x 39; |
||
|
|
0,92 |
при |
39 x 42; |
|
|
|
|
|
при |
42 x 45; |
|
|
0,96 |
|||
|
|
|
|
при |
x 45. |
|
|
1,00 |
|||
Построим график эмпирической функции распределения (рис. 2.1) и гистограмму (рис. 2.2).
14
F(x)
1,00
0,96
0,92
0,84
0,70
0,46
0 |
30 |
33 |
36 |
39 |
42 |
45 |
x |
Рис. 2.1. График эмпирической функции распределения к примеру 1
ni h
7,67
4,00
2,33
1,33
0 |
30 |
33 |
36 |
39 |
42 |
45 |
x |
Рис. 2.2. Гистограмма к примеру 1
3. Рассчитаем точечные оценки для математического ожидания и дисперсии.
Найдём точечные оценки для интервального вариационного ряда:
15
|
|
|
|
|
|
x |
l |
ni xi* |
|
|
||
|
|
|
|
|
|
i 1 |
|
|
, |
|
||
|
|
|
|
|
|
n |
|
|||||
|
|
|
|
|
|
|
|
|
|
|||
где |
x |
|
xi xi 1 |
– серединные значения вариант в интервалах; |
||||||||
|
||||||||||||
|
|
i |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
x |
1 |
23 31,5 12 34,5 7 37,5 4 40,5 4 43,5 173750 |
34,74 . |
|||||||||
50 |
||||||||||||
|
|
|
|
|
|
Dв* |
ni xi2 |
|
x 2 ; |
|
||
|
|
|
|
|
|
n |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||
Dв* 501 23 31,52 12 34,52 7 37,52 4 40,52 4 43,5234,742 61 50078,5 1206,87 14,7.
Найдём исправленную выборочную дисперсию, которая служит несмещённой оценкой генеральной дисперсии:
S*2 nn 1 Dв* 5049 14,7 15 .
4. Выберем закон распределения и по критерию Пирсона проверим гипотезу о его согласованности с результатами наблюдений при уровне значимости 0,05 .
По виду гистограммы – распределение показательное, сдвинутое на 30 единиц вправо по оси 0x (см. рис. 2.2).
Получим несдвинутое показательное распределение (прил. 1). Вычитая из всех членов вариационного ряда число 30, т.е. переходя от случайнойвеличиныX кновойвеличинеU, получаем ui xi 30 (табл. 2.2).
Таблица 2.2
Номер |
[ui; ui 1) |
ni |
n j |
Относительные частоты |
wi ni |
|
интервала i |
|
|
|
[ui ;ui 1) |
|
n |
|
|
|
|
|
|
|
1 |
[0; |
3) |
|
23 |
0,46 |
|
2 |
[3; |
6) |
|
12 |
0,24 |
|
3 |
[6; |
9) |
|
7 |
0,14 |
|
4 |
[9; 12) |
|
4 |
0,08 |
|
|
5 |
[12; |
15) |
|
4 |
0.08 |
|
16
l wi 1. i 1
Среднее выборочное несдвинутого распределения u x 30 34,74 30 4,74.
Если случайная величина распределена по показательному закону с параметром , то её математическое ожидание M U 1 . Так как
M U u 4,74 , то
|
1 |
|
1 |
0,211. |
M U |
4,74 |
Проверим гипотезу о показательном распределении с помощью критерия согласия Пирсона (критерий 2 – «хи квадрат»):
|
l |
n |
nP |
2 |
|
набл2 |
|
i |
nP |
i |
, |
|
i 1 |
|
i |
|
|
где ni эмпирические частоты;
ni nPi – теоретические частоты.
Вычислим теоретические вероятности Pi попадания случайной величины U в каждый интервал по формуле
P u b e e b ,
где и b концы интервала.
P1 P 0 u 3 e 0 e 0,211 3 1 0,5309 0, 4691;
P2 P 3 u 6 e 0,211 3 e 0,211 6 0,5309 0,2819 0,2490 ; P3 P 6 u 9 e 0,211 6 e 0,211 9 0,2819 0,1497 0,1322 ;
P4 P 9 u 12 e 0,211 9 e 0,21112 0,1497 0,0,0795 0,0702 ;
P5 P 12 u e 0,21112 e 0,0795 0 0,0795 .
Интервал 12 u 15 объединён с интервалом 15 u .
17
Для вычисления наблюдаемого значения критерия Пирсона набл2 составим расчётную табл. 2.3.
Таблица 2.3
Номер |
Интер- |
Часто- |
Гипотетиче- |
Теоретче- |
n nP |
n nP |
|
2 |
n nP |
2 |
|||
интер- |
вал |
та ni |
ская вероят- |
ская час- |
i |
i |
|
|
|||||
|
|
i |
i |
|
|
i |
i |
|
|||||
вала i |
|
|
ность Pi |
тота nPi |
|
|
|
|
|
|
|
nPi |
|
1 |
[0; 3) |
23 |
0,4691 |
23,455 |
–0,455 |
0,2070 |
|
|
0,0088 |
|
|||
2 |
[3; 6) |
12 |
0,2490 |
12,450 |
– 0,450 |
0,2025 |
|
|
0,0163 |
|
|||
3 |
[6; 9) |
7 |
0,1322 |
6,610 |
0,390 |
|
0,1521 |
|
|
0,0230 |
|
||
4 |
[9; 12) |
4 |
0,0702 |
3,510 |
0,490 |
|
0,2401 |
|
|
0,0684 |
|
||
5 |
[12; ) |
4 |
0,0795 |
3,975 |
0,025 |
|
0,0006 |
|
|
0,0002 |
|
||
|
– |
50 |
1 |
– |
– |
|
|
– |
|
|
набл2 |
0,1167 |
|
Сравним эмпирические и теоретические частоты с помощью критерия Пирсона.
|
2 |
2 |
Критическую величину кр |
1 k найдём по табл. П3 крити- |
|
ческих точек распределения |
2 |
(прил. 3). Уровень значимости задан: |
0,05 . |
. |
|
|
|
|
Подсчитаем число степеней свободы k l 1 r ,
где l – число групп (частичных интервалов) выборки;
r – число параметров предполагаемого показательного распределения.
Так как l = 5, r = 1 (подсчитан всего 1 параметр ), то k 5 1 1 3 .
Потабл. П3 (прил. 3) находимкритическуюточкуправостороннейкритической области кр2 0,952 3 7,8 . Так как набл2 0,1167 кр2 7,8 ,
то нет оснований отвергнуть нулевую гипотезу о показательном распределении генеральной совокупности.
Рассмотрим конкурирующую гипотезу H1. Предположим, что полученное эмпирическое распределение случайной величины похоже на нормальное (прил. 1). Проверим гипотезу о нормальном распреде-
лении с математическим ожиданием |
M U u 4,74 и средним |
квадратическим отклонением S |
S 2 15 3,87 . |
18
Вычислим гипотетические вероятности Pi попадания случайной
величины в интервалы ai ;bi |
по формуле |
|
|
|
|
|||
b M U |
a M |
U |
|
|||||
P a U b |
|
|
|
|
|
|
|
, |
|
|
|
|
|||||
|
|
|
|
|
|
|||
где x функция Лапласа.
Значения функции Лапласа (интеграла вероятностей) приведены в табл. П2 (прил. 2). Параметры распределения М(U) и σ заменяем на
найденные их оценки M U 4,74 , 3,87 . |
|
|
|||||
P P U |
3 3 4,74 0,45 |
||||||
1 |
|
|
|
|
|
|
|
|
3,87 |
|
|
|
|
|
|
|
0,1736 0,5 0,3264 ; |
|
|
|
|||
P P 3 U 6 |
|
6 4,74 |
|
3 |
4,74 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
3,87 |
|
|
3,87 |
|
0,3255 0,45 0,1275 0,1736 0,3011; |
|||||||
P P 6 U 9 |
|
9 4,74 |
|
6 |
4,74 |
|
|
3 |
|
|
|
|
|
|
|
|
|
|
3,87 |
|
|
3,87 |
|
1,1 0,3255 0,3643 0,1275 0,2368; |
|||||||
P P 9 |
U 12 |
|
12 4,74 |
9 4,74 |
|
||
4 |
|
|
3,87 |
|
|
3,87 |
|
|
|
|
|
|
|
||
1,876 1,1 0,4696 0,3643 0,1053;
P5 P 12 U 12 4,743,87
0,5 1,876 0,5 0,4696 0,0304.
Для определения набл2 составим табл. 2.4.
Сравним эмпирические и теоретические частоты с помощью критерия Пирсона, приняв число степеней свободы k l 1 r 5 1 2 2 , так как нормальное распределение зависит
от двух r 2 параметров M(X) и . По табл. П3 (прил. 3) находим
19
критическую |
|
точку |
правосторонней |
критической |
области |
|||||||||||
кр2 0,952 |
2 5,99 . |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2.4 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Номер |
Интер- |
Час- |
Гипотети- |
Теорети- |
n nP |
ni nPi |
2 |
|
ni nPi |
2 |
||||||
интер- |
вал |
тота |
ческая ве- |
ческая |
|
i |
i |
|
|
|
||||||
вала i |
|
|
|
ni |
роятность Pi |
частота nPi |
|
|
|
|
|
|
nPi |
|
||
1 |
[0; 3) |
23 |
0,3264 |
16,320 |
|
6,680 |
44,622 |
|
|
2,734 |
|
|||||
2 |
[3; 6) |
12 |
0,3011 |
15,055 |
|
–3,055 |
9,333 |
|
|
0,620 |
|
|||||
3 |
[6; 9) |
|
7 |
0,2368 |
11,840 |
|
–4,840 |
23,426 |
|
|
1,978 |
|
||||
4 |
[9; 12) |
|
4 |
0,1053 |
5,265 |
|
–1,265 |
1,600 |
|
|
0,304 |
|
||||
5 |
[12; ) |
|
4 |
0,0304 |
1,520 |
|
2,480 |
6,150 |
|
|
4,046 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
– |
50 |
|
1 |
– |
|
|
– |
– |
|
набл2 |
9,682 |
|||
Так как набл2 |
9,682 кр2 |
5,99 , |
то наблюдаемое значение набл2 |
|||||||||||||
превышает критическое значение кр2 . Это означает, что результаты
эксперимента не согласуются с гипотезой о нормальном распределении случайной величины.
Пример 2
Имеютсярезультаты60 измеренийнекоторойслучайнойвеличиныХ: 4,1; 3,4; 5,7; 5,1; 4,5; 6,8; 5,0; 3,8; 6,1; 6,3; 6,9; 5,6; 6,3; 5,9; 3,8; 5,6; 6,0; 4,6; 3,4; 6,9; 3,5; 6,4; 3,9; 5,8; 3,7; 3,9; 3,3; 5,9; 4,9; 6,8; 6,4; 4,7; 6,2; 3,3; 5,2; 4,4; 6,7; 3,6; 4,0; 6,6; 6,5; 5,5; 4,2; 6,1; 4,8; 5,8; 3,5; 6,7;
5,3; 4,4; 3,6; 6,6; 5,9; 5,2; 3,7; 6,5; 4,9; 4,3; 5,7; 3,8.
Принять l = 6.
1. Составим вариационный ряд:
j |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
xj |
3,3 |
3,4 |
3,5 |
3,6 |
3,7 |
3,8 |
3,9 |
4,0 |
4,1 |
4,2 |
4,3 |
4,4 |
4,5 |
4,6 |
4,7 |
4,8 |
4,9 |
5,0 |
nj |
2 |
2 |
2 |
2 |
2 |
3 |
2 |
1 |
1 |
1 |
1 |
2 |
1 |
1 |
1 |
1 |
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j |
19 |
20 |
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
xj |
5,1 |
5,2 |
5,3 |
5,5 |
5,6 |
5,7 |
5,8 |
5,9 |
6,0 |
6,1 |
6,2 |
6,3 |
6,4 |
6,5 |
6,6 |
6,7 |
6,8 |
6,9 |
nj |
1 |
2 |
1 |
1 |
2 |
2 |
2 |
3 |
1 |
2 |
1 |
2 |
2 |
2 |
2 |
2 |
2 |
2 |
nj 60 n .
Составим интервальный вариационный ряд.
20
