CC BY
177
ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2019 ЭКОНОМИКА Том 14. № 2
doi 10. 17072/1994-9960-2019-2-268-288
УДК 330.43:336.02
ББК65в661
JEL Code C1. C6. G17
МОДЕЛИРОВАНИЕ И ПРОГНОЗИРОВАНИЕ ДИНАМИКИ КУРСОВ
ФИНАНСОВЫХ ИНСТРУМЕНТОВ С ПРИМЕНЕНИЕМ ЭКОНОМЕТРИЧЕСКИХ МОДЕЛЕЙ И ФРАКТАЛЬНОГО АНАЛИЗА *
Петр Михайлович Симонов
ORCID ID: 0000-0001-6357-662Х. Researcher ID: S-5438-2016, е-mail: [email protected] Роберт Викторович Гарафутдинов
ORCID ID: 0000-0002-2130-9352. Researcher ID: A-9122-2019. e-mail: [email protected]
Пермский государственный национальный исследовательский университет (Россия. 614990. г. Пермь. ул. Букирева, 15)
Задача прогнозирования динамики изменений курсов финансовых инструментов является актуальной. так как ее решение позволило бы снизить риски и увеличить доходность от операций на финансовых рынках. Согласно классическому представлению о природе рынков. процессы ценообразования на них имеют стохастический характер и не поддаются прогнозированию. В последние годы активное развитие получила эконофизика - наука на стыке экономики и физики. применяющая к анализу экономических систем подходы. типичные для изучения физических явлений. К числу таких подходов относится фрактальный анализ. Он базируется на гипотезе фрактального рынка. согласно которой динамика изменений цен финансовых инструментов подчиняется степенным законам и ее прогнозирование возможно. Фрактальность финансовых рядов выражается в их свойстве длительное время поддерживать тенденцию изменений (долгой памяти). Существуют модификации эконометрических моделей. учитывающие фрактальные свойства временных рядов. Их применение на данных российского финансового рынка изучено недостаточно полно: единичны примеры успешных предсказаний. полученных с помощью таких моделей; отсутствуют работы. акцентирующие внимание на сравнении фрактальных и нефрактальных моделей на достаточно больших массивах ценовых данных. что не позволяет с уверенностью говорить о превосходстве моделей. учитывающих фрактальные свойства рядов. Цель исследования - проверка гипотезы о том. что учет фрактальности финансовых рядов при прочих равных условиях позволяет получать более качественные прогнозы курсов финансовых инструментов. Получены следующие результаты исследования. Гипотеза о лучшей прогностической способности фрактальных моделей не отвергается. так как в большинстве случаев точность прогнозов ARFIMA и ARFIMA-GARCH оказалась выше точности прогнозов ARIMA и GARCH. Гипотеза о переменной фрактальной структуре финансовых рядов подтверждается. о чем свидетельствуют изменчивые графики локальных фрактальных размерностей. Гипотеза о преимуществе моделей с долгой памятью перед моделями с короткой памятью не отвергается. так как модель ARIMA показала наименьшую точность прогнозов. Гипотеза об изменении фрактальных свойств рядов при трансформации подтверждается. Это доказывают различные значения показателя Херста и отличающиеся графики локальных размерностей. Не удалось подтвердить следующие гипотезы: использование локальных показателей фрактальности вместо глобальных для вычисления параметра дробного дифференцирования увеличивает точность прогнозов; при росте длины прогноза его точность снижается; участки рядов с долгой памятью лучше моделируются и прогнозируются. Сформулированы возможные направления перспективных исследований: разработка метода определения границ локальных участков финансовых рядов с разными фрактальными свойствами; формулирование рекомендаций относительно того. по какому ряду (исходному или преобразованному) следует оценивать оператор дробного дифференцирования модели ARFIMA; проверка неподтвержденных гипотез с изменением методики исследования; изучение различных фрактальных модификаций GARCH-моделей и выявление условий их применимости.
Ключевые слова: эконофизика, фрактальный анализ, долгая память, персистентность, локальные и глобальные фрактальные характеристики, показатель Херста, метод минимального покрытия, R/S-анализ, дробное дифференцирование ряда, точность прогноза, ARFIMA, ARFIMA-GARCH, финансовые рынки.
* Статья развивает положения научной школы «Функционально-дифференциальные уравнения: конструктивные методы исследования и приложения (в экономике)» кафедры информационных систем и математических методов в экономике экономического факультета ПГНИУ.
© Симонов П.М., Гарафутдинов Р.В., 2019
0 I Данная статья распространяется на условиях лицензии кНШ Creative Commons - Attribution 4.0International ICC BY 4.0)
MODELING AND FORECASTING OF FINANCIAL INSTRUMENTS DYNAMICS USING ECONOMETRICS MODELS AND FRACTAL ANALYSIS*
Pyotr M. Simonov
ORCID ID: 0000-0001-6357-662X, Researcher ID: S-5438-2016, e-mail: [email protected]
Robert V. Garafutdinov
ORCID ID: 0000-0002-2130-9352, Researcher ID: A-9122-2019, e-mail: [email protected]
Perm State University (15, Bukireva st., Perm, 614990, Russia)
The task of forecasting the dynamics of changes in the rates of financial instruments is relevant, since its solution would reduce risks and increase the profitability of operations in financial markets. According to the classical concept of the nature of markets, their pricing processes are stochastic and cannot be predicted. In recent years the general trend in scientific research is econo -physics - the science in the convergence of economics and physics that uses the approaches typical for physical phenomena investigation to analyze economic systems. One of these approaches is fractal analysis. It is based on the fractal market hypothesis, which states that the dynamics of changes in prices of financial instruments is subject to power laws and their prediction is possible. The fractal of financial series is expressed in their property to maintain the trend of changes (so called long memory) for a long time. There are modifications of econometric models that take into account the fractal properties of time series. Their use in domestic markets has not been studied enough: there are few examples of successful predictions obtained using such models; there are no works that focus on the comparison of fractal and non-fractal models on sufficiently large arrays of price data, and it does not allow us to say with certainty about the superiority of models that take into account the fractal properties of series. The main idea of the study is to compare the prediction accuracy of econometric models and their fractal extensions using the same data. The purpose of the research is to verify the hypothesis that taking into account the fractal of financial series, while all other things are equal, allows us to obtain better and more accurate forecasts of financial instruments. The following approaches and methods have been used: fractal analysis and econometric models (ARIMA, GARCH), as well as their modifications, taking into account the property of long memory. For the first time, the accuracy of a large number of similar fractal and non-fractal models has been compared with differences in the nuances of fractal parameter estimation on the same data, the problems of determining the boundaries of financial series sections with different fractal properties and selecting a series (initial or transformed) for evaluating the operator of fractional differentiation of the ARFIMA model have been formulated. We have obtained the following results. The hypothesis of the best predictive ability of fractal models is not rejected, since in most cases the forecast accuracy of ARFIMA and ARFIMA-GARCH was higher than the forecast accuracy of ARIMA and GARCH. The hypothesis of a variable fractal structure of financial series has been confirmed, as evidenced by the variable graphs of local fractal dimensions. The hypothesis of the advantage of models with long memory over models with short memory has not been rejected, since the ARIMA model has demonstrated the least accuracy of forecasts. The hypothesis of the alteration of the fractal properties of time series in the transformation has been confirmed, that has been proven with different values of the Hurst exponent and different graphics of local dimensions. We did not manage to confirm the following hypotheses: the use of local fractal indicators instead of global ones to calculate the fractional differentiation parameter increases the accuracy of forecasts; with the growth of the forecast length, its accuracy decreases; sections of series with long memory are better modeled and forecasted. Possible directions for further research have been formulated: to develop a method for determining the boundaries of local areas of financial series with different fractal properties; to formulate recommendations on which series (initial or transformed) should be evaluated by the operator of fractional differentiation of the ARFIMA model; to verify unconfirmed hypotheses with a change in the research methodology; to study various fractal modifications of GARCH models and identification of conditions for their applicability.
Keywords: econophysics, fractal analysis, long memory, persistence, local and global fractal properties, the Hurst exponent, minimum coverage method, R/S analysis, fractional differentiation of a series, forecast accuracy, ARFIMA, ARFIMA-GARCH, financial markets.
* The article improves the concepts of the scientific school "Functional-differential equations: constructive research methods and applications (in Economics)" that are suggested by the Department of Information Systems and Mathematical Methods in Economics, the Faculty of Economics, Perm State University.
Введение
Функционирование современной рыночной экономики невозможно без привлечения и перераспределения капиталов. Одним из основных механизмов для решения этих задач является финансовый рынок. Финансовые рынки как система сложны и неустойчивы, а протекающие на них процессы, в частности динамика изменений ценовых котировок финансовых инструментов, стохастичны [1]. Природа этой стохастики очевидна, ведь на курсы финансовых инструментов влияет огромное количество различных факторов - экономических, психологических (настроения и ожидания участников рынка), природных, политических и т. д. Этих факторов настолько много, что ни одна модель не в состоянии учесть их все. Срабатывает философский принцип перехода количества в качество, и создается впечатление, что динамика ценовых приращений хаотична и непредсказуема.
В то же время задача прогнозирования будущих значений курсов финансовых инструментов является актуальной, поскольку возможность формирования достаточно точных предсказаний позволила бы снизить риски и увеличить доходность от операций на рынках. Существует много подходов к решению этой задачи: фундаментальный и технический анализ, статистическое и эконометрическое моделирование, различные методы машинного обучения, такие как нейросетевой анализ [2], эконофизические подходы (фрактальный, р-адический анализ [1]) и др., но их эффективность и применимость в тех или иных условиях остается предметом дискуссии и исследований, во многом потому, что сама природа динамики финансовых рынков вызывает споры. Нет единого мнения, является ли процесс ценовых изменений полностью случайным (стохастическим) и потому принципиально непредсказуемым или же это процесс детерминированный: поведение цен зависит от конечного числа каких-то конкретных факторов и поддается моделированию.
Классические методы анализа финансовых рынков, такие как модель опти-
мального инвестиционного портфеля Г. Марковица, модель САРМ У. Шарпа и модель ценообразования опционов Блека -Шоулза, основаны на предположении Л. Башелье о том, что характеристики финансовых активов являются случайными величинами, распределенными по нормальному закону [3]. Теоретической базой, воплощающей данный постулат, является предложенная экономистом Ю. Фама гипотеза эффективного рынка (Efficient Market Hypothesis, EMH) [4], суть которой состоит в том, что вся значимая информация немедленно и в полной мере отражается в ценах, они являются практически непредсказуемыми, а их динамика описывается моделью «случайного блуждания». Соответственно, исследовать данные процессы возможно методами вероятностного подхода.
Однако крупнейшие финансовые кризисы (такие как кризис 2008 г., превратившийся в мировой экономический) не вписываются в постулаты данного подхода. В 1987 г. случился обвал американского рынка, за один день индекс Доу-Джонса упал на 29,2% [3]. Согласно кривой Гаусса такое событие должно происходить раз в миллионы лет. Однако последовали и другие кризисы, вследствие чего стало понятно, что гипотеза эффективного рынка не является универсальной. Было обнаружено, что в нестабильные периоды кривая плотности распределения изменения доходно-стей не соответствует гауссовской и риск наступления больших отклонений значительно (на 5 и более стандартных отклонений) превышает таковой при нормальном законе. Такие распределения получили названия распределений с «толстыми хвостами» (fat tails). Их форма свидетельствует о том, что динамика рыночных процессов не всегда является случайной.
Когда стало очевидным, что классический подход к пониманию принципов развития финансовых рынков не в состоянии объяснить кризисные явления, исследователи начали искать другие подходы и пытаться применять к анализу экономических процессов неортодоксальные методы, ранее для экономики не характерные. Вклад в развитие экономической науки
внесли ученые из таких областей знания. как военное дело (опционы. теория игр. логистика. исследование операций). биология (нейронные сети). информатика (информационные технологии). а также физика [5]. В физике хорошо изучено поведение неупорядоченных стохастических систем; оттуда в экономику пришли такие понятия. как «броуновское движение» и «случайное блуждание». относящиеся к поведению элементарных частиц [1]. Наука. которая использует типичные для физики методы исследований для изучения экономических систем. называется эконофизикой1.
Одним из методов эконофизики является фрактальный анализ - перспективное (о чем свидетельствует достаточно большое количество зарубежных публикаций). но недостаточно представленное в отечественной научной литературе [7] направление исследования финансового рынка. Основателем его стал американский математик Б. Мандельброт. предположивший. что динамика рынков подчиняется степенным законам [3]. Он выдвинул концепцию. альтернативную EMH, - гипотезу фрактального рынка (Fractal Market Hypothesis, FMH). согласно которой процесс ценообразования на рынках не случаен. а глобально детерминирован. прошлые цены влияют на будущие. а рынки обладают фрактальными свойствами. Наличие этих свойств объясняется тем. что для устойчивости рынка на нем должны присутствовать инвесторы с разными инвестиционными горизонтами (от нескольких часов до нескольких лет) [8]. Если рынок имеет один горизонт инвестирования. на нем возникает нехватка ликвидности и. как следствие. паника [9]. Фрактальность финансовых рядов выражается в таком их свойстве. как «долгая память» (long memory). или персистентность. Это свойство процесса поддерживать тенденцию изменения. Автокорреляционная функция персистентного ряда убывает гиперболически медленно. Персистентные процессы характеризуются независимостью от вре-
1 Термин «эконофизика» введен Р. Мантенья и Г. Стенли [6].
менного масштаба: дневные приращения зависят от прошедших дневных приращений в той же степени, что и недельные от прошедших недельных [6]. Существуют расширения авторегрессионных стохастических моделей временных рядов, учитывающие свойство долгой памяти [2], и, согласно утверждениям ряда исследователей, такие модели позволяют получать более точные прогнозы [10; 11].
Применению фрактального анализа к задаче прогнозирования динамики финансовых рынков посвящены немногочисленные исследования. Так, в работе Е.С. Остапенко и Т.А. Дунаевой [12] установлено, что на ценах акций компании Google модель с долгой памятью ARFIMA дает меньшую ошибку прогноза, чем ARMA. М.В. Прудский [10] использовал ARFIMA для моделирования курса доллара к рублю, им были сделаны выводы о том, что модель обладает способностью делать краткосрочные прогнозы и является более точной в сравнении с другими статистическими моделями из-за учета ею фрактальных свойств. И.И. Белолипцевым и С.А. Фар-хиевой [13] был предложен подход к прогнозированию финансовых рядов на основе нейросетевой модели, одним из предикторов которой явился индекс фрактальности, и на примере курсов акций ОАО «Татнефть» было показано, что его использование улучшает прогностические свойства модели. Тем не менее применение фрактальных методов прогнозирования на финансовых рынках изучено недостаточно: единичны примеры успешных предсказаний, полученных с помощью таких моделей; отсутствуют работы, акцентирующие внимание на сравнении фрактальных и нефрактальных моделей на достаточно больших массивах ценовых данных, что не позволяет с уверенностью утверждать о превосходстве моделей, учитывающих фрактальные свойства рядов. Этим обусловлена актуальность данного исследования.
Основная гипотеза исследования заключается в следующем: учет фрактально-сти финансовых рядов при прочих равных условиях позволяет получать более точные прогнозы динамики курсов финансовых
инструментов. Целью исследования является проверка этой гипотезы. В качестве методического инструментария выбран фрактальный анализ, эконометрические модели ARIMA, GARCH и их фрактальные модификации, учитывающие свойство долгой памяти анализируемых процессов. Предмет исследования - динамика курсов инструментов, обращающихся на российском финансовом рынке, а также фондового индекса (на примере индекса ММВБ).
1Сир<
3;
Приступим к достижению поставленной цели и начнем с обзора результатов ранее проведенных исследований. Теоретический обзор
адачей применения фрактального подхода к анализу финансовых рынков и прогнозированию курсов финансовых инструментов занимались отечественные и иностранные авторы. Результаты обзора некоторых актуальных научных работ в данной области представлены в табл. 1.
Таблица 1
Научные работы по фрактальному анализу финансовых рынков
Ключевые достижения Авторы Подходы (методы) Результаты
Подтвержден фрактальный характер финансовых временных рядов Ю.М. Балагула [14] Р/8-анализ, спектральный метод ^РН) Выявлено наличие долгой памяти в рядах биржевых оптовых цен на электроэнергию и фондовых индексов разных стран
А.В. Зиненко [3] Р/8-анализ
G. Caporale, M. Skare [15] Р/8-анализ, спектральный метод (ОРИ), АРРГМА-ИОАРСИ-мо делирование Выявлено наличие двойной длинной памяти (как в рядах доходности, так и в рядах волатильности) реального ВВП Великобритании
S. Zhelyazkova [16] АРЕМА-ИОАРСИ, ИУ-ОАРСИ и ИАРАРСИ-моделирование Выявлено наличие двойной длинной памяти в курсах двенадцати валют к доллару США; установлено, что динамика различных валютных пар хорошо описывается различными фрактальными модификациями ОАРСИ-моделей
Предложены новые фрактальные показатели: размерность минимального покрытия и индекс фрак-тальности Н.В. Старченко, М.М. Дубовиков [17] Метод минимального покрытия, Р/З-анализ, клеточный метод Выявлено, что для определения индекса фрактальности требуется на два порядка меньше данных, нежели для оценивания показателя Херста методом Р/8-анализа (на примере цен акций тридцати компаний, входящих в индекс Доу-Джонса)
Выполнено сопоставление различных фрактальных показателей, показано, что они согласованы, но отличаются точностью Е.К. Кривоносова и др. [18; 19] На примере цен акций Лукойл установлено, что наиболее точно фрактальная размерность определяется методом клеточного покрытия
Обоснована возможность применения фрактального анализа для предсказания кризисных ситуаций объектов разного уровня А.К. Мансуров [20] Метод ДФА Выявлено, что существенное отклонение фрактальной размерности рядов валютных курсов сигнализирует о валютных кризисах (на примере курсов двадцати трех валют к доллару США)
Е.К. Кривоносова [2] Фрактальные (Р/8-анализ, методы клеточного и минимального покрытия) и мультифрактальные (МФ-ДФА, вейвлет-преобразования, анализ с использованием показателя Гельдера) методы Показано, что наиболее точным для определения стабильности работы предприятия из фрактальных методов является клеточный, из мультифрактальных -анализ на базе показателя Гельдера (на примере финансовых показателей компаний, действующих на российском рынке производства и обслуживания нефтегазодобывающего комплекса)
Обосновано, что перси-стентные временные ряды прогнозируются точнее, чем случайное блуждание V. Kulish, V. Horak [21] Р/8-анализ На примере ряда индекса Доу-Джонса установлено, что участки ряда с показателем Херста более 0,5 могут быть предсказаны более точно, чем участки, у которых показатель близок к 0,5; сделан вывод о том, что показатель Херста может служить мерой предсказуемости
Окончание таблицы 1
Ключевые достижения Авторы Подходы (методы) Результаты
Обосновано, что фрактальные модели могут осуществлять точечные прогнозы значений котировок с более высокой точностью, чем нефрактальные Е.С. Остапенко, Т.А. Дунаева [12] К^-анализ, АКБИМА-моделирование На примере цен акций компании Google установлено, что модель ARFIMA дает более точные прогнозы, чем ARMA; сделано предположение о том, что игнорирование наличия длинной памяти приводит к большей погрешности прогноза, чем учет длинной памяти при ее отсутствии
М.В. Прудский [10] К^-анализ, АКБИМА-моделирование Сделаны выводы о том, что модель ARFIMA позволяет с более высокой точностью, чем другие модели, прогнозировать курс доллара к рублю
И.И. Белолипцев, С.А. Фархиева [13] Метод минимального покрытия, нейросетевое моделирование На примере курсов акций ОАО «Татнефть» показано, что использование индекса фрактальности в качестве предиктора улучшает прогностические свойства нейросетевой модели
На наш взгляд, общим пробелом в работах о прогнозировании будущих значений курсов финансовых инструментов является то, что в них, как правило, исследован лишь один временной ряд. Кроме того, не во всех работах [10] произведено сравнение точности прогнозов фрактальных моделей и максимально похожих на них нефрактальных.
При этом необходимо отметить, что во всех исследованиях были получены положительные результаты, продемонстрировавшие эффективность фрактального подхода к анализу финансовых рынков. Как показывает обзор, в основном фрактальный подход используется авторами для подтверждения гипотезы о фракталь-ности и персистентности финансовых временных рядов, а также для исследования возможности предсказывать критические точки (кризисы). Для определения фрактальных свойств рядов чаще всего используется исторически первый и наиболее простой метод - R/S-анализ в совокупности с показателем Херста, реже применяются альтернативные методы - ДФА, клеточный, спектральный. Перспективному мультифрактальному подходу почти не уделяется внимания, несмотря на то, что, по мнению некоторых авторов [2], реальные финансовые ряды обладают спектром фрактальных показателей на разных масштабах и лучше анализируются именно методами мультифрактального анализа. По
всей видимости. это связано с высокой сложностью применения таких методов. Также было обнаружено лишь небольшое количество работ. посвященных точечному прогнозированию значений финансовых временных рядов [10; 12; 13]. причиной чего предположительно является скептическое отношение многих исследователей к самой возможности получения таких достаточно точных прогнозов. В рассмотренных работах наибольшей популярностью пользуются эконометрические модели с долгой памятью (в частности. ARFIMA). которые. согласно выводам авторов. обладают высокой прогностической способностью. Тем не менее применение таких моделей на российском финансовом рынке изучено недостаточно хорошо. работы и. соответственно. примеры полученных данным методом успешных предсказаний единичны. поэтому требуется дальнейшее исследование их эффективности. Не было обнаружено попыток применения для прогнозирования фрактальных расширений GARCH-моделей. в то время как теоретически они должны лучше. чем ARFIMA. описывать реальные финансовые ряды. для которых характерно явление кластеризации волатильности. Кроме того. авторам не встретилось сравнения фрактальных и нефрактальных моделей на достаточно больших массивах ценовых данных. Как правило. для исследования выбирается лишь один финансовый инстру-
мент, в ряду курсов которого моделируется фиксированной участок, на основании чего делаются выводы о прогностических способностях моделей. Современные инструментальные средства и вычислительные возможности аппаратного обеспечения позволяют проводить массовое моделирование, сравнивая сотни и даже десятки тысяч моделей на множестве участков ценовых рядов. Настоящая работа не претендует на полномасштабное исследование качества различных моделей на обширном эмпирическом материале, однако, как показали его результаты, можно говорить о некоторых преимуществах фрактальных модификаций в различных условиях.
Прежде чем переходить к методике исследования, следует кратко коснуться теоретических аспектов фрактального анализа и его категориального аппарата.
Размерность - это число измерений геометрического объекта или множества. Топологической размерностью (Ох) называют минимальное количество координат, необходимое для описания объекта. может принимать только целые значения, подходит для описания идеализированных, гладких множеств, таких как прямая линия, куб, сфера [17]. Но в природе такие практически не встречаются. Для описания реальных объектов подходит размерность Хаусдорфа (Он), являющаяся обобщением понятия топологической размерности ОТ. Он идеальных объектов (гладких кривых, поверхностей) совпадает с их Ох, в случае же самоподобных структур (таких, как береговая линия) Он всегда превышает ОТ и является нецелой (дробной) величиной. Такая размерность показывает, насколько плотно и равномерно элементы множества заполняют пространство. Наконец, Б. Мандельбротом введен термин «фрактальная размерность» (О) для обозначения частного случая размерности Хаусдорфа (дробной). Множества, для которых справедливо неравенство Он>Ох, Б. Мандель-брот назвал фракталами. В дальнейшем в нашем исследовании термин «размерность» и обозначение О относятся именно к этому типу размерности.
Отметим, что не существует точного определения сущности категории «фрактал». Сам Б. Мандельброт охарактеризовал ее так: «Фракталом называется структура, состоящая из частей, которые в каком-то смысле подобны целому» [22]. Принято относить к фракталам объекты, обладающие следующими свойствами [10]: самоподобие или масштабная инвариантность (отдельные составляющие объекта подобны всему объекту целиком, вследствие чего при различных масштабах объект выглядит одинаково); масштабирование по степенному закону (изменение характеристик объекта с изменением масштаба рассмотрения имеет степенную зависимость); фрактальная размерность, которая превосходит топологическую (при увеличении соотношение между масштабом и каким-либо параметром фрактала постоянно и равно фрактальной размерности [2]).
Сущность фрактального анализа как инструмента описания различных объектов заключается в том, что производится обработка масштабированием фрактальной структуры и описывается распределение какой-либо структурной характеристики (или меры) при этом масштабировании [18]. Фрактальными свойствами могут обладать не только математические абстракции и природные сущности, но и временные ряды, являющиеся численным дискретным описанием непрерывных динамических процессов. Как правило, когда говорят об анализе и прогнозировании динамики финансовых рынков, подразумевают исследование именно ценовых рядов финансовых инструментов. Если обратить внимание на графики котировок акций или валют, можно заметить их сильную изре-занность (признак фрактальной размерности), а также то, что при рассмотрении на различных временных масштабах они выглядят похоже, нельзя с уверенностью сказать, дневные это курсы, месячные или минутные (самоподобие). На свойство масштабирования по степенному закону указывает форма распределений с «толстыми хвостами».
Если для пространственной структуры основной фрактальной характеристи-
кой является ее размерность, то в случае временного ряда говорят о наличии у процесса «долгой памяти», когда ему присуще свойство персистентности. Эти понятия тесно взаимосвязаны. Чтобы численно охарактеризовать свойство «долгой памяти», британским гидрологом Г. Херстом был введен показатель Н, позднее названный его именем [14]. Кроме того, Н.В. Стар-ченко был предложен оригинальный метод оценивания фрактальной размерности фи-
нансовых временных рядов и введен новый показатель - индекс фрактальности (ц) [8].
Финансовые временные ряды демонстрируют сложное непериодическое поведение, при котором тренды и флэты хаотическим образом сменяют случайное блуждание. Рассмотренные фрактальные показатели позволяют идентифицировать состояние процесса и делать выводы о его дальнейшем поведении. В табл. 2 приведено их сопоставление [19].
Таблица 2
Взаимосвязь показателей фрактальности временного ряда
Показатель Диапазон значений Характер ряда
Антиперсистентный (флэт) Случайный (стохастический) Персистентный (тренд)
Показатель Херста, Н [0; 1] 0 < Н < 0,5 Н = 0,5 0,5 < Н < 1
Фрактальная размерность, Б [1; 2] 1,5 < Б < 2 Б = 1,5 1 < Б < 1,5
Индекс фрактальности, ц [0; 1] 0,5 < ц < 1 ц = 0,5 0 < ц < 0,5
Из таблицы 2 видно, что все три показателя согласуются друг с другом при индикации состояний хаотического процесса. В работе [19] дана следующая характеристика этих состояний:
1) Тренд - персистентный участок, на котором поддерживается тенденция изменений. Если в текущем периоде показатель увеличивался, с большой долей вероятности он продолжит расти и в следующем периоде.
2) Флэт - антиперсистентный участок, направление изменений на котором постоянно меняется (чаще, чем при случайном характере процесса). Такой характер динамики называют «возврат к среднему». Если в текущем периоде показатель увеличивался, с большой долей вероятности в следующем периоде направление изменения показателя сменится, в результате чего его значение не может далеко уйти от среднего.
3) Случайное блуждание - промежуточное состояние между трендом и флэтом. Динамика процесса случайна, будущие значения не зависят от прошлых.
Охарактеризовав теоретические основы фрактального анализа, рассмотрим методы оценивания фрактальных характеристик финансовых рядов. Как показал обзор литературы, существует достаточно
большое количество способов вычисления как фрактальной размерности, так и показателя Херста. Рассмотрим два метода, которые и будут применены в исследовании.
Метод минимального покрытия, детально разработанный Н.В. Старченко [8], позволяет характеризовать локальную динамику процесса. Для этого сужается репрезентативный масштаб до значений, при которых временной ряд не меняет своего поведения, и определяется локальная фрактальная размерность [19]. Алгоритм метода состоит из следующих шагов. Пусть процесс характеризуется некоторой функцией у = /(^ на отрезке [ а = < < ^ < — < = Ъ], i = 1, 2, ..., m. Отрезок [а, Ь] разбивается на m сегментов равной длины 8 = = — _1. Затем грат 1 1 1 г
фик функции покрывается прямоугольниками с основанием 3 таким образом, чтобы это покрытие было минимальным по площади. Тогда высота прямоугольника на отрезке t¡\ будет соответствовать амплитуде Аг(3), равной разности максимального и минимального значения А^) на данном отрезке. Накопленная амплитуда на всем отрезке вычисляется как У^(8) = = а площадь минимального по-
крытия как 5^(5) = У^(8)8. Очевидно, что при сокращении длины отрезка 3 точность
вычисления площади покрытия будет увеличиваться. Имеет место степенная зависимость S^(ô)~ô2-D при ô — 0, где D - фрактальная размерность. Из этой формулы следует, что Vf(5)~5-^ при ô —> 0, где д = D^ — 1. Н. В. Старченко назвал показатель Dp размерностью минимального покрытия, а ц - индексом фрактальности. Показатель ц вычисляют как взятый с отрицательным знаком коэффициент при независимой переменной в уравнении линейной регрессии InVf (ô) = —^Inô + Ъ, где b - свободный член. Показано, что для приемлемо точного определения размерности минимального покрытия достаточно репрезентативных интервалов длиной 32 и иногда даже 16 дней.
Метод нормированного размаха, или R/S-анализ (Rescaled range analysis), предложен Г. Херстом в 50-х гг. XX в. и до сих пор является одним из наиболее популярных подходов в исследованиях фрактальных рядов самой различной природы, что подтверждается обзором научных работ. Можно сказать, что для вычисления H он является «родным» методом, и некоторые авторы, говоря о показателе Херста, подразумевают его оценивание именно посредством метода нормированного размаха [3; 21]. Идея R/S-анализа состоит в том, что существует степенная зависимость вида ' ÔH, где R - размах вариации (под
вариацией показателя понимается накопленное отклонение его от среднего значения), S - стандартное отклонение показателя, ô - количество значений показателя в группе, c - некоторая константа. Вычисляемая статистика R/S является размахом, нормированным стандартным отклонением, что и дало название методу. В работе [3] подробно описан его алгоритм. Пусть дан временной ряд x(t) длиной n+1. Прежде всего его приводят к так называемым «логарифмическим доходностям» y(t) по
формуле y(t) = . Далее работают с
xt-i
преобразованными данными. Ряд y(t) длиной n делится на m групп значений по ô элементов. Для каждой i-й группы (i=1, 2, ..., m) рассчитываются: среднее
значение ; накопленные отклонения от среднего = — Уд, образуя т ря-
дов Уг(0; размах R¿ = тах(У;(0) — тт(У^(0); стандартное отклонение
Si=llïôt=1(yt-yi)2;
R/S-статистика
(нормированный размах) по формуле Ri/Si. После полученные m величин R/S усредняются, и получается двойка элементов < R/S(S),S >. Такие вычисления проводят для различных значений S. Показатель Херста H определяется как коэффициент при независимой переменной в
уравнении линейной регрессии =
= H In S + с, где с - свободный член. Стоит отметить, что величина S не может принимать значения, меньшие 10. При этом очевидно, что она является собственным делителем количества элементов в ряду n, поэтому выборку нужно формировать таким образом, чтобы максимизировать число собственных делителей - это положительно скажется на точности вычисления H. Кроме того, для достижения приемлемой точности значение n должно быть достаточно высоко (несколько тысяч) [9]. В работе [17] показано, что на интервалах порядка нескольких десятков значений оценивание H не имеет смысла, так как функция R/S(S) имеет медленный выход на асимптотический режим.
Фрактальную размерность можно вычислить через показатель Херста по формуле D=2-H [17].
Далее, основываясь на методологии фрактального анализа, охарактеризуем методику исследования.
Методика исследования
Для моделирования динамики курсов финансовых инструментов нами были выбраны две популярные эконометрические модели - ARIMA и GARCH - и их фрактальные модификации ARFIMA и ARFIMA-GARCH.
ARIMA (p, d, q). Интегрированная модель авторегрессии - скользящего среднего порядков p, d, q - относится к линейным моделям и предназначена для описания нестационарных процессов. Данная
модель предложена Дж. Боксом и Г. Джен-кинсом [23]. Ее уравнение имеет следующий вид: Ф(Ь)(1 — Ь)аХ,: = ^ + &(Ь)е1 , где Хг - исследуемый процесс (временной ряд), Ь - оператор сдвига, Ф(Ь) - полином степени р от Ь, ®(Ь) - полином степени q от Ь, ё - порядок интегрирования процесса Хг. Параметры р, ё, q могут принимать целые неотрицательные значения.
ЛКИМА (р, ё, q). Дробно-интегрированная модель авторегрессии - скользящего среднего порядков р, ё, q - является обобщением модели ЛК1МЛ (р, ё, q), допускающим дробное значение параметра ё. Предложена Дж. Хоскингом [24]. Уравнение для этой модели остается неизменным, а оператор дробного дифференцирования аппроксимируется следующим рядом:
(1-Dd = T,
r(k-d) lk
Lk, где Г - гам-
'к-0 г(_й)Г(к+1)' ма-функция [11]. Забегая вперед, отметим, что для программной реализации алгоритма моделирования АКИМА нами была использована следующая формула вычисления ряда дробных разностей У1: У1 = Х1 +
+ ltk~=1i(-1)
k П?=М—+1)
k!
xt-k [25]. Каждая
разность аппроксимируется предыдущими значениями исходного ряда, при этом «глубину» памяти, т. е. количество участвующих в вычислении значений, можно ограничить. В этой особенности проявляется учет моделью свойства «долгой памяти» процесса. Для оценки параметра ё существуют специализированные методы [14], а также его можно вывести через соотношение с показателем Херста: й = Н — 0,5 [9].
ОЛКСН (р, q). Обобщенная авторегрессионная условная гетероскедастич-ность порядков р, q являет собой семейство нелинейных моделей, предназначенных для изучения волатильности в ценовых рядах и позволяющих учитывать такие ее свойства, как кластеризация («группирование» периодов с высокой и низкой во-латильностью), а также так называемый эффект «рычага» - уменьшение волатиль-ности при росте цены, и наоборот. Классическая модификация, называемая иногда БОЛКСН, предложена Т. Боллерслевом [26]. Стандартно для моделей GARCH
определяется два уравнения - условного математического ожидания доходности у] и условной дисперсии о]2 [27]. Формула классической модели имеет следующий вид: о] = а0 + Щ • У]_1 + 1%1р] • о,
J ] ut-j-
Для корректного определения условной дисперсии должны выполняться ограничения а0 > 0,ai > 0,ßj> 0. В этой модели
2
ряд у2 удовлетворяет процессу ARMA (max [p, q], q) модели. Величина служит условной дисперсией процесса ytIYt-1~N(0,ot2). На практике часто используется модель GARCH (1, 1) [28].
ARFIMA-GARCH (p, d, q). Данная модель является модификацией SGARCH-модели, в которой в качестве уравнения условного математического ожидания выступает не ARMA, а ARFIMA [29]. Таким образом, учитываются фрактальные свойства моделируемого процесса.
В процессе исследования предполагается проверить ряд гипотез, высказанных авторами научных работ по данной тематике, а также выдвинутых авторами настоящей статьи:
1) Учет моделями фрактальных свойств моделируемых процессов позволяет повысить точность прогнозов [10].
2) Финансовые временные ряды обладают переменной структурой, их фрактальные свойства на различных участках различны [19; 30].
3) Использование локальных показателей фрактальности (оцененных на моделируемом участке ряда) вместо глобальных (оцененных по всему ряду) для вычисления параметра дробного дифференцирования положительно сказывается на точности прогнозов [7].
4) Учет при моделировании «долгой памяти» процесса в том случае, если в действительности он ею не обладает, не ухудшает качества прогнозов, в то время как непринятие во внимание данной особенности процесса при ее наличии снижает их точность [12].
5) При росте длины прогноза его точность снижается.
6) При переходе к лог-доходностям фрактальные свойства ряда меняются.
Для проверки перечисленных гипотез в качестве исходных данных для моделирования следует взять ценовые ряды различных финансовых инструментов с достаточно длительной предысторией (несколько тысяч значений). Предполагается выполнить преобразование цен yt к логарифмическим доходностям по формуле хк =
= [3] и далее моделировать эти до-
yt-i
ходности. Так как предположительно ценовые ряды обладают переменной фрактальной структурой. влияющей на качество моделирования и прогнозирования. не следует аппроксимировать модели на всей истории рядов. необходимо выбрать относительно короткий локальный участок. В то же время проблема определения границ таких участков остается актуальной. в литературе не было обнаружено соответствующей методики. Было решено применить так называемый метод «скользящего окна» (sliding window) [31] и рассмотреть все возможные локальные участки рядов с их начала до конца с шагом 1. на каждом попытаться аппроксимировать модели и проверить их прогнозирующую способность.
В машинном обучении существует задача выявления «эффекта от воздействия» (treatment effect) [32]. Типичный пример ее решения - определение эффективности лекарственных препаратов. Идея метода заключается в том. чтобы оказать на объект некоторое воздействие. оценить изменение его состояния. затем не оказывать воздействий и вновь оценить изменение состояния. Разница между изменениями и будет предполагаемым эффектом от воздействия. В идеале воздействовать и не воздействовать нужно на один и тот же объект одновременно. В жизни реализовать это затруднительно. однако при анализе ретроспективных ценовых данных такой подход становится возможным. Поэтому в исследовании решено применить данную концепцию при сравнении моделей. В качестве воздействия будет выступать учет моделями свойств долгой памяти процессов. Для этого следует тестировать модели в максимально похожих условиях. т. е. на одних и тех же данных. Причем обе
ARMA-модели (фрактальная и обычная) должны пройти проверки на адекватность: оценки коэффициентов должны быть статистически значимы (так как не всегда наилучшая аппроксимированная к данным модель обладает только значимыми коэффициентами, решено принимать модели, у которых их более 50%); в остатках должна отсутствовать автокорреляция (проверка тестом Льюнга - Бокса); остатки должны быть распределены по нормальному закону (проверка тестом Лиллиефорса) [28]. Кроме того, отбрасываются модели «белого шума» ARMA (0, 0) и «случайного блуждания» ARIMA (0, 1, 0), т. е. модели, сумма порядков p, q которых меньше 1. В качестве показателя точности прогнозирования решено применять среднюю абсолютную ошибку МАЕ = ^Libt — 9tl где yt - фактическое значение наблюдения, yt - прогнозное значение наблюдения, T - длина тестового участка. В качестве локальной фрактальной характеристики рядов решено использовать размерность минимального покрытия D^, вычисляемую методом минимального покрытия. Данный метод выбран потому, что позволяет оценивать размерность на небольшом количестве наблюдений (от 16) с приемлемой точностью [17]. В качестве глобальной фрактальной характеристики (оцениваемой по всей длине ряда) выбран показатель Херста H, вычисляемый методом R/S-анализа. Конкретных рекомендаций по величине участка аппроксимации (длине скользящего окна) в литературе обнаружено не было, поэтому он составляет 128 наблюдений, такое значение близко к круглому значению 100 и при этом является степенью двойки (как того требует специфика метода минимального покрытия). Предполагается делать краткосрочные прогнозы глубиной 1, 3 и 10 шагов. Величины MAE, полученные за один проход по ряду, усредняются по каждой из моделей, такие средние значения предполагается сравнивать между собой.
Программную реализацию методики решено осуществлять на языке программирования R в среде разработки RStudio. Для
вычисления показателя Херста, реализации процедур дробного дифференцирования и метода минимального покрытия были написаны пользовательские функции в соответствии с приведенными в литературе алгоритмами. Аппроксимация моделей ARIMA и ARFIMA выполняется с использованием готовой функции auto.arima() из пакета forecast. Данная функция отбирает из нескольких конкурирующих моделей наилучшую по информационному критерию AIC. Аппроксимация моделей GARCH осуществляется через функции ugarchspec() (задание спецификации модели) и ugarchfit() (непосредственная идентификация параметров модели) из пакета rugarch. Использована базовая модификация SGARCH. При обучении модели ARFIMA-GARCH процедуру вычисления дробных разностей, как и оценивание величины d, решено доверить встроенным в пакет rugarch-алгоритмам.
Далее выполняется сравнение четырех авторегрессионных моделей: ARIMA (p,d,q); ARFIMA (p,d,q) (параметр d вычислен через локальную размерность Dp, оцененную по исходному ряду котировок); ARFIMA (p,d,q) (параметр d вычислен через локальную размерность D^, оцененную по ряду лог-доходностей); ARFIMA (p,d,q) (параметр d вычислен через глобальный показатель Херста H, оцененный по всему ряду). Параметры p, q всех моделей оценены по ряду наблюдений.
Выполняется сравнение двух GARCH (1,1)-моделей: SGARCH (модель условного математического ожидания -ARMA); ARFIMA-GARCH (модель условного математического ожидания - ARFIMA).
Прежде чем переходить непосредственно к моделированию и прогнозированию в соответствии с составленной методикой необходимо проанализировать входные данные.
Данные
В качестве моделируемых ценовых рядов выбраны котировки набора финансовых инструментов, на котором проводилось исследование в работе [10]. В набор дан-
ных внесены следующие изменения: исключены акции компании «Газпром» и биржевой индекс RTSI, добавлены валютные пары USDRUB и EURRUB (табл. 3). Использованы дневные цены закрытия.
Таблица 3
Исследуемые финансовые инструменты и индекс
№ п/п Тикер История Количество наблюдений
1 ALRS 11.29.2011 -02.05.2018 1616
2 EURRUB 01.01.2000 -29.04.2018 4637
3 GMKN 31.01.2001 -02.05.2018 4118
4 MTSS 15.10.2003 -02.05.2018 3614
5 ROSN 19.07.2006 -02.05.2018 2948
6 USDRUB 01.01.2000 -29.04.2018 4867
7 VTBR 28.05.2007 -02.05.2018 2736
8 IMOEX 05.01.2000 -02.05.2018 4574
9 LKOH 05.01.2000 -02.05.2018 4575
10 SBER 05.01.2000 -02.05.2018 4563
Источник: ФИНАМ. URL: https://www. finam.ru (дата обращения: 10.01.2019).
На рис. 1-4 приведена динамика показателей ряда на примере индекса ММВБ как индикатора состояния российского финансового рынка (характеристики других инструментов сопоставимы).
Эмпирический показатель Херста H для исходных уровней ряда составил ровно 1 с R2=0,9996, т. е. динамика IMOEX на протяжении без малого 20 лет имеет выраженный тренд, что заметно по графику. Динамика локальной размерности колеблется в диапазоне (1,05; 1,5), что, во-первых, подтверждает наличие у процесса долгой памяти, а во-вторых, укладывается в рамки гипотезы о переменной структуре исходных уровней ряда (размерность меняется в широких диапазонах).
Эмпирический показатель Херста H
для лог-доходностей равен 0,5341 с 2 „ R =0,9965, что говорит о слабо выраженной
персистентности процесса. При этом график лог-доходностей визуально напомина-
ет «белый шум» с нулевым средним, но локальные размерности в основном находятся в диапазоне (1,5; 1,75), т. е. флэта. Поэтому в данном случае глобальные и усредненные локальные фрактальные характеристики не согласуются. Наличие таких особенностей вызвало вопросы о способе расчета параметра ё модели АКБИМА: стоит ли оценивать его по моделируемому ряду, демонстрирующему устойчивую антиперси-
стентность, либо по трендовому ряду исходных уровней.
Коэффициент эксцесса распределения лог-доходностей индекса ММВБ равен 15,3592, коэффициент симметрии -0,2405. Это свидетельствует о том, что распределение не соответствует нормальному, характеризуется вытянутым пиком и «толстыми хвостами», как и ожидалось (для нормального распределения значения коэффициентов составляют соответственно 3 и 0).
Значения индекса
2500 2000 1500 1000 500 0
#> ^ ^ ^ ^ ^ ф ф 4? <$>
Отсчеты
Рис. 1. Значения индекса ММВБ Размерность минимального покрытия
1,5 1,4 1,3 1,2 1,1 1
* # ^ # ^ / / / # / ^ ^ £ ^
Отсчеты
Рис. 2. Динамика фрактальной размерности значений индекса ММВБ
Лог-доходности индекса
0,3 0,2 0,1 0 -0,1 -0,2
-0,3
* # # # ^ / / / # / ^ ^# ^ / /
Отсчеты
Рис. 3. Лог-доходности индекса ММВБ
Размерность минимального покрытия
1,8
1,7 1,6 1,5 1,4
1,3
* # # # ^ ^ ^ / / / ¿р ^ ^
Отсчеты
Рис. 4. Динамика фрактальной размерности лог-доходности индекса ММВБ
В следующем разделе статьи описаны и проанализированы полученные в ходе моделирования результаты.
Результаты и их обсуждение
езультаты моделирования и прогнозирования представлены в табл. 4-6.
гезул!
Р1
Таблица 4
Усредненные показатели MAE при глубине прогноза 1 шаг вперед*
Тикер ARIMA ARFIMA (D„) ARFIMA (Arlog) ARFIMA (H) ARIMA-GARCH ARFIMA-GARCH
ALRS 0,01922** 0,01888 0,01866 0,01826* 0,01949* 0,02044
EURRUB 0,00693** 0,00681 0,00659* 0,00677 0,00509* 0,00516
GMKN 0,01424 0,01395* 0,01408 0,01442** 0,01563 0,01534*
MTSS 0,01454 0,01472 0,01393* 0,01455** 0,01699 0,01678*
ROSN 0,01982 0,02084** 0,01898 0,01816* 0,01664 0,01636*
USDRUB 0,00536 0,00508* 0,00556** 0,00547 0,00588* 0,00595
VTBR 0,01522 0,01504 0,01433* 0,01530** 0,02363 0,02198*
IMOEX 0,01285** 0,01284 0,01270 0,01261* 0,01382 0,01342*
LKOH 0,01608** 0,01561* 0,01575 0,01598 0,01668 0,01602*
SBER 0,01590 0,01531* 0,01591** 0,01551 0,01734 0,01694*
Среднее 0,01401** 0,01391 0,01365* 0,01370 0,01512 0,01484*
* В табл. 4-6 звездочкой (*) обозначен лучший прогноз в своем классе моделей, двумя звездочками (**) -худший прогноз.
Таблица 5
Усредненные показатели MAE при глубине прогноза
3 шага вперед
Тикер ARIMA ARFIMA (D„) ARFIMA (D„-log) ARFIMA (H) ARIMA-GARCH ARFIMA-GARCH
ALRS 0,01770 0,01784** 0,01749 0,01723* 0,01809* 0,01892
EURRUB 0,00598 0,00621** 0,00585* 0,00602 0,00508 0,00506*
GMKN 0,01452 0,01446 0,01409* 0,01462** 0,01606* 0,01607
MTSS 0,01516 0,01535 0,01468* 0,01536** 0,01837* 0,01859
ROSN 0,02179 0,02222** 0,02096 0,02055* 0,01781 0,01760*
USDRUB 0,00576 0,00571* 0,00585 0,00598** 0,00587 0,00579*
VTBR 0,01158 0,01127* 0,01128 0,01165** 0,02233 0,02172*
IMOEX 0,01212 0,01218** 0,01206 0,01202* 0,01380 0,01345*
LKOH 0,01679** 0,01675 0,01652* 0,01675 0,01732 0,01679*
SBER 0,01512 0,01455* 0,01524** 0,01489 0,01745 0,01707*
Среднее 0,01365 0,01365** 0,01340* 0,01351 0,01522 0,01511*
Таблица 6
Ус
эедненные показатели MAE при глубине прогноза
Тикер ARIMA ARFIMA (D) ARFIMA (D„-log) ARFIMA (H) ARIMA-GARCH ARFIMA-GARCH
ALRS 0,01443** 0,01401 0,01411 0,01383* 0,01765* 0,01793
EURRUB 0,00639* 0,00646 0,00647 0,00648** 0,00483* 0,00484
GMKN 0,01531** 0,01514 0,01495* 0,01527 0,01564* 0,01595
MTSS 0,01297 0,01310** 0,01275* 0,01308 0,01563* 0,01566
ROSN 0,01817 0,01818** 0,01771 0,01752* 0,01730 0,01711*
USDRUB 0,00544 0,00549 0,00542* 0,00557** 0,00590 0,00582*
VTBR 0,01261 0,01242* 0,01256 0,01268** 0,02834 0,02810*
IMOEX 0,01224 0,01225** 0,01217* 0,01219 0,01453 0,01430*
LKOH 0,01683** 0,01678 0,01678* 0,01679 0,01763 0,01733*
SBER 0,01639 0,01590* 0,01639** 0,01630 0,01850 0,01804*
Среднее 0,01308** 0,01297 0,01293* 0,01297 0,01559 0,01551*
0 шагов вперед
По результатам проведенного исследования можно сделать следующие выводы.
1) Основная гипотеза исследования о лучшей прогностической способности фрактальных эконометрических моделей в сравнении с классическими не отвергается. Во всех, за исключением одного, случаях различные вариации АКБИМА оказались точнее, чем АК1МА. В большинстве случаев АКИГМА-ОАКСИ также оказались точнее, чем БОАКСИ.
2) Гипотеза о переменной фрактальной структуре финансовых временных рядов подтверждается, что доказывают изменчивые графики локальных размерностей (амплитуда колебаний составляет 0,40,5). Это закономерно и объясняется неоднородностью рыночной динамики курсов финансовых инструментов, ее хаотичностью, наличием в рядах участков как с ярко выраженным трендом, так и со случайным блужданием, наличием эффекта кластеризации волатильности. В то же время возникает задача кластеризации участков ряда с различной фрактальной размерностью. Методики ее решения в литературе обнаружено не было, однако успешная работа в этом направлении позволила бы определять, в каком состоянии процесс находится в настоящий момент, что расширило бы возможности для предсказания его будущего поведения. Кроме того, знание таких участков позволит выявить условия применимости различных моделей.
3) Гипотеза о том, что использование локальных показателей фрактальности вместо глобальных для вычисления пара-
метра дробного дифференцирования увеличивает точность прогнозов. не подтверждается. В среднем модель с глобально оцененным параметром d показала незначительно меньшую точность. чем с локальными. Предположительно это можно объяснить следующими ограничениями примененного инструментария. Во-первых. моделируемые локальные участки выбирались грубо. методом «скользящего окна» фиксированного размера. Разумно предположить. что несколько раз это окно совпало с границами истинного участка с концентрированным однотипным поведением процесса и оцененная величина его размерности также являлась истинной. Не исключено. что среди фрактальных моделей. аппроксимированных на таких участках. наибольшей точностью прогнозов обладали те. параметр d которых определялся по локальной размерности. Однако проверить это невозможно. так как ошибка прогнозов усреднялась и итоговая анализируемая величина MAE учитывает как прогнозы. полученные на истинных участках. так и прогнозы на участках. содержащих границы состояния процесса (например. тренда и случайного блуждания). а оцененная на таких «смешанных» участках величина локальной размерности адекватно описывает этот участок не в большей степени. чем величина глобальной. Во-вторых. для оценки локальной и глобальной размерностей и. следовательно. параметра d моделей были использованы разные методы (соответственно метод минимального покрытия и R/S-анализ). Данные ограничения методики
создают предпосылки для дальнейших исследований и проверок гипотезы о преимуществе применения локальных размерностей для целей прогнозирования.
4) Гипотеза о том, что учет несуществующего эффекта долгой памяти не ухудшает прогнозов, а его игнорирование существующего - ухудшает, не отвергается. Модель с короткой памятью АК1МА показала самые низкие результаты среди авторегрессионных моделей почти в каждом случае. Чтобы получить более убедительные доказательства данного утверждения, следует провести исследование с акцентом на сравнении фрактальных и нефрактальных моделей на выявленных участках с долгой памятью, а также на участках без нее. Потенциальным результатом такого исследования может стать рекомендация всегда использовать модели АКБИМА вместо АК1МА.
5) Гипотеза о том, что при росте длины прогноза его точность снижается, не подтверждается. На горизонтах прогнозирования 1, 3, 10 шагов точность остается сопоставимой. По нашему мнению, причины несогласованности результатов с этой достаточно очевидной гипотезой (проще предсказать одно ближайшее значение, чем серию, вероятность реализации такого прогноза выше) - в следующем ограничении методики: для оценки точности прогноза использовалась усредненная абсолютная величина ошибки, не учитывающая направления ошибок единичных. Возможно, с ростом глубины прогноза модели ошибаются чаще именно в направлениях динамики курсов финансовых инструментов. Поэтому данный момент требует дальнейших исследований.
6) Гипотеза об изменении фрактальных свойств при переходе к лог-доходностям подтверждается, о чем свидетельствуют различные значения показателя Херста, а также графики локальных размерностей до и после этой операции. Данный факт объясняется тем, что лог-трансформация, как и дифференцирование, устраняет в ряду тренд, поведение ряда начинает напоминать паттерн «возврат к среднему».
7) В сравнении АКИМА и ОАКСИ-моделей превосходство в точности прогнозирования осталось за АК1МА. Причиной этого может быть то, что исследуемые ценовые ряды не обладают высокой вола-тильностью, для них нехарактерен эффект ее кластеризации, поэтому ОАКСИ-модели не демонстрируют преимуществ. Однако исследование этих моментов (в частности, тестов на условную гетероскедастичность в остатках АКМА-моделей) не проводилось.
8) В среднем небольшое преимущество перед другими получила модель АКБИМА с параметром ё, вычисленным по антиперсистентному ряду лог-доходностей, а не по трендовому исходному ряду котировок. Этот факт не укладывается в рамки гипотезы о том, что участки с долгой памятью лучше моделируются [21]. Предполагаемым объяснением данного факта может послужить то, что, как подтвердили результаты исследования, при лог-трансформации ряда меняются его фрактальные свойства, в том числе и величина размерности, из которой выводится значение ё. Логичным выглядит то, что учитываться при моделировании должны фрактальные свойства непосредственно моделируемого ряда. Данный ряд имеет преимущественно антиперси-стентный характер, практически не обладает долгой памятью, и вопрос о том, насколько вообще к нему применимы модели, ее учитывающие, остается открытым. Этот момент требует дальнейших исследований.
Таким образом, можно констатировать, что полученные в результате исследования факты согласуются с четырьмя из проверяемых гипотез и противоречат двум из них. При этом причиной того, что эти гипотезы не подтвердились, предположительно стали ограничения методики, поэтому отбрасывать их нельзя, необходимы дальнейшие исследования в данном направлении.
Заключение
Основная выдвинутая авторами гипотеза о том, что при прогнозировании финансовых временных рядов фрактальные модели позволяют получать более точные прогно-
зы, чем нефрактальные, не противоречит данным исследования. Также было подтверждено, что:
1) Финансовые ряды имеют изменчивую фрактальную структуру, величина их локальной фрактальной размерности зависит от границ рассматриваемого участка.
2) Модели с долгой памятью показывают устойчивое преимущество перед моделями с короткой памятью на подавляющем большинстве данных.
3) При трансформации исходного ряда курсов финансовых инструментов его фрактальные свойства меняются.
Все эти факты согласуются с положениями теории о фрактальности финансовых рынков и результатами других исследований. Новизной исследования является то, что авторами впервые получено подтверждение этих гипотез на большом массиве данных, а также поставлена и решена проблема кластеризации участков финансовых рядов с различными фрактальными свойствами, определено, по какому ряду следует оценивать фрактальную размерность для вычисления оператора дробного дифференцирования модели ARFIMA - исходному или преобразованному (моделируемому). Это позволяет лучше определять условия применимости фрактальных эконометрических моделей и в целом развивает методологию прогнозирования динамики финансовых временных рядов, увеличивая точность прогнозов.
В ходе исследования не удалось подтвердить следующие гипотезы:
1) Применение локальных показателей фрактальности вместо глобальных
для вычисления параметра дробного дифференцирования модели АКРГМА увеличивает точность прогнозов.
2) С увеличением глубины прогноза его точность снижается.
3) Участки рядов с долгой памятью (тренды) лучше моделируются и прогнозируются. Авторами не ставилась задача проверки данной гипотезы, однако некоторые исследователи [21] ее выдвигали.
Предположительно причиной несоответствия установленных фактов этим гипотезам явились ограничения примененного авторами методического инструментария, поэтому результаты настоящего исследования не являются достаточным основанием для того, чтобы данные гипотезы отвергнуть как ошибочные.
На основании анализа полученных результатов перспективными направлениями исследований являются: разработка метода определения границ локальных участков финансовых рядов с разными фрактальными свойствами; формулирование рекомендаций относительно того, по какому ряду (исходному или преобразованному) следует оценивать оператор дробного дифференцирования модели АКИМА; проверка гипотезы о преимуществе применения локальных размерностей для целей прогнозирования, влиянии длины прогноза на его точность и о лучшей прогнозируемости рядов финансовых инструментов на их трендовых участках, исследование различных фрактальных модификаций GARCH-моделей ^ЮА^И, FIAPARCH и др.) и выявление условий их применимости.
Список литературы
1. Симонов П.М., Филимонова С.А. P-адическое моделирование динамики индекса РТС в зависимости от таймфреймов // Вестник Пермского университета. Сер. «Экономика» = Perm University Herald. Economy. 2016. № 4 (31). С. 74-85. doi: 10.17072/1994-9960-2016-4-74-85.
2. Кривоносова Е.К. Разработка методов прогнозирования и анализа кредитных и инвестиционных рисков с применением фрактальных и мультифрактальных характеристик: дисс. ... канд. экон. наук: 08.00.13. Пермь, 2015. 167 с.
3. Зиненко А.В. R/S анализ на фондовом рынке // Бизнес-информатика. 2012. № 3 (21). С. 24-30.
4. Fama E.F. Efficient capital markets: A review of theory and empirical // The Journal of Finance. 1969. Vol. 25, № 2. P. 383-417.
5. Карасева Е.И. Преимущества междисциплинарных исследований в экономике // Известия высших учебных заведений. Поволжский регион. Общественные науки. 2014. № 1 (29). С. 210-227.
6. Романовский М.Ю., Романовский Ю.М. Введение в эконофизику: статистические и динамические модели. Изд. 2-е. испр. и доп. М. Ижевск: Ин-т компьютерных исследований. 2012. 340 с.
7. Гарафутдинов Р.В. Применение моделей с долгой памятью для прогнозирования динамики фондового индекса // Математика и междисциплинарные исследования - 2018. Материалы Всеросс. науч.-практ. конф. молодых ученых с международным участием. 14-19 мая 2018 г. Пермь: Пермский государственный национальный исследовательский университет. 2018. С. 158-161.
8. Старченко Н.В. Индекс фрактальности и локальный анализ хаотических временных рядов: дисс. ... канд. физ.-мат. наук: 05.13.18. 01.01.03. М.. 2005. 122 с.
9. Петерс Э. Фрактальный анализ финансовых рынков: применение теории хаоса в инвестициях и экономике. М.: Интернет-трейдинг. 2004. 304 с.
10. Прудский М.В. Фрактальный анализ финансовых рынков // Информационные системы и математические методы в экономике: сб. науч. тр. / общ. ред М.В. Радионовой. Пермь: Перм. гос. нац. исслед. ун-т. 2012. Вып. 5. С. 109-120.
11. Балагула Ю.М., Абакумова Ю.А. Длинная память на рынке нефти: спектральный подход. Препринт Ес-01/11. СПб.: ЕУСПб. 2011. 40 с.
12. Остапенко Е.С., Дунаева Т.А. Прогнозирование временных рядов с долговременной памятью с помощью моделей класса ARFIMA // Вюник Кшвського нацюнального ушверситету технологш та дизайну. 2010. № 7. C. 270-273.
13. Белолипцев И.И., Фархиева С.А. Предсказание финансовых временных рядов на основе индекса фрактальности // Мир Науки. 2014. Вып. 3. URL: https://mir-nauki.com/PDF/01EMN314.pdf (дата обращения: 10.11.2018).
14. Балагула Ю.М. Фрактальные характеристики длинной памяти в ценах на электроэнергию. Препринт Ec-03/16. СПб.: ЕУСПб. 2016. 13 с.
15. Caporale G.M., Skare M. Long memory in UK real GDP, 1851-2013: An ARFIMA-FIGARCH analysis // DIW Berlin Discussion Paper. 2014. № 1395. doi: 10.2139/ssrn.2459806.
16. Zhelyazkova S. ARFIMA-FIGARCH, HYGARCH and FIAPARCH Models of Exchange Rates // Izvestia Journal of the Union of Scientists - Varna. Economic Sciences Series. 2018. Vol. 7 (2). P. 142-153.
17. ДубовиковМ.М., Старченко Н.В. Эконофизика и фрактальный анализ финансовых временных рядов // Успехи физических наук. 2011. Т. 181. № 7. С. 779-786.
18. Кривоносова Е.К., Первадчук В.П. Использование фрактального подхода для анализа стабильности многоуровневых структур // Вестник ПНИПУ. Машиностроение. материаловедение. 2013. № 1. С. 63-69.
19. Кривоносова Е.К., Первадчук В.П., Кривоносова Е.А. Сравнение фрактальных характеристик временных рядов экономических показателей // Современные проблемы науки и образования (электронный научный журнал). 2014. № 6. URL: https://www.science-education.ru/ru/article/view?id=15974 (дата обращения: 10.01.2019).
20. Мансуров А.К. Прогнозирование валютных кризисов с помощью методов фрактального анализа // Проблемы прогнозирования. 2008. № 1 (106). С. 145-158.
21. Kulish V, Horak V. Forecasting the behavior of fractal time series: Hurst exponent as a measure of predictability // Review of the Air Force Academy. 2016. № 2 (32). P. 61-68. doi: 10.19062/18429238.2016.14.2.8.
22. Мандельброт Б. Фрактальная геометрия природы. М.: Ин-т компьютерных исследований. 2002. 656 с.
23. Box G., Jenkins G., Reinsel G. Time series analysis: Forecasting and control. 4th ed. N. Y.: Wiley, 2008. 784 p.
24. Hosking J. Fractional differencing // Biometrika. 1981. Vol. 68. № 1. P. 165-176.
25. Haubrich J.G. Consumption and fractional differencing: Old and new anomalies // The Review of Economics and Statistics. 1993. Vol. 75. № 4. P. 767-772. doi: 10.2307/2110038.
26. Bollerslev T. Generalized autoregressive conditional heteroscedasticity // Journal of Econometrics. 1986. Vol. 31. P. 307-327.
27. Аганин А.Д. Сравнение GARCH и HAR-RV моделей для прогноза реализованной волатильности на российском рынке // Прикладная эконометрика. 2017. Т. 48. С. 63-84.
28. Ахуньянова С.А., Симонов П.М. Моделирование и прогнозирование на финансовых рынках с помощью эконометрики и эконофизики: монография. Пермь: Перм. гос. нац. исслед. ун-т, 2017. 203 с. URL: https://elis.psu.ru/node/486405 (дата обращения: 06.04.2019).
29. Ling S., Li W. On fractionally integrated autoregressive moving-average time series models with conditional heteroscedasticity // Journal of the American Statistical Association. 1997. Vol. 92, № 439. P. 1184-1194.
30. СизовА.А. Модели, способы и программные средства поддержки принятия решений на основе прогнозирования временных рядов с переменной структурой: диса ... канд. техн. наук: 05.13.17. Смоленск, 2014. 139 с.
31. Chou J.S., NgoN.T. Time series analytics using sliding window metaheuristic optimization-based machine learning system for identifying building energy consumption patterns // Applied Energy. 2016. Vol. 177. P. 751-770.
32. Jorda O., Taylor A.T. The time for austerity: Estimating the average treatment effect of fiscal policy // The Economic Journal. 2016. Vol. 126, Iss. 590. P. 219-255. doi: 10.1111/ecoj.12332.
Статья поступила в редакцию 21.01.2019, принята к печати 24.04.2019
Сведения об авторах
Симонов Петр Михайлович - доктор физико-математических наук, профессор, профессор кафедры информационных систем и математических методов в экономике, Пермский государственный национальный исследовательский университет (Россия, 614990, г. Пермь, ул. Букирева, 15; e-mail: [email protected]).
Гарафутдинов Роберт Викторович - аспирант кафедры информационных систем и математических методов в экономике, Пермский государственный национальный исследовательский университет (Россия, 614990, г. Пермь, ул. Букирева, 15; e-mail: [email protected]).
References
1. Simonov P.M., Filimonova S.A. P-adicheskoe modelirovanie dinamiki indeksa RTS v zavisimosti ot taimfreimov [P-adic modeling of the RTS index dynamics depending on the timeframes]. Vestnik Permskogo universiteta. Seriya "Ekonomika" [Perm University Herald. ECONOMY], 2016, no. 4 (31), pp. 74-85. (In Russian). doi: 10.17072/1994-9960-2016-4-74-85.
2. Krivonosova E.K. Razrabotka metodov prognozirovaniya i analiza kreditnykh investitsionnykh riskov s primeneniem fraktal'nykh i mul'tifraktal'nykh kharakteristik. Diss. kand. ekon. nauk [Development of methods for forecasting and analyzing credit and investment risks using fractal and multifractal characteristics. Cand. econ. sci. diss.]. Perm, 2015. 167 p. (In Russian).
3. Zinenko A.V. R/S analiz na fondovom ryinke [R/S analysis of stock market]. Biznes-informatika [Business Informatics], 2012, no. 3 (21), pp. 24-30. (In Russian).
4. Fama E.F. Efficient capital markets: A review of theory and empirical. The Journal of Finance, 1969, vol. 25, no. 2, pp. 383-417.
5. Karaseva E.I. Preimushchestva mezhdistsiplinarnykh issledovanii v ekonomike [Advantages of interdisciplinary research in economics]. Izvestiya vysshikh uchebnykh zavedenii. Povolzhskii region. Obshchestvennye nauki [University Proceedings. Volga region. Social Sciences], 2014, no. 1 (29), pp. 210-227. (In Russian).
6. Romanovskii M.Yu., Romanovskii Yu.M. Vvedenie v ekonofiziku: statisticheskie i dinamicheskie modeli. Izd. 2-e, ispr. i dop. [Introduction to econophysics: statistical and dynamic models. 2nd ed., rev. and add.]. Moscow, Izhevsk, Institut komp'yuternykh issledovanii Publ., 2012. 340 p. (In Russian).
7. Garafutdinov R.V. Primenenie modelei s dolgoi pamyat'yu dlya prognozirovaniya dinamiki fondovogo indeksa [Stock exchange index forecasting using long memory models]. Matematika i mezhdistsiplinarnye issledovaniya - 2018. Materialy Vserossiiskoi nauchno-prakticheskoi konferentsii molodykh uchenykh s mezhdunarodnym uchastiem, 14-19 maya 2018 g. [Mathematics and interdisciplinary research - 2018. Proceedings of Russian Scientific and Practical Conference of Young
Scientists with International Participation, May 14-19, 2018]. Perm, Permskii gosudarstvennyi natsional'nyi issledovatel'skii universitet Publ., 2018, pp. 158-161. (In Russian).
8. Starchenko N.V. Indeks fraktal'nosti i lokal'nyi analiz khaoticheskikh vremennykh ryadov. Diss. kand. fiz.-mat. nauk [Fractality index and local analysis of chaotic time series. Cand. phys.-math. sci diss.]. Moscow, 2005. 122 p. (In Russian).
9. Peters E. Fraktal'nyi analiz finansovykh rynkov: primenenie teorii khaosa v investitsiyakh i ekonomike [Fractal analysis of financial markets: Application of chaos theory in investments and the economy]. Moscow, Internet-treiding, 2004. 304 p. (In Russian).
10. Prudskii M.V. Fraktal'nyi analiz finansovykh rynkov [Fractal analysis of financial markets]. Informatsionnye sistemy i matematicheskie metody v ekonomike. Sb. nauch. tr. Obshch. red. M.V. Radionovoi. Vyp. 5 [Information systems and mathematical methods in economics. Coll. scientific works. Ed. by M.V. Radionova. Iss. 5]. Perm, Perm. gos. nats. issl. un-t Publ., 2012, pp. 109-120. (In Russian).
11. Balagula Yu.M., Abakumova Yu.A. Dlinnaya pamyat' na rynke nefti: spektral'nyi podkhod. Preprint Es-01/11 [Long memory in the oil market: A spectral approach. Preprint Es-01/11]. St. Petersburg, EUSPb Publ., 2011. 40 p. (In Russian).
12. Ostapenko E.S., Dunaeva T.A. Prognozirovanie vremennykh ryadov s dolgovremennoi pamyat'yu s pomoshch'yu modelei klassa ARFIMA [Forecasting time series with long memory using ARFIMA models]. Visnik Kiivs'kogo natsional'nogo universitetu tekhnologii ta dizainu [Bulletin of Kiev National University of Technology and Design], 2010, no. 7, pp. 270-273. (In Russian).
13. Beloliptsev I.I., Farkhieva S.A. Predskazanie finansovykh vremennykh ryadov na osnove indeksa fraktal'nosti [Forecasting financial time series based on the fractal index]. Mir nauki [World of Science], 2014, vol. 3. (In Russian) Available at: https://mir-nauki.com/PDF/01EMN314.pdf (accessed 10.11.2018).
14. Balagula Yu.M. Fraktal'nye kharakteristiki dlinnoi pamyati v tsenakh na elektroenergiyu. Preprint Es-03/16 [Fractal characteristics of long memory in electricity prices. Preprint Es-03/16]. St. Petersburg., EUSPb Publ., 2016. 13 p. (In Russian).
15. Caporale G.M., Skare M. Long memory in UK real GDP, 1851-2013: An ARFIMA-FIGARCH analysis. DIWBerlin Discussion Paper, 2014, no. 1395. doi: 10.2139/ssrn.2459806.
16. Zhelyazkova S. ARFIMA-FIGARCH, HYGARCH and FIAPARCH Models of Exchange Rates. Izvestia Journal of the Union of Scientists - Varna. Economic Sciences Series, 2018, vol. 7 (2), pp.142-153.
17. Dubovikov M.M., Starchenko N.V. Ekonofizika i fraktal'nyi analiz finansovykh vremennykh ryadov [Econophysics and fractal analysis of financial time series]. Uspekhi fizicheskikh nauk [Advances in Physical Sciences], 2011, vol. 181, no. 7, pp. 779-786. (In Russian).
18. Krivonosova E.K., Pervadchuk V.P. Ispol'zovanie fraktal'nogo podkhoda dlya analiza stabil'nosti mnogourovnevykh struktur [Fractal approach to the multilevel structure stability analyze]. VestnikPNIPU. Mashinostroenie, materialovedenie [Bulletin PNRPU. Mechanical Engineering, Materials Science], 2013, no. 1, pp. 63-69. (In Russian).
19. Krivonosova E.K., Pervadchuk V.P., Krivonosova E.A. Sravnenie fraktal'nykh kharakteristik vremennykh ryadov ekonomicheskikh pokazatelei [Comparison of the fractal characteristics of economic indicators time series]. Sovremennye problemy nauki i obrazovaniya (elektronnyi nauchnyi zhurnal). [Modern Problems of Science and Education (scientific journal)], 2014, no. 6. (In Russian) Available at: https://www.science-education.ru/ru/article/view?id=15974 (accessed 10.01.2019).
20. Mansurov A.K. Prognozirovanie valyutnykh krizisov s pomoshch'yu metodov fraktal'nogo analiza [Forecasting currency crises by fractal analysis techniques]. Problemy prognozirovaniya [Forecasting Problems], 2008, no. 1 (106), pp. 145-158. (In Russian).
21. Kulish V, Horak V. Forecasting the behavior of fractal time series: Hurst exponent as a measure of predictability. Review of the Air Force Academy, 2016, no. 2 (32), pp. 61-68. doi: 10.19062/18429238.2016.14.2.8.
22. Mandel'brot B. Fraktal'naya geometriyaprirody [Fractal geometry of nature]. Moscow, Institut komp'yuternykh issledovanii Publ., 2002. 656 p. (In Russian).
23. Box G., Jenkins G., Reinsel G. Time series analysis: Forecasting and control. 4th Ed. New York, Wiley, 2008. 784 p.
24. Hosking J. Fractional differencing. Biometrika, 1981, vol. 68, no. 1, pp. 165-176.
25. Haubrich J.G. Consumption and fractional differencing: Old and new anomalies. The Review of Economics and Statistics, 1993, vol. 75, no. 4, pp. 767-772. doi: 10.2307/2110038.
26. Bollerslev T. Generalized autoregressive conditional heteroscedasticity. Journal of Econometrics, 1986, vol. 31, pp. 307-327.
27. Aganin A.D. Sravnenie GARCH i HAR-RV modelei dlya prognoza realizovannoi volatil'nosti na rossiiskom rynke [Forecast comparison of volatility models on Russian stock market]. Prikladnaya ekonometrika [Applied Econometrics], 2017, vol. 48, pp. 63-84. (In Russian).
28. Akhun'yanova S.A., Simonov P.M. Modelirovanie i prognozirovanie na finansovykh rynkakh s pomoshch'yu ekonometriki i ekonofiziki: monografiya. [Modeling and prediction of financial markets by means of econometrics and econophysics: Monograph]. Perm, Perm. gos. nats. issled. un-t Publ., 2017. 203 p. (In Russian) Available at: https://elis.psu.ru/node/486405 (accessed 10.01.2019).
29. Ling S., Li W. On fractionally integrated autoregressive moving-average time series models with conditional heteroscedasticity. Journal of the American Statistical Association, 1997, vol. 92, no. 439, pp. 1184-1194.
30. Sizov A.A. Modeli, sposoby i programmnye sredstva podderzhki prinyatiya reshenii na osnove prognozirovaniya vremennykh ryadov s peremennoi strukturoi. Diss. kand. tekhn. nauk [Models, methods and software for decision support based on the prediction of time series with a variable structure. Cand. tech. sci. diss.]. Smolensk, 2014. 139 p. (In Russian).
31. Chou J.S., Ngo N.T. Time series analytics using sliding window metaheuristic optimization-based machine learning system for identifying building energy consumption patterns. Applied Energy, 2016, vol. 177, pp. 751-770.
32. Jorda O., Taylor A.M. The time for austerity: Estimating the average treatment effect of fiscal policy. The Economic Journal, 2016, vol. 126, iss. 590, pp. 219-255. doi: 10.1111/ecoj.12332.
Simonov Pyotr Mikhailovich - Doctor of Physical and Mathematical Sciences, Professor, Professor at the Department of Information Systems and Mathematical Methods in Economics, Perm State University (15, Bukireva st., Perm, 614990, Russia; e-mail: [email protected]).
Garafutdinov Robert Viktorovich - Postgraduate Student at the Department of Information Systems and Mathematical Methods in Economics, Perm State University (15, Bukireva st., Perm, 614990, Russia; e -mail : rvgarafutdinov@gmail .com).
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:
Симонов П.М., Гарафутдинов Р.В. Моделирование и прогнозирование динамики курсов финансовых инструментов с применением эконометрических моделей и фрактального анализа // Вестник Пермского университета. Сер. «Экономика» = Perm University Herald. Economy. 2019. Том 14. № 2. С. 268-288. doi: 10.17072/1994-9960-2019-2-268-288
Please cite this article in English as:
Simonov P.M., Garafutdinov R.V. Modeling and forecasting of financial instruments dynamics using econometrics models and fractal analysis. Vestnik Permskogo universiteta. Seria Ekonomika = Perm
University Herald. Economy, 2019, vol. 14, no. 2, pp. 268-288. doi: 10.17072/1994-9960-2019-2-268-288
Received January 21, 2019; accepted April 24, 2019
Information about the Authors
