книги / Прикладная статистика
..pdfзрения принятия выдвинутой гипотезы о законе распределения генеральной совокупности X произошло достоверное событие. Гипотеза считается не противоречащей опытным данным и принимается. Если же оказывается, что > Х2^* то выдви нутая гипотеза отвергается, считается, что она противоречит опытным данным.
В пашем случае к = 11, г - 11-1 -1 = 9 (по выборке был опреде лен один параметр— X). Если положить р =0,95 ( а =0,05 — наибо лее употребительное значение уровня значимости), то по табли це распределения х1накодим, что = 16,92. Между тем = = 2,45 < х7,р- Так что мы можем считать, что случайная величина X имеет показательное распределение с параметром X = 0,00115. Если бы мы объединили три последних интервала в одци, то имс-
либы:г = 9 - 1 |
- 1 - 7;у* |
= 14,07;х' = 2,33< х2 . |
Случайная |
величина х |
называется критерием х • Критерии |
X2был предложен Карлом Пирсоном в 1900 г. До этого времени совпадение экспериментальных результатов с теоретический оценивалось по тому, как онн выглядят на графике.
Мам осталось ответить на вопросы, поставленные в пункте 6.1. Мы считаем справедливым показательный закон с парамет ром X = 0,00] 15. Следовательно, ЩХ) = /А » 870 (ч).
р{Х> 1000) = е-х"|™ - е <- = <г,-,3»0,32; р(Х< 200) - е° - = I - «0,21.
6 . 3 . Д Р У Г И Е П Р И М Е Р Ы
6.3.1. Проверка гипотезы о нормальном законе распределения
Заказчику необходимы валы с допустимым отклонением диаметра от номинального размера ±0,1 мкм. Прежде чем по купать партию из 1000 валов, он приобрел партию из 200 валов, чтобы оценить ожидаемую долю неподходящих ему изделий. Результаты измерений представлены в табл. 6.4.
111
Таблица 6.4
200откл<тстЬ| диаметрадала от номинального размера (мкм)
Середина интервала |
-0,14 |
-0.12 |
-0,10 |
-0,08 |
-0,06 -0,04 -0,02 |
|||
Частота |
3 |
8 |
и |
20 |
27 |
36 |
29 |
|
Середина шпервала |
0,00 |
0,02 |
0,04 |
0.06 |
0,08 |
0,10 0.12 |
||
Частота |
18 |
17 |
17 |
8 |
А |
1 |
1 |
|
Здесь Ь - |
0,02 мкм; п ® 200; мД = 4. |
|
|
|
|
|||
Гистограмма показана на рссс.6.3. Высоты гистограммы та |
||||||||
ковы: |
2; ЬЛ- 2,75; А, = 5; /г, = 5,75; А,= 9; А, = 7,25; А,= 4,5; |
|||||||
А,= 0,75; |
||||||||
К ш !,и = 4 '251А/г= 2‘> |
= 1: К = К |
= °>25- |
|
|||||
Числовые характеристики: л = -0,028 (мкм); 5= 0,05 (мкм). |
|
|||||||
Судя по гистограмме, можно заключить, что случатыя вели- |
||||||||
чина X — отклонение диаметра вала от номинального размера |
— |
имеет нормальное распределение. Функция плотности нормально го закона зависит от двух параметров — аист:
|
]__ |
■(Х-Д)3 |
о х |
в |
2-а* |
л/2гс |
|
Как известно, М(Х) = а9а(Х) = с . Для определения а и о по
ложим, что а = х9о = Отсюда а = -0,03; с = |
0,05 (значение " |
х округлено, исходя из соображений здравого |
смысла). Тогда |
112
|
1 |
-(.г+0,03)* |
|
/ г о - |
е г-о,оо25 = 8 х е - 2001*+0-03)2 |
||
0,05 х -Л л |
|||
|
Значения функции плотности вероятности на границах ин тервалов таковы (табл. 1.5):
|
|
|
|
|
|
|
Таблица 6.5 |
|
** |
-0.15 |
-0,13 |
-о,м |
-0,09 |
-0,07 |
-0,05 |
-0,03 |
-0,01 |
№ |
0,45 |
1.08 |
2,22 |
3,89 |
5,81 |
7,38 |
8,00 |
7,38 |
\ |
0,01 |
0,03 |
0,05 |
0,07 |
0,09 |
0,11 |
0,13 |
- |
м |
5,81 |
3,89 |
2,22 |
1,08 |
0,45 |
0,16 |
0,05 |
- |
График функции плотности вероятности показан на рис. 6.3. Вычислим теоретические вероятности попадания в интер валы. Формула вычисления вероятности попадания в интервал л,) нормально распределенной случайной вел™ ош X такова:
/>(т,_1 < Х < л ,) = ~<»^А/ д - " ]-
где Ф(х) — функция Лапласа.
Значения функции Лапласа приведены в приложении 1. От сюда:
/>|(-0,15<*< "
=-0,477 -(-0,492)^0,0] 5;
^ |
0,05 |
) |
I. |
0,05 |
] |
= -0,445-(-0,477) = 0,032
Дальнейшие вычисления приведены в табл. 6.6.
Из-за того, что значение параметра я случайно совпало
113
с одной из границ» значения вероятностейр4оказались симметрич ныотносительно интервала (-0,05; -0,01). Двапоследних интервала [0,09; 0,11) н [0,11; 0,13) объединены ввиду их малочисленности.
Положим а = 0,05. Число степеней свободы г - 13 - 2 - I - = 10, х2^ = 18,3 > х 2^ = 8,06. Нет оснований отвергнуть выдви нутую нами гипотезу о нормальном законе распределения от клонений диаметра вала от номинального значения.
|
|
|
|
|
|
|
|
|
Та&мца 6.6 |
К г х> |
|
( |
-<Л |
Р* |
|
|
пГ »рй |
(л, -ирг)2 |
|
и |
ф - — |
|
|
* |
"Р. |
||||
|
\ |
У |
/ |
|
|
|
|
||
|
|
|
|
|
|
|
|||
- |
-2.4 |
-0.492 |
|
- |
- |
- |
- |
- |
|
{-0,15;-0,13) |
-2 |
-0,477 |
|
0,015 |
3 |
3 |
0 |
0 |
|
[-0,13;*0,11) |
-1,6 |
-0,445 |
|
0,032 |
6,4 |
8 |
1.6 |
0.4 |
|
[-0,11; -.09) |
-1.2 |
-0.387 |
|
0,058 |
П.6 |
11 |
-0,6 |
0,03 |
|
[0,09;-0,07) |
-0,8 |
-0,288 |
|
0,099 |
19.8 |
20 |
0,2 |
0,002 |
|
(-0,07.-0,05) |
-0,4 |
41,155 |
|
0,133 |
26,6 |
27 |
0,4 |
0,006 |
|
[-0,05;-0,03) |
0 |
0,000 |
|
0,155 |
31 |
36 |
5 |
0,81 |
|
(-0,03,-0,01) |
0,4 |
0,155 |
|
0,155 |
31 |
29 |
-2 |
0ДЗ |
|
[-0,01 ;0,0!) |
0,8 |
0,288 |
|
0,133 |
26,6 |
18 |
-8,6 |
2,78 |
|
[0,01ДОЗ) |
1.2 |
0,387 |
|
0,099 |
19,8 |
17 |
-2,8 |
0,4 |
|
[0,03;0,05) |
1,6 |
0,445 |
|
0,058 |
И.6 |
17 |
5,4 |
2.51 |
|
[0,05;0,07) |
2 |
0,477 |
|
0,032 |
6.4 |
8 |
1,6 |
0,4 |
|
[0,07;0,09) |
2,4 |
0,492 |
|
0,015 |
3 |
4 |
1 |
0.33 |
|
[0.09Д10 |
2.8 |
0,497 |
|
0,005 |
‘■Ч |
ч |
0 |
0 |
|
[0,11,0,13) |
3.2 |
0,499 |
|
0,002 |
0,4/ |
1/ |
0,6 |
0,26 |
|
X |
— |
— |
|
0,991 |
198,2 |
200 |
— |
|
Оценим долю валов, подходящих заказчику. Вероятность того, что диаметр вала соответствует требованиям заказчика равна/?(—0,1 <Х<0,1) =
= Ф{2,6)+Ф(1,4)=0,495+0.419=0,914.
114
В среднем около 9% валов окажутся непригодными для за казчика.
6,3.2. Проверка гипотезы о равномерном законе распределения
В точение 10 часов регистрировали время прибытия машин к бензоколонке (табл. 6.7).
Таблица 6.7
В ремя |
[ 8 - |
1 9 - |
[ 1 0 - |
[ 1 1 - |
[1 2 |
- |
[ 1 3 - |
[ 1 4 - |
[ 1 5 - |
[ 1 6 - |
[ 1 7 - |
пр и б ы ти я |
9 ) |
10) |
П ) |
12) |
13) |
14) |
15) |
16) |
17) |
18) |
|
(часы ) |
|
|
|
|
|
|
|
|
|
|
|
", |
22 |
30 |
22 |
16 |
2 8 |
|
13 |
17 |
2 0 |
17 |
15 |
При уровне значимости а = 0,05 проверить гипотезу о том, что время прибытия машин — случайная величина, имеющая равнодIсрнос распределенне.
Построим гистограмму. Так как п = 200, /г = 1, то высоты гистограммы таковы:
*, = 2?0 =0,П;/'2=^ =0Д5;Ал =0,11;/' ^ 0,С>8;
Л3 = 0,14; Л, = 0,065; А, - 0,085; ЬЙ= 0,1; Н9= 0,085; = 0,075.
Гистограмма приведена на рис. 6.4.
Если мы считаем, что время прибытия машин имеет равно мерное распределение, мы должны определить два параметра (я и Ь) равномерного закона. Как известно, функция плотности вероятности/(х) равномерного закона такова:
1
(■Ь - а )
0,дг«(<х,6)
И 5
Рис. 6.4.
При этап щ х ) = Е ± Ь ; 0 ( х ) = ! ^ А - 1 ф ) = ~ ^ -
Так что для определения а и Ь можно записать два уравне
ния:
я + 6 —
2 = а*;
Ь- а
2 Л
откуда а = х - 5л/3; Ь = х + з4 Ъ .
Но мы поступил! проще и разумнее. Ноша выборка расположе на на интервале (в,18), поэтому положим: а ~ 8, Ь = 18,Дх) = 0,1
(*е(8,18)).
График функции плотности вероятности /(х) также показан на рис. 6.4. Все теоретические вероятностир, одинаковы и равны
— - — = 0,1. Дальнейшие расчеты представлены в табл. 6.8.
|
|
|
|
|
Таблица 6,8 |
|
Л |
Щ |
|
п - щ |
(п,-пр,)'- |
|
|
Щ |
|||
1 |
|
3 |
|
|
|
2 |
4 |
5 |
6 |
||
[8;« |
0.1 |
20 |
22 |
2 |
0,2 |
Р:10) |
0.1 |
20 |
30 |
10 |
5 |
[Ю;Щ |
0.1 |
20 |
22 |
2 |
0.2 |
_ ГЧ:»2) |
Г А 1 ____ |
20 |
16 |
- А |
0.8 |
116
1 |
2 |
3 |
П2.13) |
0.1 |
20 |
[ 13; 1-4) |
0.1 |
20 |
[И :15) |
ОД |
20 |
115,16) |
ОД |
20 |
Ц6;17> |
ОД |
20 |
| 17;I8) |
о д |
20 |
— |
5 р. - I |
1ир = 200 |
4
28
13
17
20
17
15
1 II
Окончание табл. 6.8
5 |
6 |
8 |
5.2 |
- 7 |
2,45 |
- 3 |
0.45 |
0 |
0 |
- з |
0,45 |
-5 |
*.25 |
—Х * = 14 *■1М
Итак, = 14. Найдем ^ . Мы ис определяли по выбор ке параметров закона — ирсмя работы бензоколонки задано за
ранее. Поэтому число степенен свободы г = 10 - |
1 = 9. Тогда |
||||||
Х ^= 16,9 > х \ нхп’ Выдвинутую гипотезу можно принять. |
|
||||||
6.3.3. |
П роверка гипотезы о биномиальном |
|
|||||
|
законе распределения |
|
|
||||
Семь, монет подбрасывались |
1536 раз. Каждый раз отмеча |
||||||
лось число X выпавших гербов (табл. 6.9). |
|
|
|
||||
|
|
|
|
|
|
Таблица6.9 |
|
0 |
1 |
2 |
3 |
4 |
5 |
б |
7 |
12 |
78 |
270 |
456 |
386 |
252 |
69 |
13 |
При уровне значимости а = 0,05 проверить гипотезу о том» что монеты правильные.
Если вес монеты правильные, то вероятность выпадения гер ба длл каждой из них равнар = 0,5. Тогда случайная величина X — число выпавших гербов при бросании семи монет — име ет биномиальное распределение с параметрами п = 1 н р = 0,5. Биномиальное распределение дискретно, поэтому нужно вы числить теоретические вероятности р, каждого из 8 возможных значении случайной величины X. Эти вероятности считают по формуле Бернулли:
117
ЯЙГ= 0)^С 7У ^ = 0,57= 0,0078; Р{Х=\) = С *р> < ?= 7x0,57 = 0,055; Р(Х= 2)= су< ? = 21x0,5'=0,164;
7 ^ = 3 ) = С7У<? =35x0,57=0,273; ДАГ= 4 )= С /рУ = 35x0,5* = 0,273; Р(Х= 5) = С /р У =21x0,57=0,164;
7x0,5'= 0,055; /ф Г = 7)= С /дУ = 0,5’= 0,0078. Теперь можно вычислить математические ожидания чисел появлении каждого из значений случайной величины Л’при. 1536 бросаниях семи монет, сравшггь их с экспериментальными дан
ными и вычислить |
Результаты свсдеЕсы в табл. 6.10. |
|||||
|
|
|
|
|
|
Таблица б. 10 |
л! |
Р/ |
|
>*р< |
иг |
пр, |
0». ~"1>У |
|
пр, |
|||||
0 |
0,4078 |
|
12 |
12 |
0 |
0 |
1 |
0,055 |
|
84 |
78 |
-6 |
0,43 |
2 |
0,164 |
|
252 |
270 |
18 |
1,29 |
3 |
0.273 |
|
420 |
456 |
36 |
3,09 |
4 |
0,273 |
|
420 |
386 |
-34 |
2,75 |
5 |
0,164 |
|
252 |
252 |
0 |
0 |
б |
0,055 |
|
84 |
69 |
-15 |
2,68 |
7 |
0,0078 |
|
12 |
13 |
1 |
0,08 |
— |
Э у 1 |
Тлр= 1536 |
Е » ,-1536 |
— |
Х '„ - 10,32 |
|
Найдем |
В случае дискретной случайной величины при |
подсчете г вместо числа интервалов берут число различных зна чений лгг В нашем случае г - 8 - 1 = 7 , так как ни одного пара метра по выборке мы не находим. Тогда %^ = 14,1> х*пгн= 10,32. Ист оснований опровергнуть гипотезу о правильности монет.
6.3.4. Проверка гипотезы о законе распределения Пуассона
В таблице приведены числа п} участков равной площади (0,25 км2) южной части Лондона, на каждый из которых при ходилось по х{ попадании самолетов-снарядов во время второй мировой войны (табл. 6.11).
118
|
|
|
|
|
|
Таблица 6.Н |
|
0 |
1 |
2 |
3 |
4 |
5 к больше |
". |
229 |
2И |
93 |
35 . |
7 |
1 |
Всего и - 576 участков. При уровне значимости а =* 0,05 про верить гипотезу о том, что случайная величина А"— число са молетов-снарядов, попавших на участок, имеет раслрсдслеЕше Пуассона.
Вероятность того, что случайная величина Х> имеющая рас
пределение Пуассона, примет значение |
равна: |
|
||||||
|
|
|
Р{Х=1) =^ е - \ |
|
|
|
||
где |
А.> 0 — параметр закона, |
1 = 0 ,1 ,2 ....... |
|
М(Х) - X, |
||||
Оценим значение параметра А. но выборке. Так как |
||||||||
то положим X =х, |
|
|
|
|
|
|||
.V= — |
(0*229 + 1x211 + 2x93 + 3к35 + 4x7 + 5х 1)=0,93. |
|||||||
|
576 |
|
|
|
|
|
|
|
Положим к - |
0,93. Теперь можно нейти вероятности р1 = |
|||||||
~ р (Х -г), I = 0,1,2,3,4,5. |
|
|
|
|
|
|||
Р0 = Р (Х = 0) = |
=0,395; Р\ = Н Х = 0 |
= = |
9^67; |
|||||
р г = Р ( Х |
=2) = |
= 0,170’ |
Р з = Р ( Х =3) = |
= 0.053; |
||||
|
|
|
2' р А - |
Р ( Х |
= 4) = 0,01 Ъ |
|
|
|
|
Р $ |
- Р ( Х |
2 5) = 1 - |
р „ - |
р, - р 2 - |
Р з - |
Р а = 0,003. |
|
Остальные вычисления сведены в табл. 6.12. |
|
|||||||
|
|
|
|
|
|
|
Табгшца6Л2 |
|
|
|
Р. |
|
|
|
|
(1Г|-мр.)а |
|
|
|
|
|
|
'ГЧР, |
"Р. |
||
1 |
|
2 |
3 |
|
|
|
|
|
|
|
4 |
5 |
6 |
||||
0 |
0,395 |
227,5 |
|
229 |
1.5 |
0,01 |
119
|
|
|
|
Окхтчаиие табл, б.}2 |
|
1 |
2 |
3 |
4 |
5 |
6 |
1 |
0,367 |
211,4 |
211 |
|
0,001 |
2 |
0,170 |
97,9 |
93 |
-4*9 |
0,25 |
3 |
0,053 |
30,5 |
35 |
4,5 |
0,66 |
4 |
ОД]2 |
6,9] |
71 |
-0,6 |
0,04 |
55 |
0,003 |
1.7 |
>1 |
||
— |
1р г * |
Т п р ,= 576 |
Е»,= 576 |
— |
х'«, =0>96 |
Два последних значения п4и пр пр4и пр5 объединены, что бы обеспечить выполнение условия лр. > 5. Таким образом, ос талось 5 разных значении случайной величины: 0,1, 2, 3 и вес, «по больше или равно 4. Число степеней свобода равно г = 5 -1 - I = 3, так как по выборке было определено значение параметра X. Тог да = 7,8 > хажп “ 0,96. И в этом случае можно считать спра ведливой выдвинутую гипотезу.
6.3.5. Последний пример
Согласно закону Геллпна, предложенному им в 1655 г., веро ятности рождения двоен, троен и четверней есть соответственно р, р \ р \ где р ■— число, постоянное для данной группы населе ния. На основании приведенных ниже данных проверить, вы полняется ли закон Геллина для многоплодных рождений среди японцев н белого населения США. В табл. 6.13 через у, обозначены относительные частоты рождений двоен, троен и четверней соответственно за указанные периоды.
|
|
|
|
Таблица 6.13 |
|
Годы |
Население |
Числорож |
Ъ |
*4 |
|
дений |
|||||
|
|
||||
1922-1936 |
Белое население |
27939615 |
0,0И29 0,0001088 0,00000177 |
||
1926-1931 |
Лловцы |
1226106 |
0,00697 0,0000473 |
— |
Прежде всего оценим ло нашим выборкам неизвестные зна чения р . Положим, что сумма частот у2 + у, + V, равна сумме
120