Экспонента Херста рассчитывается с помощью следующих формул:
где H - экспонента Херста; n - длина временного ряда- размах временного ряда; S - стандартное отклонение;
Кроме того, найдя экспоненту Херста, можно вычислить предполагаемое
значение дробной разности модели ARFIMA по формуле
Таблица 13
Значения экспоненты Херста для кластеров остатков денежных средств
|
Кластер |
Экспонента Херста |
Предполагаемое значение разности d |
|
Черный |
0,94 |
0,44 |
|
Красный |
0,97 |
0,47 |
|
Зеленый |
0,97 |
0,47 |
|
Синий |
0,94 |
0,44 |
Можно видеть, что значения экспоненты имеет схожие значения для всех
кластеров. Значение экспоненты больше 0,5, что говорит о персистентности
процесса. Также данные ряды называются трендоустойчивыми, то есть, если в данных
на определенном промежутке времени наблюдался восходящий или нисходящий тренд,
то данные ряды будут сохранять этот тренд еще какое-то время. Наблюдения ряда
не являются независимыми, они сохраняют данные о том, что было на предыдущих
интервалах. Также можно предположить, что расчетное значение дробной разности в
модели будет близко к значениям, полученным с помощью экспоненты Херста.
Экспонента Херста подтверждает возможность использования модели ARFIMA. Для
всех кластеров была оценена модель ARFIMA. Результаты оценивания приведены в
табл. 14.
Таблица 14
Результаты оценивания модели ARFIMA для четырех кластеров
|
|
Черный |
Красный |
Зеленый |
Синий |
|
Разность d |
0,49*** |
0,29*** |
0,41*** |
0,27*** |
|
Константа |
154860*** |
234480*** |
110162*** |
14721*** |
|
AR(1) |
-0,22*** |
-0,28*** |
-0,89*** |
-0,23*** |
|
MA(1) |
0,67*** |
0,65*** |
0,92*** |
0,27*** |
*** - Значимость на уровне 0,001
Недостатком данной модели является сложность вычисления оценок параметров: в модель невозможно добавить сезонный фактор, вычисления не производятся для больших лагов. Поэтому в модель были включены лаги AR(1) и MA(1), так как коррелограмма показала значимость на первом лаге. Все полученные оценки модели значимы на уровне 0,001. Полученная разность модели наиболее близка к расчетной для зеленого и черного кластеров. Несмотря на схожесть поведения рядов, оценки имеют существенно различающиеся значения. Оценка AR(1) значительно меньше у зеленого кластера, в то время как MA(1) наоборот больше, чем у остальных. Константа отражает среднее значение ряда.
Еще одним недостатком данной модели является сложность прогнозирования
будущих наблюдений, а также краткосрочный интервал прогнозирования. Для
упрощения вычислений модель ARFIMA может аппроксимироваться с помощью формулы:
Расчет страховой суммы представлен на примере красного кластера. Были
предсказаны значения остатков денежных средств в точках в 2017 году.
Рис. 31. Прогноз остатков денежных средств модели ARFIMA в торговых точках красного кластера на 2017 год
Наибольшая страховая сумма была спрогнозирована для красного кластера:
среднее значение ряда равно 252 тысячам рублей, а верхняя граница
доверительного интервала не превосходит 400 тысяч. Страховые суммы для зеленого
и синего кластера отличаются незначительно. Во всех кластерах, кроме красного,
наблюдается небольшой спад остатков денежных средств.
Таблица 15
Характеристики прогноза ARFIMA на 2017 год для четырех кластеров
|
Кластер |
Красный |
Черный |
Зеленый |
Синий |
|
Среднее значение ряда |
252 005 |
116 421 |
167 607 |
185 137 |
|
Верхняя граница 95% интервала страхового покрытия |
371 318 |
236 276 |
197 384 |
268 970 |
|
Динамика остатков денежных средств |
1% |
-1% |
-3% |
-1% |
Итак, нами была построена адекватная модель для прогнозирования остатков денежных средств без выбросов. Оптимальной моделью для прогнозирования является модель ARFIMA. Были получены значимые оценки модели ARFIMA для лагов AR(1) и MA(1), также в модель вошла константа, отражающая среднее значение временного ряда. Был построен прогноз для остатков денежных средств на 2017 год и рассчитаны страховые суммы, требуемые для страхования исследуемых рисков.
В данном параграфе моделируются временные ряды остатков денежных средств без исключения выбросов с помощью модели ARIMA с интервенцией. Данная модель учитывает наиболее высокие непостоянные скачки, резко меняющие среднее значение. Рассмотрим графическое изображение временных рядов без исключения выбросов.
Рис. 32. Остатки денежных средств с выбросами в торговых точках в четырех
кластерах
Можно видеть, что стабильное поведение ряда сменяется резкими скачками. Данные скачки могут сильно влиять на оценку моделей, поэтому важно экзогенно ввести в модель даты, в которые интервенции (выбросы) происходят, и учесть их влияние. Также необходимо понимать природу интервенций: влияют ли они на общий тренд модели или происходят и исчезают внезапно. В случае изучаемых данных, изменения временные и на общее поведение ряда не влияют. Это объясняется особенностью исследуемой предметной области - поведением покупателей электроники. Большая часть продаж приходится на праздничные дни и выходные, после чего спрос возвращается к обычным значениям.
Для построения модели необходимо знать номера наблюдений, в которых
происходят скачки. Наиболее простым способом обнаружения интервенций является
графический анализ. Как правило, для данных остатков денежных средств
интервенции в несколько раз превышают другие значения. Далее оценки модели
подтвердят или опровергнут значимость интервенции. Можно видеть (табл. 16), что
некоторые значения интервенций в кластерах повторяются.
Таблица 16
Номера дат (наблюдений) предполагаемых интервенций в четырех кластерах
|
|
Черный |
Красный |
Зеленый |
Синий |
|
Интервенция 1 |
121 |
480 |
360 |
121 |
|
Интервенция 2 |
480 |
700 |
519 |
359 |
|
Интервенция 3 |
828 |
828 |
952 |
707 |
|
Интервенция 4 |
|
|
986 |
828 |
Для определения лагов в модели ARIMA с интервенцией применяются те же
методы, что и для обычной ARIMA модели. В предыдущем параграфе были рассмотрены
характеристики данных без выбросов и выявлена недельная сезонность. Поэтому в
модель ARIMA с интервенцией будет включен сезонный AR(1)7 лаг. Также
были подобраны обычные AR и МА лаги путем исключения незначимых значений.
Результат оценки ARIMA с интервенцией показан на примере зеленого кластера
(табл. 17), остальные оценки приведены в приложении (прил.1).
Таблица 17
|
|
Оценка параметра |
Значимость |
Номер интервенции |
|
AR(1) |
0,8 |
0,00 |
|
|
AR(1)7 |
0,5 |
0,00 |
|
|
Omega(1) |
96311 |
0,00 |
360 |
|
Delta(1) |
1,0 |
0,00 |
360 |
|
Omega(2) |
60420 |
0,00 |
519 |
|
Delta(2) |
1,0 |
0,00 |
519 |
|
Omega(3) |
54253 |
0,00 |
952 |
|
Delta(3) |
-0,3 |
0,13 |
952 |
|
Omega(4) |
113358 |
0,00 |
986 |
|
Delta(4) |
0,1 |
0,55 |
986 |
Параметр Delta отвечает за значимость интервенции. Так как в данной
модели интервенции задаются с помощью бинарных фиктивных переменных, значимые
значения Delta должны быть близки к единице. Показатель Omega отражает
коэффициент при бинарной переменной. Для зеленого кластера наблюдения 952 и 986
не являются интервенциями, несмотря на то, что оценки их параметров значимы.
Наибольшее влияние оказывает интервенция с номером 360. Также в модели значимы
AR лаги, как сезонный, так и обычный.
Рис. 33. Прогноз остатков денежных средств зеленого кластера в 2017 году
Прогнозирование для данной модели имеет те же недостатки, что и
прогнозирование модели ARFIMA. Основным результатом прогнозирования является
95% доверительный интервал для остатков денежных средств в 2017 году (табл.
18).
Таблица 18
Характеристики прогноза на 2017 год модели ARIMA с интервенцией по четырем кластерам
|
|
Красный |
Черный |
Зеленый |
Синий |
|
Среднее значение ряда |
110 749 |
198 159 |
133 712 |
86 753 |
|
Верхняя граница 95% интервала страхового покрытия |
407 281 |
283 068 |
222 267 |
291 512 |
|
Динамика остатков денежных средств |
-82% |
8% |
11% |
-76% |
Согласно данной модели, в красном и синем кластере в 2017 году произойдет большой спад. Также в этих кластерах верхняя граница доверительного интервала почти в 4 раза больше среднего значения. Черный и зеленый кластеры спрогнозированы более достоверно: в них наблюдается небольшой рост, а страховые суммы равны 198 и 133 тысячи соответственно.
Выводы по главе 3
Точки были разбиты на четыре кластера методом k-средних. При кластеризации учитывались не только значения ряда и его первой разности, но и характеристики выбросов, а именно даты, в которые они происходили. Значимость различий между средними хотя бы двух кластеров подтверждена дисперсионным анализом.
Для всех кластеров, кроме черного, страховые суммы, рассчитанные моделью ARFIMA, превышают значения прогноза с интервенцией. Это может быть объяснено тем, что для модели ARFIMA были удалены лишь некоторые выбросы и влияние на среднее значение не было полностью исключено. Однако прогнозы модели ARFIMA кажутся более достоверными, так как модель с интервенцией предсказывает значительный рост и спад остатков денежных средств, хотя до 2017 года тренды не наблюдаются.
Таким образом, в данном параграфе были получены значимые оценки модели ARIMA с интервенцией. Данная модель определяет значимость скачков остатков денежных средств и их влияние на оценки модели. Кроме интервенций в модель вошли лаги AR(1) и AR(1)7. Страховые суммы, полученные с помощью модели c интервенцией, ниже оценок модели ARFIMA и показывают значительный спад и рост остатков денежных средств.
Обе модели могут быть использованы для прогнозирования страховых сумм
остатков денежных средств в точках. Модель ARFIMA хорошо прогнозирует временные ряды с малым числом
выбросов. В данном случае ритейлер может исключать выбросы и отдельно
страховать денежные средства в праздничные дни с более высокой суммой. Это
поможет снизить общий уровень страховых сумм для обычных дней. Модель ARIMA с интервенцией позволяет учесть
влияние выбросов на среднее значение ряда. Она проверяет значимость экзогенно
значимых выбросов и более точно предсказывает страховые суммы.
Заключение
В данной работе рассчитываются страховые покрытия для запасов товаров и денежных средств в магазинах одного из лидеров на рынке продаж электроники. Ритейлер мобильной техники выступает в роли страхователя и вынужден прогнозировать страховые суммы для волатильных объектов страхования. Компания представляет собой ритейловую сеть, состоящую из множества торговых точек розничной продажи. Данный вид торговли появился еще 1990-х годах и продолжает развиваться и претерпевать изменения и по сей день. Текущая экономическая ситуация в России значительно повлияла на ритейловые компании: ослабление рубля привело к снижению покупательной активности населения, расходы людей на электронику резко упали, выручка ритейлеров снизилась. Кроме того, ритейл-продажи быстрыми темпами переходят в онлайн-сегмент, что негативно сказывается на продажах торговых точек. Это вынуждает ритейлеров разрабатывать стратегии по снижению издержек. Кроме того, ритейловые компании могут получить существенные потери вследствие различных несчастных случаев и противоправных действий. Одним из способов предотвращения убытков от грабежей и несчастных случаев является страхование запасов товаров и остатков денежных средств. Так как данные активы страхуются на год вперед, необходимо спрогнозировать страховую сумму, на основе которой будут рассчитываться страховые премии.
Для практического использования, обобщения информации и построения небольшого числа универсальных моделей торговые точки ритейловой компании были разбиты на несколько значимых кластеров. Визуально разбить точки на группы невозможно, поэтому были применены методы кластерного анализа - иерархический и метод k-средних. Для запасов товаров было выделено четыре однородных кластера по следующим переменным:
. Среднее значение логарифма запасов товаров в данной торговой точке
. Стандартное отклонение логарифма запасов товаров;
. Размах временного ряда логарифма запасов товаров (max-min);
. Максимальный размах между соседними наблюдениями;
. Среднее первой разности логарифма запасов товаров;
. Стандартное отклонение первой разности логарифма запасов товаров.
Дальнейший анализ был проведен на усредненных временных рядах. Исследование показало, что поведение рядов запасов товаров становится стационарным после взятия первой разности: тест на единичный корень Дики-Фуллера отвергает гипотезу о наличии единичных корней. Автокорреляционная функция показала значимость 3, 4 и 12 лагов, периодограмма не выявила наличия сезонности, поэтому путем исключения незначимых лагов была построена модель ARIMA для первой разности. В модель вошли лаги AR(12) и MA(12). Модель позволила спрогнозировать с заданной надежностью страховые суммы для страхования запасов товаров на 2017 год.