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

книги из ГПНТБ / Рыжов, П. А. Математическая статистика в горном деле учеб. пособие

.pdf
Скачиваний:
25
Добавлен:
19.10.2023
Размер:
7.54 Mб
Скачать

Примяв

пятипроцентный уровень

значимости,

по

табл. 32 находят значение критической

области t„. При

k = 488 —2 = 486, ^=1,96.

 

 

Среднее содержание цинка по блоку оказывается

в

критической

области (4,2 значительно

больше 1,96), по­

этому гипотеза об однородности средних по блокам от­

вергается. Следовательно, средние по блокам

имеют су­

щественное отличие от общей средней.

 

 

 

 

Проверка гипотезы о равенстве дисперсий. F-pacnpe-

деление. При проверке гипотезы о равенстве

дисперсий

используется ^-распределение.

 

 

 

 

Пусть имеются две независимых выборки

объемами

пх и Пц из двух нормально распределенных

совокупно­

стей с дисперсиями а \ х ) и

а2о(и).

а2х

 

а2ѵ

 

Отношение выборочных

дисперсий

и

имеет

/•"-распределение со степенями свободы

к — пх1

и ко —

= / і „ — 1 .

 

 

 

 

 

В табл. 34 приведены пятипроцентные пределы откло­

нения величины F в зависимости от степеней

свободы ki

и ко. При проверке гипотезы о равенстве

дисперсии зна­

чение, вычисленное по выборкам, сравнивают с границей допустимых значений.

Табл. 34 составлена для отношения большей выбороч­ ной дисперсии к меньшей, поэтому F будет всегда боль­ ше единицы.

Если определенное при заданном уровне значимости F больше вычисленного по данным опыта (по данным вы­ борок), то считается, что опытные данные не противоре­ чат гипотезе о равенстве дисперсии. В противном случае различие в дисперсиях признается существенным.

Пример 4.5. На руднике производилось эксперимен­ тальное сравнение данных бороздового и шпурового ме­ тодов опробования.

При проходке горизонтальной выработки вкрест про­ стирания мощной залежи отбирались пары проб (бороз­ довая и шпуровая), которые брали одну рядом с другой.

Расстояние

между

парами проб 2 м. Всего было взято

26 пар проб. Длина

борозды и шпура составляла 1 м. Ре­

зультаты

химических анализов проб приведены в

табл. 35.

 

 

Проверку гипотезы о равенстве дисперсий выполняют

вследующей последовательности.

1.Вычисляют несмещенные значения выборочных дис­ персий для каждого метода опробования.

120 —

 

 

 

 

Т а б л и ц a â5

 

Бороздовое

 

Шиуроное

 

>Ь опыта

опробома-

( г , - С , ) 2

онробопа-

( г , . - С , ) 2

 

II не

 

ннс f.,

 

 

 

 

% "

 

1

1,65

+0,92

0,846

2,50

2

1,20

+0,47

0,221

1,10

3

0,24

—0,49

0,240

0,31

4

0,27

—0,46

0,212

0,64

5

1,15

+0,42

0,176

1,30

6

0,30

—0,43

0,185

0,25

7

0,15

—0,58

0,336

0,30

8

1,30

+0,57

0,325

1,20

9

0,45

—0,28

0,078

0,25

10

1,00

+0,27

0,075

0,75

11

0,75

+0,02

0,000

0,75

12

0,80

+0,07

0,005

1,10

13

0,20

—0,53

0,281

0,10

14

1,05

+0,32

0,102

1,20

15

0,65

—0,08

0,006

1,10

16

0,17

—0,56

0,314

0,20

17

1,20

+0,47

0,221

1,15

18

0,50

—0,23

0,053

0,40

19

0,95

+0,22

0,048

0,80

20

0,75

+0,02

0,000

0,45

21

0,14

—0,59

0,348

0,65

22

0,73

—0,00

0,000

1,00

23

0,65

—0,08

0,006

0,35

24

1,50

+0,77

0,593

1,10

25

1,20

+0,47

0,221

1,50

26

Следы

—0,73

0,533

0,10

 

18,96

—0,03

5,423

20,55

 

Сі=0,73;

 

 

 

+

1,71

2,924

+

1,31

0,096

—0,48

0,230

—0,15

0,022

+0,51

0,260

—0,54

0,292

—0,49

0,240

+0,41

0,168

—0,54

0,292

—0,04

0,002

—0,04

0,002

+0,31

0,096

—0,69

0,476

+0,41

0,168

+0,31

0,096

—0,59

0,348

+0,36

0,130

—0,39

0,152

+0,01

0,000

—0,34

0,116

—0,14

0,020

+0,21

0,044

—0,44

0,194

+0,31

0,096

+0,71

0,504

—0,69

0,476

+0,01

7,444

С 2 = 0 , 7 9

По данным бороздового опробования

2

id-С)*

5,423

 

п—

1

= 0,227.

25

 

По данным шпурового опробования

г 7,444 Ог= — = 0,297.

121 —

2. Вычисляют опытное значение F:

 

 

 

 

 

о\

0,217

 

 

3.

При &і=/г2 = 26—1=25 и уровне

значимости

а =

= 5%;

по табл. 34 определяют границу критической

обла­

сти: F = 1,96.

 

 

 

 

F =

Полученное

по

данным

опробования значение

= 1,37

попадает

в

область

допустимых

значений крите­

рия F, поэтому

данные опыта не противоречат гипотезе

о равенстве дисперсий.

 

 

 

Следовательно,

данные

бороздового

опробования в

отношении случайной ошибки не являются более точными по сравнению с данными шпурового опробования.

Отметим, что F не зависит от центров сравниваемых распределений, поэтому проверка гипотезой о равенстве дисперсий не устанавливает наличия систематических ошибок метода.

Г л а в а 5 ДИСПЕРСИОННЫЙ А Н А Л И З

§ 1. СУЩНОСТЬ МЕТОДА

Данные, полученные в результате исследования того или иного явления, часто зависят от ряда независимых, изменяющихся в течение опыта -параметров (факторов). Количественная оценка влияния отдельных факторов на изучаемое явление при некоторых условиях может быть произведена с помощью дисперсионного анализа.

Рассмотрим сущность этого метода на примере ана­ лиза данных разведки никелевого месторождения.

Пример 5.1. Разведка крупного месторождения нике­ левых руд производится с помощью скважин, располо­ женных по определенной сетке. При опробовании руды и вмещающих пород керны отбирались через каждый метр.

Контроль разведки производился с помощью шур­ фов, закладываемых по центрам ранее пройденных сква­ жин. Для опробования руды в шурфах отбирались вер­ тикальные бороздовые пробы метровой длины. Содержа­ ние никеля в пробах, взятых по шурфам, в большинстве случаев незначительно превышало содержание никеля в пробах, взятых по кернам в тех же интервалах.

Проверка бороздового ' опробования производилась валовым опробованием. Вся руда, извлеченная при углуб­ лении шурфа на один метр, отбиралась как одна валовая проба. Сравнение валового опробования с бороздовым показало отсутствие систематических ошибок при бороз­ довом опробовании шурфов. Поэтому контроль разведки

— 123 —

месторождения с помощью шурфов в данном случае до­ статочно обоснован.

Задача исследования состоит в установлении причин, обусловливающих ошибки при разведке скважинами.

В зависимости от условий проведения скважин, отбо­ ра керна и свойств полезного ископаемого причинами расхождений в значениях содержания металла в пробах, отобранных по кернам, и в пробах, отобранных в шур­ фах, могут быть:

1)перемещение с поверхности керна минералов, со­ держащих никель. В этом случае ошибки содержания никеля по скважинам зависят от глубины расположения проб;

2)разрушение хрупких рудных минералов и вдавли­ вание их в стенки скважины. В этом случае ошибки зави­ сят от величины содержания никеля в пробах.

Исходный материал для дисперсионного анализа за­ писывают в формуляр (табл. 36).

 

 

 

 

Т а б л и ц а 36

 

 

 

Содержание

никеля

Интервалы

опро­

 

 

 

бования,

м

0 - 6

6 - 2 *

(r-\)b-:b

 

 

0—1

/—2/

(k—l) l—kl

Для каждой пары сравниваемых содержаний опре­ деляют расхождение:

А с =

сС І І В

Сщур.

Полученные значения

Ас

заносят в клетки табл. 36

в соответствии с глубиной расположения проб и содер­

жанием никеля по шурфу. Табл. 36 позволяет

проследить

влияние каждого из отмеченных

факторов

на

ошибки

определения

содержания по

скважинам.

Рассмотрим,

например, изменение ошибки с увеличением

глубины.

В первой

строке приведены

Ас

для глубины

0 I м

при всех возможных значениях содержания никеля. Вто­ рая строка содержит ошибки для интервала опробования /—21 м также при всех значениях содержания никеля и т. д.

— 124 —

Таким образом, различие в средних значениях, вычис­ ленных по строкам табл. 36, зависит от глубины взятия пробы. Аналогично этому изменение средних значений содержания металла по графам обусловлено только значением содержания никеля в пробах.

Предположим, что ни один из факторов не влияет на возникновение ошибки. Тогда средние значения, вычис­ ленные по строкам и графам табл. 36, будут иметь толь­ ко случайные отклонения. Сущность дисперсионного ана­ лиза заключается в проверке гипотезы о равенстве средних. Дисперсионный анализ разработан для нор­ мального распределения изучаемой величины X. Практи­ чески он может быть применен и при отклонениях этих распределений от нормального закона,

§ 2, ОДНОФАКТОРНЫЙ ДИСПЕРСИОННЫЙ АНАЛИ З

Рассмотрим наиболее простой случай дисперсионного

анализа при проверке

гипотезы о влиянии

фактора В на

результаты

испытаний.

 

 

 

 

 

Пусть фактор В может встретиться в г вариантах. Для

каждого из возможных вариантов В\, В2,

ß r

получены

выборки объемами;

п2,

пг.

 

в

 

формуляр

Результаты

опыта

записывают

 

(табл. 37).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 37

 

 

Наблюдения

л

 

 

 

Средние зна ­

Варианты в

 

,Чнсло

наблюдений

чения но

 

 

 

 

 

 

 

 

вариантам

 

х \\< х \% •

• • i Х\п

 

 

х \

в 2

Х2\,

Х22

 

А'2м2

 

п 2

 

х 2

Вг

х г \ ,

хг2

 

хгпг

 

п г

 

х г

Общее число наблюдений

N = 2 пи

— 125 —

Среднее значение для всех наблюдений

Предположим, что фактор В не оказывает

влияния

на изучаемую величину X. Тогда различие в выборочных

средних по вариантам будет носить случайный

характер.

Пусть Хй

= а. Проверяемая

 

гипотеза заключается

в вы­

полнении

равенства

 

 

 

 

 

 

 

 

 

ai = а2

=

... =

а,

 

 

 

 

 

где а —- средние значения

по вариантам фактора

В.

 

Для построения доверительной

области

необходимо

знать дисперсии средних значений ХІ П О вариантам.

 

Возьмем один из вариантов. Для него

можно

соста­

вить равенство

 

 

 

 

 

 

 

 

 

2 ( * < - а ) * = 2 ( x t - x + x - a ) *

 

=

 

 

 

 

= 2 (хо-х)г+

 

 

 

 

(5-1)

где п — число наблюдений по варианту.

 

 

 

 

 

Левая

часть выражения

(5.1)

имеет

распределение

а2, X 2 с п

степенями свободы. Первое из слагаемых

пра­

вой части подчиняется распределению б2Х

с

(п—\)

сте­

пенями свободы,так как

 

 

 

 

 

 

 

 

n

Последнее слагаемое имеет одну степень свободы (как квадрат одной нормально распределенной вели­ чины) .

Таким образом, сумма п независимых квадратов раз­ бивается на две части, одна из которых имеет (п1) степеней свободы, а другая — одну степень свободы.

Несмещенная оценка дисперсии по варианту может быть получена по формуле:

— 126 —

или

üb = п (х — а)2.

В первом случае а2 не зависит от а и, следовательно, среднее по .варианту не будет равно а, во втором —оцен­ ка db2 справедлива только в том случае, если среднее по

варианту

равно а.

 

 

 

 

Если

а в .варианте отличается от

среднего

значения

во всей совокупности, то отношение дисперсий

F =

6b2lo2,

подчиненное

/•'-распределению,

будет

больше

единицы.

Величина

F

служит критерием

проверки гипотезы

при

установлении влияния фактора В на изучаемое явление.

Распределение ХІ

П О

вариантам предполагается нор­

мальным с одинаковыми

дисперсиями

 

2

2

г

01

=

02 = • • • =

0*г.

Примем гипотезу

Ü2 = ... = аг = а.

Согласно выражению (5.1), по каждому варианту фактора В можно написать

п

п

2

(ХІ - а)2 = 2 (Хі - х)2 + ru (х - а)2.

І1

Суммируем выражения по всем вариантам

г

и

г

п

г

2 2 ( * * - t f ) 2 = 2 2

( ^ - * ) 2 + 2 M * - ß ) 2 - (5-2)

i

l

i

l

i

 

r

 

n

 

 

2 S !*< - *>

 

Величина tri? =

 

-

представляет собой не-

 

 

 

N

r

смещенную оценку дисперсий, обусловленную

колебания­

ми X внутри вариантов. Очевидно, оу2 не зависит от сред­

них значений по вариантам.

 

 

Второе значение дисперсии

 

 

г

 

 

2 "»•(*»•-я)2

 

о-6 = —

-

(5.3)

— 127 —

отображает колебания средних в отдельных вариантах относительно общей средней X.

Отношение дисперсий а2ь/ог2 имеет /-"-распределение с (/•—1) и (N—/') степенями свободы.

Задавшись уровнем значимости при известных степе­ нях свободы, по табл. 34 определяем границу допустимых значений F. Если вычисленное по данным опыта F пре­ вышает установленный предел, то нулевая гипотеза от­

вергается и считается,

что на изучаемое

явление

влияет

фактор В.

 

 

 

 

 

 

 

 

 

Пример

5.2. В табл. 38 приведены

скорости бурения,

замеренные

на различных

участках

шпура

диаметра

34 мм.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

38

 

 

 

Скорость

бурении,

смімин

 

 

 

Глубина

шпура, .и

 

 

 

 

 

 

 

 

 

 

опыт № 1

опыт № 2

опыт N> 3

опыт

опыт а

5

0 - 0 , 2

13,0

15,0

13,7

 

17,6

 

14,8

 

0 , 2 - 0 , 4

15,2

9

, 3

11,1

 

9 , 5

 

12,4

 

0 , 4 — 0 , 6

12,1

1 2

, 9

13, 7

 

12,0

 

1 5 , 2

 

0 , 6

— 0 , 8

9 , 8

8

, 9

13, 4

 

14,7

 

11,8

 

0 , 8

— 1 , 0

11, 3

10,8

14,1

 

1 2 , 2

 

10, 0

 

Вмикрокварцитах опытное бурение производилось молотком ОМ-506 при постоянном давлении воздуха. Требуется установить, влияет ли глубина шпура (в пре­ делах 1 иг) на скорость бурения? Иными словами, нужно проверить гипотезу о равенстве скоростей бурения на различных участках шпура. Для решения этой задачи применяем методику однофакторного дисперсионного анализа.

Втабл. 39 и 40 приведены числовые значения сумм и квадратов сумм скоростей, полученные по наблюдаемым значениям скорости бурения для каждого шпура.

По формуле

г п

2

2

2

< * - * ' > •

i

l

/Г-

о,- =

•найдем оценку дисперсии,

 

обусловленную колебаниями

х внутри вариантов.

 

 

 

128 —

 

 

 

 

 

 

Т а б л и ц а 39

 

 

 

Скорость бурения,

смінин

 

Глубина

 

 

 

 

 

 

6

шнура,

 

-г*

 

-Ѵ.|

-Г 5

 

м

 

 

 

?*'

 

 

 

 

 

 

 

0,1

13,0

15,0

13,7

17,6

14,8

74,1

1110,5

0,3

15,2

9,3 •

11,1

9,5

12,4

57,5

684,7

0,5

12,1

12,9

13,7

12,0

15,2

65,9

875,5

0,7

9,8

8,9

13,4

14,7

11,8

58,6

710,1

0,9

11,3

10,8

14,1

12,2

10,0

58,4

691,9

 

61,4

56,9

66,0

66,0

61,2

314,5

4072,7

 

 

 

 

 

Т а б л и ц а 40

Глубина

 

п

H

л

л

я

 

 

шпура,

п і

f i

S

(ЪхГ-

 

Цх —.г)2

м

 

 

1

1

я

 

 

 

 

 

 

 

1

2

3

4

5

6

7 = ( 4 ) -

 

 

 

 

 

 

- ( 6 )

0,1

5

74,1

1110,5

5490,8

1098,2

12,3

0,3

5

57,5

684,7

3306,2

661,2

23,5

0,5

5

65,9

875,5

4342,8

868,6

6,9

0,7

5

58,6

710,1

3434,0

686,8

23,3

0,9

5

58,4

691,9

3410,6

682,1

9,8

 

25

314,5

4072,7

19984,4

3996,9

75,8

Значение 22 найдем в графе 7 табл. 40; общее число

наблюдений п—25; число вариантов фактора п = 5 , сле­ довательно,

г75,8

Дисперсию вь2, выражающую колебания средних в вариантах относительно общей средней, определяем по

5-1924

— 129 —

Соседние файлы в папке книги из ГПНТБ