Курсовая работа (т): Анализ и построение зависимостей

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

7. Проверка, что выборка осуществлялась по случайному закону

Будем использовать критерий согласия Пирсона

. (1)

Если, как это обычно имеет место,

,

то указанную формулу можно преобразовать к виду более удобному для вычислений

. (2)

Проверяемая гипотеза состоит в том, что выбор 150 изделий из нормальной генеральной совокупности с μ = 60 и σ = 10 был произведён по случайному закону. Так как известны параметры нормальной генеральной совокупности, то критерий χ2 может только проверить случайный характер выборки.

В качестве уровня значимости выберем α = 0,05.

Для того, чтобы получить визуальное представление о степени соответствия нашей выборки нормальной кривой, воспользуемся гистограммой частот 150 изделий (см. рис. Б) с наложенной на неё нормальной кривой с параметрами μ = 60 и σ = 10. Вычисленные значения ординат нормальной кривой приведены в таблице 7.

Таблица 7

Ординаты нормальной кривой: μ=60; σ=10; n=150; ∆Х=8

Х

Z (x/σ)

T (z2/2)

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

± zs

Теоретический процент наблюдений

Фактический процент наблюдений

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

Источник: https://www.bibliofond.ru/detail.aspx?id=877693