2013
|
2,203
|
0,070
|
0,529
|
46,544
4.1 Исключение
ложных данных
Ложные результаты проникают в исходную
совокупность случайным путем из-за действия достаточно большого числа случайных
причин, которые не удается заранее предусмотреть при проведении измерений или
при сборе статистических сведений (описки, опечатки). Не исключен случай, когда
ошибки в статистические данные вносятся злонамеренно, к примеру, с целью
искажения или сокрытия истинного положения вещей.
Процедура исключения ложных данных позволяет
выявить такие данные в статистической совокупности с большой степенью
вероятности.
К источникам ложных данных относятся:
грубые (недостаточно точные) измерения;
нарушение условий эксперимента;
использование неисправного оборудования;
ошибки при обработке информации и др.
Наличие ложных данных приводит к необоснованным
выводам, поэтому ложные данные должны быть исключены из выборочной совокупности
в процессе первичной обработки опытных данных [26,30,33,35].
Процедура исключения ложных данных для
представленных совокупностей:
. Имеющиеся значения совокупностей
значений располагаются в порядке возрастания: х(1) ≤ х(2)
≤ х(m),
и определяются крайние элементы представленной совокупности:
х(1) = хmin,
х(n)
= х(max).
Таблица 7 - Упорядоченные по возрастанию
совокупности значений
L,
чел.
|
m1,
чел.
|
m2,
чел.
|
m3,
чел.
|
N,
дни
|
|
278200
|
625
|
10
|
18
|
29090
|
|
283730
|
725
|
12
|
18
|
29482
|
|
290000
|
858
|
13
|
32
|
35398
|
|
291008
|
899
|
17
|
36
|
35476
|
|
292200
|
912
|
18
|
36
|
35604
|
|
324500
|
1133
|
20
|
38
|
40278
|
|
339100
|
1313
|
28
|
40
|
40689
|
|
343600
|
1542
|
30
|
47
|
42046
|
|
351900
|
1650
|
31
|
51
|
47039
|
|
357900
|
1722
|
33
|
58
|
48838
|
|
390500
|
2059
|
44
|
63
|
52710
|
|
411700
|
2387
|
46
|
63
|
56094
|
|
427300
|
2540
|
51
|
65
|
58420
|
|
2,203
|
0,034
|
0,529
|
23,000
|
1527,6
|
2,500
|
0,041
|
0,621
|
23,500
|
1566,1
|
3,084
|
0,047
|
1,062
|
24,417
|
1778,3
|
3,089
|
0,052
|
1,150
|
24,660
|
2221,3
|
3,121
|
0,062
|
1,171
|
25,600
|
2519,4
|
3,492
|
0,070
|
1,237
|
26,121
|
3155,6
|
3,872
|
0,081
|
1,313
|
35,826
|
4194,5
|
4,308
|
0,085
|
1,369
|
38,899
|
4992,5
|
4,689
|
0,087
|
1,449
|
39,604
|
5488,4
|
5,012
|
0,097
|
1,474
|
40,665
|
7007,6
|
5,273
|
0,103
|
1,579
|
41,256
|
7306,8
|
5,586
|
0,118
|
1,613
|
43,105
|
8339,3
|
6,170
|
0,124
|
1,688
|
46,544
|
11161,5
|
17,8
|
10360
|
11,5
|
65,6
|
41400
|
18
|
10886
|
13,1
|
78,3
|
49800
|
18,2
|
11690
|
13,6
|
89
|
56900
|
18,5
|
11729
|
15,5
|
100,8
|
64800
|
20,5
|
11811
|
26,9
|
140
|
90400
|
20,8
|
11929
|
34,3
|
164,8
|
107000
|
23,3
|
12124
|
41,3
|
205,8
|
134100
|
25,8
|
12314
|
42,3
|
226,1
|
148200
|
26,5
|
12500
|
44,5
|
240,3
|
157200
|
26,5
|
12550
|
49,9
|
264,5
|
173800
|
30,9
|
12975
|
57,2
|
318,8
|
187300
|
33,3
|
13690
|
60,9
|
335,4
|
220800
|
33,5
|
14383
|
71,8
|
404,8
|
244700
|
2. Вычисляются значения ,
отклонений
крайних значений случайной величины хmin,
хmax
от ее среднего значения с учетом разброса значений
случайной величины x в выборке
по следующим формулам:
, . (5)
Среднее значение опытных данных вычисляется
по исходной совокупности объемом m
по формуле:
(6)
Среднее значение для величины L:
Lср==337049,1
По исходной совокупности опытных данных
вычисляется величина , характеризующая
разброс опытных значений величины x
вокруг среднего
Разброс опытных данных для величины L:
. (7)
1,18
1,80
. Рассчитанные значения относительно
отклонения (α;m)
сравниваются с его критическими (теоретическими) значениями .
Так как выборка небольшая (m≤25)
для исключения ложных данных используется таблица квантилей распределения
максимального относительного отклонения .
Значение квантиля максимального относительного отклонения для
уровня значимости α=0,05 принимается
равным 2,43. [30,35]
Таблица 8 -Квантили максимального относительного
отклонения
α,
%
|
Объем
выборочной совокупности m
|
|
8
|
9
|
10
|
11
|
12
|
13
|
14
|
15
|
16
|
1
|
2,37
|
2,46
|
2,54
|
2,61
|
2,66
|
2,71
|
2,76
|
2,8
|
2,84
|
5
|
2,17
|
2,24
|
2,29
|
2,34
|
2,39
|
2,43
|
2,46
|
2,49
|
2,52
|
Сравниваются расчетное и критическое значения по
формуле:
≤.
Оба значения и
1,80
меньше =2,43,
следовательно, оба неравенства выполнены и оба значения не считаются ложными и
остаются в выборке.
Таким образом, рассчитываются оставшиеся
выборочные совокупности. Результаты расчетов представлены в таблице
Таблица 9 - Результаты расчетов ,
отклонений
крайних значений случайных величин хmin,
хmax
от среднего значения
|
хmin
|
хmax
|
|
|
|
|
L
|
278200
|
427300
|
337049,1
|
50071,51
|
1,18
|
1,80
|
m1
|
625,00
|
2540,00
|
1412,69
|
633,14
|
1,24
|
1,78
|
m2
|
10,00
|
51,00
|
27,15
|
13,63
|
1,26
|
1,75
|
m3
|
15,00
|
65,00
|
43,23
|
16,46
|
1,72
|
1,32
|
N
|
29090,00
|
58420,00
|
42397,23
|
9609,77
|
1,38
|
1,67
|
Kч
|
2,203
|
6,170
|
4,031
|
1,197
|
1,787
|
0,8
|
Kл
|
0,034
|
0,124
|
0,077
|
0,028
|
1,685
|
0,02
|
Kп.з.
|
0,529
|
1,688
|
1,250
|
0,341
|
1,285
|
0,25
|
Kт
|
23,000
|
46,544
|
33,323
|
8,472
|
1,561
|
5,0
|
S
|
1527,60
|
11161,50
|
4712,3
|
3024,73
|
1,05
|
2,13
|
D
|
17,80
|
33,50
|
24,12
|
5,78
|
1,09
|
1,62
|
E
|
10360,00
|
14383,00
|
12226,23
|
1,75
|
2,02
|
I
|
11,50
|
71,80
|
37,14
|
20,00
|
1,28
|
1,73
|
V
|
65,50
|
404,80
|
202,70
|
108,25
|
1,27
|
1,87
|
Vp
|
41400,00
|
244700,00
|
128953,8
|
66864,25
|
1,31
|
1,73
|
Как видно из таблицы все значения и
выборочных
совокупностей удовлетворяют неравенству ≤.
В результате проведенной процедуры можно утверждать, что в таблице значений
случайных величин ложных данных нет.
4.2 Проверка опытных
данных на их случайность и независимость
При проведении статистического исследования
возможны случаи нарушения условий проведения эксперимента, не относящиеся к
наличию в выборочной совокупности ложных данных. Случайность и независимость
опытных данных - необходимое условие репрезентативности выборочной
совокупности.
Наблюдение считается статистически независимым,
если результаты, полученные в результате отдельного наблюдения, не связаны с
данными предыдущих и последующих наблюдений. Необходимы критерии, которые
позволяют установить случайность и независимость данных в выборочной
совокупности.
Для статистической проверки случайности и
независимости результатов наблюдения применяются:
) критерий серий, основанный на использовании
медианы выборки;
) критерий «восходящих» и «нисходящих»
серий.
В данном случае выбрана проверка опытных данных
с помощью критерия, основанного на использовании медианы выборки при уровне
значимости α=0,05.
Критерий, основанный на использовании медианы
выборки, позволяет заметить монотонное смещение среднего выборочного значения в
ходе эксперимента. [26]
1. Имеются 13 совокупностей значений
случайных величин, прошедшие процедуру исключения ложных данных. В данном
случае n=m,
так как ни одного значения исключено не было (m
- количество значений в выборке, n
- число наблюдений, оставшихся после исключения ложных результатов).
Совокупности располагаются в порядке возрастания.
(Таблица
7)
. Находится выборочное значение медианы (n)
по следующей формуле:
Так как число n=13
является нечетным, то согласно формуле получаем:
Так для каждого значения выписывается седьмое
значение из упорядоченной выборочной совокупности и заносится в таблицу:
Таблица 10 - Значения медианы
Показатели
|
|
L
|
339100
|
m1
|
1313
|
m2
|
28
|
m3
|
40
|
N
|
40689
|
Kч
|
3,872
|
Kл
|
0,081
|
Kп.з.
|
1,313
|
Kт
|
35,826
|
S
|
4194,5
|
d
|
23,3
|
e
|
12314
|
I
|
44,5
|
V
|
205,8
|
Vp
|
134100
|
3. В исходной (неупорядоченной) выборке x1,
x2,
…,xn вместо каждого
числа xi
ставится
«+» (плюс), если xi
> xmed(n),
и «-» (минус), если
xi<
xmed(n).
Значениям xi=
xmed(n)
никакого знака не присваивается.
Таблица 11 - Последовательность
№
|
L
|
m1
|
m2
|
m3
|
N
|
Kч
|
Kл
|
Kп.з.
|
1
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
2
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
3
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
4
|
+
|
+
|
|
+
|
+
|
+
|
|
+
|
5
|
+
|
+
|
+
|
+
|
|
+
|
+
|
+
|
6
|
+
|
+
|
+
|
+
|
-
|
+
|
+
|
|
7
|
|
|
+
|
-
|
+
|
|
+
|
-
|
8
|
-
|
-
|
-
|
-
|
+
|
-
|
-
|
-
|
9
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
10
|
-
|
-
|
-
|
|
-
|
-
|
-
|
+
|
11
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
12
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
13
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
№
|
Kт
|
S
|
D
|
E
|
I
|
V
|
Vp
|
1
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
2
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
3
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
4
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
5
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
6
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
7
|
|
|
|
|
|
|
|
8
|
+
|
+
|
+
|
+
|
-
|
+
|
+
|
9
|
+
|
+
|
+
|
-
|
-
|
+
|
+
|
10
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
11
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
12
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
13
|
+
|
+
|
+
|
+
|
+
|
+
|
+
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Полученные последовательности плюсов и минусов
характеризуются числом серий ν(α;n)
и длиной (α;n)
самой длинной серии. Под серией понимается последовательность идущих подряд
плюсов или минусов.
Производится подсчет числа νрасч(n)
подряд идущих знаков «+» и подряд идущих знаков «-», а также длина расч(α;n)
самой длиной серии плюсов или минусов. Пустые клетки таблицы при подсчетах не
учитываются.
Таблица 12 - Значения νрасч
и расч
Показатель
|
νрасч(n)
|
расч(n)
|
L
|
2
|
6
|
m1
|
2
|
6
|
m2
|
2
|
6
|
m3
|
2
|
6
|
N
|
4
|
5
|
Kч
|
2
|
6
|
Kл
|
2
|
6
|
Kп.з.
|
4
|
5
|
Kт
|
2
|
6
|
S
|
2
|
6
|
D
|
2
|
6
|
E
|
4
|
6
|
I
|
4
|
6
|
V
|
2
|
6
|
Vp
|
2
|
6
|
4. Рассматривается гипотеза о случайности и
независимости данных в рассматриваемых выборочных совокупностях (при уровне значимости
α=0,05).
Если хотя бы одно из условий системы окажется
невыполненным, то предположение о независимости результатов наблюдения
отвергается с вероятностью α=0,05 совершить
ошибку первого рода.
Если ≤,
то есть если расчетное число серий не будет
превосходить критическую величину числа
серий, то данные исследуемой выборочной совокупности следует признать
неслучайными и зависимыми при заданном уровне значимости α=0,05.
Это
же верно и в случае, если > ,
то есть если расчетная длина самой длинной из
серий превосходит критическую величину ,
вычисленную по формуле, либо равна ей.
Сравнив полученные результаты с
расчетными для всех выборок, напрашивается вывод, что с вероятностью 1-α=1-0,05=0,95
гипотеза
о случайности и независимости совокупности исследуемых выборочных значений
случайных величин не должна быть отвергнута.
5. Статистический
анализ выборочных совокупностей
После процедур первичной обработки исходная
совокупность опытных результатов представляется выборочной совокупностью, для которой
выполнены важнейшие требования, обеспечивающие ее репрезентативность.
5.1 Проверка
гипотезы о нормальном распределении генеральной совокупности
Использование регрессионного анализа с целью
прогноза опирается на гипотезу о нормальности распределения исследуемой
совокупности опытных данных. Подбор вида любого из предполагаемых законов
распределения рекомендуется начинать с построения гистограммы опытных значений.
Для построения гистограммы сначала нужно найти
наименьшее xmin
и
наибольшее xmax
из выборочных значений, а исходные данные xi
(i=1,2,…,n),то
есть промежуток [xmin,
xmax], разбить
на ряд частичных интервалов. При разбиении исходного промежутка [xmin,
xmax] на ряд
частичных интервалов применяется следующая формула для вычисления искомой длины
частичных интервалов:
(8)
Число l
таких частичных интервалов в дальнейшем может быть уменьшено, если в отдельные
интервалы не попадает ни одного значения (или попадает слишком мало значений).
[26,33]
Таблица 13 - Интервалы значений
|
хmin
|
хmax
|
|
|
L
|
278200
|
427300
|
31886
|
|
m1
|
625,00
|
2540,00
|
410
|
|
m2
|
10,00
|
51,00
|
8,8
|
|
m3
|
15,00
|
65,00
|
10,7
|
|
N
|
29090,00
|
58420,00
|
62,72
|
Kч
|
2,203
|
6,170
|
0,8
|
Kл
|
0,034
|
0,124
|
0,02
|
Kп.з.
|
0,529
|
1,688
|
0,25
|
Kт
|
23,000
|
46,544
|
5,0
|
S
|
1527,60
|
11161,50
|
2060
|
d
|
17,80
|
33,50
|
3,4
|
e
|
10360,00
|
14383,00
|
860
|
I
|
11,50
|
71,80
|
12,9
|
V
|
65,50
|
404,80
|
72,6
|
Vp
|
41400,00
|
244700,00
|
43477
|
|
|
|
|
|
|
|
Гистограммы:
Рис.2
Рис.3
Рис.4
Рис.5
Рис.6 Рис.7
Рис.8
Рис.9
Рис.10 Рис.11
Рис.12
Рис.13
Рис.14
Рис.15
Рис. 16
Так как объем рассмотренных выборочных
совокупностей невелик, можно было не выполнять проверку на нормальность их
распределения.
.2 Проверка
гипотезы о значимости коэффициента корреляции
Проверку значимости коэффициента корреляции двух
совокупностей опытных данных можно выполнить двумя способами: с помощью
критерия Стьюдента и с помощью критерия Фишера-Снедекора.
Формула для нахождения коэффициента парной
корреляции:
(9)
Таблица 14 - Вычисленные значения коэффициентов
корреляции
|
|
|
L
|
S
|
D
|
e
|
I
|
V
|
Vp
|
m1
|
0,976
|
-0,848
|
-0,890
|
-0,599
|
-0,920
|
-0,935
|
-0,948
|
m2
|
0,940
|
-0,835
|
-0,800
|
-0,524
|
-0,806
|
-0,836
|
-0,857
|
m3
|
0,874
|
-0,780
|
-0,936
|
-0,721
|
-0,968
|
-0,972
|
-0,979
|
N
|
0,899
|
-0,740
|
-0,775
|
-0,545
|
-0,772
|
-0,833
|
-0,838
|
Kч
|
0,941
|
-0,861
|
-0,930
|
-0,650
|
-0,943
|
-0,965
|
-0,977
|
Kл
|
0,873
|
-0,824
|
-0,772
|
-0,505
|
-0,745
|
-0,787
|
-0,812
|
Kп.з.
|
0,663
|
-0,631
|
-0,876
|
-0,776
|
-0,890
|
-0,893
|
-0,897
|
Kт
|
-0,865
|
0,838
|
0,923
|
0,613
|
0,933
|
0,931
|
0,937
|
1. Проверка значимости коэффициента
корреляции с помощью критерия Фишера-Снедекора.
Формула для расчета опытного значения F-критерия
Фишера-Снедекора для проверки значимости коэффициента rxy(x,y)
парной линейной корреляции при заданном объеме n
выборочной совокупности:
(10)
Таблица 15 - Вычисленные опытные значения
критерия Фишера
|
L
|
S
|
D
|
|
|
|
|
|
|
|
|
|
|
m1
|
0,95
|
0,05
|
220,81
|
0,72
|
0,28
|
28,26
|
0,79
|
0,21
|
41,89
|
m2
|
0,88
|
0,12
|
83,94
|
0,70
|
0,30
|
25,25
|
0,64
|
0,36
|
19,49
|
m3
|
0,76
|
0,24
|
35,44
|
0,61
|
0,39
|
17,11
|
0,88
|
0,12
|
77,91
|
N
|
0,81
|
0,19
|
46,51
|
0,55
|
0,45
|
13,35
|
0,60
|
0,40
|
16,52
|
Kч
|
0,89
|
0,11
|
84,86
|
0,74
|
0,26
|
31,41
|
0,87
|
0,13
|
70,98
|
|
L
|
S
|
D
|
|
|
|
|
|
|
|
|
|
|
Kл
|
0,76
|
0,24
|
35,35
|
0,68
|
0,32
|
23,21
|
0,60
|
0,40
|
16,26
|
Kп.з.
|
0,44
|
0,56
|
8,61
|
0,40
|
0,60
|
7,29
|
0,77
|
0,23
|
36,14
|
Kт
|
0,75
|
0,25
|
32,59
|
0,70
|
0,30
|
25,85
|
0,85
|
0,15
|
63,63
|
|
E
|
I
|
|
|
|
|
|
|
|
m1
|
0,36
|
0,64
|
6,15
|
0,85
|
0,15
|
60,90
|
m2
|
0,27
|
0,73
|
4,15
|
0,65
|
0,35
|
20,41
|
m3
|
0,52
|
0,48
|
11,94
|
0,94
|
0,06
|
162,59
|
N
|
0,30
|
0,70
|
4,64
|
0,60
|
0,40
|
16,18
|
Kч
|
0,42
|
0,58
|
8,03
|
0,89
|
0,11
|
89,17
|
Kл
|
0,26
|
0,74
|
3,77
|
0,56
|
0,44
|
13,75
|
Kп.з
|
0,60
|
0,40
|
16,63
|
0,79
|
0,21
|
42,06
|
Kт
|
0,38
|
0,62
|
6,62
|
0,87
|
0,13
|
73,99
|
|
V
|
Vp
|
|
|
|
|
|
|
m1
|
0,87
|
0,13
|
76,70
|
0,90
|
0,10
|
97,91
|
m2
|
0,70
|
0,30
|
25,53
|
0,73
|
0,27
|
30,44
|
m3
|
0,94
|
0,06
|
184,87
|
0,96
|
0,04
|
259,79
|
N
|
0,69
|
0,31
|
24,95
|
0,70
|
0,30
|
26,05
|
Kч
|
0,93
|
0,07
|
147,82
|
0,95
|
0,05
|
232,57
|
|
V
|
Vp
|
|
|
|
|
|
|
|
Kл
|
0,62
|
0,38
|
17,88
|
0,66
|
0,34
|
21,33
|
Kп.з
|
0,80
|
0,20
|
43,20
|
0,80
|
0,20
|
45,22
|
Kт
|
0,87
|
0,13
|
71,84
|
0,88
|
0,12
|
79,02
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вычисленные опытные значения сравниваются
с критическим (табличным) значением =F(α;k1;k2)
критерия Фишера-Снедекора. Критические значения критерия Фишера-Снедекора
содержатся в статистических таблицах и зависят:
от выбранного уровня значимости α;
- от показателей k1
и k2 степеней
свободы критерия F.
В рассматриваемом случае парной линейной
корреляции k1=m-2,
k2=1,
то есть
=F(α;k1;k2)=
F(α;m-2;1)
Для данного случая k1=11,
k2=1
при уровне значимости α=0,05 и
критическое значение критерия Фишера-Снедекора =4,84.
Если выполняется неравенство >,
то с вероятностью α=0,05 будет
неверно отвергать гипотезу о значимости коэффициента парной линейной
корреляции. Если <,
делается вывод о том, что с вероятностью p=1-α=1-0,05=0,95
исследуемая
корреляционная связь случайных величин незначима и ею можно пренебречь.
Проанализировав полученные критерии Фишера,
следующие корреляционные взаимосвязи можно считать незначимыми, так как не
удовлетворяют неравенству >,:
между коэффициентом несчастных случаев (Кс.и)
и количеством обученных по охране труда в организациях Удмуртской Республики
(е);
между числом дней нетрудоспособности (N)
и количеством обученных по охране труда в организациях Удмуртской Республики
(е);
между числом пострадавших со смертельным исходом
(m2)
и количеством обученных по охране труда в организациях Удмуртской Республики
(е).
2. Проверка значимости коэффициента
корреляции с помощью критерия Стьюдента.
Критические значения tкр(α;k)
при заданном уровне значимости α
и заданном числе степеней свободы k=n-2
находятся в таблице распределения Стьюдента.
Для данного случая k=n-2=13-2=11
при заданном уровне значимости α=0,05
критическое значение tкр=2,201.
Опытное значение критерия Стьюдента вычисляется
по следующей формуле:
. (11)
Таблица 16 - Вычисленные значения критерия
Стьюдента
|
|
|
L
|
S
|
d
|
е
|
I
|
V
|
Vp
|
m1
|
14,87
|
-5,32
|
-6,48
|
-2,48
|
-7,81
|
-8,77
|
-9,91
|
m2
|
9,17
|
-5,03
|
-4,42
|
-2,04
|
-4,52
|
-5,06
|
-5,52
|
m3
|
5,96
|
-4,14
|
-8,84
|
-3,46
|
-12,76
|
-13,61
|
-16,13
|
N
|
6,83
|
-3,66
|
-4,07
|
-2,16
|
-4,03
|
-5,00
|
-5,11
|
Kч
|
9,22
|
-5,61
|
-8,43
|
-2,84
|
-9,45
|
-12,17
|
-15,27
|
Kл
|
5,95
|
-4,82
|
-4,04
|
-1,94
|
-3,71
|
-4,23
|
-4,62
|
Kп.з.
|
2,94
|
-2,70
|
-6,02
|
-4,08
|
-6,49
|
-6,58
|
-6,73
|
Kт
|
-5,71
|
5,09
|
7,99
|
2,58
|
8,61
|
8,48
|
8,90
|
Если , случайные величины X и Y с
вероятностью p=1-α следует
считать независимыми.
Если , случайные величины X и Y связаны
линейной корреляционной зависимостью.
Следующие величины следует считать
независимыми, так как не удовлетворяют неравенству :
число пострадавших со смертельным
исходом (m2) и
количество обученных по охране труда в организациях Удмуртской Республики (е);
коэффициент частоты смертельных
исходов (Кс.и.) и количество обученных по охране труда в организациях
Удмуртской Республики (е);
число дней нетрудоспособности (N) и
количество обученных по охране труда в организациях Удмуртской Республики (е).
6. Корреляционно-регрессионный
анализ
Принятие решений на основе анализа опытных
данных опирается на использование вероятностных законов распределения и
корреляционно-регрессионного анализа. В отдельных случаях анализ опытных данных
дает основания предположить, что между некоторыми из выборочных совокупностей
существует определенная зависимость.
Отличие от нуля коэффициента корреляции rxy
двух случайных величин X
и Y означает, что эти
две величины связаны линейной зависимостью и можно ставить вопрос о поиске вида
этой зависимости. Соответствующее уравнение, описывающее функциональную
зависимость величины Y
от величины X, называется
уравнением регрессии величины X
на случайную величину Y.
Уравнение регрессии, построенное по выборочным
совокупностям опытных данных, позволяет с определенной вероятностью
прогнозировать поведение генеральных совокупностей исследуемых величин в рамках
некоторого горизонта прогноза и может быть использовано для расчетов с целью
принятия решений на основе установленных закономерностей.[26]
6.1 Уравнения
множественной регрессии
В условиях данной задачи будет рассматриваться
уравнение множественной регрессии. Для обнаружения значимой статистической
зависимости между случайными величинами Y,
X1,
X2,
…, Xn
ставится
задача отыскания вида этой зависимости.
В общем случае зависимость ищется в виде функции
n переменных: y=f(x1,
x2,
…, xn). Здесь =(x1,
x2,
…, xn) - n-мерная
случайная величина, y - значение
функции f(x1,
x2,
…, xn). Функцию y=f(x1,
x2,
…, xn) требуется
определить так, чтобы при каждом из значений аргумента =(x1,
x2,
…, xn) значение функции f(x1,
x2,
…, xn) было максимально
приближено к соответствующему значению случайной величины Y.
Функция f
предполагается линейно зависящей от своих аргументов и уравнение регрессии
ищется в виде:
y = a0 + a1x1
+…+ajxj +…+ anxn. (12)
Для нахождения неизвестных параметров a0,
a1,…,
an функции необходимо
решить следующую систему уравнений:
……………………………………………………………………………..
………………………………………………………………………………..
В условиях данной задачи была
выбрана система уравнений с тремя неизвестными:
где у - значения показателей производственного
травматизма и профессиональной заболеваемости (Кч, Кл., Кп.з.,
Кт);
m - объем
выборки, для данного случая m=13;
x1, x2, x3 -
показатели, которые были выбраны, исходя из результатов проверки значимости
парных коэффициентов корреляции.
В уравнении для коэффициента частоты
травматизма Кч:
x1 - средства,
израсходованные на мероприятия по охране труда в расчете на одного работающего
(S);
x2 -
инвестиции в основной капитал в фактически действовавших ценах (I);
x3 - валовой
региональный продукт (V).
В уравнении для коэффициента частоты
смертельных исходов Кл.:
x1 - средства,
израсходованные на мероприятия по охране труда в расчете на одного работающего
(S);
x2 -
инвестиции в основной капитал в фактически действовавших ценах (I);
x3 - валовой
региональный продукт (V).
В уравнении для коэффициента частоты
профессиональных заболеваний Кп.з.:
x1 - валовой
региональный продукт (V);
x2 - средства,
израсходованные на мероприятия по охране труда в расчете на одного работающего
(S);
x3 -
инвестиции в основной капитал в фактически действовавших ценах (I).
В уравнении для коэффициента тяжести
Кт:
x1 - доля
работников, занятых в условиях, не отвечающих санитарно-гигиеническим нормам,
от общего числа работающих (d);
x2 - средства,
израсходованные на мероприятия по охране труда в расчете на одного работающего
(S);
x3 - валовой
региональный продукт (V).
Система уравнений для коэффициента
частоты производственного травматизма (Кч):
52,4 = 13a0
+ 61259,9a1 + 482,8a2 + 2635,1a3
,3= 61259,9a0
+ 398462769,9a1 + 2854899,2a2 + 15804692,9a3
,9= 482,8a0
+ 2854899,2a1 + 22730,3a2 + 123312,3a3
9060,1= 2635,1a0 + 15804692,9a1 + 123312,3a2 + 674761,3a3
Система уравнений для коэффициента
частоты смертельных исходов (Кл):
1,00 = 13a0
+ 61259,9a1 + 482,8a2 + 2635,1a3
,3 =
61259,9a0 ++ 398462769,9a1 + 2854899,2a2 +
15804692,9a3
,1 = 482,8a0
+ 2854899,2a1 + 22730,3a2 + 123312,3a3
173,6= 2635,1a0 +15804692,9a1 + 123312,3a2 + 674761,3a3
Система уравнений для коэффициента
частоты профессиональных заболеваний (Кп.з.):
16 = 13a0
+ 2635,1a1 + 61259,9a2 + 482,8a3
,8= 2635,1a0
+674761,3 a1 +15804692,9 a2 + 123312,3a3
,3= 61259,9a0
+ 15804692,9 a1 + 398462769,9a2 + 2854899,2a3
527,9= 482,8a0 + 123312,3a1 + 2854899,2a2 + 22730,3a3
Система уравнений для коэффициента
тяжести (Кт):
433,0 = 13a0
+ 313,6a1 + 61259,9a2 + 2635,1a3
= 313,6a0
+ 7966,2a1 + 1650687,6a2 + 70811,3a3
,9= 61259,9a0
+ 1650687,6a1 +398462769,9 a2 + 15804692,9 a3
98476,2= 2635,1a0 +70811,3 a1 +15804692,9
a2 +674761,3 a3
Значения неизвестных параметров
системы находятся методом Крамера.
1. Решение системы уравнений для коэффициента
частоты травматизма (Кч):
Искомые значения коэффициентов
уравнения находятся по формуле:
(13)
По найденным параметрам составляется
уравнение регрессии:
Кч = 6,28 - 5,5*S - 732,3*I -845,2*V
2. Решение системы уравнений для коэффициента
частоты смертельных исходов Кс.и.:
Уравнение регрессии:
Кл = 0,1179 - 0,54*S - 1,8*I - 7,2*V
3. Решение системы уравнений для
коэффициента частоты профессиональных заболеваний Кп.з.:
Уравнение регрессии:
Кп.з. = 1,76 - 507*V + 6,74*S + 493,5*I
4. Решение системы уравнений для коэффициента
тяжести Кт:
Уравнение регрессии:
Кт = 8,757 + 0,6718d + 0,000431S + 0,03114V
6.2
Коэффициент множественной корреляции
Если частные коэффициенты корреляции модели
множественной регрессии оказались значимыми, т. е. между результативной
переменной и факторными модельными переменными действительно существует
корреляционная взаимосвязь, то в этом случае построение множественного коэффициента
корреляции считается целесообразным.
С помощью множественного коэффициента корреляции
характеризуется совокупное влияние всех факторных переменных на результативную
переменную в модели множественной регрессии.[36]
Формула для определения коэффициента корреляции
уравнения множественной регрессии через матрицу парных коэффициентов
корреляции:
. (14)
где - определитель матрицы парных
коэффициентов корреляции;
- определитель матрицы межфакторной
корреляции.
Как видно из формул, величина множественного
коэффициента корреляции зависит не только от корреляции результата с каждым из
факторов, но и от межфакторной корреляции. Рассмотренная формула позволяет
определять совокупный коэффициент корреляции, не обращаясь при этом к уравнению
множественной регрессии, а используя лишь парные коэффициенты корреляции.
Таблица 17 - Результаты расчетов множественного
коэффициента корреляции
Показатель
травмирования
|
r
|
|
|
Кч
|
0,0005
|
0,0084
|
0,967
|
Кл
|
0,0026
|
0,0084
|
0,830
|
Кп.з.
|
0,0011
|
0,0084
|
0,935
|
Кт
|
0,0021
|
0,0177
|
0,939
|
6.3 Оценка качества
построенной модели
Коэффициентом множественной детерминации R2
называется квадрат множественного коэффициента корреляции.
Коэффициент множественной детерминации
характеризует, на сколько процентов построенная модель регрессии объясняет
вариацию значений результативной переменной относительно своего среднего
уровня, т. е. показывает долю общей дисперсии результативной переменной,
объяснённой вариацией факторных переменных, включённых в модель регрессии. Чем
больше значение коэффициента множественной детерминации, тем лучше построенная
модель регрессии характеризует взаимосвязь между переменными.
Для коэффициента множественной детерминации
всегда выполняется неравенство вида:
. (15)
Следовательно, включение в линейную модель
регрессии дополнительной факторной переменной не снижает значения коэффициента
множественной детерминации.[26]
Таблица 18 - Рассчитанные коэффициенты
детерминации
Показатель
травмирования
|
Кч
|
Кл
|
Кп.з.
|
Кт
|
R2
|
0,935
|
0,690
|
0,874
|
0,881
|
Для того чтобы не допустить преувеличения
тесноты связи, применяется скорректированный индекс множественной детерминации,
который содержит поправку на число степеней свободы и рассчитывается по
формуле:
. (16)
где n - объем выборки, m - число переменных в
уравнении множественной регрессии. При небольшом числе наблюдений
нескорректированная величина коэффициента множественной детерминации R2
имеет тенденцию переоценивать долю вариации результативного признака, связанную
с влиянием факторов, включенных в регрессионную модель.
Таблица 19 - Скорректированный индекс
множественной детерминации
Показатель
травмирования
|
Кч
|
Кл
|
Кп.з.
|
Кт
|
0,9130,5870,8320,842
|
|
|
|
|
Высокие величины коэффициентов детерминации R2
указывают на то, что модели регрессии хорошо аппроксимируют исходные данные и
такими регрессионными моделями можно воспользоваться для прогноза значений
результативного показателя.
Проверить значимость (качество) уравнения
регрессии - значит установить, соответствует ли математическая модель,
выражающая зависимость между переменными, экспериментальным данным, достаточно
ли включенных в уравнение объясняющих переменных для описания зависимой
переменной. Чтобы иметь общее суждение о качестве модели, по каждому наблюдению
из относительных отклонений определяют среднюю ошибку аппроксимации. Проверка
адекватности уравнения регрессии (модели) осуществляется с помощью средней
ошибки аппроксимации, величина которой не должна превышать 12-15% (максимально
допустимое значение).
Формула для расчета средней ошибки
аппроксимации:
, (17)
где n - число
переменных в уравнении множественной регрессии; f(xi1, xi2, …, xin) - i-е расчетное
значение переменной y; - i-е опытное
значение переменной y.
Таблица 20 - Средняя ошибка аппроксимации
Показатель
травмирования
|
Кч
|
Кл
|
Кп.з.
|
Кт
|
,
%6,310,56,55,6
|
|
|
|
|
Как видно из результата расчетов, средние ошибки
аппроксимации не превышают допустимые значения в 12-15%, что говорит об
адекватности полученных моделей.
Проверка значимости коэффициентов линейного уравнения
множественной регрессии.
Проверка значимости отдельных коэффициентов
уравнения означает, что если коэффициент при некоторой переменной незначим, то
доверять влиянию этой переменной на значения результирующей функции y
нельзя. Незначимый коэффициент следует положить равным нулю, т.е.
соответствующую переменную следует исключить из дальнейшего рассмотрения.
Для проверки значимости каждого из коэффициентов
a0,
a1,…,
an используется t-статистика
Стьюдента, опытное значение которой вычисляется по формуле:
, (i=0,1,…,n), (18)
где ai -
коэффициент при переменной xi, -
среднеквадратическая ошибка этого коэффициента,
, (19)
где - среднее квадратичное отклонение
для значений переменной y; - среднее квадратичное отклонение
для значений xi; -
коэффициент множественной детерминации для уравнения регрессии в целом; -
коэффициент множественной детерминации, характеризующий зависимость между
фактором xi и
остальными факторами (x1, x2,…, xi-1, xi+1,…, xn)уравнения
регрессии.
Каждое из опытных значений
статистики сравнивают
с критическим значением (i=1,2,…,n), которое
ищется по таблице распределения Стьюдента при заданном уровне значимости α и числе
степеней свободы k, равном k=m-n-1. В данном
случае при уровне значимости α=0,05 и k=13-3-1=9 =2,26.[26]
Таблица 21 - Рассчитанные опытные
значения t -
статистики Стьюдента
Показатель
травмирования
|
|
|
a1
|
a2
|
a3
|
Кч
|
7,32
|
7,35
|
6,45
|
Кл
|
8,17
|
8,67
|
7,23
|
Кп.з.
|
6,17
|
6,39
|
2,29
|
Кт
|
7,72
|
3,46
|
Если > , то гипотеза о значимости
коэффициента ai не
отвергается, и соответствующая переменная xi остается в
уравнении. В противном случае коэффициент ai считается
незначимым и соответствующую ему переменную следует исключить из уравнения
регрессии. Таким образом, сравнив полученные опытные значения с
критическим , можно
сделать вывод, что незначимых коэффициентов во всех четырех уравнениях нет.
6.4 Проверка
значимости линейного уравнения множественной регрессии в целом
Если окажется, что при заданном уровне
значимости α уравнение
незначимо, то пользоваться им нельзя, а найденной зависимостью следует
пренебречь.
Для проверки значимости уравнения регрессии
используется опытная F-статистика
Фишера:
, (20)
где m
- объем выборки; n - число
переменных в уравнении множественной регрессии; f(xi1,
xi2,
…, xin) - i-е
расчетное значение переменной y;
-
среднее опытных значений случайной величины Y.[26]
Полученные опытные значения критерия
Фишера сравниваются с критическими значениями =F(α;k1;k2)
при выбранном уровне значимости α. Число
степеней свободы k1
=
m - n
- 1, k2
= n.
При выбранном уровне значимости α=0,05
и
числе степеней свободы k1
=
13 - 3 - 1= 9, k2
= 3 =8,81
Таблица 22 - Рассчитанные опытные значения
критерия Фишера
Показатель
травмирования
|
Fоп
|
Кч
|
42,74
|
Кл
|
21,24
|
Кп.з.
|
20,8
|
Кт
|
22,4
|
При сравнении опытных значений критериев Фишера
с критическим (при уровне значимости α=0,05 Fкр=8,81),
все они удовлетворяют неравенству Fоп
> Fкр
и делается вывод, что с вероятностью p=1-α=0,95
все
уравнения значимы, и мы получаем определенные основания доверять построенным
уравнениям регрессии.
6.5 Оценка точности
линейного уравнения множественной регрессии
Заключительная статистическая процедура - оценка
точности построенных уравнений регрессии.
Оценка близости опытных значений yi случайной
величины Y и ее
расчетных значений f(xi),
получаемых с помощью уравнения линейной регрессии, выполняется с помощью
среднеквадратической погрешности по следующей формуле:
. (21)
Таблица 23 - Результаты расчета
среднеквадратичной погрешности уравнений
|
|
|
Кч
|
1,223
|
0,318
|
Кл
|
0,003
|
0,016
|
Кп.з.
|
10,984
|
0,955
|
Кт
|
111,746
|
3,045
|
7. Расчет
прогнозных значений
Для вычисления прогнозных значений воспользуемся
данными из Программы социально-экономического развития республики.
Целесообразно сравнить планируемые значения показателей с фактическими.
Результаты сравнения представлены в таблицах №24,25
Таблица 24- Выполнение основных показателей
Программы социально-экономического развития Удмуртской Республики на 2005 -
2013 годы - Валовой региональный продукт
Год
|
Факт,
млрд.руб.
|
План,
млрд.руб.
|
2005
|
140
|
120,7
|
2006
|
164,8
|
132,9
|
2007
|
205,8
|
145,9
|
2008
|
240,3
|
159,4
|
2009
|
226,1
|
174,1
|
2010
|
264,5
|
243,8
|
2011
|
318,8
|
264,5
|
2012
|
335,4
|
346,6
|
2013
|
404,8
|
414,5
|
Таблица 25 - Выполнение основных показателей
Программы социально-экономического развития Удмуртской Республики на 2005 -
2013 годы - Инвестиции в основной капитал
Год
|
Факт,
млрд.руб.
|
План,
млрд.руб.
|
2005
|
26,9
|
17,7
|
2006
|
34,3
|
20,1
|
2007
|
44,5
|
22,9
|
2008
|
49,9
|
26,0
|
2009
|
41,3
|
29,1
|
2010
|
42,3
|
40,1
|
2011
|
60,9
|
50,2
|
Таблица 25 - Выполнение основных показателей
Программы социально-экономического развития Удмуртской Республики на 2005 -
2013 годы - Инвестиции в основной капитал(продолжение)
Год
|
Факт,
млрд.руб.
|
План,
млрд.руб.
|
2012
|
57,2
|
65,6
|
2013
|
71,8
|
72,3
|
Сравнение показало, что данными из Программы
социально-экономического развития можно воспользоваться, так как плановые
значения чаще ниже фактических.
Статистические, полученные расчетные и
прогнозные данные представлены в таблице.
Таблица 26 - Статистические, расчетные и
прогнозные показатели травматизма
Год
|
Кч
|
Кч
прог.
|
Кл
|
Кл
прог
|
Кп.з.
|
Кп.з.прог
|
Кт
|
Кт
прог
|
2001
|
6,170
|
5,543
|
0,124
|
0,107
|
1,579
|
1,598
|
23,000
|
23,309
|
2002
|
5,586
|
5,448
|
0,103
|
0,106
|
1,474
|
1,525
|
23,500
|
23,978
|
2003
|
5,273
|
5,334
|
0,118
|
0,104
|
1,613
|
1,493
|
25,600
|
25,931
|
2004
|
5,012
|
5,192
|
0,081
|
0,100
|
1,688
|
1,475
|
24,417
|
25,136
|
2005
|
4,689
|
4,761
|
0,085
|
0,096
|
1,449
|
1,353
|
24,660
|
26,106
|
2006
|
4,308
|
4,454
|
0,087
|
0,090
|
1,313
|
1,301
|
26,121
|
29,03
|
2007
|
3,872
|
3,984
|
0,103
|
0,082
|
1,062
|
1,219
|
35,826
|
32,353
|
2008
|
3,492
|
3,582
|
0,052
|
0,072
|
1,171
|
1,158
|
43,105
|
35,621
|
2009
|
3,089
|
3,792
|
0,062
|
0,076
|
1,237
|
1,154
|
39,604
|
35,45
|
2010
|
3,121
|
3,121
|
0,041
|
0,040
|
1,369
|
1,380
|
38,899
|
39,261
|
2011
|
3,084
|
2,754
|
0,047
|
0,058
|
1,150
|
0,917
|
41,300
|
42,047
|
2012
|
2,500
|
2,624
|
0,034
|
0,055
|
0,621
|
0,834
|
40,665
|
44,417
|
2013
|
2,203
|
1,874
|
0,070
|
0,044
|
0,529
|
0,624
|
46,544
|
46,809
|
2014
|
-
|
1,812
|
-
|
0,033
|
-
|
0,714
|
-
|
47,940
|
2015
|
-
|
1,496
|
-
|
0,027
|
-
|
0,637
|
-
|
50,028
|
2016
|
-
|
1,179
|
-
|
0,020
|
-
|
0,561
|
-
|
52,116
|
2017
|
-
|
0,862
|
-
|
0,014
|
-
|
0,484
|
-
|
54,204
|
8. Расчет риска
травмирования
В настоящее время многие авторы отмечают, что
поток несчастных случаев на производстве распределяется по дискретному
пуассоновскому закону. Этот закон предполагает, что указанный поток обладает
свойствами стационарности - интенсивность несчастных случаев не зависит от
времени работы, ординарности - вероятность возникновения двух несчастных
случаев и более на малом отрезке времени является величиной более высокого
порядка малости по сравнению с вероятностью одного случая травмирования,
отсутствием последействия - на любых двух не перекрещивающихся отрезках времени
числа проявлений несчастных случаев независимы. Стационарность доказывается
тем, что число случаев травмирования по годам меняется незначительно,
ординарность - тем, что групповые несчастные случаи происходят крайне редко.
Отсутствие последействия подтверждается тем, что число несчастных случаев,
например, в декабре не зависит от того, сколько их было в предыдущем месяце,
т.е. несчастные случаи появляются на оси времени по причинам, не зависимым от
самих этих случаев. Вероятность несчастных случаев, по данным В.М. Минько,
может быть представлена формулой
, (22)
где Р(к) - вероятность k несчастных случаев, k =
0, 1, 2, 3, ...; N - число работающих; t - продолжительность работы, лет; β
- повышающий
коэффициент, использующийся тогда, когда есть основания считать данные о
несчастных случаях заниженными. Имеются результаты исследований, из которых
вытекает, что 1 < β < 5.
Выражение позволяет получать прогностические
оценки различных событий, связанных с производственным травматизмом. Может быть
получен и риск R травмирования
, (23)
где Р(о) вычисляется по выражению. При k = 0
выражение получает вид
, (24)
В выражение вместо коэффициента частоты можно
ввести коэффициент смертности Кс.и., что позволяет записать
, (25)
где - вероятность Кс.и. (Кс.и. = 0, 1,
2, ...) несчастных случаев со смертельным исходом.
Таблица 27 - Риски травмирования
Год
|
Риск
травмирования R
|
|
Кч
|
Кл
|
Кп.з.
|
2014
|
1812*10-6
|
32,897*10-6
|
714*10-6
|
2015
|
1496*10-6
|
26,582*10-6
|
637*10-6
|
2016
|
1179*10-6
|
20,267*10-6
|
561*10-6
|
2017
|
862*10-6
|
13,952*10-6
|
484*10-6
|
9. Апробация
полученных моделей
Апробация полученных моделей на основе
ретроспективных данных, полученные прогнозы и риски представлены на графиках.
Рис.17. Результаты прогнозирования коэффициента
частоты несчастных случаев Кч и риска травмирования на 2014-2017
г.г.
Рис. 18. Результаты прогнозирования коэффициента
частоты смертельных исходов Кл. и риска смертельного исхода на
2014-2017 г.г.
Рис. 19. Результаты прогнозирования коэффициента
частоты профзаболеваемости Кп.з. и риска профзаболеваемости на 2014
- 2017 г.г.
Рис. 20. Результаты прогнозирования коэффициента
тяжести травматизма Кт на 2014 - 2017 г.г
Заключение
1. В результате проведенного исследования
были собраны статистические данные показателей травматизма в Удмуртской
Республике за временной промежуток с 2001 по 2013 годы.
. Был проведен анализ статистических
данных производственного травматизма в Удмуртской Республике за временной
промежуток с 2001 по 2013 годы
Проведенный корреляционный анализ показал, что
коэффициент частоты травматизма (Кч) имеет обратную корреляционную
зависимость с валовым региональным продуктом (V)
с коэффициентом корреляции r=-0,965,
с количеством инвестиций в основной капитал (I)
с коэффициентом корреляции r=-0,943
и количеством средств, израсходованных на мероприятия по охране труда (S)
с коэффициентом корреляции r=-0,861,
т.е. с увеличением социально-экономических показателей коэффициент частоты
травматизма будет уменьшаться.
Коэффициент частоты смертельных исходов (Кл.)
имеет аналогичную корреляционную зависимость: с валовым региональным продуктом
(V) коэффициент
корреляции r=-0,787, с
инвестициями в основной капитал (I)
r=-0,745 и со
средствами, израсходованными на мероприятия по охране труда (S)
r=-0,824.
Коэффициент частоты профессиональных заболеваний
(Кп.з.) взаимосвязан с инвестициями в основной капитал (I)
с r=-0,890, с валовым
региональным продуктом (V)
с коэффициентом корреляции r=-0,893
и с количеством средств, израсходованных на мероприятия по охране труда (S)
с коэффициентом корреляции r=-0,631.
Коэффициент тяжести (Кт) имеет
положительную корреляционную зависимость с долей работников, занятых в
условиях, не отвечающих санитарно-гигиеническим нормам (d)
с коэффициентом корреляции r=0,923,
со средствами, израсходованными на мероприятия по охране труда в расчете на
одного работающего (S) c
коэффициентом корреляции r=0,838
и валовым региональным продуктом (V)
c коэффициентом
корреляции r=0,931
. Разработаны математические модели,
основанные на статистике травматизма за 2001 - 2013 годы в Удмуртской
Республике, предназначенные для вычисления прогнозных показателей травматизма и
профессиональной заболеваемости, показавших удовлетворительную сходимость
полученных прогнозов с фактическими значениями:
Кч
= 6,28 - 5,5*S - 732,3*I -845,2*V
Кл
= 0,1179 - 0,54*S - 1,8*I - 7,2*V
Кт
= 8,757 + 0,6718d + 0,000431S + 0,03114V
Кп.з. = 1,76 - 507*V + 6,74*S + 493,5*I
Результаты прогнозирования показателей
травматизма и профзаболеваемости на 2014-2017 годы на основе данных из
Программы социально-экономического развития Удмуртской Республики:
Кч : 2014 -1,812 ±0,318, 2015 -
1,496±0,318, 2016 -1,179±0,318, 2017-0,862±0,318
Кл.: 2014 - 0,033±0,016, 2015 - 0,027
±0,016, 2016 - 0,020±0,016, 2017-0,014±0,016
Кп.з:. 2014 - 0,714±0,955,
2015 - 0,637±0,955, 2016 -0,561±0,955, 2017-0,484±0,955
Кт: 2014 - 47,940±3,045, 2015 -
50,028±3,045, 2016 - 52,116±3,045, 2017-54,204±3,045
. Для полученных прогнозов вычислен риск
травмирования, смертельного исхода и регистрации профессионального заболевания R
для одного человека в прогнозируемый год.
Для коэффициента частоты травматизма (Кч)
риск (R) составит: 2014г.
- 1812*10-6, 2015 г. - 1496*10-6, 2016 г. - 1179*10-6,
2017-862*10-6
Для коэффициента частоты смертельных исходов (Кл)
риск (R) составит: 2014 г.
- 32,897*10-6, 2015 г. - 26,582*10-6, 2016 г. - 20,267*10-6,
2017-13,952*10-6
Для коэффициента частоты профессиональных
заболеваний (Кп.з). риск (R)
составит: 2014г. -714*10-6, 2015 г. - 637*10-6,2016г.-561*10-6,2017-
484*10-6
Литература
1. Конституция Российской
Федерации (принята всенародным голосованием 12.12.1993 г. ) // «Российская газета» от 25 декабря 1993 г.
2. Федеральный закон от 30.12.2001
г. № 197-ФЗ «Трудовой Кодекс Российской Федерации».
3. Федеральный закон от 18 июля
2011 г. № 238-ФЗ «О внесении изменений в Трудовой кодекс РФ».
4. Об утверждении положения о
расследовании и учете профессиональных заболеваний: Постановление Правительства
РФ от 15 декабря 2000 г. № 967.
5. Об утверждении положения об
оплате дополнительных расходов на медицинскую, социальную и профессиональную
реабилитацию застрахованных лиц, получивших повреждение здоровья вследствие
несчастных случаев на производстве и профессиональных заболеваний :
Постановление Правительства РФ от 15 мая 2006 г. № 286 : [утв. Постановлением
Правительства РФ от 15 мая 2006 г. № 286 (с изм. от 27 октября 2008 г.)].
6. Об утверждении форм документов,
необходимых для расследования и учета несчастных случаев на производстве, и
Положения об особенностях расследования несчастных случаев на производстве в
отдельных отраслях и организациях : Постановление Минтруда РФ от 24 октября
2002 г. № 73 : [зарег. в Минюсте РФ 05.12.2002 г. Рег. № 3999].
7. Об утверждении инструкции о
порядке заполнения формы программы реабилитации пострадавшего в результате
несчастного случая на производстве и профессионального заболевания, утвержденной
Постановлением Министерства труда и социального развития РФ от 18 августа 2001 г. № 56 : Постановление Минтруда
РФ от 30 января 2002 г. № 5 : [зарег. в Минюсте РФ 14 февраля 2002 г. Рег. №
3246].
8. Об утверждении форм документов
о результатах установления Федеральными государственными учреждениями
медико-социальной экспертизы степени утраты профессиональной трудоспособности в
процентах и рекомендациях по их заполнению: приказ Минздравсоцразвития РФ от 20
октября 2005 г. № 643.
9. Об утверждении типового положения
о комитете (комиссии) по охране труда: приказ Минздравсоцразвития РФ от 29 мая
2006 г. № 413.
10. О форме документов, необходимых
для расследования несчастных случаев на производстве: приказ
Минздравсоцразвития РФ от 15 апреля 2005 г. № 275.
11. О совершенствовании системы
расследования и учета профессиональных заболеваний в Российской Федерации:
приказ Минздрава РФ от 28 мая 2001 г. № 176.
12. ГОСТ 12.0.230-2007 ССБТ.
Системы управления безопасностью и гигиеной
труда. Общие требования : [введ. 01.07.2009 г.]. - М. :
Стандартинформ, 2007.
13. ГОСТ Р 12.0.007-2009 ССБТ.
Система управления охраной труда в организации. Общие требования по разработке,
применению, оценке и совершенствованию.
15. Декларация Российского научного
общества анализа риска «Об экономической оценке жизни среднестатистического
человека» / Проблемы анализа риска - М., 2007. - т. 4 - № 2.
16. Декларация Российского научного
общества анализа риска «О предельно-допустимых уровнях риска» / Проблемы
анализа риска - М., 2006. - т. 3 - № 2.
17. Генеральное соглашение между
общероссийскими объединениями профсоюзов, общероссийскими объединениями
работодателей и Правительством Российской Федерации на 2011 - 2013 годы.
18. Регионы России. Социально-экономические
показатели // Федеральная служба государственной статистики Российской
Федерации. [Электронный ресурс] - URL:
<http://www.gks.ru>
19. Труд и занятость в России // Федеральная служба
государственной статистики Российской Федерации. [Электронный ресурс] - URL:
<http://www.gks.ru>
20. Срочная информация по актуальным вопросам //
Территориальный орган федеральной службы государственной статистики по
Удмуртской Республике. [Электронный ресурс] - URL:
<http://udmstat.gks.ru/public/>
21. Программа социально-экономического развития
Удмуртской Республики на 2010 - 2014 годы «Приложение к Закону Удмуртской
республики от 18 декабря 2009 года №68-РЗ» // Правительство Удмуртской
Республики. [Электронный ресурс]. -
URL:
<http://www.udmurt.ru/region/economic/program_2014/>
22. Доклад «Состояние условий и
охраны труда в Удмуртской Республике в 2009 году и меры по их улучшению» //
Министерство труда Удмуртской Республики. [Электронный ресурс] - URL:
http://mintrud.udmurt.ru/
23. Доклад «Состояние условий и
охраны труда в Удмуртской Республике в 2010 году и меры по их улучшению» //
Министерство труда Удмуртской Республики. [Электронный ресурс] - URL:
http://mintrud.udmurt.ru/
24. Доклад «Состояние условий и
охраны труда в Удмуртской Республике в 2011 году и меры по их улучшению» //
Министерство труда Удмуртской Республики. [Электронный ресурс] - URL:
http://mintrud.udmurt.ru/
25. Доклад «Состояние условий и
охраны труда в Удмуртской Республике в 2012 году и меры по их улучшению» //
Министерство труда Удмуртской Республики. [Электронный ресурс] - URL:
http://mintrud.udmurt.ru/
26. Лялькина Г.Б. Математическая
обработка результатов эксперимента: учеб. Пособие / Г.Б. Лялькина, О.В.
Бердышев - Пермь: Изд-во Перм.нац.исслед.политехн.ун-та, 2013. - 78 с.
27. Руководство по системам
управления охраной труда. - МОТ - СУОТ, 2001.
28. Вентцель, В.Н. Теория
вероятностей: Учеб. Для вузов/Е.С. Вентцель. - 8-е изд., стер.- М.: Высш. Шк.,
2002. - 575 с.: ил.
29. Кобзарь, А.П. Прикладная
математическая статистика. - М.: Физматлит, 2006. - 816 с.
30. Айвазян С.А., Мхитарян В.С
Прикладная статистика в задачах и упражнениях. - М.:ЮНИТИ, 2001.- 207с.
31. Орлов, А.И. Эконометрика:
Учебное пособие для ВУЗов. - М.: Издательство «Экзамен», 2002. - 576 с.
32. Андерсон, Т. Статистический
анализ временных рядов. - М.: Мир,
1976. - 756 с.
33. Гмурман В.Е Теория вероятностей
и математическая статистика: учеб. пособие.- М.: Высшее образование, 2004-
479с.
34. Бокс, Дж., Дженкинс, Г. Анализ
временных рядов, прогноз и управление. - М.: Мир, 1974. - Ч. 1. - 404 с.
35. Коваленко И.Н., Филипова А.А
Теория вероятностей и математическая статистика: учеб. пособие.- М.: Высшее
образование, 1982,- 256с.
36. Эконометрика: учеб. под. ред.
И.И. Елисеевой - М.: Финансы и статистика, 2001- 344с.
37. Лукашин, Ю.П. Адаптивные методы
краткосрочного прогнозирования временных рядов. - М.: Финансы и статистика,
2003. - 416 с.
38. Каримов Р.Н. Обработка
экспериментальной информации. - Саратов: Изд-во СарГУ, 2001. - Ч.4. - 103 с.
39. Колесниченко В.И. Обработка и
представление результатов эксперимента / Перм. гос. техн. ун-т. - Пермь, 2000.
- 74 с.
40. Маркин Н.С. Основы теории
обработки результатов измерений. - М.: Изд-во стандартов, 1991. - 173 с.
41. Гмурман, В.Е. Теория
вероятностей и математическая статистика. - М.: Высш. школа, 1977. - 479 с.
42. Бородич, С.А. Вводный курс
эконометрики: Учебное пособие - Мн.: БГУ, 2000. - 354 с.
43. Севастьянов, Б. В.
Информирование о риске повреждения здоровья - обязанность работодателя / Б. В.
Севастьянов, Е. Б. Лисина // Библиотека инженера по охране труда. - 2007. - №
2. - С. 70-83.
44. Салтыков, А. М. Система
управления охраной труда в Удмуртской Республике // Справочник специалиста по
охране труда. - 2002. - № 11. - С. 6-14.
45. Севастьянов, Б. В. Система
менеджмента в области охраны труда и предупреждения профессиональных
заболеваний - составная часть общей системы менеджмента организации //
Примышленная и экологическая безопасность. - 2007. - № 9. - С. 13-17.
46. Севастьянов, Б. В. Опасность
травмирования и защита от него / Б. В. Севастьянов, Е. Б. Лисина, Н. В.
Селюнина, Р. М. Хазеев // Промышленная и экологическая безопасность. - 2009. -
№ 11. - С. 34-38.
47. Севастьянов, Б.В. Организация и
управление охраной труда в Удмуртской Республике : монография /
Б.В.Севастьянов, Р.О.Шадрин - Ижевск : Изд-во ИжГТУ, 2011. - 132 с.
48. Севастьянов, Б.В. Управление
безопасностью труда : учеб. пособие для студ. высш. учеб. заведений. В 2 ч. /
Б.В.Севастьянов, Е.Б.Лисина,
И.Г.Тюрикова ; под общ. ред. проф. Б.В.Севастьянова. - Ч. I.
Государственное управление охраной труда. - Ижевск : Изд-во ИжГТУ, 2010. - 296
с. - (Безопасность технологических процессов и производств).
49. Горлов, А. Г. Охрана труда в
муниципальных образованиях // Справочник специалиста по охране труда. - 2002. -
№ 4. - С. 12-15.
50. Соловьев, А. П. Система
управления охраной труда в Российской Федерации // Справочник специалиста по
охране труда. - Октябрь, 2001.
51. Охрана труда : Правовое
регулирование. Практика. Основные документы / под ред. Ю. Л. Фадеева. - М. :
Эксмо, 2008. - 224 с.
52. Севастьянов, Б. В. Исследование
профессиональных рисков возникновения производственных травм в организациях
Удмуртской Республики. Отчет о НИР по договору от 20.06.2006 г. № БЖД-1-06/11
(заключ.) / Б. В. Севастьянов, А. П. Тюрин,
Э. А. Поликарпов. - Ижевск : ИжГТУ, 2006.
53. Севастьянов, Б. В. Анализ
причин возникновения несчастных случаев с причинением тяжкого вреда здоровью и
смертельными исходами в отраслях экономики Удмуртской Республики (по данным
организационно-технических судебных экспертиз) / Б. В. Севастьянов, В. А.
Никешкин, П. С. Васильева, Е. А. Здобяхина // Примышленная и экологическая
безопасность. - 2007. - № 2. - С. 56-59.
54. Отчет по НИР по контракту с
Министерством труда УР от 23 августа 2010 № 28/МТ-10 на тему «Разработка модели
прогнозирования и управления рисками повреждения здоровья работающими»,
исполнители: Б.В.Севастьянов, А.П.Тюрин, Р.О.Шадрин, И.Г.Русяк, В.Г. Суфиянов и
И.В.Васильева.
55. Ожегов, С. И. Словарь русского
языка. - Екатеринбург : Изд-во «Урал-Советы», 1994. - 800 с. \
Похожие работы на - Анализ и прогнозирование производственного травматизма
|