Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Прикладная статистика

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
9.53 Mб
Скачать

зрения принятия выдвинутой гипотезы о законе распределения генеральной совокупности 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

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