МЕТОД ОЦЕНКИ ВЛИЯНИЯ ПРИЧИН ВОЗГОРАНИЯ НА ОБСТАНОВКУ С ПОЖАРАМИ НА ПРОМЫШЛЕННЫХ ПРЕДПРИЯТИЯХ РЕГИОНА В ИНТЕРЕСАХ ГОСУДАРСТВЕННОГО ПОЖАРНОГО НАДЗОРА
М.В. Довиденко; А.П. Корольков; А.Г. Филиппов. Санкт-Петербургский университет ГПС МЧС России
Рассматривается метод оценки влияния причин возгорания на число пожаров, основанный на интервальных статистических оценках и корреляционном анализе.
Ключевые слова: пожарная опасность, математическая статистика, интервальные статистические оценки, выборочная дисперсия, среднеквадратическое отклонение, корреляционный анализ
METHOD OF AN ESTIMATION OF INFLUENCE OF THE REASONS OF IGNITION ON CONDITIONS WITH FIRES AT THE INDUSTRIAL ENTERPRISES OF REGION IN INTERESTS SFS
M.V. Dovidente; A.P. Korolkov; A.G. Filippov.
Saint-Petersburg university of State fire service of EMERCOM of Russia
In article the method of an estimation of influence of the reasons of ignition on number of the fires, based on interval statistical estimations and the correlation analysis is considered.
Key words: fire danger, mathematical statistics, interval statistical estimations, selective dispersion, correlation analysis
Для выработки полноценных предложений по повышению эффективности органов Государственной противопожарной службы в целом и Государственного пожарного надзора, в частности, необходимо знать причины, оказывающие существенное влияние на показатели эффективности.
Зависимость степени пожарной опасности от причин носит вероятностный характер. Например, появление одной и той же причины в некоторых случаях ведет к возникновению пожара, а в других - нет. Поэтому для количественной оценки степени влияния различных причин на показатели пожарной опасности региона нужны методы математической статистики.
Авторами рассматривается метод оценки влияния причин возгорания на число пожаров, основанный на интервальных статистических оценках и корреляционном анализе. Пусть в нашем распоряжении имеются данные по причинам пожаров. В табл. 1. приведены частоты причин возникновения пожаров по годам.
В этой таблице: m - количество лет; n - количество причин; xy - количество возгораний в i-м году по j-й причине; yi - количество возгораний в i-м году. Найдем: - среднее значение количества возгораний по каждой причине
- выборочную дисперсию количества возгораний по каждой причине
47
- среднее количество возгорании по всем причинам
Уср=1ш1=1шу1;
- выборочную дисперсию количества возгораний по всем причинам
;=1
- коэффициент корреляции между количеством возгорания по каждой причине и общим количеством возгораний.
■ Л. -л- 1 .
т — 1
5 5 -'.Г "Ъ
. Таблица 1
Годы Причины Всего пожаров
П1 П2 Пз П4 Пп-1 Пп
1 х1,1 Х1,2 Х1,з Х1,4 Х1,п-1 Х1,п У1
2 Х2,1 Х2,2 Х2,3 Х2,4 Х2,п-1 Х2,п У2
3 Х3,1 х3,2 Хзз Х3,4 Х3,п-1 Х3,п Уз
ш Хп,1 хп,2 Хп,3 Хп,4 Хш,п-1 Хш,п4 Уш
Первичные описательные статистики причин возгораний представлены в табл. 2. Результаты вычислений, отображенные в таблице, являются величинами случайными, так как они рассчитаны на основании данных, выбранных из генеральной совокупности случайным образом.
............................................................................................................................................. Таблица 2
Статистики Причины По всем причинам
П1 П2 Пз П4 Пп-1 Пп
Среднее выборочное '' "СР1 *срж ** *СР/П Уср
Выборочное среднеквадратическое отклонение Ву
Коэффициент корреляции у) К«а-у) 'Ога-у) —
Очевидно, что если из этой генеральной совокупности взять другие данные, то и статистики будут несколько другие. Поэтому, после вычисления статистик необходимо оценить вероятности того, что истинные, неизвестные нам коэффициенты корреляции отличаются от нуля и вычислить размер интервалов, внутри которых находятся неизвестные нам математические ожидания генеральных совокупностей.
Наличие связи между пожаром и причиной оценивается с помощью выборочного коэффициента корреляции, рассчитываемого по формуле (1). Однако мы установили, что этот коэффициент является величиной случайной. Поэтому в результате совпадения ряда случайных обстоятельств выборочный коэффициент корреляции может оказаться не равным нулю, в то время как истинный, неизвестный нам коэффициент равен нулю. Для контроля этого предположения мы используем методы теории проверки статистических гипотез.
Используем нулевую гипотезу (Но) о том, что на самом деле связи между событием, которое мы назвали причиной и пожаром нет. В формальном виде эта гипотеза записывается выражением
В математической статистике [1, 2] показано, что в качестве критерия проверки нулевой гипотезы в данном случае целесообразно использовать случайную величину
Эта случайная величина при справедливости нулевой гипотезы, распределена по закону Стьюдента с к= (т - 2) степенями свободы, где т - объем выборки. В [2, 3, 4] показано, что принятые гипотезы проверяются с помощью критерия Фишера. Вначале определяется величина расчетного критерия Фишера по формулам:
Затем определяются степени свободы Vi = n-1, v2 = m-1, где n - объем выборки значений числителя, m - объем выборки значений знаменателя в формуле (3).
Затем задается уровень значимости, например а=0,05. Используя вычисленные значения v1, v2 и а определяется критическое значение критерия Фишера (Ркрит) по специальным таблицам или с помощью математических пакетов прикладных программ. Например, в Excel для этих целей используется функция КОБР.ПХ(а, v1, v2).
Если критическое значение критерия Фишера окажется меньше расчетного, то основная гипотеза отвергается, и считается, что дисперсии случайных величин x и y различны. Тогда расчетный критерий Стьюдента определяется по формуле:
Альтернативная гипотеза запишется в виде
2
3
1расч=хв-увБх2п+8у2т ,
а число степеней свободы равно максимально возможному целому, меньшему, чем
Если критическое значение критерия Фишера окажется больше расчетного, то основная гипотеза принимается, и считается, что дисперсии случайных величин х и у одинаковы. Тогда расчетный критерий Стьюдента определяется по формуле:
1расч=(хв-ув)пш(п+ш-2)(п+ш)п-18х2+ш- 1Бу2 , а число степеней свободы по формуле
V = п + т -г
Вычисление критического значения критерия Стьюдента осуществляется с помощью соответствующих статистических таблиц или математических пакетов. Таким образом последовательность расчетов при проверке гипотезы в отношении коэффициентов корреляции между причинами пожаров и пожарами (табл. 3) включает следующие этапы: в начале проверки выбирается уровень значимости а, затем для этого уровня и числа степеней свободы к вычисляется критическое значение критерия Стьюдента 1:КРИТ, затем это критическое значение сравнивается с расчетным, вычисленным по формуле (2).
...............................................................................................................Таблица 3
Причины пожаров
П1
П2
Пз
Пд
Пп-1
Пт
Вычисление Т расчетных
Т1
Т2
Тз
Тд
Тп-1
Вычисление Ткрит для а=0,1
Вычисление Ткрит для а=0,05
Вычисление Ткрит ддя а=0,01
Сравнение Т1 и Ткрит, формулировка выводов для а=0,1
Сравнение Т1 и Ткрит, формулировка выводов для а=0,05
Сравнение Т1 и Ткрит, формулировка выводов для а=0,01
1
п
Если критическое значение критерия Стьюдента окажется больше расчетного, то принимается гипотеза Н0, если меньше, то принимается альтернативная гипотеза. О величине интервала можно судить по величине выборочного среднего и выборочного среднеквадратического отклонения. Но поскольку обе эти величины есть величины случайные, то и сам интервал - величина случайная.
Поэтому нельзя говорить о величине интервала, внутри которого находится неизвестное нам математическое ожидание. Можно только говорить о вероятности накрытия неизвестного математического ожидания данным интервалом.
Таким образом, математическое ожидание оценивается не числом, а интервалом. Такие оценки математического ожидания называются интервальными [3]. Вероятность накрытия математического ожидания интервалом называется доверительной вероятностью или надежностью интервальной оценки (у), а величина интервала - точностью (1у). Величина ^ находится из статистических таблиц или компьютерных программ по известным величинам ш и у (объем выборки и надежности соответственно). В математической статистике доказано, что для оценки точности и надежности интервальной оценки математического ожидания целесообразно использовать случайную величину 1 которая имеет распределение Стьюдента с к = (ш-1) степенями свободы:
- а
где Хср - выборочное среднее; Б - «исправленное» выборочное среднеквадратическое отклонение; ш - объем выборки; а - неизвестное математическое ожидание генеральной совокупности.
Поскольку плотность распределения Стьюдента 8(1;, ш) - четная функция, то вероятность выполнения неравенства
| - а
■ < *■,
вычисляется по формуле
Отсюда, выполнив элементарные преобразования, перейдем к двойному неравенству
Таким образом, с помощью распределения Стьюдента мы нашли доверительный
интервал ( ^ V''! ' Р V"?), покрывающий неизвестное нам математическое ожидание с вероятностью (надежностью) у.
Задавшись надежностью у и вычислив интервалы для каждой причины, мы для каждой причины возгорания найдем статистически обоснованный интервал математического ожидания числа пожаров. Сведем вычисления интервальных оценок в таблицу. В табл. 4. приведены данные о точности и надежности оценок математического ожидания. Теперь мы можем оценить минимальное и максимальное число пожаров произошедших при наличии каждой причины. Однако, вычислять интервалы мы будем только для тех причин, для которых коэффициент корреляции (1) статистически значимо отличается от нуля.
............................................................................................................................................. Таблица 4
Статистики Причины По всем
П1 П2 Пз П4 . • Пп-1 Пп причинам
Надежность у
Параметр ^ Стьюдента 11,у 12у 1з,у 14,у . ;п-1,у 1у,у
Нижняя граница Х1шт Х2шт Х3шт Х4шт • • Хп-1шт Хйшш Ушт
интервала
Верхняя граница Х1шаХ Х2шаХ Х3шаХ Х4шаХ • • Хп-1 шаХ ХпшаХ УшаХ
интервала
Рассмотренный метод оценки влияния причин возгорания на число пожаров, основанный на интервальных статистических оценках и корреляционном анализе может быть использован и для оценки влияния причин на число погибших, на число травмированных и на величину материального ущерба.
Литература
1. Айвазян С. А., Бухштабер В.М., Енюков И.С. Прикладная статистика: Классификация и снижение размерности. М.: Финансы и статистика, 1989.
2. Гмурман В.Е. Теория вероятностей и математическая статистика. М.: Высшая школа, 2001.
3. Куватов В.И., Величко Г.А. Исследование операций. СПб.: ВМИРЭ, 2005.