Мультипликативная модель временного ряда

Простейший подход к моделированию сезонных колебаний – это расчет значений сезонной компоненты методом скользящей средней и построение аддитивной или мультипликативной модели временного ряда.
Общий вид мультипликативной модели выглядит так:
Y = T x S x E
где T - трендовая компонента, S - сезонная компонента и E - случайная компонента.
Назначение. С помощью данного сервиса производится построение мультипликативной модели временного ряда.
Инструкция. Укажите количество данных (количество строк), нажмите Далее. На втором шаге выберите диапазон моделирования. Полученное решение сохраняется в файле MS Word и MS Excel (см. пример решения).
Количество строк (исходных данных)

Здесь будет отображаться решение.

см. также Мультипликативная индексная двухфакторная модель.

Алгоритм построения мультипликативной модели

Построение мультипликативной моделей сводится к расчету значений T, S и E для каждого уровня ряда.
Процесс построения модели включает в себя следующие шаги.
  1. Выравнивание исходного ряда методом скользящей средней.
  2. Расчет значений сезонной компоненты S.
  3. Устранение сезонной компоненты из исходных уровней ряда и получение выровненных данных (T x E).
  4. Аналитическое выравнивание уровней (T x E) с использованием полученного уравнения тренда.
  5. Расчет полученных по модели значений (T x E).
  6. Расчет абсолютных и/или относительных ошибок. Если полученные значения ошибок не содержат автокорреляции, ими можно заменить исходные уровни ряда и в дальнейшем использовать временной ряд ошибок E для анализа взаимосвязи исходного ряда и других временных рядов.

Пример. Построить аддитивную и мультипликативную модель временного ряда, характеризующую зависимость уровней ряда от времени.
Решение. Построение мультипликативной модели временного ряда.
Общий вид мультипликативной модели следующий:
Y = T x S x E
Эта модель предполагает, что каждый уровень временного ряда может быть представлен как сумма трендовой (T), сезонной (S) и случайной (E) компонент.
Рассчитаем компоненты мультипликативной модели временного ряда.
Шаг 1. Проведем выравнивание исходных уровней ряда методом скользящей средней. Для этого:
1.1. Найдем скользящие средние (гр. 3 таблицы). Полученные таким образом выровненные значения уже не содержат сезонной компоненты.
1.2. Приведем эти значения в соответствие с фактическими моментами времени, для чего найдем средние значения из двух последовательных скользящих средних – центрированные скользящие средние (гр. 4 табл.).

tytСкользящая средняяЦентрированная скользящая средняяОценка сезонной компоненты
1898---
27941183.25--
314411200.51191.881.21
416001313.512571.27
59671317.751315.630.74
612461270.751294.250.96
714581251.751261.251.16
814121205.51228.631.15
98911162.751184.130.75
1010611218.51190.630.89
111287---
121635---
Шаг 2. Найдем оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние (гр. 5 табл.). Эти оценки используются для расчета сезонной компоненты S. Для этого найдем средние за каждый период оценки сезонной компоненты Sj. Сезонные воздействия за период взаимопогашаются. В мультипликативной модели это выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна числу периодов в цикле. В нашем случае число периодов одного цикла равно 4.
Показатели1234
1--1.211.27
20.740.961.161.15
30.750.89--
Всего за период1.491.852.372.42
Средняя оценка сезонной компоненты0.740.931.181.21
Скорректированная сезонная компонента, Si0.730.911.161.19
Для данной модели имеем:
0.744 + 0.927 + 1.183 + 1.211 = 4.064
Корректирующий коэффициент: k=4/4.064 = 0.984
Рассчитываем скорректированные значения сезонной компоненты Si и заносим полученные данные в таблицу.
Шаг 3. Разделим каждый уровень исходного ряда на соответствующие значения сезонной компоненты. В результате получим величины T x E = Y/S (гр. 4 табл.), которые содержат только тенденцию и случайную компоненту.
Находим параметры уравнения методом наименьших квадратов.
Система уравнений МНК:
a0n + a1∑t = ∑y
a0∑t + a1∑t2 = ∑y•t
Для наших данных система уравнений имеет вид:
12a0 + 78a1 = 14659.84
78a0 + 650a1 = 96308.75
Из первого уравнения выражаем а0 и подставим во второе уравнение
Получаем a1 = 7.13, a0 = 1175.3
Среднее значения
tyt2y2t•yy(t)(y-ycp)2(y-y(t))2
11226.8111505062.021226.811182.4326.591969.62
2870.354757510.321740.71189.56123413.31101895.13
31238.1691533048.663714.491196.69272.591719.84
41342.37161801951.565369.471203.8214572.0919194.4
51321.07251745238.056605.371210.969884.6512126.19
61365.81361865450.098194.891218.0920782.6321823.45
71252.77491569433.898769.391225.22968.3759.1
81184.64641403371.149477.121232.351369.992276.31
91217.25811481689.2610955.221239.4819.42494.41
101163.031001352627.8211630.251246.613437.216987
111105.841211222883.4712164.251253.7513412.5121875.75
121371.731441881649.2116460.791260.8822523.7712288.93
7814659.8465018119915.4996308.7514659.84210683.05203410.13
Шаг 4. Определим компоненту T данной модели. Для этого проведем аналитическое выравнивание ряда (T + E) с помощью линейного тренда. Результаты аналитического выравнивания следующие:
T = 1175.298 + 7.132t
Подставляя в это уравнение значения t = 1,...,12, найдем уровни T для каждого момента времени (гр. 5 табл.).
tytSiyt/SiTTxSiE = yt / (T x Si)(yt - T*S)2
18980.731226.811182.43865.511.041055.31
27940.91870.351189.561085.210.7384801.95
314411.161238.161196.691392.741.032329.49
416001.191342.371203.821434.871.1227269.14
59670.731321.071210.96886.41.096497.14
612460.911365.811218.091111.231.1218162.51
714581.161252.771225.221425.931.021028.18
814121.191184.641232.351468.870.963233.92
98910.731217.251239.48907.280.98264.9
1010610.911163.031246.611137.260.935814.91
1112871.161105.841253.751459.130.8829630.23
1216351.191371.731260.881502.871.0917458.67
Шаг 5. Найдем уровни ряда, умножив значения T на соответствующие значения сезонной компоненты (гр. 6 табл.).
Расчет ошибки в мультипликативной модели производится по формуле:
E = Y/(T * S) = 12
Для сравнения мультипликативной модели и других моделей временного ряда можно использовать сумму квадратов абсолютных ошибок:
Среднее значения
ty(y-ycp)2
1898106384.69
2794185043.36
3144147016.69
41600141250.69
596766134.69
61246476.69
7145854678.03
8141235281.36
9891111000.03
10106126623.36
1112873948.03
121635168784.03
7814690946621.67


Следовательно, можно сказать, что мультипликативная модель объясняет 79% общей вариации уровней временного ряда.
Проверка адекватности модели данным наблюдения.

где m - количество факторов в уравнении тренда (m=1).
Fkp = 4.96
Поскольку F> Fkp, то уравнение статистически значимо
Шаг 6. Прогнозирование по мультипликативной модели. Прогнозное значение Ft уровня временного ряда в мультипликативной модели есть сумма трендовой и сезонной компонент. Для определения трендовой компоненты воспользуемся уравнением тренда:T = 1175.298 + 7.132t
Получим
T13 = 1175.298 + 7.132*13 = 1268.008
Значение сезонного компонента за соответствующий период равно: S1 = 0.732
Таким образом, F13 = T13 + S1 = 1268.008 + 0.732 = 1268.74
T14 = 1175.298 + 7.132*14 = 1275.14
Значение сезонного компонента за соответствующий период равно: S2 = 0.912
Таким образом, F14 = T14 + S2 = 1275.14 + 0.912 = 1276.052
T15 = 1175.298 + 7.132*15 = 1282.271
Значение сезонного компонента за соответствующий период равно: S3 = 1.164
Таким образом, F15 = T15 + S3 = 1282.271 + 1.164 = 1283.435
T16 = 1175.298 + 7.132*16 = 1289.403
Значение сезонного компонента за соответствующий период равно: S4 = 1.192
Таким образом, F16 = T16 + S4 = 1289.403 + 1.192 = 1290.595

Пример построения мультипликативной модели

Задание. На основе данных, скорректированных на инфляцию, о прибыли компании за 12 кварталов (табл.) построить мультипликативной модель тренда и сезонности для прогнозирования прибыли компании на следующие два квартала. Дать общую характеристику точности модели и сделать выводы.

Решение. Общий вид мультипликативной модели следующий:
Y = T x S x E
Эта модель предполагает, что каждый уровень временного ряда может быть представлен как сумма трендовой (T), сезонной (S) и случайной (E) компонент.
Рассчитаем компоненты мультипликативной модели временного ряда.
Шаг 1. Проведем выравнивание исходных уровней ряда методом скользящей средней. Для этого:
1.1. Найдем скользящие средние (гр. 3 таблицы). Полученные таким образом выровненные значения уже не содержат сезонной компоненты.
1.2. Приведем эти значения в соответствие с фактическими моментами времени, для чего найдем средние значения из двух последовательных скользящих средних – центрированные скользящие средние (гр. 4 табл.).

t yt Скользящая средняя Центрированная скользящая средняя Оценка сезонной компоненты
1 375 - - -
2 371 657.5 - -
3 869 653 655.25 1.33
4 1015 678 665.5 1.53
5 357 708.75 693.38 0.51
6 471 710 709.38 0.66
7 992 718.25 714.13 1.39
8 1020 689.25 703.75 1.45
9 390 689.25 689.25 0.57
10 355 660.5 674.88 0.53
11 992 678.25 669.38 1.48
12 905 703 690.63 1.31
13 461 685 694 0.66
14 454 690.5 687.75 0.66
15 920 - - -
16 927 - - -

Шаг 2. Найдем оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние (гр. 5 табл.). Эти оценки используются для расчета сезонной компоненты S. Для этого найдем средние за каждый период оценки сезонной компоненты Sj. Сезонные воздействия за период взаимопогашаются. В мультипликативной модели это выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна числу периодов в цикле. В нашем случае число периодов одного цикла равно 4.
Показатели 1 2 3 4
1 - - 1.33 1.53
2 0.51 0.66 1.39 1.45
3 0.57 0.53 1.48 1.31
4 0.66 0.66 - -
Всего за период 1.74 1.85 4.2 4.28
Средняя оценка сезонной компоненты 0.58 0.62 1.4 1.43
Скорректированная сезонная компонента, Si 0.58 0.61 1.39 1.42
Для данной модели имеем:
0.582 + 0.617 + 1.399 + 1.428 = 4.026
Корректирующий коэффициент: k=4/4.026 = 0.994
Рассчитываем скорректированные значения сезонной компоненты Si и заносим полученные данные в таблицу.
Шаг 3. Разделим каждый уровень исходного ряда на соответствующие значения сезонной компоненты. В результате получим величины T x E = Y/S (гр. 4 табл.), которые содержат только тенденцию и случайную компоненту.
Находим параметры уравнения методом наименьших квадратов.
Система уравнений МНК:
a0n + a1∑t = ∑y
a0∑t + a1∑t2 = ∑y•t
Для наших данных система уравнений имеет вид:
16a0 + 136a1 = 10872.41
136a0 + 1496a1 = 93531.1
Из первого уравнения выражаем а0 и подставим во второе уравнение
Получаем a0 = 3.28, a1 = 651.63
Среднее значения
y = ∑yi/n = 10872.41/16 = 679.53
t y t2 y2 t•y y(t) (y-ycp)2 (y-y(t))2
1 648.87 1 421026.09 648.87 654.92 940.05 36.61
2 605.46 4 366584.89 1210.93 658.2 5485.32 2780.93
3 625.12 9 390770.21 1875.35 661.48 2960.37 1322.21
4 715.21 16 511519.56 2860.82 664.76 1273.1 2544.83
5 617.72 25 381577.63 3088.6 668.04 3819.95 2532.22
6 768.66 36 590838.18 4611.96 671.32 7944.97 9474.64
7 713.6 49 509219.75 4995.17 674.6 1160.83 1520.44
8 718.73 64 516571.58 5749.83 677.88 1536.93 1668.26
9 674.82 81 455381.82 6073.38 681.17 22.14 40.28
10 579.35 100 335647.52 5793.51 684.45 10034.93 11045.26
11 713.6 121 509219.75 7849.56 687.73 1160.83 669.14
12 637.7 144 406656.13 7652.35 691.01 1749.71 2842.39
13 797.67 169 636280.07 10369.73 694.29 13958.53 10687.5
14 740.92 196 548957.15 10372.83 697.57 3768.85 1878.69
15 661.8 225 437983.3 9927.05 700.85 314.08 1524.97
16 653.2 256 426667.57 10451.17 704.14 693.14 2594.6
136 10872.41 1496 7444901.2 93531.1 10872.41 56823.71 53162.96

Шаг 4. Определим компоненту T данной модели. Для этого проведем аналитическое выравнивание ряда (T + E) с помощью линейного тренда. Результаты аналитического выравнивания следующие:
T = 651.634 + 3.281t
Подставляя в это уравнение значения t = 1,...,16, найдем уровни T для каждого момента времени (гр. 5 табл.).
t yt Si yt/Si T TxSi E = yt / (T x Si) (yt - T*S)2
1 375 0.58 648.87 654.92 378.5 0.99 12.23
2 371 0.61 605.46 658.2 403.31 0.92 1044.15
3 869 1.39 625.12 661.48 919.55 0.95 2555.16
4 1015 1.42 715.21 664.76 943.41 1.08 5125.42
5 357 0.58 617.72 668.04 386.08 0.92 845.78
6 471 0.61 768.66 671.32 411.36 1.14 3557.43
7 992 1.39 713.6 674.6 937.79 1.06 2938.24
8 1020 1.42 718.73 677.88 962.03 1.06 3359.96
9 390 0.58 674.82 681.17 393.67 0.99 13.45
10 355 0.61 579.35 684.45 419.4 0.85 4147.15
11 992 1.39 713.6 687.73 956.04 1.04 1293.1
12 905 1.42 637.7 691.01 980.66 0.92 5724.7
13 461 0.58 797.67 694.29 401.25 1.15 3569.68
14 454 0.61 740.92 697.57 427.44 1.06 705.39
15 920 1.39 661.8 700.85 974.29 0.94 2946.99
16 927 1.42 653.2 704.14 999.29 0.93 5225.65

Шаг 5. Найдем уровни ряда, умножив значения T на соответствующие значения сезонной компоненты (гр. 6 табл.).
Расчет ошибки в мультипликативной модели производится по формуле:
E = Y/(T * S) = 16
Для сравнения мультипликативной модели и других моделей временного ряда можно использовать сумму квадратов абсолютных ошибок:
Среднее значения
y = ∑yi/n = 10874/16 = 679.63
t y (y-ycp)2
1 375 92796.39
2 371 95249.39
3 869 35862.89
4 1015 112476.39
5 357 104086.89
6 471 43524.39
7 992 97578.14
8 1020 115855.14
9 390 83882.64
10 355 105381.39
11 992 97578.14
12 905 50793.89
13 461 47796.89
14 454 50906.64
15 920 57780.14
16 927 61194.39
136 10874 1252743.75

R2 = 1 - 43064.467/1252743.75 = 0.97
Следовательно, можно сказать, что мультипликативная модель объясняет 97% общей вариации уровней временного ряда.
Проверка адекватности модели данным наблюдения.
F = R2/(1 - R2)(n - m -1)/m = 0.972/(1 - 0.972)(16-1-1) = 393.26
где m - количество факторов в уравнении тренда (m=1).
Fkp = 4.6
Поскольку F > Fkp, то уравнение статистически значимо
Шаг 6. Прогнозирование по мультипликативной модели. Прогнозное значение Ft уровня временного ряда в мультипликативной модели есть сумма трендовой и сезонной компонент. Для определения трендовой компоненты воспользуемся уравнением тренда:T = 651.634 + 3.281t
Получим
T17 = 651.634 + 3.281*17 = 707.416
Значение сезонного компонента за соответствующий период равно: S1 = 0.578
Таким образом, F17 = T17 + S1 = 707.416 + 0.578 = 707.994
T18 = 651.634 + 3.281*18 = 710.698
Значение сезонного компонента за соответствующий период равно: S2 = 0.613
Таким образом, F18 = T18 + S2 = 710.698 + 0.613 = 711.311
T19 = 651.634 + 3.281*19 = 713.979
Значение сезонного компонента за соответствующий период равно: S3 = 1.39
Таким образом, F19 = T19 + S3 = 713.979 + 1.39 = 715.369
T20 = 651.634 + 3.281*20 = 717.26
Значение сезонного компонента за соответствующий период равно: S4 = 1.419
Таким образом, F20 = T20 + S4 = 717.26 + 1.419 = 718.68

Пример. На основе поквартальных данных построена мультипликативная модель временного ряда. Скорректированные значения сезонной компоненты за первые три квартала равны: 0,8 — I квартал, 1,2 — II квартал и 1,3 — III квартал. Определите значение сезонной компоненты за IV квартал.
Решение. Поскольку сезонные воздействия за период (4 квартала) взаимопогашаются, то имеем равенство: s1 + s2 + s3 + s4 = 4. Для наших данных: s4 = 4 - 0.8 - 1.2 - 1.3 = 0.7.
Ответ: Сезонная компонента за IV квартал равна 0.7.

Евгения
Помогу написать студенческую работу
×
Здравствуйте, меня зовут Евгения. Хотите мы оценим стоимость вашей работы
Меньше минуты назад
Это абсолютно бесплатно и ни к чему вас не обязывает
Меньше минуты назад