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