
книги из ГПНТБ / Рыжов, П. А. Математическая статистика в горном деле учеб. пособие
.pdfПримяв |
пятипроцентный уровень |
значимости, |
по |
табл. 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 Х\п |
|
1Ц |
|
х \ |
||
в 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 — |