УДК 338.27
Г орбатков Станислав Анатольевич
Уфимский филиал ФГОБУ ВПО «Финансового университета при Правительстве Российской Федерации»
Россия, Уфа1 Профессор, доктор технических наук E_Mail: sgorbatkov@mail.ru
Белолипцев Илья Игоревич
Уфимский филиал ФГОБУ ВПО «Финансового университета при Правительстве Российской Федерации»
Россия, Уфа Преподаватель E-Mail: red7315@gmail.com
Фархиева Светлана Анатольевна
Уфимский филиал ФГОБУ ВПО «Финансового университета при Правительстве Российской Федерации»
Россия, Уфа Доцент, кандидат технических наук E-Mail: ok-xi@yandex.ru
Логистический нейросетевой метод построения динамических моделей банкротств при неполных данных
Аннотация: До настоящего времени не были исследованы нейросетевые
динамические модели оценки вероятности риска банкротств, которые представляют большой практический интерес. Знание динамики развивающего процесса банкротства позволяет в каждый момент времени определить стадию банкротства и применить упреждающие управляющие воздействия. Однако на пути создания динамических моделей риска банкротства возникают труднопреодолимые препятствия в виде неполноты данных. В данной статье предложен оригинальный нейросетевой итерационный гибридный метод оценки вероятности развивающегося во времени риска банкротства предприятий. Метод ориентирован на применение в сложных условиях моделирования: неполноты данных и сильного их зашумления. Вычислительные эксперименты подтвердили эффективность основной идеи разработки гибридного нейросетевого метода построения динамической модели банкротств с непрерывным временем. Разработанный метод может быть использован в широком спектре прикладных задач: оценки кредитоспособности заемщика, обеспечения экономической безопасности предприятий путем постоянного мониторинга их финансового состояния, оценка инвесторами объектов инвестиций, оценки финансовоэкономической состоятельности контрагентов в сделках и др. Предложенные инструментарии выходят за рамки рассмотренных выше задач и могут быть использованы в других динамических задачах с неполными данными при моделировании экономических объектов.
Ключевые слова: Диагностика банкротств; неполные данные; нейросеть;
динамическая модель; риск; логистический метод; гибридный метод; оценка.
Идентификационный номер статьи в журнале 136EVN214
1 450015, г. Уфа, ул. Мустая Карима 69/1, Уфимский филиал Финансового университета при Правительстве Российской Федерации
Введение
Математические модели оценки вероятности риска банкротств в настоящее время актуальны в связи со многими приложениями в реальном и финансовом секторах экономики. К ним относятся задачи оценки банкротств корпораций; оценки финансово-экономической состоятельности контрагентов в сделках; оценки инвестиционной привлекательности предприятий; обеспечения экономической безопасности предприятий путем постоянного мониторинга их финансового состояния; оценки банком кредитоспособности заемщиков; при выработке управленческих решений по налоговому регулированию (управлению) на местном уровне и во многих других экономических приложениях.
Методы моделирования оценки риска банкротств развивались начиная с пионерских работ Альтмана (1968г.) [1] и его последователей, где использовались различные версии линейного метода множественного дискриминантного анализа (MDA). Эти методы оказались практически неприменимыми для российских предприятий по трем причинам:
1. различия в законодательстве (инвестиции на Западе более защищены законодательством, чем в России);
2. различия в стандартах бухгалтерской отчетности;
3. отсутствие глубокой адаптации моделей банкротств для различных стран и отраслей экономики.
В связи с указанными причинами в России были предприняты многочисленные попытки совершенствования MDA - линейных моделей, а также появившихся в 1980 г. логистических моделей банкротств [2]. Бурное развитие в России в последнее 20-летие интеллектуальных информационных технологий, в частности нейросетевых [3], открыло принципиально новые возможности в аспектах:
• Учета нелинейной связи моделируемого показателя У с вектором экономических показателей (факторов) анализируемого объекта и окружающей среды
У = ф\ х = (х1,х2,..., х^,..., хп), (1)
где р(-) - нелинейная функция.
• Создания моделей банкротств в условиях сильного зашумления данных, отягченного их дефицитом.
Однако до настоящего времени не исследованы нейросетевые динамические модели оценки вероятности риска банкротств, которые представляют наибольший практический интерес. Знание динамики развивающего процесса банкротства позволяет в каждый момент времени определить стадию банкротства и применить упреждающие управляющие воздействия.
На пути создания динамических моделей риска банкротства предприятий возникают труднопреодолимые препятствия в виде неполноты данных. Суть проблемы заключается в следующем. Пусть имеются данные наблюдения вектора х в нескольких временных срезах {^}, к = 1,2,..., N, т.е. известны значения всех факторов во временных срезах {х]к},
7 = 1, п; к = 1, N. Пусть обучающая выборка Б содержит наблюдения для нескольких
предприятий (у8к, х&к g = 1, G, где < , > - знак кортежа; g - номер предприятия. При этом в
последнем временном срезе (к=Щ, как правило, известны векторы значения факторов {х ^ }, а
также значения эндогенной (выходной) переменной {у ^}. Это позволяет оценить вероятность
риска банкротства в последнем временном срезе по формуле логистической трансформации вероятности риска банкротства:
Р& = 1/1 + еХР(-Yg (XgN ))] . (2)
При этом показатель экспоненты У (х^) в нашем исследовании является нелинейной функцией от вектора факторов х м и аппроксимируется нейросетевой моделью [4].
Однако, в других временных срезах {їк}, к Ф N известны только векторы факторов
{хяі;}, а значения эндогенной переменной {у^} и, соответственно, вероятности р }
неизвестны, поскольку часть предприятий имеет промежуточную (неизвестную) стадию банкротства, которую по статистической модели вида (2), связывающей х и Р, определить не представляется возможным. Причиной этого является деформация связи вида (2) с течением времени во временных срезах {їк }.
В итоге столбцы (либо строки) данных во временных срезах оказываются неполными и обучение нейросетевой динамической модели оценки риска банкротства невозможно.
В данной работе предлагается гибридный нейросетевой итерационный метод, суть которого основана на следующее идее. Будем считать, что нелинейная функциональная связь между У и вектором х в последнем временном срезе ї = N приближенно сохраняется и в других временных срезах {^}, к Ф N. При таком допущении эндогенная переменная
у к, к = 1,2,..., N -1 может быть восстановлена поочередно в каждом временном срезе подстановкой вектора х в нейросетевую модель вида (4) (см. ниже). После этого данные становятся полными (комплектными), несущими в себе динамическую информационную через функции х8 (їк) и у8(їк), где временной аргумент їк, к = 1,2,..., N -1 является дискретным.
Естественно, допущение о сохранении нелинейной зависимости У (х^) во всех
временных срезах вызывает определенную погрешность в оценке динамики риска банкротства, поэтому для уменьшения этой погрешности вводится специальный
итерационный процесс исправления функции У (х (ї)), в итоге которой мы получим
динамическую модель оценки риска банкротства с непрерывным временем ї є [о, т], что и
является целью исследования в работе.
1. Постановка задачи
Пусть известны неполные данные Б' ^ Б, где термин «неполноты» понимается в смысле отсутствия наблюдаемых значений эндогенной переменной {у } в части временных срезов для некоторых предприятий. Здесь Б гипотетически полные данные.
Требуется построить логистическую динамическую нейросетевую модель оценки вероятности риска банкротства вида
р(ї) = і/[1 + ЄХР(-У(х(ї), і))], (3)
где показатель экспоненты У(х(X), X) восстанавливается по неполным данным О с помощью нейросетевого отображения вида
У (х, X) = ¥ (х,Ж, X), (4)
Здесь оператор нейросетевого отображения
¥ : х е Я(п) ^ У е Я(1); (5)
где Я(п+1), Я(1) - пространства действительных чисел размерности п и 1
соответственно; W - матрица параметров (синаптических весов) нейросети.
Поскольку требуемая нейросетевая модель (3) имеет четко выраженную прикладную
направленность, то постановку задачи приблизим к реальным сложным условиям
моделирования:
1) помимо оговоренной неполноты данных и возникающей при этом
неопределенности в оценке динамической ситуации риска банкротства, будем считать, что данные О' ^ О сильно зашумлены (вплоть до сознательного искажения части данных, что может иметь место в бухгалтерской отчетности);
2) априорный вид закона распределения шумов неизвестен;
3) количество наблюдений в выделенном кластере предприятий ограничено
несколькими сотнями.
В литературе такие сложные условия моделирования относят классу задач с триадой «НЕ-факторов» (неполнота, неопределенность, неточность данных).
2. Описание алгоритма предлагаемого гибридного метода
Алгоритм метода (Рис. 1) содержит два итерационных цикла:
• внутренних итераций с индексом к, к = 1,2,..., N -1, целью которых является восстановление эндогенной переменной {у } во временных срезах {їк };
• внешних итераций с индексом S = 1,2,..., целью которых является введение в нейросетевую модель вида (4) информации о времени и, соответственно, получение нейросетевой динамической модели с непрерывным временем.
Построение промежуточной квазистатической нейросетевой модели оценки
риска банкротства (КСНСМ)
Опишем алгоритм предлагаемого итерационного метода по шагам.
Шаг 1. (к=1; 8=1) Используя данные последнего временного среза (і=Щ, где известны вектор факторов х ы и метки «банкрот» (РgN=1), «небанкрот» (РgN=0) для всех объектов
1, G
и, соответственно значения аргумента экспоненты = -7 и = 7, в блоке III
строим статическую нейросетевую модель (СНСМ) III, например, типа многослойного персептрона (MLP) с алгоритмом обучения типа «обратного распространения ошибки» (BP)
[4].
На первом шаге внешних итераций (£=1) важно восстановить наиболее точно функциональную связь У (х), скрытую в данных, в сложных условиях триады «НЕ-факторов», оговоренных во введении к статье. От качества НСМ на первом шаге будет на последующих шагах зависеть качество восстановления динамических закономерностей P(х, ц) во
временных срезах ц, k = 1,2,..., N -1. Поэтому применим для построения КСНСМ
модифицированные метод вложенных математических моделей (МВММ) подробно описанный в [5, 6], который обеспечивает регуляризацию НСМ на основе байесовского подхода.
Процедуры 11-17, входящие в МВММ, по сути являются процедурами предрегуляризации исходных «сырых» данных, сформированных в блоке 11. по стандартным данным бухгалтерской отчетности. Эти процедуры, кроме \.7, совпадают с аналогичными процедурами схемы МВММ из [5, Рис. 3.5]. Общей концепцией предобработки (предрегуляризации) данных в блоках \.1 -\.7 является последовательное (сверху - вниз) сужение пространства переменных [X х У ] и, соответственно, повышение качества (однородности и информативности) данных.
Сделаем некоторые комментарии к оригинальным процедурам \.4 и \.7, которые детально описаны в [5, 6]. Формирование систем показателей факторов в процедуре \.4 - это по существу формализованный оптимальный выбор по минимума критерия Кт взвешенной
суммы ошибок \ и \\ рода идентификации объектов («банкрот» - «не банкрот»). {Ик} - это множество анализируемых систем показателей, играющих роль априорных байесовских гипотез в ансамбле согласно квазибайесовскому методу регуляризации нейросетевой модели из [7]. Критерий Кт определяется формулой:
зафиксированной байесовской мета-гипотезы N с номером к; Ы1т, - количество ошибок
первого и второго рода соответственно при идентификации; г , г , г - удельные веса показателей в (6), характеризующие их предпочтения, т.е. кофициенты Фишберна [8]:
п - количество ранжируемых показателей, входящих в критерий (6). Формулу (7) можно использовать в случае, если множество показателей N*m, , N^ упорядочено в
порядке убывания их предпочтения (> N^11) .
В блоке 1.7 применяется нейросетевой аналог известного метода «включения» факторов в регрессионную модель [8]. Для выбранной в блоке 1.4 системы факторов (мы использовали 16 показателей А. О. Недосекина [7, 9], разбитых на 4 группы:
«рентабельность», «ликвидность и платежеспособность», «деловая активность», «финансовая устойчивость»). Вначале выбрано по одному фактору из каждой группы и построены вспомогательные НСМ вида (2), (5). Затем на каждом шаге итераций добавлялось по одному новому фактору, выбранному случайно. После полного перебора факторов по индексу р = 1,16 фиксируется тот набор р*, который соответствует при тестировании модели условию
(6)
где Nнk - количество верно идентифицированных предприятий для данной
2(п - і +1) (п + 1)п ’
(7)
p* = ш^тKнк,Р\Н*) .
P
(8)
Нам удалось сократить размерность факторного пространства с 16 до 6 без заметного снижения критерия (8), т.е. процента правильной идентификации объектов по банкротству.
В блоке III подмодели II проводится кластеризация предприятий в пространстве кортежей (у8к, Х^ У g = 1, G по усеченному в блоке I.? вектору факторов х'^ ^ х^. Целью
кластеризации является дальнейшее углубление степени однородности данных в аспекте учета различия динамики вероятности банкротства Р(г) в разных кластерах (см. ниже, Рис. 2).
«Различие» состоит в том, что значения факторов {х} в последнем временном срезе ^ = Т предопределяют восстанавливаемые значения (у } в КСНСМ IV и динамической модели V.
Для кластеризации можно использовать нейросети Кохонена [10], либо более простой метод к-средних [11].
Мы выделяли 3 кластера и дальнейшее построение моделей КСНСМ с дискретным временем }, к = 1, N -1 и ДНСМ с непрерывным временем г є [0;Т] осуществлялось отдельно для каждого кластера.
В блоке Ш.1 осуществляется регуляризация СНСМ оценки риска банкротства на байесовском ансамбле нейросетей на основе МВММ, подробно описанного в [5, 6]. Качество регуляризированной модели определяется критерием
тах[кнкр IН*] ^ т =1,3 > (9)
т
где т - номер кластера; р* - номер наилучшего набора факторов; Н - наилучшая система факторов (мета-гипотеза); % - экспертно задаваемое число.
В итоге шага 1 в блоке Ш.2 для фиксированного времени ік = Т в последнем
временном срезе мы получаем статическую логистическую нейросетевую модель оценки риска банкротства (СЛНСМ) (2), (5):
У (х) = ^ (х,Ж), (10)
где ^ (•) - оператор нейросетевого отображения [10].
Квазистатический нейросетевой метод оценки риска банкротства (КСНСМ) с
дискретным временем
Шаг 2. (к=0,1,2,...,Ы-1; 8=1). Далее в блоке IV будем восстанавливать эндогенную переменную (у к} в каждом временном срезе во внутреннем итерационном цикле
(к=0,1,2,...,Ы-1), используя статическую нейросетевую модель (2), (10), полученную на шаге
1, для каждого §-го предприятия. Итогом этой процедуры является множество вероятностей риска банкротства во временных срезах:
РУ (*(<к))]=!/[!+ ехР(Ч(Х('к)))]; к = 1,2,..., N -1; (11)
к (гк )}к:г=Pgh^ рл-.,; g=уаг • (12)
Замечание 1. Модель (4, 5), (11), (12) - это уже квазистатическая нейросетевая модель оценки риска банкротства с логистической трансформацией вероятности (ЛКСНСМ) и дискретным временем (гк }. Здесь уже можно приближенно оценивать стадии развивающегося
процесса банкротства по вероятностям во временных срезах. Однако, динамическая информация в моменты времени между срезами отсутствует.
На языке функционального анализа модель (5), (11), (12) представляет собой композицию двух нелинейных операторов:
х е Я(п) ——^Р{Хк) еМ(Ю; р = А о Е; к = 1,N -1, (13)
где Г - оператор нейросетевого отображения вида (4)-(5); А - оператор логистической трансформации вероятности вида (11); М(^ - конечномерное пространство дискретных функций размерности N.
Динамический нейросетевой метод (ДНСМ) оценки вероятности риска банкротства с непрерывным временем
Шаг 3. (к=1,2,...,Ы-1; £=2). На шаге 3 в блоке V алгоритма мы переходим к внешним итерациям введения в динамическую модель непрерывного времени г е [0;Т]. Для этого в каждой кластерной НСМ в состав вектора факторов х введем в качестве независимого нового фактора - время Г.
X ^ х2 ,...;> х} ,.. х„ , х„+1 = Х) . (14)
Используя восстановленные данные, т.е. кортежи
(У8к, хяк, *к), к = 1,2,..., N - 1N , (15)
обучаем заново и тестируем нейросеть:
У (х(г), г) = Е (х(г), г,Ж), (16)
где в нейросетевой модели время I входит явно и является непрерывным. Значит модель (3), (16) можно в полной мере считать динамической, поскольку при обучении НСМ учитывается, что нелинейная зависимость (3) «плывет» во времени. Причем можно получить информацию о вероятности риска банкротства, как во временных срезах, так и в любой точке между ними. Следовательно, появляется возможность уточнения восстанавливаемых значений эндогенных переменных в кортежах:
У2} = Е(х(г), Хк,Ж(1)), к = 1,2,..., N -1, (17)
где синаптические веса Ж(1) берутся из предыдущего шага внешних итераций (£=1).
Шаг 4. (к=1,2,...,Ы-1; £=3). Заново обучаем ДНСМ с учетом новых восстановленных данных, т.е. реализуем шаги 2, 3 и т.д.
Критерий сходимости внешнего итерационного процесса по индексу £ восстановления эндогенных переменных у к и уточнения ДНСМ выбираем в виде:
тах
к
у2} - у2_1)
<£; я = 1,0; к = 1, N -1; £ = 3,4,..., (18)
Замечание 2. Вычислительный процесс построения ДНСМ является довольно трудоемким, поскольку восстанавливать приходится эндогенные переменные {у(£)} на шагах
3,4,... для всех предприятий, входящих в кластер да и во всех временных срезах. Это естественная плата за переход от квазистатической нейросетевой модели к динамической для
данного сектора экономики и данного региона (страны). Зато последующее использование на практике ДНСМ довольно просто - оно сводится к подстановке на входе обученной и верифицированной нейросети, единой для всего кластера модели вектора х(ї) факторов для анализируемого £-го объекта в интересующих момент времени Ї.
Замечание 3. Основная ценность модели банкротства, получаемой с помощью ДНСМ
- это возможность получения качественного прогноза банкротства с информацией о стадиях его развития.
3. Модификация динамического нейросетевого метода оценки риска банкротства (МДНСМ)
На практике встречаются экономические объекты с большой инерцией. В таких объектах при исследовании динамики риска банкротства важно учесть «медленное время» с помощью лаговых переменных разного типа.
Авторы статьи предлагают, аналогично [12], ввести в состав вектора факторов х* также лаговые переменные вида:
р-і; ^ р-2; ^ р-і /р-2 • (19)
Эффективность такого приема можно оценить только путем вычислительных экспериментов по критерию скорости сходимости итерационного процесса (18).
4. Количественные оценки построения динамической модели банкротств
Использовались реальные данные для строительных организаций из [7].
Этап 1 (Итерация 0). Был обучен ансамбль из десяти сетей на выборке из 136 наблюдений для 1=5 (последний временной срез). Результаты тестирования и структура сетей представлены в таблицах 1 и 2.
Таблица 1
Структура байесовского ансамбля
№ сети Количество скрытых слоев Оптимальное количество нейронов в скрытых слоях Вид активационной функции в скрытых слоях
НСМ1 1 7 Г иперболический тангенс
НСМ2 1 6 Сигмоид
НСМ3 2 14, 9 Г иперболический тангенс в обоих слоях
НСМ4 2 14, 11 Г иперболический тангенс в первом слое, сигмоид - во втором
НСМ5 2 8, 13 Сигмоид в первом слое, гиперболический тангенс - во втором
НСМ6 2 13, 12 Сигмоид в обоих слоях
НСМ7 3 14, 13, 14 Г иперболический тангенс в первом и третьем слое, сигмоид - во втором
НСМ8 3 8, 14, 12 Г иперболический тангенс во всех слоях
НСМ9 3 12, 14, 10 Г иперболический тангенс во всех слоях
НСМ10 3 14, 13, 15 Г иперболический тангенс во всех слоях
Рис. 1. Логическая схема алгоритма динамического метода оценки вероятности риска
банкротства предприятий (начало)
Рис. 1. Логическая схема алгоритма динамического метода оценки вероятности риска
банкротства предприятий (продолжение)
Таблица 2
Результаты тестирования нейросетей, обучавшихся на данных за последний год (1=5)
№ сети Доля верно распознанных предприятий Всего ошибок идентификации Количество ошибок 1-го рода Количество ошибок 2-го рода Кн,
НСМ1 65,85% 14 5 9 10,33
НСМ2 80,49% 8 5 3 14,33
НСМ3 80,49% 8 3 5 14,66
НСМ4 78,05% 9 5 4 13,66
НСМ5 78,05% 9 5 4 13,66
НСМ6 85,37% 6 3 3 16,00
НСМ7 75,61% 10 6 4 12,83
НСМ8 82,93% 7 3 4 15,33
НСМ9 80,49% 8 3 5 14,66
НСМ10 78,05% 9 3 6 14,00
Среднее на отфильтрованном ансамбле 85,37% 6 3 3 16,00
После фильтрации, обученным сетям были предъявлены данные за предшествующие 4 года. Получены осредненные оценки величины Y для всех временных срезов t=1,2,3,4.
Была сформирована общая выборка, состоящая из 680 наблюдений (136*5=680).
Этап 2 (Итерация 1). В общую выборку был добавлен фактор времени (относительное время). Данные были разбиты на обучающее и тестовое множество следующим образом: 610 наблюдений в обучающем множестве и 70 - в тестовом. На этих данных был переобучен ансамбль нейросетей (оптимальное количество нейронов в скрытых слоях подбиралось заново). Результаты тестирования показаны в таблице 3.
Таблица 3
Результаты тестирования байесовского ансамбля нейросетей,
____________обучавшихся на данных за все 5 лет_________________________
№ сети Доля верно распознанных предприятий Всего ошибок Количество ошибок 1-го рода Количество ошибок 2-го рода Критерий КИ, Средняя на отфильтрованном ансамбле ошибка обобщения 5
НСМ1 84,29% 11 5 6 26,83 163%
НСМ2 85,71% 10 1 9 28,16 173%
НСМ3 84,29% 11 6 5 26,66 170%
НСМ4 88,57% 8 4 4 29 172%
НСМ5 85,71% 10 1 9 28,16 177%
НСМ6 71,43% 20 0 20 21,66 171%
НСМ7 87,14 9 4 5 28,33 179%
НСМ8 80,0% 14 3 11 25,16 191%
НСМ9 81,43% 13 6 7 25,33 174%
НСМ10 85,71% 10 3 7 27,83 178%
Среднее на отфильт- рованном ансамбле 87,14% 9 3 6 28,5
Критерий фильтрации КНк > 26 .
Оставшимся после фильтрации сетям (№№ 1, 2, 3, 4, 5, 7, 10) были предъявлены данные и получены значения У .
Дополнительно для каждой нейросети была рассчитана средняя относительная погрешность аппроксимации по формуле
N..
5=1
і=1
Уі
^ NІеагп , (20)
где Ыкат - количество наблюдений в обучающем множестве.
Этап 3 (Итерация 2). В данных значения У (показатель экспоненты) были заменены на значения У, полученные на предыдущей итерации. Цикл обучения повторяется заново. Результаты в таблице 4.
Таблица 4
Результаты тестирования байесовского ансамбля нейросетей, обучавшихся на данных за все 5 лет (итерация 8=2)
№ сети Доля верно распознанных предприятий Всего ошибок Количество ошибок 1-го рода Количество ошибок 2-го рода КН, 5
НСМ1 85,71% 10 1 9 28,16 68.37%
НСМ2 85,71% 10 1 9 28,16 65.98%
НСМ3 85,71% 10 1 9 28,16 78.92%
НСМ4 85,71% 10 1 9 28,16 69.61%
НСМ5 85,71% 10 1 9 28,16 66.38%
НСМ7 85,71% 10 1 9 28,16 70.84%
НСМ10 87,14% 9 1 8 28,83 70.03%
Среднее на отфильтро ванном ансамбле 85,91% 10 1 9 28
Результат: практически все сети (кроме №10) показали один и тот же результат. При этом ошибки идентификации наблюдаются в одних и тех же точках тестового множества. Этот результат можно объяснить тем, что при переходе к следующей итерации нейросетевое отображение «сжимает» невязку для у , т.е. отображение (3) является сжимающим [13].
Относительная погрешность также уменьшилась примерно в 2,5 раза.
Полученная динамическая модель с непрерывным временем позволяет достаточно точно отслеживать характерную динамику вероятности банкротства и своевременно принимать соответствующие управляющие воздействия. Для иллюстрации на рис. 2 приведены графики динамики вероятности банкротств трех предприятий из выборки 136 предприятий строительных организаций из [7].
Рисунки а), б), и в) соответствуют объектам из трех различных кластеров, где кластерообразующие признаки х входят в качестве аргумента в показатель экспоненты У(х,ї) в логистическую функцию (3). Видно, что различие в признаках х трансформируется в
динамических моделях в различие динамики риска банкротства. На рис. а) предприятие, близкое к банкротству в период с 2005 по 2008 годы, благодаря принятым антикризисным мерам в 2009 году улучшило свое финансово-экономическое состояние, и вероятность р иска банкротства снизилась до 0,42.
На рис. б) предпринятые антикризисные меры оказались не достаточно эффективными: вероятность риска банкротства осталась к 2010 году не высоком уровне (Р=0,8).
Рис. 2. Динамика вероятности банкротств
На рис. в) видно принимаемые антикризисные меры в период с 2008года по 2009 год оказались не достаточными и к 2010 году предприятие обанкротилось.
Таким образом, вычислительные эксперименты подтвердили эффективность основной идеи разработки гибридного нейросетевого метода построения динамической модели банкротств с непрерывным временем, т.е. идеи восстановления непрерывной динамической зависимости вероятности риска банкротств предприятий по неполным данным.
5. Ограничения для предложенных модификаций нейросетевого метода оценки риска банкротства
По сути, имеются три ограничения:
1. В последнем временном срезе ^ = Т должны быть известны метки «банкрот-небанкрот» для эндогенной переменной (Р=1 либо Р=0).
2. Выборка должна содержать несколько сотен наблюдений.
3. Среди наблюдений должны быть как предприятия «банкроты», так и
предприятия «небанкроты» (примерно одинаковое количество в последнем временном срезе
= Т ).
Как указывалось во Введении, допускается сильное зашумление данных, т.е. триада «НЕ-фактров».
Выводы
1. Разработан метод, который может быть использован в широком спектре динамических прикладных задач:
• оценки кредитоспособности заемщика;
• обеспечение экономической безопасности корпораций или компаний путем
непрерывного мониторинга их кредитоспособности;
• оценка инвесторами объектов инвестиций;
• оценка в сделках контрагентов;
• оценка инвестиционной привлекательности регионов при стратегическом
планировании развитии территории России;
и др.
2. Цели и задачи исследования достигнуты: предложен и апробирован на реальных данных строительных организаций [7] гибридный метод построения динамических моделей диагностики банкротств при неполных данных, которые подтвердили его правильность и применимость.
3. Разработанные инструментарии и методы выходят за рамки рассмотренных
выше задач и могут быть использованы в других динамических задачах с неполными
данными при моделировании экономических объектов.
ЛИТЕРАТУРА
1. Altman E. I. Financial ratios, discriminant analysis and the predition of corporate bankruptcy. Journal of Finance, 1968, vol. 23, №9, pp. 589-609.
2. Ohlson J. A. Financial Ratios and the Probabilistic Prediction of Bankruptcy. Journal of Accounting Research, 1980,№18(1), p. 109-113.
3. Галушкин А. И. Теория нейронных сетей. Кн. 1: Учебн. пособие для вузов. - М.: ИПРЖ, 2000. - 416 с.
4. Горбатков С. А., Фархиева С. А., Белолипцев И. И. Обобщенная нейросетевая модель риска банкротства предприятия // XVI Всероссийская научнотехническая конференция «Нейроинформатика - 2014» 20: Сборник научных трудов в 3-х частях. Ч.3. - М.: НИЯУ МИФИ. 2013, с. 210-220.
5. Горбатков С. А., Полупанов Д. В., Макеева Е. Ю., Бирюков А. Н. Методологические основы разработки нейросетевых моделей экономических объектов в условиях неопределенности Монография/ Под ред. С.А. Горбаткова.
- М.: Издательский дом «Экономическая газета», 2012. - 494 с.
6. Горбатков С. А., Белолипцев И. И., Фархиева С. А. Приближенный метод
байесовской регуляризации и двухступенчатая оценка адекватности гибридной нейросетевой модели налогового контроля // Научная сессия НИЯУ МИФИ -2011: XIII Всероссийская научно-техническая конференция
"Нейроинформатика-2011":Сб. научных трудов - М.:МИФИ, 2011, с. 144-154.
7. Горбатков С. А., Белолипцев И. И., Макеева Е. Ю. Выбор системы экономических показателей для диагностики и прогнозирования банкротств на основе байесовского подхода и интеллектуальных информационных технологий // Вестник Финансового Университета при Правительстве РФ. - 2013. - №4(76).
- с. 50-62.
8. Горбатков С. А., Полупанов Д. В., Фархиева С. А., Коротнева М. В. Эконометрика: Учеб. пособие / Под ред. С. А. Горбаткова. - Уфа: Изд-во РИЦ БашГУ, 2012. - 204 с.
9. Недосекин А. О. Комплексная оценка риска банкротства корпораций на основе нечетких описаний // http://sedok.narod.ru/sc_group.html
10. Хайкин С. Нейронные сети: полный курс: Учебник / 2-е издание: Пер. с англ. -М. : Издательский дом «Вильямс», 2006. - 1104 с.
11. Айвазян С. А., Мхитарян В. С. Прикладная статистика и основы эконометрики:Учебник. - М.: ЮНИТИ, 1998. - 1023с.
12. Алексеева Ю. А. Оценка финансового состояния и прогнозирование банкротства предприятия: Автореферат дис. ... канд. экон. наук. - М.: «Высшая школа экономики»ю - 31с.
13. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука. Физмалит, 1986-288с.
Рецензент: Черняховская Лилия Рашитовна, доктор технических наук, профессор кафедры «Техническая кибернетика», ФГБОУ ВПО «Уфимский государственный авиационный технический университет».
Stanislav Gorbatkow
The Ufa branch of the Finance Academy under the Government of the Russian Federation
Russia, Ufa E-Mail: sgorbatkov@mail.ru
Ilya Beloliptsev
The Ufa branch of the Finance Academy under the Government of the Russian Federation
Russia, Ufa E-Mail: red7315@gmail.com
Svetlana Farkhieva
The Ufa branch of the Finance Academy under the Government of the Russian Federation
Russia, Ufa E-Mail: ok-xi@yandex.ru
Logistics neural network method for constructing dynamic models of bankruptcies with incomplete data
Abstract: Hitherto not been investigated Neural network dynamic models estimating the probability of the risk of bankruptcy are of great practical interest today. But they have not been sufficiently investigated yet. Knowledge of developmental dynamics of the bankruptcy process allows the manager to determine the stage of bankruptcy and to apply proactive control actions. However, the creation of dynamic models of the bankruptcy risk is complicated by the presence of incomplete data. In this paper, we proposed an original hybrid neural network iterative method for estimating the probability of developing in time the risk of bankruptcy of enterprises. The method is focused on the using in difficult conditions: the incompleteness of the data and their strong noising. Computational experiments have confirmed the effectiveness of the basic idea of creating a dynamic model of bankruptcies with continuous time. The developed method can be used in a wide range of applications: assessing the creditworthiness of the borrower, economic security by continually monitoring their financial condition, investment attractiveness assessment, evaluation of financial and economic viability of the counterparties, etc. The proposed tools are beyond the scope of the above problems and can be used in other dynamic problems with incomplete data for modeling economic objects.
Keywords: Diagnosis failures; incomplete data; neural networks; dynamic model; risk; logistics method; hybrid method; the evaluation.
Identification number of article 136EVN214
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
9.
10. 11. 12.
13.
Altman E. I. Financial ratios, discriminant analysis and the predition of corporate bankruptcy. Journal of Finance, 1968, vol. 23, №9, pp. 589-609.
Ohlson J. A. Financial Ratios and the Probabilistic Prediction of Bankruptcy. Journal of Accounting Research, 1980,№18(1), p. 109-113.
Galushkin A. I. Teorija nejronnyh setej. Kn. 1: Uchebn. posobie dlja vuzov. - M.: IPRZh, 2000. - 416 s.
Gorbatkov S. A., Farhieva S. A., Belolipcev I. I. Obobshhennaja nejrosetevaja model' riska bankrotstva predprijatija // XVI Vserossijskaja nauchno-tehnicheskaja konferencija «Nejroinformatika - 2014» 20: Sbornik nauchnyh trudov v 3-h chastjah. Ch.3. - M.: NIJaUMIFI. 2013, s. 210-220.
Gorbatkov S. A., Polupanov D. V., Makeeva E. Ju., Birjukov A. N. Metodologicheskie osnovy razrabotki nejrosetevyh modelej jekonomicheskih ob#ektov v uslovijah neopredelennosti :Monografija/ Pod red. S.A. Gorbatkova. - M.: Izdatel'skij dom «Jekonomicheskaja gazeta», 2012. - 494 s.
Gorbatkov S. A., Belolipcev I. I., Farhieva S. A. Priblizhennyj metod bajesovskoj reguljarizacii i dvuhstupenchataja ocenka adekvatnosti gibridnoj nejrosetevoj modeli nalogovogo kontrolja // Nauchnaja sessija NIJaU MIFI - 2011: XIII Vserossijskaja nauchno-tehnicheskaja konferencija "Nejroinformatika-2011":Sb. nauchnyh trudov -M. MIFI, 2011, s. 144-154.
Gorbatkov S. A., Belolipcev I. I., Makeeva E. Ju. Vybor sistemy jekonomicheskih pokazatelej dlja diagnostiki i prognozirovanija bankrotstv na osnove bajesovskogo podhoda i intellektual'nyh informacionnyh tehnologij // Vestnik Finansovogo Universiteta pri Pravitel'stve RF. - 2013. - №4(76). - s. 50-62.
Gorbatkov S. A., Polupanov D. V., Farhieva S. A., Korotneva M. V. Jekonometrika: Ucheb. posobie / Pod red. S. A. Gorbatkova. - Ufa: Izd-vo RIC BashGU, 2012. - 204
Nedosekin A. O. Kompleksnaja ocenka riska bankrotstva korporacij na osnove
nechetkih opisanij // http://sedok.narod.ru/sc_group.html
Hajkin S. Nejronnye seti: polnyj kurs: Uchebnik / 2-e izdanie: Per. s angl. - M.:
Izdatel'skij dom «Vil'jams», 2006. - 1104 s.
Ajvazjan S. A., Mhitarjan V. S. Prikladnaja statistika i osnovy
jekonometriki:Uchebnik. - M.: JuNITI, 1998. - 1023s.
Alekseeva Ju. A. Ocenka finansovogo sostojanija i prognozirovanie bankrotstva predprijatija: Avtoreferat dis. ... kand. jekon. nauk. - M.: «Vysshaja shkola jekonomiki»ju - 31s.
Tihonov A.N., Arsenin V.Ja. Metody reshenija nekorrektnyh zadach. M.: Nauka. Fizmalit, 1986-288s.