В табл. 1представлены квартальные данные об изменении прибыли фирмы за последние четыре года, тыс. долл.
Таблица 1
Год |
2013 |
2014 |
2015 |
2016 |
|||||||||||||||
Квартал |
I |
II |
III |
IV |
I |
II |
III |
IV |
I |
II |
III |
IV |
I |
II |
III |
IV |
|||
Прибыль, yt |
85,1 |
52,0 |
44,0 |
115,5 |
98,0 |
68,5 |
53,0 |
134,1 |
129,0 |
83,1 |
69,2 |
159,6 |
147,5 |
98,0 |
77,9 |
174,0 |
Построим коррелограмму временного ряда прибыли (рис. 1):
Рисунок 1 – Коррелограмма временного ряда
Анализ графика данного временного ряда (рис. 2) и коррелограммы (рис. 1) позволяет сделать вывод о наличии общей возрастающей тенденции уровней ряда и сезонных колебаний периодичностью в четыре квартала.
Рисунок 2 – График временного ряда прибыли
Поскольку амплитуда сезонных колебаний увеличивается, можно предположить наличие мультипликативной модели: Y=T·S·E. Определим ее компоненты.
Шаг 1. Проведем выравнивание исходных уровней ряда методом скользящей средней (табл. 2).
Таблица 2 – Расчет оценок сезонной компоненты
№ квартала, t |
Прибыль, yt |
Итого по 4 кварталам |
Скользящая средняя по кварталам |
Центрированная скользящая средняя |
Оценка сезонной компоненты |
1 |
2 |
3 |
4 |
5 |
6 |
1 |
85,1 |
- |
- |
- |
- |
2 |
52,0 |
296,6 |
74,2 |
- |
- |
3 |
44,0 |
309,5 |
77,4 |
75,76 |
0,581 |
4 |
115,5 |
326,0 |
81,5 |
79,44 |
1,454 |
5 |
98,0 |
335,0 |
83,8 |
82,63 |
1,186 |
6 |
68,5 |
353,6 |
88,4 |
86,08 |
0,796 |
7 |
53,0 |
384,6 |
96,2 |
92,28 |
0,574 |
8 |
134,1 |
399,2 |
99,8 |
97,98 |
1,369 |
9 |
129,0 |
415,4 |
103,9 |
101,83 |
1,267 |
10 |
83,1 |
440,9 |
110,2 |
107,04 |
0,776 |
11 |
69,2 |
459,4 |
114,9 |
112,54 |
0,615 |
12 |
159,6 |
474,3 |
118,6 |
116,71 |
1,367 |
13 |
147,5 |
483,0 |
120,8 |
119,66 |
1,233 |
14 |
98,0 |
497,4 |
124,4 |
122,55 |
0,800 |
15 |
77,9 |
- |
- |
- |
- |
16 |
174,0 |
- |
- |
- |
- |
Шаг 2. Найдем оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние (табл. 2).
Найдем (табл. 3) средние за каждый квартал (по всем годам) оценки сезонной компоненты Si.
Таблица 3 – Расчет значений сезонной компоненты в мультипликативной модели
Показатели |
№ квартала, i |
|||
I |
II |
III |
IV |
|
Сезонная компонента 1-го года |
- |
- |
0,581 |
1,454 |
Сезонная компонента 2-го года |
1,186 |
0,796 |
0,574 |
1,369 |
Сезонная компонента 3-го года |
1,267 |
0,776 |
0,615 |
1,367 |
Сезонная компонента 4-го года |
1,233 |
0,800 |
||
Итого за i-й квартал (за все годы) |
3,686 |
2,372 |
1,770 |
4,190 |
Средняя оценка сезонной компоненты для i-гo квартала, |
1,229 |
0,791 |
0,590 |
1,397 |
Скорректированная сезонная компонента, Si |
1,227 |
0,789 |
0,589 |
1,395 |
Взаимопогашаемость сезонных воздействий в мультипликативной модели выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна числу периодов в цикле, т.е. четырем. Для данной модели имеем: .
Определим корректирующий коэффициент:
.
Рассчитаем скорректированные значения сезонной компоненты, умножив ее средние оценки на корректирующий коэффициент k .
Получены следующие значения сезонной компоненты по кварталам года: S1 = 1,263; S2 = 0,751; S3 = 0,586; S4 = 1,399.
Занесем полученные значения в табл. 4 для соответствующих кварталов каждого года.
Таблица 4 – Расчет выравненных значений тренда и ошибок в мультипликативной модели
t |
yt |
Si |
T·E= = yt / Si |
T |
T·S |
E=yt: :(T·S) |
E=yt – – (T·S) |
E2 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
1 |
85,1 |
1,227 |
69,37 |
65,97 |
80,92 |
1,052 |
4,18 |
17,46 |
2 |
52,0 |
0,789 |
65,87 |
70,41 |
55,59 |
0,935 |
-3,59 |
12,88 |
3 |
44,0 |
0,589 |
74,68 |
74,86 |
44,10 |
0,998 |
-0,10 |
0,01 |
4 |
115,5 |
1,395 |
82,82 |
79,31 |
110,61 |
1,044 |
4,89 |
23,92 |
5 |
98,0 |
1,227 |
79,89 |
83,76 |
102,75 |
0,954 |
-4,75 |
22,53 |
6 |
68,5 |
0,789 |
86,77 |
88,20 |
69,63 |
0,984 |
-1,13 |
1,28 |
7 |
53,0 |
0,589 |
89,96 |
92,65 |
54,59 |
0,971 |
-1,59 |
2,51 |
8 |
134,1 |
1,395 |
96,15 |
97,10 |
135,42 |
0,990 |
-1,32 |
1,75 |
9 |
129,0 |
1,227 |
105,16 |
101,55 |
124,57 |
1,036 |
4,43 |
19,62 |
10 |
83,1 |
0,789 |
105,26 |
106,00 |
83,68 |
0,993 |
-0,58 |
0,33 |
11 |
69,2 |
0,589 |
117,46 |
110,44 |
65,07 |
1,064 |
4,13 |
17,08 |
12 |
159,6 |
1,395 |
114,44 |
114,89 |
160,23 |
0,996 |
-0,63 |
0,40 |
13 |
147,5 |
1,227 |
120,24 |
119,34 |
146,40 |
1,008 |
1,10 |
1,22 |
14 |
98,0 |
0,789 |
124,14 |
123,79 |
97,72 |
1,003 |
0,28 |
0,08 |
15 |
77,9 |
0,589 |
132,23 |
128,23 |
75,55 |
1,031 |
2,35 |
5,53 |
16 |
174,0 |
1,395 |
124,76 |
132,68 |
185,05 |
0,940 |
-11,05 |
122,02 |
Шаг 3. Разделим каждый уровень исходного ряда на соответствующие значения сезонной компоненты. Получим: Т · Е = Y / S, которые содержат только тенденцию и случайную компоненту.
Шаг 4. Определим компоненту Т в мультипликативной модели модели. Для этого рассчитаем параметры линейного тренда, используя уровни (Т·Е). Результаты аналитического выравнивания следующие:
.
Подставив в это уравнение значения t = 1, ..., 16, найдем уровни Т для каждого момента времени.
График уравнения тренда приведен на рис. 3.
Рисунок 3 – Моделирование динамики прибыли компании
Шаг 5. Найдем уровни ряда по мультипликативной модели, умножив уровни Т на значения сезонной компоненты для соответствующих кварталов. Графически значения (Т·S) представлены на рис. 3.
Шаг 6. Расчет ошибок в мультипликативной модели проводится по формуле . Если временной ряд ошибок не содержит автокорреляции, его можно использовать вместо исходного ряда для изучения взаимосвязи с другими временными рядами.
Библиографический список:
Дайитбегов, Д.М. Компьютерные технологии анализа данных в эконометрике: Монография / Д.М. Дайитбегов. - 3-e изд., испр. и доп. - М.: Вузовский учебник: НИЦ Инфра-М, 2013. - XIV, 587 с.
Тимофеев, В.С. Эконометрика [Электронный ресурс]: учебник/ В.С. Тимофеев, А.В. Фаддеенков, В.Ю. Щеколдин. – Электрон. текстовые дан. - 3-e изд., перераб. и доп. - Новосибирск: Изд-во НГТУ, 2013. – 340 с. – Режим доступа: http://znanium.com/bookread2.php?book=546264
Федеральная служба государственной статистики http://www.gks.ru/