Экономические прогнозы по временным рядам. Учебное пособие
.pdfНомер |
Цена акции, yt |
Индекс фондового |
|
наблюдения, t |
рынка, xt |
||
|
|||
13 |
73 |
61 |
|
14 |
67 |
61 |
|
15 |
70 |
59 |
|
16 |
80 |
63 |
|
17 |
81 |
65 |
|
18 |
90 |
67 |
|
19 |
88 |
68 |
Результаты расчета сведены в следующую таблицу: Таблица 5
Расчет параметров регрессионной модели
|
t |
х |
t |
y |
t |
х – х |
(х – х )2 |
(y – y ) |
(х – хср) |
a |
0 |
+ a x |
t |
Е |
t |
|
|
|
|
ср |
ср |
ср |
(y – y ) |
|
1 |
|
|||||
|
|
|
|
|
|
|
|
|
ср |
|
|
|
|
|
|
|
1 |
45 |
25 |
–10 |
100 |
–31,7 |
317 |
|
30,3 |
|
–5,3 |
||||
|
2 |
45 |
25 |
–10 |
100 |
–31,7 |
317 |
|
|
Зоз |
|
–5,3 |
|||
|
3 |
47 |
34 |
–8 |
64 |
–22,7 |
181,6 |
|
35,6 |
|
–1,6 |
||||
|
4 |
46 |
39 |
–9 |
81 |
–17,7 |
159,3 |
|
32,9 |
|
6,1 |
||||
|
5 |
50 |
42 |
–5 |
25 |
–14,7 |
73,5 |
|
43,5 |
|
–1,5 |
||||
|
6 |
48 |
37 |
–7 |
49 |
–19,7 |
137,9 |
|
38,2 |
|
–1,2 |
||||
|
7 |
48 |
45 |
–7 |
49 |
–11,7 |
81,9 |
|
38,2 |
|
6,8 |
||||
|
8 |
50 |
44 |
–5 |
25 |
–12,7 |
63,5 |
|
43,5 |
|
0,5 |
||||
|
9 |
54 |
55 |
–1 |
1 |
–1,7 |
1,7 |
|
54,1 |
|
0,9 |
||||
|
10 |
55 |
52 |
0 |
0 |
–4,7 |
0 |
|
56,7 |
|
47 |
||||
|
11 |
57 |
67 |
2 |
4 |
10,3 |
20,6 |
|
|
62 |
|
5 |
|
||
|
12 |
56 |
63 |
1 |
1 |
6,3 |
6,3 |
|
59,3 |
|
3,7 |
||||
|
13 |
61 |
73 |
6 |
36 |
16,3 |
97,8 |
|
72,5 |
|
0,5 |
||||
|
14 |
61 |
67 |
6 |
36 |
10,3 |
61,8 |
|
72,5 |
|
–5,5 |
||||
|
15 |
59 |
70 |
4 |
16 |
13,3 |
53,2 |
|
67,3 |
|
2,7 |
||||
|
16 |
63 |
80 |
8 |
64 |
23,3 |
186,4 |
|
77,8 |
|
2Д |
||||
|
17 |
65 |
81 |
10 |
100 |
24,3 |
243 |
|
83,1 |
|
–2,1 |
||||
|
18 |
67 |
90 |
12 |
144 |
33,3 |
399,6 |
|
88,4 |
|
1,6 |
||||
|
19 |
68 |
88 |
13 |
169 |
31,3 |
406,9 |
|
|
91 |
|
3 |
|
||
Итого |
190 |
1045 |
1077 |
0 |
1064 |
–0,3 |
2809 |
1077,2 |
|
–0,2 |
|||||
Средние |
10 |
55 |
56,7 |
|
56 |
–0,02 |
147,84 |
56,69 |
|
–0,01 |
|||||
Расчеты параметров, полученного прогноза приведены в таблице 6.
21
22
Таблица 6
Проверка адекватности модели линейной регрессии
|
t |
хt |
yt |
Еt |
Точки поворота |
Еt2 |
Еt – Еt – 1 |
(Еt – Еt – 1)2 |
Еt *Еt – 1 |
|Еt / y| |
|
1 |
45 |
25 |
–5,3 |
– |
28,09 |
– |
– |
– |
0,21 |
|
2 |
45 |
25 |
–5,3 |
0 |
28,09 |
0 |
0 |
28,09 |
0,21 |
|
3 |
47 |
34 |
–1,6 |
0 |
2,56 |
3,7 |
13,69 |
8,48 |
0,05 |
|
4 |
46 |
39 |
6,1 |
1 |
37,21 |
7,7 |
59,29 |
–9,76 |
0,16 |
|
5 |
50 |
42 |
–1,5 |
1 |
2,25 |
–7,6 |
57,76 |
–9,15 |
0,04 |
|
6 |
48 |
37 |
–1,2 |
0 |
1,44 |
0,3 |
0,09 |
1,8 |
0,03 |
|
7 |
48 |
45 |
6,8 |
1 |
46,24 |
8 |
64 |
–8,16 |
0,15 |
|
8 |
50 |
44 |
0,5 |
1 |
0,25 |
–6,3 |
39,69 |
3,4 |
0,01 |
|
9 |
54 |
55 |
0,9 |
1 |
0,81 |
0,4 |
0,16 |
0,45 |
0,02 |
|
10 |
55 |
52 |
–4,7 |
1 |
22,09 |
–5,6 |
31,36 |
–4,23 |
0,09 |
|
11 |
57 |
67 |
5 |
1 |
25 |
9,7 |
94,09 |
–23,5 |
0,07 |
|
12 |
56 |
63 |
3,7 |
0 |
13,69 |
–1,3 |
1,69 |
18,5 |
0,06 |
|
13 |
61 |
73 |
0,5 |
0 |
0,25 |
–3,2 |
10,24 |
1,85 |
0,01 |
|
14 |
61 |
67 |
–5,5 |
1 |
30,25 |
–6 |
36 |
–2,75 |
0,08 |
|
15 |
59 |
70 |
2,7 |
1 |
7,29 |
8,2 |
67,24 |
–14,85 |
0,04 |
|
16 |
63 |
80 |
2,2 |
0 |
4,84 |
–0,5 |
0,25 |
5,94 |
0,03 |
|
17 |
65 |
81 |
–2,1 |
1 |
4,41 |
–4,3 |
18,49 |
–4,62 |
0,03 |
|
18 |
67 |
90 |
1,6 |
1 |
2,56 |
3,7 |
13,69 |
–3,36 |
0,02 |
|
19 |
68 |
88 |
–3 |
– |
9 |
–4,6 |
21,16 |
–4,8 |
0,03 |
|
|
|
|
|
|
|
|
|
|
|
Итого |
190 |
1045 |
1077 |
–0,2 |
11 |
266,32 |
2,3 |
528,89 |
–16,67 |
1,34 |
Средние |
10 |
55 |
56,68 |
–0,01 |
0,58 |
14,02 |
0,12 |
27,84 |
–0,88 |
0,07 |
Проверим адекватность полученной модели и найдем точечный и интервальный прогноз для t = 20 (k = 1, n = 19) .
Случайность Да |
g = 11 > 9 |
|
|
|
|
Независимость Да |
d = 2 > 1,4 |
|
|
|
|
Нормальность Да |
S = 3,8 |
Emin = –5,5 |
RS = 3,198 |
||
Emax |
= 6,8 |
3,14 < RS < 4,43 |
|||
|
|
||||
Точность модели: MSE = 14,02; МАРЕ= 7,05; ME = 3,75
Доверительный интервал: u (1) = 9,3 при р = 0,05; R2 = 0,96 Результат прогноза yТ (t) = a0 + a1t будет:
t |
yT (t) |
Нижняя гр. |
Верхняя гр. |
20 |
93,7 |
84,4 |
103,0 |
Сравнивая точечные прогнозные оценки модели регрессии с оценками по линейной временной модели, можно отметить их явную близость, однако доверительный интервал регрессионной модели заметно шире, что снижает ее практическую значимость.
Адаптивная модель статистически полностью адекватна и имеет достаточно высокие точностные характеристики. Ее результаты можно взять в качестве прогноза.
Результаты сопоставления трех методов прогнозирования рассматриваемого ряда y (t) можно представить в следующей таблице:
|
|
|
|
|
|
Таблица 7 |
|
Прогнозы, полученные различными методами |
|||||||
|
|
и адекватность моделей |
|
|
|
||
|
|
|
|
|
|
|
|
|
Точеч- |
Нижняя |
Верхняя |
Свойства остатков |
|||
Метод |
ный |
граница |
граница |
Случай- |
Нормаль- |
Незави- |
|
|
прогноз |
|
|
ность |
ность |
|
симость |
Экстраполя- |
92,8 |
83,6 |
102,0 |
да |
да |
|
да |
ция тренда |
|
|
|
|
|
|
|
Сглажива- |
93,3 |
88,4 |
98,2 |
да |
да |
|
да |
ние Брауна |
|
||||||
|
|
|
|
|
|
|
|
Линейная |
93,7 |
84,4 |
103,0 |
да |
да |
|
да |
регрессия |
|
||||||
|
|
|
|
|
|
|
|
23
Близость прогнозных оценок, полученных различными методами, свидетельствует о достаточно высоком качестве выполненного прогноза. Вместе с тем, практика показывает, что для оперативных и краткосрочных прогнозов предпочтительнее модельадаптивногосглаживания,адлясреднесрочныхидолгосрочных — модель экстраполяции тренда. Модель регрессионного анализа носит универсальный характер, но точность прогноза по ней обычно ниже первых двух.
§ 3. Временные ряды с сезонной компонентой
Временные ряды характеризуются не только трендом, но и подвержены систематическим колебаниям. Периодичность таких колебаний может варьировать в весьма широких пределах.
Если продолжительность цикла превышает год, то говорят о циклических колебаниях. При колебаниях с периодом длительностью не более одного года обычно используют термин «сезонные колебания». В зависимости от характера изучаемых процессов длительность того или иного сезона может варьировать в самых широких пределах: от полугодия до одного дня или даже часа.
Для выявления и последующего выделения систематических колебаний необходимо, чтобы временной ряд заключал в себе несколько полных циклов. В данной работе мы не будем непосредственнорассматриватьпримерыссобственноциклическимиколебаниями, хотя изложенные ниже подходы применимы и во многих подобных случаях. Остановимся подробнее на прогнозировании временных рядов, содержащих сезонную компоненту [4].
Запишеммодельвременногорядасаддитивнойсезоннойкомпонентой. Представим фактическое значение временного ряда yt в виде суммы трех компонент:
y(t) = T(t) +S(t) +E(t) |
(20) |
где T (t) — систематическая компонента, или значение тренда; S (t) — сезонная компонента;
E (t) — случайная компонента.
Процедура прогнозирования будущих значений переменной с учетом сезонной компоненты состоит из нескольких этапов.
24
На первом этапе проводится расчет значений сезонной компоненты. Для этого рассчитывается скользящая средняя по всей совокупности значений временного ряда. Сущность метода скользящей средней заключается в исчислении среднего значения для определенного числа, называемого параметром сглаживания, первых по счету элементов временного ряда, затем для такого же числа элементов, начиная со второго по счету, и так далее. Таким образом, средняя как бы скользит по ряду, передвигаясь на один интервал или момент времени.
В зависимости от особенностей (характера, типа) динамического ряда можно усреднять и другими методами, например, — двойнойскользящейсредней,экспоненциальновзвешеннойсредней и т.п.
Следуетиметьввиду,чтопричетномзначениипараметрасглаживаниязначенияскользящейсреднейотносятсякнесколькоиным моментам времени, чем фактические данные, располагаясь в промежутке между моментами или на пересечении интервалов. Иначе говоря, при отнесении скользящей средней к интервалам фактических данных ее значения окажутся опережающими (у+) или запаздывающими (у–) на половину интервала. В этом случае проводится дополнительная синхронизация путем расчета центрированной скользящейсредней.Этаоперациясводитсякповторениюпроцедуры сглаживания, но уже для ряда полученных значений скользящей средней и при значении параметра сглаживания, равном двум. При нечетном значении параметра сглаживания полученные оценки скользящей средней относятся к середине выбранного интервала.
Рассмотримнапримереособенностипрогнозированиядинамического ряда, характеризующегося наличием сезонных колебаний.
Пример 3. Известны данные о продажах товара за 11 последних кварталов:
Год |
Квартал |
Продано, тыс. шт. |
1 |
I |
40 |
1 |
II |
76 |
1 |
III |
80 |
1 |
IV |
60 |
25
Год |
Квартал |
Продано, тыс. шт. |
2 |
I |
80 |
2 |
II |
122 |
2 |
III |
120 |
2 |
IV |
100 |
3 |
I |
120 |
3 |
II |
160 |
3 |
III |
160 |
3 |
IV |
? |
4 |
I |
? |
Требуется дать прогноз продаж товара на 12 и 13 квартал. При построении модели принимается предположение, что амплитуда и фаза сезонного цикла не подвержена случайным изменениям. Это предположение достаточно грубое и при накоплении данных о рассматриваемом явлении должно уточняться. Однако, на начальном этапе исследования экономического явления, в условиях дефицита данных, такое предположение обычно
применяется.
Построение модели и прогноз включает этапы:
–ввод данных;
–выбор типа модели;
–выравнивание исходного ряда для выделения сезонного
хода;
–вычисление сезонной компоненты;
–удаление сезонной компоненты из исходного ряда;
–вычисление тренда и оценка уровня ошибок;
–прогноз значений ряда с учетом тренда, сезонного цикла и уровня ошибок.
Процесс разработки модели сопровождается построением графиков, по виду которых принимаются эвристические решения на различных этапах.
Если по поведению исходных данных не удаётся сделать однозначный выбор в пользу аддитивной или мультипликативной модели, расчеты выполняются для обоих вариантов. Обе модели тестируются на обучающих примерах и выбирается та, которая дает более точные результаты для типовых ситуаций.
26
Приведём ход построения прогноза в среде Excel по данным
ссезонным циклом.
1.Ввод данных и построение графика продаж по кварталам.
2.Выбор типа модели с использованием вида графика.
•Аддитивная T + S + E;
•Мультипликативная T · S · E,
где T — значение тренда, S — сезонный ход, E — ошибка.
Внашей модели амплитуда сезонных колебаний меняется слабо, а следовательно аддитивная модель предпочтительней.
3. Выравнивание исходного ряда методом скользящей средней.
Внашем случае — скользящее среднее строится по четырем точкам.
Поскольку среднее за четыре точки попадает между вторым и третьим кварталами, необходимо их привести к точкам наблюдения по схеме двойного осреднения, как показано в таблице:
Таблица 8
Схема двойного осреднения для чётного числа отсчётов в сезонном цикле
Номер кв |
Продано |
Скользящее |
Приведённые |
Сезонная |
|
среднее |
значения |
компонента |
|||
|
|
||||
1 |
40 |
|
|
|
|
2 |
76 |
|
|
|
|
3 |
80 |
64 |
69 |
11 |
|
4 |
60 |
74 |
79,75 |
–19,8 |
|
5 |
80 |
85,5 |
90,5 |
–10,5 |
|
6 |
122 |
95,5 |
100,5 |
21,5 |
|
7 |
120 |
105,5 |
110,5 |
9,5 |
|
8 |
100 |
115,5 |
120,25 |
–20,3 |
|
9 |
120 |
125 |
130 |
–10 |
|
10 |
160 |
135 |
|
|
|
11 |
160 |
|
|
|
|
12 |
|
|
|
|
|
13 |
|
|
|
|
В нашем случае z’2 = (40 + 76 + 80 + 60) / 4 = 64;
z’3 = (76 + 80 + 60 + 80) / 4 = 74; z3 = (z’2 + z’3) / 2 = 69 и т.д.
27
4. Оценка поквартальных значений сезонной компоненты S; В силу предположения о неслучайности сезонной компоненты, общая сумма всех сезонных амплитуд за год должна быть нулевой. С этой целью значения амплитуд сезонного хода корректи-
руются как показано в следующей таблице.
Таблица 9
Расчет коррекции сезонного цикла
Год |
Оценка сезонной компоненты |
сумма |
||||
кв1 |
кв2 |
кв3 |
кв4 |
|||
|
|
|||||
1 |
|
|
11 |
–19,8 |
|
|
2 |
–10,5 |
21,5 |
9,5 |
–20,3 |
|
|
3 |
–10 |
|
|
|
|
|
Амплитуда |
–10,3 |
21,5 |
10,3 |
–20,1 |
1,4 |
|
Коррекция |
–10,65 |
21,15 |
9,95 |
–20,45 |
0 |
|
Амплитуда сезонного хода считается как среднее многолетнее по соответствующим кварталам. Коррекция вносится в каждый квартал, чтобы суммарный сезонный ход был нулевым.
5. Устранение сезонной компоненты из исходных уровней ряда (получение выровненных данных T + E) путем вычета из исходного объема продаж откорректированных сезонных данных:
Квартал |
Сезонный ход |
Десезональные данные |
|
y |
оценка: T + E |
||
1 |
–10,65 |
40 |
50,65 |
2 |
21,15 |
76 |
54,85 |
3 |
9,95 |
80 |
70,05 |
4 |
–20,45 |
60 |
80,45 |
5 |
–10,65 |
80 |
90,65 |
6 |
21,15 |
122 |
100,85 |
7 |
9,95 |
120 |
110,05 |
8 |
–20,45 |
100 |
120,45 |
9 |
–10,65 |
120 |
130,65 |
10 |
21,15 |
160 |
138,85 |
11 |
9,95 |
160 |
150,05 |
12 |
|
|
|
13 |
|
|
|
28
6. Получение уравнения тренда путем построения регрессии номера квартала (i) и оценки (T + E):
T = 39,003 + 10,129 · i
Среднеквадратическая ошибка аппроксимации тренда 1,76.
Модель описывает 99% изменчивости тренда. Уровень значимости 10–12.
7. Модель дает следующий прогноз на 2 квартала:
Номер квартала |
Тренд (T) |
Сезон (S) |
T + S |
12 |
160,55 |
–20,45 |
140,10 |
13 |
170,68 |
–10,65 |
160,03 |
Согласно модели в 4 квартале текущего года ожидается снижение объема продаж на 20 тыс. единиц, а в 1 квартале следующего года — увеличение объема продаж до уровня 3 квартала текущего года.
Расчеты по мультипликативной модели могут быть выполнены аналогично, но при этом, средние арифметические заменяются на средние геометрические, суммы — на умножение, а разности — на деление.
Литература
1.Льюис К.Д. Методы прогнозирования экономических показателей. – М.: Финансы и статистика, 1986.
2.Тюрин Ю.Н., Макаров А.А. Анализ данных на компьюте-
ре. – М.: Инфра-М, 2003.
3.Федосеев В.В. Математическое моделирование в экономике и социологии труда. – М.: ЮНИТИ-ДАНА, 2007.
4.Эконометрика: Учебник для вузов / Под ред. И.И. Елисеевой. – М.: Финансы и статистика, 2007.
29
Приложения
Приложение 1
Критические значения для проверки гипотез на уровне значимости 0,05 для критерия Дарбина–Уотсона d1 и d2 для модели с одним (m = 1) и двумя (m = 2) параметрами:
Количество |
r (1) |
m = 1 |
m = 2 |
||
наблюдений |
d1 |
d2 |
d1 |
d2 |
|
10 |
0,360 |
– |
– |
– |
– |
15 |
0,328 |
1,08 |
1,36 |
0,95 |
1,54 |
20 |
0,300 |
1,20 |
1,41 |
1,10 |
1,54 |
25 |
0,276 |
1,28 |
1,45 |
1,20 |
1,55 |
30 |
0,257 |
1,35 |
1,49 |
1,28 |
1,57 |
30
