Санкт-Петербургский Государственный Университет
Кафедра вычислительных методов механики деформируемого
тела
Рединских Надежда Дмитриевна
Магистерская диссертация
Статистическое оценивание стратегии развития
фирмы в условиях коррупции
Направление 010400
Прикладная математика и информатика
Магистерская программа «Надежность и безопасность сложных систем»
Научный руководитель,
доктор физ.-мат.наук,
профессор
Малафеев О.А.
Санкт-Петербург
2016 г.
Оглавление
Введение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
Постановка задачи . . . . . . . . . . . . . . . . . . . . . . . . . . .
4
Обзор литературы . . . . . . . . . . . . . . . . . . . . . . . . . . .
5
Глава 1. Статистическое оценивание стратегии развития фирмы в
условиях коррупции . . . . . . . . . . . . . . . . . . . . . . .
7
1.1. Связь между характеристиками функционирующей фирмы 8
1.2. Распределения, используемые для вычисления вероятности функционирования фирмы . . . . . . . . . . . . .
10
1.2. Возобновление деятельности фирмы, прекратившей функционирование, вследствие негативного воздействия коррупционных факторов . . . . . . . . . . . . . . . . . .
15
Глава 2. Статистическое оценивание и прогноз эффективности стратегии развития конгломерата фирм в условиях коррупционного воздействия . . . . . . . . . . . . . . . . . . . . . . . . .
20
Глава 3. Экспериментальное определение характеристик функционирующей фирмы в условиях негативного коррупционного
воздействия . . . . . . . . . . . . . . . . . . . . . . . . . . . .
23
Глава 4. Примеры . . . . . . . . . . . . . . . . . . . . . . . . . . .
26
Выводы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
31
Заключение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
32
Список литературы . . . . . . . . . . . . . . . . . . . . . . . . . .
32
Приложение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
35
2
Введение
Исследование вопроса коррупции в настоящее время имеет особую
актуальность. Все чаще в средствах массовой информации появляется информация о проявлениях коррупции в различных сферах. За период с января по октябрь 2015 года в России возбудили 5,8 тыс. уголовных дел по
статье 290 УК РФ «Получение взятки». Материальный ущерб по оконченным уголовным делам коррупционной направленности за приведенный
период по официальным сведениям составил 39,1 миллиарда рублей [20].
Средний размер взятки по России за 2015 год составил 172,9 тыс.руб.[22].
Расхищение государственных средств лицами, причастными к коррупционной системе, и их присвоение ради собственной выгоды приводит к снижению уровня жизни населения, сокращению нормального функционирования предпринимательства, к повышению криминализации общества. Ежегодно неправительственная международная организация по борьбе с коррупцией Transparency International составляет индекс, отражающий уровень восприятия коррупции аналитиками и предпринимателями различных стран. Согласно оценке за 2015 год Россия в рейтинге стран занимает
119 место из 168 возможных.
На сегодняшний момент принимаются различные меры по снижению
уровня коррумпированности в России. 1 апреля 2016 года Президентом
Российской Федерации принят указ №147 «О Национальном плане противодействия коррупции на 2016 - 2017 годы»[23]. Данный указ предусматривает проведение мероприятий направленных на решение следующих задач:
усиление контроля за сведениями о доходах/расходах чиновников; проведение процедуры по предотвращению хищений в сфере госзакупок, муниципального и государственного имущества; обнаружение фактов подкупа
иностранных должностных лиц. Предполагается проведение социологических исследований для оценки уровня коррупции.
3
Постановка задачи
В качестве задачи была поставлена проблема формализации и анализа стратегии развития фирмы в условиях негативного коррупционного
воздействия. Необходимо было определить основные понятия и определения функционирования фирмы в условиях негативного коррупционного
воздействия и исследовать связь между введенными понятиями. Требовалось рассмотреть введенные характеристики функционирования фирмы в условиях негативного коррпуционного воздействия для различных
распределений (распределение Вейбулла, нормальное распределение, логарифмическое распределение, усеченно-нормальное распределение, Гаммараспределение, распределение Рэлея). После прекращения функционирования фирмы в условиях негативного коррупционного воздействия, она
может возобновить свое функционирование, так что процесс ее развития
будет продолжаться. Для такого типа функционирования необходимо было рассмотреть возобновление деятельности фирмы, прекратившей свое
функционирование вследствие негативного коррупционного воздействия.
Также необходимо было формализовать и исследовать процесс функционирования конгломерата (действующих связанным образом между собой)
фирм в условиях негативного коррупционного воздействия. Требовалось
ввести основные понятия и определения для случая функционирования
конгломерата фирм в условиях негативного коррупционного воздействия.
В качестве подзадачи, необходимо было проанализировать зависимость прекращения функционирования конгломерата фирм от вероятности прекращения функционирования составляющих его подразделений, для случая,
если прекращение функционирования подразделений происходит независимо. Необходимо было сформулировать экспериментальное определение
характеристик вероятности функционирования фирмы в условиях негативного коррупционного воздействия, рассмотреть численные примеры для
случаев функционирования фирмы и конгломерата фирм в условиях негативного коррупционного воздействия.
4
Обзор литературы
В работе [6] рассматриваются причины появления коррупции в современной Российской экономике. Коррупция представлена как сложный
процесс, охватывающий такие сферы экономики, как функционирование
бизнеса, государственного аппарата и жизнедеятельность домохозяйств.
Также в статье проводится анализ причин распространения коррупции в
различных странах. На основе проведенных исследований предложен ряд
мер, направленных на противодействие коррупции в России. Также рассмотрены возможности формирования общественных и государственных
институтов по борьбе с коррупцией.
Различные подходы к исследованию коррупции представлены в работе [13]. Рассмотрена классификация подходов к изучению коррупции в
экономике – построение экономико-математических моделей, проведение
эконометрических исследований, и социологии – проведение социологических опросов и экспертных оценок. Практика противодействия коррупции
в России и направления борьбы с коррупцией исследованы в статье [3].
В работе вводится понятие коррупции, рассматриваются законы о противодействии коррупции, принятые, начиная с 1992 года. Проводится сравнительный анализ имеющихся антикоррупционных законодательств в России. Предлагается методология улучшения имеющегося антикоррупционного законодательства.
Различные способы оценки уровня коррупции представлены в [17].
Рассматриваются причины, по которым необходимо исследовать коррупцию, методики составления индекса восприятия коррупции и исследуется
точность приведенных методик. На основе полученных данных ставится
вопрос о построении оптимальной антикорупционной стратегии. Рассмотрена методика проведения региональных исследований уровня коррупции.
Использование математических моделей различного рода для оценки коррупции рассмотрено в работах [18],[19]. Влияние теневой экономики
на уровень коррупции в различных странах рассматривается в работе [16].
Автор предполагает, что коррупция и теневая экономика неразрывно связаны в странах с развитой экономикой и дополняют друг друга в странах
5
с развивающейся экономикой. Сделанное предположение протестировано
на данных о 98 странах. Результаты исследования показали, что в большинстве случаев предположение подтверждается исследованием. В статье
[15] исследуется вопрос становления коррупции в современной России. Автор рассматривает разветвленную коррумпированную структуру, которая
имеется в настоящий момент в России, оценивает ее негативное влияние на
различные сферы развития экономики и общества.
Вопросы коррупции в сетевой торговле рассмотрены в работе [4].
Исследуются теоретические положения и факты появления коррупции в
сетевой торговле. Выявлены причины и факторы по которым коррупция
возникает в торговле: поставка низкокачественного товара по завышенной цене, поставка товаров торговым организациям, которые не покупают
товар без взяток. Предложены методы по противодействию коррупции в
сетевой торговле. Методология оценки масштабов коррупционной деятельности представлена в работе [11]. На основе проведенного опроса было выявлено отсутствие зависимости между средним размером взятки, который
изменяется от региона к региону, и интенсивностью получения взяток. Выявлен резкий рост коррупции в период 2001-2005 гг. Рассмотрена методология противодействия коррупции. Результаты социологических исследований уровня коррупции в России приведены в статье [12]. Представленные
исследования проведены фондом ИНДЕМ с 2001 по 2010 годы. В книге
анализируются различные проявления коррупции с помощью качественных и количественных эмпирических методов. Проведен анализ динамики
коррупции за исследуемый период. Изложена методология исследования,
практика деловой и бытовой коррупции в России.
6
Глава 1.Статистическое оценивание стратегии
развития фирмы в условиях коррупции
Формализуется и исследуется задача оценивания качества стратегии
φ развития фирмы в условиях негативного коррупционного воздействия с
помощью аппарата, методов и результатов теории надежности.
Введем основные понятия и определения используемые в работе.
Момент T прекращения функционирования фирмы при начальном
моменте времени (начале процесса функционирования фирмы) t = 0 в
условиях негативного коррупционного воздействия – случайная величина.
Вероятность функционирования фирмы p(t) – вероятность функционирования фирмы P (T ) в течение заданного времени t в условиях негативного коррупционного воздействия, при условии, что фирма начала функционировать в момент времени t0 = 0, или вероятность того, что время функционирования фирмы до момента прекращения T функционирования окажется больше заданного времени функционирования t, т.е. p(t) = P (T ≥ t).
Вероятность прекращения функционирования фирмы q(t) в условиях негативного коррупционного воздействия – вероятность того, что
фирма прекратит свое функционирование в течение заданного времени
функционирования t, начав функционировать в начальный момент времени t0 = 0, q(t) = 1 − p(t).
Плотность распределения времени функционирования фирмы в условиях негативного коррупционного воздействия f (t) – функция, определяемая по формуле f (t) =
d
dt q(t)
= − dtd p(t). Смысл плотности распределения
времени функционирования фирмы в условиях негативного коррупционного воздействия состоит в том, что она показывает, частоту появления
случайной величины в некоторой окрестности точки t при повторении наблюдений.
Интенсивность прекращения функционирования фирмы в момент
времени t – плотность распределения времени функционирования фирмы
к моменту времени t при условии, что до этого момента прекращение функционирования фирмы не произошло: λ(t) =
7
f (t)
p(t)
Среднее время прекращения функционирования фирмы в условиях
коррупционного воздействия tavg – математическое ожидание случайной
+∞
∫∞
∫
tf (t)dt = tf (t)dt, (f (t) = 0, при t >
величины T есть tavg = M (T ) =
−∞
0
0).
1.1. Связь между характеристиками функционирующей
фирмы
Введенные величины связаны между собой и по одной из них можно
найти другие. После прекращения функционирования фирмы в условиях негативного коррупционного воздействия она может возобновить свое
функционирование так, что процесс ее развития будет продолжаться (возможно, под иным названием). В этом случае, из определения функции распределения и ее плотности стандартным образом получаем:
f (t) = F ′ (t) = q ′ (t).
(1)
Все эти величины равны нулю при t < 0. Т.к. q(t) = 1 − p(t), то
f (t) = −p′ (t).
(2)
f (t) = p(t)λ(t).
(3)
Для вывода (3) воспользуемся следующим: пусть А – событие, означающее функционирование фирмы к моменту времени t, B – событие, означающее функционирование фирмы к моменту времени t+∆t. Событие AB
– прекращение функционирования фирмы в интервале (t, t+∆t), т.е. время
прекращения функционирования лежит в интервале (t, t + ∆t). Тогда по
теореме умножения P (AB) = P (A) · P (B|A), но по определению события
A и формуле p(t) = P (T ≥ t), получаем P (A) = P (T ≥ t) = p(t). Согласно
определению события B и формуле λ(t)∆t ≃ P (t ≤ T ≤ t + ∆t|T ≥ t)
следует, что P (B|A) = P (t ≤ T ≤ t + ∆t|T ≥ t) = λ(t)∆t. В соответствии с
формулой f (t)∆t ≃ P (t ≤ T ≤ t + ∆t) получаем, что P (AB) = P (t ≤ T ≤
t + ∆t) = f (t)∆t. Следовательно, f (t)∆t = p(t)λ(t)∆t. После сокращения
8
на ∆t получается выражение (3). Учитывая (1), (2) λ(t) =
f (t)
1−F (t)
=
F ′(t)
1−F (t)
=
q′(t)
1−q(t)
=
f (t)
p(t)
=
f (t)
1−q(t)
=
p′(t)
p(t) .
Функция распределения определяется через плотность распределе∫t
ния: F (t) =
f (x)dx. Вероятность прекращения функционирования: q(t) =
F (t) =
∫t
−∞
f (x)dx, следовательно p(t) = 1 − q(t) = 1 −
0
+∞
∫
По свойству плотности распределения
при t ≤ 0,
∫∞
∫∞
f (x)dx =
∫∞
f (x)dx.
t
0
f (x)dx = 1. Поскольку f (t) = 0
−∞
∫t
f (x)dx = 1. Тогда p(t) = 1− f (x)dx =
0
∫t
0
∫∞
∫t
f (x)dx− f (x)dx =
0
0
f (x)dx.
t
Можно также выразить p(t) через λ(t), решив дифференциальное
уравнение
λ(t)dt = −
dp(t)
.
p(t)
Заменяя t на x и интегрируя обе части от 0 до t (p(0) = 1), получаем
∫t
∫t
λ(x)dx = −
0
Следовательно,
dp(t)
= − ln p(t) |t0 = −l ln p(t).
p
0
∫t
− λ(x)dx
p(t) = e
0
∫t
− λ(x)dx
=e
0
.
(4)
Интегральные характеристики можно выразить через любые из функций f (t), q(t), p(t), λ(t). Если интенсивность прекращения функционирования фирмы λ(t) постоянна (что можно интерпретировать как наличие
системной коррупции, порождаемой внешними факторами) при которой
вероятность прекращения функционирования фирмы за любой короткий
промежуток времени постоянна, то в этом случае
∫∞
tavg = M (T ) =
tλe−λt dt = 1/λ,
0
σ = M (T ) = 1/λ.
9
1.2. Распределения, используемые для вычисления вероятности функционирования фирмы
Экспоненциальное распределение
Распределение вероятности функционирования фирмы в условиях
негативного коррупционного воздействия называется экспоненциальным,
если
p(t) = e−λt ; t ≥ 0.
Оно характеризуется параметром λ = const. Вероятность прекращения
функционирования фирмы в условиях негативного коррупционного воздействия определяется по формуле:
0, t ≤ 0
q(t) = F (t) =
1 − e−λt , t > 0
Плотность распределения времени функционирования фирмы в условиях негативного коррупционного воздействия: f (t) = F ′ (t) = λe−λt , интенсивность прекращения функционирования фирмы в таком случае λ(t) =
p(t)
f (t)
= λ = const. При экпоненциальном законе распределения интенсив-
ность прекращения функционирования постоянна, что означает, что вероятность прекращения функционирования фирмы в интервале ∆t не зависит от предыстории фирмы и равна λ∆t.
Распределение Вейбулла
При экспоненциальном распределении интенсивность прекращения
функционирования фирмы λ(t) остается постоянной. Однако во многих
случаях интенсивность прекращения функционирования фирмы является или убывающей или возрастающей функцией времени. В этом случае
для описания процесса функционирования фирмы удобно воспользоваться распределением Вейбулла, при котором вероятность функционирования
фирмы есть
p(t) = e−λt ,
α
(5)
где λ, α – суть положительные константы. Естественным образом интер-
10
претируются формулы
f (t) = −p′ (t) = αtα−1 λe−λt .
α
(6)
λ(t) = αλtα−1 .
При α = 1 распределение Вейбулла совпадает с экспоненциальным распределением, при α < 1 дает убывание интенсивности прекращения функционирования с временем, а при α > 1 – возрастание. Отсюда следует, что
для распределения Вейбулла λ(t) не равна постоянной λ. Среднее время
функционирования определяется по формуле:
∫∞
tavg = M (T ) =
∫∞
0
e−λt dt.
α
f (t)dt =
0
Этот интеграл может быть выражен через гамма-функцию Γ(s) =
( )1
1
1
Используя подстановку x = λtα , откуда t = λx α = λ α x α ;
dt =
∫∞
xs−1 e−x dx.
0
1 − 1 1 −1
λ α x α dx.
α
Так как λ и α > 0 (иначе распределение не имело бы смысла при t = 0 или
t = ∞), пределы интегрирования не изменяются
∫∞
λ− α
dt =
α
1
e−λt
α
∫∞
x
1
α −1
(
)
− α1
λ
1
· e−x dx =
Γ −
.
α
α
0
0
Следовательно
tavg = M (T ) =
Γ
(1)
α
1
αλ α
.
Распределение Рэлея
Если интенсивность прекращения функционирования фирмы линейно растет со временем, то распределение Вейбулла в таком случае называется распределением Рэлея. Это соответствует значению α = 2. Параметр
λ в этом случае равен 1/2σ 2 . При данных значениях получаем следующие
11
формулы для вероятности функционирования фирмы, плотности среднего времени прекращения функционирования, интенсивности прекращения
функционирования:
p(t) = exp− 2σ2
t
tavg
(7)
t2
t
− 2σ
f (t) = 2 · exp 2
σ
t
λ(t) = 2
σ
(1)
√
√ √
Γ 2
π· 2
π
= ( )1 =
=
σ
1 2
2
2
2
(8)
(9)
(10)
2σ 2
Гамма-распределение
В случае, когда анализируемая фирма имеет сложную, разветвленную структуру, отдельные подразделения которой могут в результате коррупционного воздействия прекратить функционирование и тем самым обусловить прекращение функционирования всей фирмы, вероятность функционирования фирмы выражается с помощью гамма-функции более сложным образом. В таком случае
(λt)α−1
f (t) = λ
exp−λt , (t ≥ 0)
Γ(α)
(11)
Используя замену λx = s, λdx = ds, получаем следующее:
∫∞
p(t) =
α−1
(λx)
e−λx dx =
λ
Γ(α)
t
∫∞
sα−1 −s
e ds
Γ(α)
(12)
λt
При α = k интеграл может быть выражен через элементарные функции с помощью интегрирования по частям, тогда (12) имеет вид
p(t) =
k−1
∑
(λt)i
i=0
12
i!
e−λt
(13)
Среднее время функционирования фирмы определяется по формуле:
∫∞
tavg =
∫∞
tf (t)dt =
0
Γ(α + 1) α
(λt)α−1 λ −λt
e dt =
= .
Γ(α)
λΓ(α)
λ
(14)
0
Нормальное распределение
Допустим теперь, что прекращение функционирования фирмы происходит в результате целенаправленного действия другой заинтересованной фирмой (например, рейдером, или связанным с ним коррумпированным чиновником), стремящейся прекратить функционирование фирмыконкурента в процессе взаимодействия и накопления дефектов управления у рассматриваемой фирмы. В этом случае вероятность прекращения
функционирования фирмы согласно центральной предельной теореме будет иметь распределение, близкое к нормальному, с тем, однако, уточнением, что время функционирования фирмы всегда положительно. Поэтому
здесь вводится «усечение» нормального распределения с плотностью
2
1
− (t−a)
2
2σ
,
f (t) = √
e
2πσb
(t ≥ 0).
(15)
Коэффициент b введен для нормировки распределения так, чтобы
∫+∞
f (t)dt = 1.
(16)
−∞
Вероятность функционирования фирмы при усеченном нормальном распределении имеет вид
∫∞
p(t) =
t
1
f (x)dx =
b
∫∞
√
t
(x−a)1
1
e− 2σ2 dx.
2π · σ
(17)
Среднее время функционирования фирмы
a2
2σ
( a )] · e− 2σ2
tavg = M (T ) = a + √ [
2π 1 + Φ σ
13
(18)
при a ≫ 3σ tavg ≈ a.
Логарифмически-нормальное распределение
При учете фактора времени в процессе функционирования фирмы
подверженной воздействию безнаказанной коррупции, с течением которого приобретается опыт коррупционной деятельности лучшее, чем усеченное
нормальное распределение доставляет логарифмически-нормальное распределение, при котором нормально распределен логарифм времени функционирования фирмы. В этом случае
t−a)2
1
− (ln2σ
2
f (t) = √
,
e
σ 2π · t
(t ≥ 0),
(19)
Интенсивность прекращения функционирования фирмы λ(t) равна
отношению f (t)/p(t). Заметим, что величины a и σ не являются математическим ожиданием и среднеквадратическим отклонением для логарифмически нормального распределения. Для среднего времени функционирования фирмы tavg , используя подстановку s =
∫∞
tavg = M (T ) =
∫∞
tf (t)dt =
0
0
1
·e
t√
2π · σ · t
ln t−a
σ ,
t−1)2
− (ln2σ
2
получаем
1
=√
2π
∫+∞
s2
ea+σs · e− 2 ds
−∞
(20)
(так как при t = 0, s = −∞, а для t = +∞, s = ∞). Вводя подстановку
s − σ = u, ds = du получаем
tavg
1 a+ σ2
=√ e 2
2π
∫+∞ 2
2
− u2
a+ σ2
e du = e
,
(21)
−∞
где a – есть математическое ожидание ln T , но ea не есть математическое
ожидание T .
14
1.2. Возобновление деятельности фирмы, прекратившей
функционирование, вследствие негативного воздействия
коррупционных факторов
До сих пор рассматривался случай функционирования отдельной фирмы. Однако, часто сторона, управляющая функционированием фирмы, в
случае прекращения ее функционирования открывает аналогичную фирму
под другим названием. В этом случае требуется, чтобы фирма функционировала время t. В случае прекращения функционирования фирмы она
вновь открывается под другой вывеской до тех пор, пока не наступит конечный момент T .
Рассмотрим зависимость вероятности функционирования фирмы от
числа ее закрытий. Полагая, что вероятность функционирования фирмы
при одном ее открытии есть p, вычислим вероятность функционирования
фирмы после (k −1) ее открытий. На всех циклах функционирования фирмы величину p считаем одной и той же. Полагая циклы независимыми и
обозначив событие A – фирма функционирует в течение всего цикла, а
событие Ā – противоположное событие, имеем
pk = P (ĀĀ . . . ĀA) = [P (Ā)]k−1 p(A).
(22)
Так как P (A) = p, P (A) = 1 − p = q, то
Pk = q k−1 p, k = 1, 2, . . .
(23)
Данное распределение числа закрытий и возобновлений функционирования фирмы называется геометрическим. Для вычисления математического
ожидания и среднеквадратического отклонения используем производящую
функцию
γx (t) =
∞
∑
k=1
k
pk t =
∞
∑
k=1
pq k−1 tk =
pt
.
1 − qt
(24)
Здесь использованы формулы суммы бесконечно убывающей геометриче-
15
ской прогрессии
a1 + a1 q + . . . + a1 q n + . . . =
a1
.
1−q
(25)
Знаменателем прогрессии является произведение qt, меньшее 1, так как
параметр t сколь угодно мал, а a1 = pt. В результате выкладок получаем
γx′ (t) =
p(1 − qt) + qpt
p
=
.
(1 − qt)2
(1 − qt)2
Откуда
m = M (X) = γx′ (1) =
p
p
1
=
=
.
(1 − q)2
p2
p
(26)
(27)
Вычислим вторую производную:
γx′′ (t) =
2pq
.
(1 − qt)3
Поэтому
M (X · (X − 1)) = γx′′ (1) =
2pq
2q
=
.
(1 − q)3
p2
(28)
(29)
Отсюда
M (X 2 ) =
2q
2q
2q 1 p + 2q
1+q
=
+
M
(X)
=
+
=
=
,
p2
p2
p2 p
p2
p2
1+q
1
q
−
=
,
p2
p2
p2
√
√
q
σ = D(X) =
.
p
D(X) = M (X 2 ) − [M (X)]2 =
(30)
(31)
(32)
Рассмотрим теперь простую задачу о возобновлении деятельности
фирмы, прекратившей функционирование, вследствие негативного воздействия коррупционных факторов. Пусть для фирмы имеет место экспоненциальное распределение вероятности функционирования с интенсивностью
прекращения функционирования λ. В случае прекращения функционирования фирма возобновляет свое функционирование, причем время возобновления функционирования также случайная величина с экспоненциальным распределением с параметром µ. Пусть p(t) – вероятность того, что
16
фирма в момент t функционирует. Эту величину можно назвать функцией жизнеспособности или нестационарным коэффициентом жизнеспособности. Для определения p(t) составим дифференциальное уравнение. Рассмотрим изменение состояния фирмы за интервал ∆t. Пусть событие A
– функционирование фирмы в момент t + ∆t, вероятность этого события
P (A) = p(t + ∆t). Рассмотрим гипотезы H1 и H2 , составляющие полную
группу. Гипотеза H1 – функционирование фирмы в момент t, P (H1 ) = p(t).
Гипотеза H2 – в момент t функционирование фирмы прекратилось. Тогда P (H2 ) = 1 − P (H1 ) = 1 − p(t). Найдем теперь условную вероятность
P (A|H1 ) того, что фирма, функционирующая в момент t, остается функционирующей и в момент t + ∆t. Для этого необходимо, чтобы фирма функционировала на интервале [t, ∆t]. Но вероятность прекращения функционирования фирмы на этом интервале, если она функционировала в момент
t, есть λ · ∆t, где λ – интенсивность воздействия негативных факторов. Тогда имеем, что P (A|H1 ) = 1 − λ∆t. При этом мы полагаем, что за малый
промежуток времени не может быть более одного прекращения функционирования фирмы, а также невозможно возобновление ее деятельности.
Условная вероятность P (A|H2 ) – есть вероятность того, что фирма,
возобновляемая в момент t, возобновит свое функционирование за время
∆t. Для распределения времени возобновления фирмы после ее прекращения функционирования параметр µ имеет то же значение и смысл, что и
параметр λ для распределения времени функционирования фирмы, поэтому P (A|H2 ) = µ∆t. По формуле полной вероятности
P (A) = P (H1 ) · P (A|H1 ) + P (H2 ) · P (A|H2 ).
Следовательно
p(t + ∆t) = p(t)(1 − λ∆t) + [1 − p(t)]µ∆t.
Перенесем p(t) влево и разделим на ∆t, получаем
p(t + ∆t) − p(t) = −λp(t)∆t + µ[1 − p(t)] · ∆t,
17
(33)
p(t + ∆t) − p(t)
= −λp(t)∆t + µ[1 − p(t)].
∆t
Устремляя t к нулю, получаем дифференциальное уравнения для p(t):
dp(t)
= µ − (λ + µ)p(t).
dt
(34)
Уравнение (34) – неоднородное линейное дифференциальное уравнение первого порядка с постоянными коэффициентами. Его решение – сумма общего решения однородного и частного решения неоднородного уравнения. Характеристическое уравнение – r+(λ+µ) = 0, корень которого r1 = −(λ+µ).
Поэтому общее решение однородного уравнения
p(t) = Ce−(λ+µ)t .
Так как µ =const, то частное решение неоднородного уравнения ищем в
виде p1 (t) = A =const и подставляем в (34). Тогда
0 = µ − (λ + µ)A,
A=
µ
.
λ+µ
Общее решение уравнения (34) имеет вид
p(t) =
µ
+ Ce−(λ+µ)t .
λ+µ
Для определения постоянной C используем налагаемое условие. Если в
момент t = 0 фирма функционирует, то p(0) = 1. Тогда
p(0) =
µ
+ Ce0 = 1.
λ+µ
Отсюда
C =1−
так что
p(t) =
µ
λ
=
,
λ+µ λ+µ
µ
λ
+
· e−(λ+µ)t .
λ+µ λ+µ
18
(35)
Это есть вероятность нахождения фирмы в состоянии функционирования
в момент t, если в начальный момент она функционировала. Считаем, что
если фирма прекратила функционирование, то она заменяется аналогичной, при этом на возобновление требуется случайное время с экспоненциальным распределением с параметром µ. На практике часто представляет
интерес величина p(t) при t → ∞. Такой предел можно назвать коэффициентом эффективной стабильности фирмы и обозначить ke = µ/(λ + µ).
При длительном времени функционирования фирмы коэффициент ke получает статистическую интерпретацию. Допустим, t1 – общее время нахождения фирмы в состоянии функционирования, а t2 – общее время, в
течение которого фирма находится в состоянии обновления и обновления
функционирования. Отсюда длина всех интервалов t1 + t2 = t̃. Тогда
p(t) ≈
t1
,
t1 + t2
e
= tlim
→∞
1
t2 →∞
t1
.
t1 + t2
(36)
Для частного случая экспоненциальных распределений для времени функционирования, состоянии обновления и обновления интерпретация коэффициентов λ, µ очевидна. Пусть T1 – случайное время функционирования
фирмы, а T2 – случайное время ее нахождения в состоянии обновления.
Тогда по свойствам показательного распределения
tavg = M (T1 ) =
1
,
λ
λ=
1
(37)
tavg
для среднего времени нахождения в состоянии функционирования, и для
среднего времени возобновления функционирования
tavgrec = M (T2 ) =
1
,
µ
µ=
1
tavgrec
.
(38)
Подставляя выражения (37), (38) в (34), получаем
ke =
tavg
.
tavg + tavgrec
(39)
Поэтому коэффициент эффективной стабильности фирмы равен отношению среднего времени пребывания фирмы в состоянии функционирования
19
к среднему промежутку времени между двумя соседними выходами фирмы в состояние функционирования (равному сумме среднего промежутка времени нахождения фирмы в состоянии функционирования и времени
прекращения функционирования). Величина tavg +tavgrec есть средний промежуток времени между последовательными прекращениями функционирования фирмы. Обратная величина доставляет среднее число закрытий
фирмы в единицу времени, она называется плотностью возобновления и
обозначается h(t)
h=
1
λµ
=
.
tavg + tavgrec
λ+µ
(40)
По величине h находится среднее число прекращений функционирования
фирмы за промежуток t. Эта величина называется функцией возобновления и обозначается H(t). При больших t
H = ht = t ·
t
.
tavg + tavgrec
(41)
Глава 2. Статистическое оценивание и прогноз
эффективности стратегии развития конгломерата фирм в условиях негативного коррупционного воздействия
Введенные ранее понятия и определения для стохастического оценивания и прогноза эффективности стратегии развития фирмы в условиях
негативного коррупционного воздействия равным образом относятся и к
процессу стохастического оценивания и прогноза эффективности стратегии
развития конгломерата фирм, действующих связанным между собой образом в условиях негативного коррупционного воздействия. Для них также
вводятся следующие понятия. Вероятность функционирования конгломерата фирм в условиях негативного коррупционного воздействия характеризует отношение числа функционирующих фирм N (t) в момент времени t
к числу фирм N (t0 ) функционирующих момент времени t0 = 0 в условиях
20
негативного коррупционного воздействия.
P (t) = N (t)/N (t0 )
(42)
Вероятность прекращения функционирования характеризующая отношение числа фирм N (t + ∆t) функционирующих до момента времени
t + ∆t к числу фирм N (t), функционирующих в момент времени t в условиях негативного коррупционного воздействия
Q(t) = N (t + ∆t)/N (t),
(43)
Плотность вероятности функционирования конгломерата фирм в условиях негативного коррупционного воздействия характеризующая отношение
числа фирм N (t + ∆t) − N (t), прекративших свое функционирование в
интервале t + ∆t к произведению числа функционирующих фирм N (0) в
начальный момент времени t0 = 0 на длительность интервала времени ∆t.
f (t) =
N (t + ∆t) − N (t)
.
N (0)∆t
(44)
Интенсивность прекращения функционирования в условиях негативного
коррупционного воздействия характеризующая отношение числа прекративших функционирование фирм в интервале времени t + ∆t к произведению числа функционирующих фирм в момент времени t на длительность
интервала времени ∆t, при условии, что прекратившие свое функционирование фирмы не восстанавливаются и вместо них не открывают новые
λ(t) =
N (t + ∆t) − N (t)
.
N (t)∆t
(45)
Среднее время прекращения функционирования группы фирм в условиях
негативного коррупционного воздействия определяется как
N
1 ∑
tcp =
ti
N i=0
(46)
где N –количество фирм, функционирующих в условиях негативного кор21
рупционного воздействия, ti – время функционирования i-ой фирмы.
Вероятность прекращения функционирования конгломерата фирм зависит сложным образом от вероятности прекращения функционирования
составляющих ее подразделений. Считая вероятность продолжения деятельности конгломерата фирм сложным событием, а вероятность продолжения деятельности составляющих подразделений элементарным событием, получаем, что каждому набору состояний подразделений соответствует вполне определенное состояние функционирования конгломерата фирм.
Для простоты можно считать вначале, что как подразделение, так и конгломерат фирм могут находиться в двух состояниях – либо функционирования, либо нефункционирования. Тогда с каждым состоянием свяжем
соответствующую булеву переменную. Булева переменная конгломерата
фирм будет функцией булевых переменных составляющих ее подразделений. Выписывая эту функцию с помощью логических сложения и умножения, можно связать вероятность функционирования конгломерата фирм с
вероятностью функционирования отдельных подразделений. Считаем для
простоты, что прекращение функционирования отдельного подразделения
происходит независимо, так что условные вероятности не будут необходимы. Рассмотрим для примера два случая.
1. Прекращение функционирования двух подразделений означает прекращение функционирования конгломерата фирм. Обозначим A событие, состоящее в функционировании одного, а B — событие, состоящее в функционировании другого подразделения. Тогда по теореме об
умножении независимых событий получаем событие AB – (конъюнкцию)
P (AB) = P (A) · P (B).
(47)
2. Пусть теперь функционирование конгломерата фирм из двух подразделений требует функционирования хотя бы одного подразделения из
двух. Тогда функционирование конгломерата фирм соответствует событию A + B и по теореме сложения для совместных событий
P (A + B) = P (A) + P (B) − P (AB).
22
(48)
Если A и B – независимые события, то
P (A + B) = P (A) + P (B) − P (A) · P (B).
(49)
Если группа состоит более чем из двух подразделений, то следует использовать формулы несколько раз.
Глава 3. Экспериментальное определение характеристик функционирующей фирмы в условиях негативного коррупционного воздействия
На практике величины определенные для вычисления характеристик
функционирования фирмы в условиях негативного коррупционного воздействия нужно определять опытным путем. Для этого требуется использовать статистическое истолкование вероятности. Пусть имеется N0 одинаковых фирм функционирующих в условиях коррупционного воздействия.
Если одновременно начать их исследовать на наличие коррупционной составляющей, то к моменту времени t будет N (t) фирм, функционирующих
в условиях коррупционного воздействия и n(t) фирм прекративших свое
функционирование. Очевидно,
N (t) + n(t) = N0 .
(50)
Отношение количества фирм, функционирующих в условиях коррупционного воздействия к первоначальному числу будет
p∗ (t) =
N (t)
N0
(51)
и будет эмпирической вероятностью функционирования фирм в условиях
негативного коррупционного воздействия. В соответствии с теоремой Бернулли, при достаточном числе N0 она будет приближаться к величине p(t),
23
определенной выражением p(t) = P (T ≥ t). Аналогично
n(t)
,
N0
q ∗ (t) =
(52)
будет приближением к q(t), определенной формулой q(t) = 1 − p(t) =
P (T < t). Если построить график величины q ∗ (t), то получится ступенчатая линия (величина q ∗ (t) изменяется на
1
N
когда прекращает функциони-
рование очередная фирма). Полученный график эмпирической функции
распределения представляет собой эмпирическое приближение к графику
теоретической функции распределения q(t) = F (t). С помощью критериев математической статистики можно установить, имеется ли отличие от
теоретической кривой в пределах допустимого.
Для получения эмпирической оценки плотности вероятности прекращения функционирования фирмы в условиях коррупционного воздействия
и интенсивности прекращения функционирования необходимо ввести величину ∆n – число фирм, отказавших в интервале (t; t+∆t), или для лучшего
приближения выбрать интервал (t −
∆t
2 ;t
+
∆t
2 ).
Тогда согласно формуле
f (t)∆t ≃ P (t ≤ T ≤ t + ∆t) и формулы λ(t)∆t ≃ P (t ≤ T ≤ t + ∆t/T ≥ t)
величины
f ∗ (t) =
∆n
,
N0 ∆t
(53)
λ∗ (t) =
∆n
N (t)∆t
(54)
будут соответственно приближениями к величинам f (t) и λ(t). Действительно, величина
∆n
N0
представляет собой относительную частоту прекра-
щения функционирования в интервале ∆t поэтому будет оценкой плотности вероятности. Если нанести величины
∆n
N0
на график, то получается
гистограмма, являющаяся приближением к кривой распределения. Аналогично, величина
∆n
N (t)
представляет собой относительную частоту прекра-
щения функционирования фирмы в условиях коррупцинного воздействия
в интервале ∆t для тех элементов, которые были исправны к моменту t
(поскольку число фирм, прекративших свое функционирование в условиях коррупционного воздействия ∆n делится на число фирм, функционирующих в условиях коррупцинного воздействия N (t)), т.е. является стати24
стической оценкой условной вероятности прекращения функционирования
фирм в условиях коррупционного воздействия в интервале ∆t. Отношение этой величины к ∆t будет поэтому оценкой условной плотности, т.е.
интенсивности коррупционного воздействия.
Оценкой среднего времени функционирования фирмы в условиях коррупционного воздействия tavg , в соответствии с теоремой Чебышева, является эмпирическое среднее
N
∑
t∗avg
tk
k=1
,
(55)
N0
где tk – время функционирования k−ой фирмы в условиях коррупцион=
ного воздействия. Эта величина характеризует среднее время функционирования фирмы до прекращения ее функционирования. Практическое
определение величины t∗avg по формуле (55) связано с трудностями, так
как требует проведения опыта прекращения функционирования в условиях коррупционного воздействия для всех подразделений. Поэтому данную
величину описывают косвенным путем по другим параметрам, либо берут
приближенно. Например, проводят опыт до момента t0 , а затем считают
n
∑
t∗avg =
tk + (N0 − n)t0
k=1
N0
,
(56)
т.е. принимают, что в момент окончания эксперимента все оставшиеся фирмы, прекращают свое функционирование.
В случае фирм, возобновляющих функционирование, можно вместо
испытаний нескольких фирм использовать одну фирму. Если считать, что
свойства фирм после восстановления не меняются, то за T можно принять время между последовательными прекращениями функционирования
фирмы в условиях коррупционного воздействия. Тогда вместо выражения
(55) получим
m
∑
t∗avg
k=1
tk
,
(57)
m
где tk – полученный на опыте промежуток времени между последователь=
25
ными прекращениями функционирования, m – число фирм, прекративших
свое функционирование.
Глава 4. Примеры
Пример 1.
Z-score model (анализ Альтмана) - математическая формула, характеризующая вероятность банкротства отдельной фирмы (предприятия),
сформулированная американским экономистом Эдвардом Альтманом в 1968.
Модель представляет собой функцию от показателей, характеризующих
экономический потенциал предприятия и результаты его работы за истекший период. На данный момент имеется 5 моделей для различных типов компаний. Рассмотрим пятифакторную модель Альтмана для частных компаний, которая формализована следующим образом: Z ∗ = 0.717 ·
X1 + 0.847 · X2 + 3.107 · X3 + 0.420 · X4 + 0.998 · X5 , где X1 = Оборотный капитал/Активы, X2 = Нераспределенная прибыль/Активы, X3 =
Операционная прибыль/Активы, X4 = Собственный капитал/ Обязательства, X5 = Выручка/Активы. В зависимости от полученного значения
Z ∗ компанию можно отнести к той или иной группе: если Z ∗ > 2, 9 зона финансовой устойчивости ("зеленая"зона), если 1, 23 < Z ∗ < 2, 9
- зона неопределенности ("серая"зона), если Z ∗ < 1, 23 - зона финансового риска ("красная"зона). В данном контексте будем рассматривать
банкротство отдельной фирмы как прекращение функционирования вследствие негативного коррупционного воздействия. Для вычисления значения Z ∗ воспользуемся бухгалтерской отчетностью с поквартальным периодом компании «A» за последние 18 месяцев до наступления банкротства.
Вычислим значения X1 = −0.2452, X2 = 0.0076, X3 = 0.0108, X4 = 0.0266,
X5 = 0.8725. Используя пятифакторную модель Альтмана вычислим значения Z ∗ -функции для каждого квартала, что будет означать вероятность
банкротства в заданный период времени. Для каждого вычисленного значения Z ∗ , можно сделать вывод, что выбранная компания находится на
протяжении полугода в зоне финансового риска, т.к. вероятность функционирования достаточно низкая.
26
Таблица 1: Показатели финансовой отчетности
Название показателя
Оборотный капитал
Активы
Собственный
капитал
Обязательства
Выручка
Операционная
прибыль
Нераспределенная
прибыль
(убыток)
12мес. 2014
9мес. 2014
6мес. 2014
3мес. 2014
12мес. 2013
9мес. 2013
63657793
23631697
23414708
25160955
22355487
23522501
128862567
14713745
122544792
2757373
117414312
2809307
120635676
2892196
104996142
2772921
95897664
2489464
114148821
117313131
-19558501
119787418
86057737
307473
114605004
48815467
285476
117743477
20904467
170146
102223221
110150759
1599358
93408199
83677109
1037971
-19322469
129771
176341
120274
1029986
736380
Таблица 2: Вероятность прекращения функционирования фирмы
12мес. 2014
0,013374105
9мес. 2014
0,582294073
6мес. 2014
0,292456342
3мес. 2014
0,040398132
12мес. 2013
0,962853319
9мес. 2013
0,747994339
Пример 2.
Для вычисления финансовой устойчивости фирмы также можно воспользоваться коэффициентом соотношения заемных и собственных средств,
который показывает количество заемных средств на единицу собственного
капитала. Коэффициент обозначается Kз/c и вычисляется по формуле
Kз/c =
Заемный капитал предприятия
=
Собственный капитал
Краткосрочные обязательства + Долгосрочные обязательства
.
Собственный капитал
В зависимости от значения коэффициента рассматриваемое предприятие
=
(фирму) можно отнести к различным степеням устойчивости / неустойчивости. Оптимальное рекомендуемое значение коэффициента – меньше
1. Чем ниже значение коэффициента, тем выше финансовая устойчивость
фирмы, соответственно, при увеличении значения коэффициента говорят
о повышении финансовой неустойчивости фирмы. Рассмотрим фирму, для
которой известны значения коэффициента соотношения заемных и собственных средств за последние 24 месяца (таблица 3).
Согласно таблице 3, можно сделать вывод о том, что фирма близка
27
Таблица 3: Значения коэффициента соотношения заемных и собственных
средств за последние 24 месяца
Месяцы
Значения коэффициента
Месяцы
Значения коэффициента
1
0.5
13
1
2
0.7
14
1.2
3
0.9
15
0.7
4
1.2
16
0.8
5
0.7
17
0.6
6
0.6
18
1.1
7
1
19
1
8
1.2
20
1.2
9
1.4
21
0.9
10
1
22
0.6
11
1.1
23
0.9
12
0.9
24
0.6
к прекращению функционирования. Для определения вероятности прекращения функционирования фирмы с заданными коэффициентами необходимо проверить гипотезу, что выборка коэффициентов согласуется с нормальным законом распределения при уровне значимости α = 0.01. Пусть H0 –
различия между эмпирическим и теоретическим законом несущественны,
H1 – различия между эмпирическим и теоретическим законом существенны. Вычислим выборочное среднее xn = 0.66 и выборочную дисперсию
σв = 0.819. Вычислим теоретические частоты, учитывая, что n = 24,
h = 0.1, σв = 0.819 по формуле n′i = σnh · φ(ui ). Сравним эмпирические
в
Таблица 4
i
1
2
3
4
5
6
7
8
9
xi
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
x −x
ui = iσ в
в
4,06779661
2,847457627
2,847457627
1,627118644
1,627118644
1,627118644
1,627118644
0,406779661
0,406779661
ϕ(ui )
0,397
0,0396
0,1354
0,397
0,0387
0,0387
0,2943
0,1354
0,397
n′i
1,162738983
0,115981017
0,396561356
1,162738983
0,113345085
0,113345085
0,861949831
0,396561356
1,162738983
и теоретические частоты. Составим расчетные таблицы 4, 5 по которым
∑
найдем наблюдаемое значение χ2набл = (ni − n′i )2 /n′i . Из таблицы 5 найдем χ2набл = 16.79. По таблице (см.приложение) χ2крит (0.01, 6) = 16.8.
Так как χ2набл < χ2крит , то гипотезу о нормальном распределении
генеральной совокупности принимаем. Далее вычислим вероятность функционирования фирмы для заданного периода времени t = 24 мес.: p(t) =
28
Таблица 5
i ni
1 4
2 3
3 3
4 2
5 2
6 2
7 2
8 1
9 1
1
2 [1
n′i
1,162738983
0,115981017
0,396561356
1,162738983
0,113345085
0,113345085
0,861949831
0,396561356
1,162738983
ni − n′i
0,662738983
-0,484018983
-0,303438644
0,362738983
-0,786654915
-0,886654915
-0,238050169
-0,803438644
-0,137261017
(ni − n′i )2 (ni − n′i )2 /n′i
0,43922296 0,377748546
0,234274376 2,019937246
0,092075011 0,232183518
0,13157957 0,113163463
0,618825956 5,459662914
0,786156939 6,935959689
0,056667883 0,06574383
0,645513655 1,627777506
0,018840587 0,016203625
1
24−0.66
− Φ( t−a
σ )] = 2 [1 − Φ( 0.81 )] = 0.5.
Пример 3.
Группа фирм состоит из трех подразделений, прекращение функционирования каждого из которых означает прекращение функционирования всей группы. Интенсивность прекращения функционирования подразделений λ1 , λ2 , λ3 . Распределение в режиме функционирования экспоненциально. Требуется найти вероятность функционирования группы в
течение 2000 дней, если λ1 = 1 · 10−3 прекращ.функц./дней, λ2 = 0.2 ·
10−3 прекращ.функц./дней, λ3 = 1 · 10−3 прекращ.функц./дней. Вероятности функционирования каждого подразделения группы суть p1 = e−λ1 t =
e−2·10
3
·10−5
= e−0.02 = 0.99, p2 = e−λ2 t = e−2·10
p3 = e−λ3 t = e−2·10
3
3
·10−3
·0.2·10−3
= e−0.4 = 0.67,
= e−2 = 0.1353, p = p1 · p2 · p3 = 0.0889.
Пример 4.
Пусть группа фирм состоит из двух подразделений, прекращение
функционирования каждого из которых означает прекращение функционирования группы. Среднее время функционирования первого подразделения – 300 суток, второго – 400 суток. Вероятность нормального функционирования подразделений распределена по Рэлею. Найти среднее время функционирования группы к моменту 1800 суток. Пусть вероятность
функционирования подразделений обозначена через p1 (t) и p2 (t) соответственно. Тогда вероятность функционирования группы P (t) = p1 (t) · p2 (t).
При этом ln P (t) = ln p1 (t) + ln p2 (t). Дифференцируя, получаем λ(t) =
29
′
′
(t)
− PP (t)
= − pp11 (t)
(t) −
p′2 (t)
p2 (t)
= λ1 (t) + λ2 (t). Для распределения Релея интенсив-
+ σt2 , где σ1 , σ2 – парамет2
√
2
ры соответствующих распределений. Тогда при σ1 =
σ2 =
π t1avg ,
√
2
π t2avg . Подставляя численные значения, получаем λ(t) = λ(8000) =
(
)
π
t
2
= 9.75 · 10−2 Для вычисления среднего времени функци2 t2 +!t2avg
ность прекращения функционирования λ(t) =
t
σ12
1avg
онирования проинтегрируем в соответствии с определением вероятности
функционирования P (t). Используя выражение для распределения
( Рэлея
)
−
p1 (t) = e
t2
2
2σ1
−
,
p2 (t) = e
t2
2
2σ2
2
− t2
, получаем P (t) = p1 (t) · p2 (t) = e
1
1
2 + σ2
σ1
2
.
Следовательно, вероятность функционирования группы имеет распреде1
ление Рэлея с параметром σ, которое вычисляется√из выражений
σ2 =
√
1
+ σ12 , σ = √σ12σ2 2 . И тогда из определения tavg = π2 · σ = π2 · √σ12σ2 2 .
σ2
1
2
Иначе tavg =
√ 300·400
3002 +4002
σ1 +σ2
·t2avg
√t1avg
.
2
t1avg +t22avg
σ1 +σ2
Подставляя численные значения, имеем tavg =
= 240суток.
Пример 5.
Группа фирм состоит из десяти одинаковых подразделений, прекращение функционирования каждого из которых означает прекращение функционирования группы. Вероятность функционирования каждого подразделения имеет нормальное распределение с параметрами a = 200 суток и σ =
40 суток. Найти вероятность функционирования за 150 суток. Здесь a > 3σ.
[
( )]
Вероятность функционирования каждого подразделения p = 21 1 − Φ t−a
=
σ
[
(
)]
1
150−200
≈ 0.854. Вероятность функционирования группы P =
2 1−Φ
40
P 10 = 0.89410 = 0.306.
30
Выводы
В работе введены основные характеристики процессов функционирования фирмы в условиях негативного коррупционнго воздействия и исследована связь между введенными характеристиками. Введенные характеристики рассмотрены для различных распределений. Рассмотрен случай
возобновления деятельности фирмы, прекратившей свое функционирование вследствие негативного коррупционного воздействия. Введенные понятия и характеристики для функционирования отдельной фирмы обобщены
для процесса функционирования конгломерата (действующих связанным
образом между собой) фирм в условиях негативного коррупционного воздействия. Исследована зависимость вероятности прекращения функционирования конгломерата фирм от вероятности прекращения функционирования составляющих его подразделений. Рассмотрено прекращение функционирования конгломерата фирм состоящего из двух подразделений, для
случая, если прекращение функционирования подразделений происходит
независимо. Формализовано условие для функционирования конгломерата
фирм, состоящего больше чем из двух фирм. Формализовано экспериментальное определение характеристик вероятности функционирования фирмы в условиях негативного коррупционного воздействия. Для различных
случаев функционирования фирмы и конгломерата фирм в условиях негативного коррупционного воздействия рассмотрены численные примеры.
31
Заключение
Коррупция в России по прежнему остается одной из самых значимых проблем страны. Потери от коррупции не ограничиваются размерами
прибыли от взяток, выплачиваемых ежегодно для подкупа должностных
лиц в различных отраслях экономики. Негативные последствия коррупции
препятствуют прогрессивному развитию общества и экономики, а также
представляют угрозу национальной безопасности страны.
В данной работе исследовалась проблема внутренней коррупции фирмы. К внутренней коррупционной деятельности фирмы можно отнести
следующее: мошеннические операции с активами компании, вывод прибыли в специальные организации, мошенничество с финансовой отчетностью компании. Внутренняя коррупция достигает максимальной концентрации в компаниях с государственным участием. Такие компании, как
правило, имеют сложную организацию структуры бизнеса, специфические
условия заключаемых договоров, предполагающие отсрочку платежей или
осуществление убыточных или сверхприбыльных операций, иными словами, работа таких компаний направлена на расхищение внутренних активов.
В формализованной в работе модели исследуется статистическое оценивание стратегии развития фирмы в условиях негативного воздействия
коррупционных факторов. Также в работе формализуется и исследуется статистическое оценивание стратегии развития конгломерата фирм в
условиях негативного коррупционного воздействия. Полученные в работе
результаты могут иметь прикладную направленность – они могут быть
использованы для построения прогноза динамики прекращения функционирования фирмы и конгломерата фирм в условиях негативного коррупционного воздействия. Кроме того, результаты данной работы могут быть
использованы в качестве рекомендаций государственным органам управления различного уровня.
32
Список литературы
[1] Барлоу Р., Прошан Ф. Статистическая теория надежности и испытания
на безотказность. М.: Наука, 1984, 488 C.
[2] Гмурман В.Е., Теория веротяности и математическая статистика:
Учеб.пособие для вузов, М.:Высш.шк., 2003, 479 С.
[3] Долгова А.И. Реагирование на коррупцию: практика противодействия
России и направления оптимизации борьбы// В сборнике: Коррупция:
состояние противодейтсвия и направления оптимизации борьбы Российская криминологическая ассоциация, 2015, C.7-24.
[4] Жигун Л.А., Е.М.Фиоктистова, Дискус о противодействии коррупции
в сетевой торговле, Актуальные проблемы экономики и права, № 4 (36)
/ 2015, С.50-59.
[5] Зубова А.Ф., Малафеев О.А., Моделирование социально-экномических
систем, СПб.: СПбГУ, 2006, – 907 С.
[6] Клейнер В.Г. Коррупция в России, Россия в коррупции: есть ли выход?
// Вопросы экономики, 2014, № 6, С. 81-96.
[7] Кокс Д., Смит В. Теория восстановления. М.:Советское радио, 1967, 299
С.
[8] Кузовков Ю.В. Мировая история коррупции. — М.: Анима-Пресс, 2010,
137 C.
[9] Кун В.В. Математические основы теории надежности. Л., 1977, 48 C.
33
[10] Петросян Л. А. Теория игр: учебник / Л. А. Петросян, Н. А. Зенкевич,
Е. В. Шевкопляс. — 2-е изд., перераб. и доп. — СПб.: БХВ-Петербург,
2012 — 432 C.
[11] Г. Сатаров, Как измерять и контролировать коррупцию. //Вопросы
экономики, № 1, 2007, С. 5-11.
[12] Г.Сатаров, Российская коррупция: уровень, структура, динамика.
Опыты социоло- гического анализа / Под ред. Г.А. Сатарова. Москва :
Фонд "Либеральная Миссия 2013, 752 С.
[13] Шуманов И.В. Анализ подходов к исследованию феномена коррупции
и современные исследовательские практики коррупции, Вестник Калининградского филиала Санкт-Петербургского университета МВД России, 2011, № 4, С. 26-31.
[14] Airson, Handbook of reliability edited by Airson, t. 1,2,3, Moscow, 1969.
[15] Cheloukhine S., King J. Corruption networks as a sphere of investment
activities in modern Russia // Communist and Post-Communist Studies,
Vol. 40, Iss. 1, March 2007, P. 107-122.
[16] Dreher A., Schneider S. Corruption and the shadow economy: an empirical
analysis // Public Choice, Vol. 144, No. 1/2 (July 2010), P. 215-238
[17] Galtung F., Shacklock A., Connors C., Sampford C., Measuring corruption,
Ashgate Publishing, Ltd., P. 318.
[18] Pak Hung Mo, Corruption and Economic Growth // Journal of
Comparative Economics, Vol. 29, Iss. 1, March 2001, P. 66–79
[19] Okadaa K.,Samreth S., The effect of foreign aid on corruption: A quantile
regression approach // Economics Letters, Vol. 115, Iss. 2, May 2012, P.
240–243
[20] В
взятки
МВД
в
подсчитали
России
34
средний
размер
[Электронный
ресурс].
http://www.rbc.ru/economics/09/12/2015/5667f4069a79473b1ae0b836
(Дата обращения 23.03.2016)
[21] Коррупция
ся
со
без
границ:
как
за
рубежом
взятками
Россия
будет
бороть-
[Электронный
ресурс].
http://www.rbc.ru/economics/09/10/2015/55db18319a7947ed6c110d4a
(Дата обращения 23.03.2016)
[22] Средняя сумма взятки в Москве увеличилась вдвое [Электронный ресурс]. http://www.rbc.ru/politics/28/01/2016/56aa10349a7947e20092f544
(Дата обращения 23.03.2016)
[23] Указ Президента РФ от 1 апреля 2016 г. N 147 "О Национальном
плане противодействия коррупции на 2016 - 2017 годы"[Электронный
ресурс]. http://www.garant.ru/hotlaw/federal/706953/ (Дата обращения
5.04.2016)
Приложение
Теорема Чебышева
Если X1 , X2 , . . . , Xn , . . . – попарно независимые случайные величины,
причем дисперсии их равномерно ограничены (не превышают постоянного числа C), то, как бы мало ни было положительное число ε, вероятность
неравенства | X1 +X2n+···+Xn − M (X1 )+M (Xn2 )+···+M (Xn ) | < ε будет как угодно близка к единице, если число случайных величин достаточно велико.
Центральная предельная теорема
Пусть X1 , . . . , Xn , . . . есть бесконечная последовательность независимых одинаково распределенных случайных величин, имеющих конечное
математическое ожидание и дисперсию. Обозначим математическое ожиn
∑
дание и дисперсию через µ и σ 2 , соответственно. Пусть также Sn =
Xi .
Тогда
Sn √
−µn
σ n
i=1
→ N (0, 1) при n → ∞, где N (0, 1) – нормальное распределе-
ние с ну левым математическим ожиданием и стандартным отклонением,
равным единице. Обозначив символом X n – выборочное среднее первых
n
∑
1
n величин, т.е. X n = n
Xi результат центральной предельной теореi=1
35
мы можно переписать в виде:
√ X n −µ
n σ → N (0, 1) по распределению при
n → ∞.
Статистическая гипотеза. Нулевая и конкурирующая гипотезы.
Часто необходимо знать закон распределения генеральной совокупности. Если закон распределения неизвестен, то имеются основания предположить, что он имеет определенный вид. Выдвигают гипотезу: генеральная
совокупность распределена по закону «А». Таким образом, в этой гипотезе
идет речь о виде предполагаемого распределения. Наравне с выдвинутой
гипотезой рассматривают и противоречащую ей гипотезу. Если выдвинутая гипотеза будет отвергнута, то имеет место противоречащая гипотеза.
Нулевой гипотезой называют выдвинутую гипотезу H0 . Конкурирующей гипотезой называют гипотезу H1 , противоречащую нулевой.
Критическая область. Область принятия гипотезы. Критические точки.
После выбора критерия множество всех его возможных значений разбивают на два непересекающихся подмножества: одно из них содержит значения критерия, при которых гипотеза отвергается, а другая – при которых
гипотеза принимается. Критической областью называют совокупность значений критерия, при которых нулевую гипотезу отвергают. Областью принятия гипотезы называют совокупность значений критерия, при которых
гипотезу принимают. Основной принцип проверки статистических гипотез:
если наблюдаемое значение критерия принадлежит критической области,
гипотезу отвергают, если наблюдаемое значение принадлежит области принятия гипотезы, гипотезу принимают. Критические точки – точки, отделяющие критическую область от области принятия гипотезы.
Правосторонней называют критическую область, определяемую неравенством K > kкр , где kкр – положительное число. Левосторонней называют критическую область, определяемую неравенством K < kкр , где kкр
– отрицательное число. Односторонней называют правостороннюю или левостороннюю критическую область.
Двусторонней называют критическую область, определяемую неравенствами K < k1 , K > k2 , где k2 > k1
36
Эмпирические и теоретические частоты.
Дискретное распределение.
Рассматривается дискретная случайная величина X, закон распределения которой неизвестен. Пусть произведено n испытаний, в результате
которых величина X приняла n1 раз значение x1 , n2 раз значение x2 , nk
∑
раз значение xk , причем
ni = n. Эмпирическими частотами называют
фактически наблюдаемые частоты ni . Пусть имеются основания предположить, что величина X распределена по некоторому определенному закону.
Чтобы проверить, согласуется ли предположение с наблюдением, вычисляют теоретические частоты наблюдаемых значений. Теоретическими частотами называют частоты n′i , найденные теоретически (вычислением). Теоретические частоты находят с помощью равенства n′i = nPi , где n – число
испытаний, Pi – вероятность наблюдаемого значения xi , вычисленная при
допущении, что X имеет предполагаемое распределение. Таким образом,
теоретическая частота наблюдаемого значения xi дискретного распределения равна произведению числа испытаний на вероятность этого наблюдаемого значения.
Непрерывное распределение.
В случае непрерывного распределения, вероятности отдельных возможных значений равны нулю (согласно следствию: вероятность того,
что непрерывная случайная величина X примет одно определенное значение равна нулю см. [2]). В таком случае весь интервал возможных значений делят на k непересекающихся интервалов и вычисляют вероятности Pi попадания X в i-й частичный интервал, а затем, как и для дискретного распределения, умножают число испытаний на вероятности. Теоретические частоты непрерывного распределения находят по равенству
n′i = nPi , где n – число испытаний, Pi – вероятность попадания X в i-й
частичный интервал, вычисленная при допущении, что X имеет предполагаемое распределение. Если предполагается, что случайная величина X
распределена нормально, то теоретические частоты могут быть найдены
по формуле n′i =
nh
σ
ϕ(ui ), где n – число испытаний (объем выборки), h
В
– длина частичного интервала, σВ – выборочное среднее квадратическое
отклонение. ui = (xi − xВ )/σВ (xi – середина i-го частичного интервала),
37
ϕ(u) =
2
√1 e−u /2 .
2π
Проверка статистических гипотез.
Если закон распределения генеральной совокупности не известен, но
есть предположение, что генеральная совокупность распределена по закону «A», в этом случае проверка гипотезы о предполагаемом законе распределения производится при помощи специально подобранной величины
– критерия согласия.
Критерий согласия – критерий проверки гипотезы о предполагаемом
законе распределения генеральной совокупности.
Воспользуемся критерием согласия Пирсона для проверки гипотезы
о нормальном распределении совокупности (критерий может быть использован для проверки других распределений). Критерий Пирсона устанавливает на принятом уровне значимости согласие гипотезы с данными наблюдений.
Пусть по выборке объема n получено эмпирическое распределение:
варианты
xi x1 x2 . . .
эмпирические частоты ni n1 n2 . . .
xs
ns
Пусть для нормального распределения генеральной совокупности вычислены теоретические частоты n′i . При уровне значимости α требуется
проверить нулевую гипотезу: генеральная совокупность распределена нормально. В качестве критерия проверки нулевой гипотезы примем случайную величину
χ2 =
∑
(ni − n′i )2 /n′i .
(58)
Величина (58) – случайная, так как в различных опытах она принимает
различные, заранее не известные значения. Чем меньше различаются эмпирические и теоретические частоты, тем меньше величина критерия, характеризующего близость эмпирического и теоретического распределений.
Доказано, что при n → ∞ закон распределения случайной величины (58)
независимо от того, какому закону распределения подчинена генеральная
совокупность, стремится к закону распределения χ2 с k степенями свободы. Поэтому случайная величина (58) обозначена через χ2 , а сам крите38
рий называют критерием согласия «хи квадрат». Число степеней свободы
находят по равенству k = s − 1 − r, где s – число групп (частичных интервалов) выборки; где s – число групп (частичных интервалов) выборки;
r – число параметров предполагаемого распределения, которые оценены
по данным выборки. Поскольку односторонний критерий более «жестко»
отвергает нулевую гипотезу, чем двусторонний, построим правостороннюю
критическую область, исходя из требования, чтобы вероятность попадания
критерия в эту область в предположении справедливости нулевой гипотезы была равна принятому уровню значимости α: P [χ2 > χ2кр (α; k)] = α.
Таким образом правосторонняя критическая область определяется неравенством χ2 > χ2кр (α; k), а область принятия нулевой гипотезы – неравенством χ2 < χ2кр (α; k). Пусть значение критерия, вычисленное по данным
наблюдений обозначено через χ2набл . Для проверки нулевой гипотезы при
заданном уровне значимости, необходимо, для начала, вычислить теоретические частоты, а затем наблюдаемое значение критерия:
χ2набл =
∑
(ni − n′i )2 /ni .
(59)
По таблице критических точек распределения χ2 по заданному уровню значимости α (уровень значимости – достаточно малая вероятность, при которой событие можно считать практически невозможным.) и числу степеней свободы k = s − 3 найти критическую точку χ2кр (α, k). Если χ2набл <
χ2кр – нет оснований отвергнуть нулевую гипотезу, если χ2набл > χ2кр –
нулевая гипотеза отвергается.
39
Отзывы:
Авторизуйтесь, чтобы оставить отзыв