Экономические прогнозы по временным рядам. Учебное пособие
.pdf
Для описания тенденций наиболее часто используют следую-
щие функции: |
f (t) = a + a t; |
|
||||
– |
линейная |
|
||||
– |
параболическая |
f (t) = a0 |
+ a1 t + a t2; |
|||
– |
степенная |
f (t) = a0 t a 1;1 |
2 |
|||
– |
экспоненциальная |
f (t) = a |
0 |
exp (a t); |
||
0 |
||||||
– |
логарифмическая |
|
+a1 ln |
1 |
||
f (t) = a0 |
(t); |
|||||
– |
кривая Торнквиста |
f (t) = a0 |
+ t / (t + a1); |
|||
– |
кривая Перла |
f (t) = exp (a0 – a1 exp(t)); |
||||
и другие. |
|
|
|
|
||
Покажем, как выполняется расчет тренда на примере линей- |
||||||
ной кривой роста: |
|
|
|
|
||
f (t) = a0 + a1 t |
|
|
|
(9) |
||
Параметры a0 и a1 оцениваются методом наименьших квадратов (МНК), т.е. подбираются таким образом, чтобы график
функции кривой роста располагался на минимальном удалении от данных наблюдений. Математически критерий оценки параметров линейной модели записывается так:
∑(y(t) −a0 −a1t)2 →min . |
(10) |
Приравняв нулю частные производные по а0 и а1, можно по- |
|
лучить систему двух уравнений с двумя неизвестными:
(11)
Решив эту систему, мы получим параметры линейного роста:
(12)
,
11
где
.
После расчета коэффициентов линейного тренда следует оценить ошибки модели и проверить её адекватность.
Для оценки качества аппроксимации полученного уравнения тренда применяется можно использовать коэффициент детерминации R2, показывающий долю вариации прогнозируемой результативной переменной у, объясняемую влиянием включенных в модель факторов, в данном случае фактором времени t.
Коэффициент детерминации рассчитывается по формуле:
(13)
Интервальнаяоценкапрогнозаполучаетсяисходяизточечной оценки, стандартной ошибки прогноза и значения t-статистики Стьюдента tp для заданной вероятности ошибки прогноза, p с числом степеней свободы, n – 1.
Стандартная ошибка прогноза определяется по формуле:
(14)
,
где k — период упреждения, s — стандартная среднеквадратическая ошибка уравнения регрессии, определяемая следующим образом:
s = |
|
∑n Et2 |
|
(15) |
t=1 |
||||
|
|
n −2 |
|
|
Абсолютная величина доверительного интервала прогноза U (k) при вероятности ошибки p составит:
12
U (k) = t p S(k) |
(16) |
Доверительный интервал прогноза будет иметь следующие границы:
верхняя граница: уT (n + k) + U (k); нижняя граница: уT (n + k) – U (k),
где уT(n + k) — прогнозное (теоретическое) значение тренда в момент времени t = n + k .
Пример 1.
Еженедельный индекс курса акций представлен следующим рядом:
Таблица 1
Цена обыкновенной корпоративной акции, в среднем за квартал (у.е.)
Номер наблюдений, t |
Цена акции, yt |
1 |
25 |
2 |
25 |
3 |
34 |
4 |
39 |
5 |
42 |
6 |
37 |
7 |
45 |
8 |
44 |
9 |
55 |
10 |
52 |
11 |
67 |
12 |
63 |
13 |
73 |
14 |
67 |
15 |
70 |
16 |
80 |
17 |
81 |
18 |
90 |
19 |
88 |
13
Рассчитаем коэффициенты тренда, сформировав следующую таблицу:
Таблица 2
Расчеты для оценки коэффициентов и проверки адекватности линейного тренда
|
|
|
|
|
|
|
|
|
|
|
t |
y |
E |
Точки |
E 2 |
E –E |
(E –E )2 |
|
|E /y| |
||
пово- |
E ·E |
|||||||||
|
t |
t |
рота |
t |
t |
t–1 |
t |
t–1 |
t t–1 |
t t |
1 |
25 |
0,8 |
– |
0,6 |
|
– |
|
– |
– |
0,03 |
2 |
25 |
–2,8 |
1 |
7,8 |
|
–3,6 |
|
13,0 |
–2,24 |
0,11 |
3 |
34 |
2,6 |
0 |
6,8 |
|
5,4 |
|
29,2 |
–7,28 |
0,08 |
4 |
39 |
4,0 |
1 |
16,0 |
|
1,4 |
|
2,0 |
10,4 |
0,1 |
5 |
42 |
3,3 |
0 |
10,9 |
|
–0,7 |
|
0,5 |
13,2 |
0,08 |
6 |
37 |
–5,3 |
1 |
28,1 |
|
–8,6 |
|
74,0 |
–17,49 |
0,14 |
7 |
45 |
–0,9 |
1 |
0,8 |
|
4,4 |
|
19,4 |
4,77 |
0,02 |
8 |
44 |
–5,5 |
1 |
30,3 |
|
–4,6 |
|
21,2 |
4,95 |
0,13 |
9 |
55 |
1,9 |
1 |
3,6 |
|
7,4 |
|
54,8 |
–10,45 |
0,03 |
10 |
52 |
–4,7 |
1 |
22,1 |
|
–6,6 |
|
43,6 |
–8,93 |
0,09 |
11 |
67 |
6,7 |
1 |
44,9 |
|
11,4 |
|
130,0 |
–31,49 |
0,1 |
12 |
63 |
–0,9 |
1 |
0,8 |
|
–7,6 |
|
57,8 |
–6,03 |
0,01 |
13 |
73 |
5,5 |
1 |
30,3 |
|
6,4 |
|
41,0 |
–4,95 |
0,08 |
14 |
67 |
–4,1 |
0 |
16,8 |
|
–9,6 |
|
92,2 |
–22,55 |
0,06 |
15 |
70 |
–4,8 |
1 |
23,0 |
|
–0,7 |
|
0,5 |
19,68 |
0,07 |
|
|
|
|
|
|
|
|
|
|
|
14
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
y |
E |
Точки |
E 2 |
E –E |
(E –E )2 |
|
|E /y| |
||
|
пово- |
E ·E |
|||||||||
|
|
t |
t |
рота |
t |
t |
t–1 |
t |
t–1 |
t t–1 |
t t |
|
16 |
80 |
1,6 |
1 |
2,6 |
|
6,4 |
|
41,0 |
–7,68 |
0,02 |
|
17 |
81 |
–1,0 |
1 |
1,0 |
|
–2,6 |
|
6,8 |
–1,6 |
0,01 |
|
18 |
90 |
4,4 |
1 |
19,4 |
|
5,4 |
|
29,2 |
–4,4 |
0,05 |
|
19 |
88 |
–1,2 |
– |
1,4 |
|
–5,6 |
|
31,4 |
–5,28 |
0,01 |
|
|
|
|
|
|
|
|
|
|
|
|
Итого |
190 |
1077 |
–0,4 |
14 |
267,1 |
|
–2,0 |
686,9 |
–77,37 |
1,22 |
|
|
|
|
|
|
|
|
|
|
|
||
Средние |
10 |
56,68 |
–0,02 |
0,74 |
14,06 |
–0,11 |
36,15 |
–4,07 |
0,06 |
||
|
|
|
|
|
|
|
|
|
|
|
|
Проверим адекватность полученной модели (табл. 2), и найдем точечный и интервальный прогноз для t = 20 (k = 1, n = 19).
Случайность Да |
g = 14 >7 |
|
|
|
Независимость Да |
d = 2,57 –d = 1,429 > 1,40 |
|
||
Нормальность Да |
S = 3,85 |
Emin = –5,5 |
RS = 3,167 |
|
Emax = 6,7 |
3,14 < RS < 4,43 |
|||
|
|
|||
Точность модели: MSE = 14,06; МАРЕ = 6,42 Доверительный интервал: u (1) = 9,2 при р = 0,05; R2 = 0,965 Результат прогноза yТ (t) = a0 + a1t будет:
t |
yT (t) |
Нижняя гр. |
Верхняя гр. |
20 |
92,8 |
83,6 |
102,0 |
2.2. Адаптивные методы прогнозирования
В практических задачах часто бывает важно учитывать не тенденцию, которая сложилась в среднем на всем интервале развития процесса, а последние значения прогнозируемого показателя. В силу инерционности значений уровней временного ряда, краткосрочный прогноз может быть выполнен на основании учета последних значений или построения прогнозируемой функции таким образом, чтобы последним отсчетам придавались большие
15
веса, чем более ранним. Опыт показывает, что при краткосрочном прогнозировании более эффективны адаптивные методы, которыеучитываютнеравноценность информации о недавних и более раннихзначенияхвременногорядаибыстроприспосабливающие свою структуру к изменяющимся условиям. Выделяют два класса адаптивных методов: скользящего среднего и авторегрессии.
Посхемескользящегосреднегооценкойтекущегоуровняявляется взвешенное значение всех предшествующих уровней, причем веса наблюдений убывают по мере удаления от последнего уровня. Такие модели хорошо реагируют на изменения в тенденции.
Вавторегрессионных моделях оценка текущего уровня является взвешенная сумма не всех, а только нескольких предшествующих уровней, причем учет уровня в оценке определяется не близостью к последнему отсчету, а теснотой связи между ними.
Впрактике социально-экономического прогнозирования наиболее часто используются схемы скользящего среднего. Самый простой метод заключается в осреднении m последних значений ряда. Метод носит название скользящих или движущихся средних. Величина m подбирается эмпирически. Скользящее среднее вычисляется по формуле:
yˆ = 1 ∑y(t) . |
(17) |
N
m t=N −m+1
Этот метод, наряду с другими, применяется при прогнозе котировок акций.
Наиболее часто для краткосрочного прогнозирования используются методы экспоненциального сглаживания, в частности, методы Брауна и метод Холта. Познакомимся с методами экспоненциального сглаживания на примере линейного метода Брауна. Этот метод не вырабатывает ошибок при прогнозировании значений линейно изменяющихся рядов. Оценка текущего значения y(t) осуществляется неявно, с помощью двойного сглаживания исходного ряда по формулам:
16
Для начала работы метода требуется значение u (0). Поскольку метод адаптируется к текущим значениям ряда, через некоторое количество шагов метод «забывает» начальные условия и можно выбрать u (0) = y (1). Значение константы 0 < α < 1 подбирается эмпирически.Еслиαблизкокединице,сглаживаниеслабоеивоценке большиевесаимеютпоследниеданные.Еслиαблизкокнулю,сглаживание сильное и оценка приближается к обычному среднему.
Выполним прогноз значений того же самого ряда линейным методом Брауна.
Промежуточные расчеты по модели линейного экспоненциального сглаживания Брауна при значении параметра α = 0,4 (соответственно, 1 – α = 0,6) представлены в таблице 2, где
Et = yt − yt — ошибка (остаток) аппроксимации.
Таблица 3
Расчеты экспоненциально сглаженных средних по линейному методу Брауна
t |
y |
t |
u |
t |
u' |
t |
2u |
t |
– u' |
E |
t |
E 2 |
u |
t |
– u' |
t |
a |
1 |
(t) |
|
|
|
|
|
t |
|
t |
|
|
|
|
||||||||
– |
– |
25 |
25 |
|
25 |
– |
|
– |
|
|
– |
|
|
– |
|||||
1 |
25 |
25 |
25 |
|
25 |
0 |
|
0 |
|
|
0 |
|
|
0 |
|||||
2 |
25 |
25 |
25 |
|
25 |
0 |
|
0 |
|
|
0 |
|
|
0 |
|||||
3 |
34 |
28,6 |
26,4 |
30,8 |
3,2 |
10,24 |
|
|
2,2 |
|
1,47 |
||||||||
4 |
39 |
32,8 |
29 |
36,6 |
2,4 |
5,76 |
|
|
3,8 |
|
2,53 |
||||||||
5 |
42 |
36,5 |
32 |
|
41 |
1 |
|
1 |
|
|
4,5 |
|
|
3 |
|||||
6 |
37 |
36,7 |
33,9 |
39,5 |
–2,5 |
6,25 |
|
|
2,8 |
|
1,87 |
||||||||
7 |
45 |
40 |
36,3 |
43,7 |
1,3 |
1,69 |
|
|
3,7 |
|
2,47 |
||||||||
8 |
44 |
41,6 |
38,4 |
44,8 |
–0,8 |
0,64 |
|
|
3,2 |
|
2,13 |
||||||||
9 |
55 |
47 |
41,8 |
52,2 |
2,8 |
7,84 |
|
|
5,2 |
|
3,47 |
||||||||
10 |
52 |
49 |
44,7 |
53,3 |
–1,3 |
1,69 |
|
|
4,3 |
|
2,87 |
||||||||
17
|
t |
y |
t |
u |
t |
u' |
t |
2u |
t |
– u' |
E |
t |
E 2 |
u |
t |
– u' |
t |
a |
1 |
(t) |
|
|
|
|
|
|
t |
|
t |
|
|
|
|
||||||||
|
11 |
67 |
56,2 |
49,3 |
63,1 |
3,9 |
15,21 |
|
|
6,9 |
|
4,6 |
||||||||
|
12 |
63 |
58,9 |
53,1 |
64,7 |
–1,7 |
2,89 |
|
|
5,8 |
|
3,87 |
||||||||
|
13 |
73 |
64,5 |
57,7 |
71,3 |
1,7 |
2,89 |
|
|
6,8 |
|
4,53 |
||||||||
|
|
|
|
|||||||||||||||||
|
14 |
67 |
65,5 |
60,8 |
70,2 |
–3,2 |
10,24 |
|
|
4,7 |
|
3,13 |
||||||||
|
|
|
|
|||||||||||||||||
|
15 |
70 |
67,3 |
63,4 |
71,2 |
–1,2 |
1,44 |
|
|
3,9 |
|
2,6 |
||||||||
|
|
|
|
|||||||||||||||||
|
16 |
80 |
72,4 |
67 |
77,8 |
2,2 |
4,84 |
|
|
5,4 |
|
3,6 |
||||||||
|
|
|
|
|||||||||||||||||
|
17 |
81 |
75,8 |
70,5 |
81,1 |
–0,1 |
0,01 |
|
|
5,3 |
|
3,53 |
||||||||
|
|
|
|
|||||||||||||||||
|
18 |
90 |
81,5 |
74,9 |
88,1 |
1,9 |
3,61 |
|
|
6,6 |
|
4,4 |
||||||||
|
|
|
|
|||||||||||||||||
|
19 |
88 |
84,1 |
78,6 |
89,6 |
–1,6 |
2,56 |
|
|
5,5 |
|
3,67 |
||||||||
|
|
|
|
|||||||||||||||||
Итого |
190 |
1077 |
988,4 |
907,8 |
1069 |
8 |
|
78,8 |
80,6 |
|
53,74 |
|||||||||
Средние |
10 |
56,7 |
52,02 |
47,78 |
56,26 |
0,42 |
4,15 |
4,24 |
|
2,83 |
||||||||||
Проверим адекватность адаптивной модели.
Расчеты для проверки адекватности представлены в таблице 3. В результате расчетов получены следующие оценки.
Случайность Да |
g = 13 > 9 |
|
|
|
Независимость Да |
d = 2,6 |
4 – d = 1,401 > 1,40 |
|
|
Нормальность Да |
S = 2,09 |
E |
= –3,2 |
RS = 3,394 |
Emin |
= 3,9 |
3,14 < RS < 4,43 |
||
|
|
max |
|
|
Модель является несмещенной, поскольку
,
где 2,1 — значение t-статистики для вероятности ошибки p = 0,05 и числе степеней свободы: n – 1 = 18.
Точность модели: MSE = 4,15; МАРЕ = 3,21. Доверительный интервал U (1) = 4,9 при р = 0,05.
Результат прогноза yТ (20) = y (19) + (α / l – α) · (u (19) u/ (19)) · l
будет:
t |
yT (t) |
Нижняя гр. |
Верхняя гр. |
20 |
93,3 |
88,4 |
98,2 |
18
2.3. Регрессионная модель прогнозирования
В экономическом моделировании и прогнозировании часто требуется давать прогноз с учетом ведущих факторов. Эти задачи решают эконометрические методы и, в частности, методы регрессионного анализа. Пусть на интересующий нас процесс оказывает влияние один фактор. В этом случае между независимой (факторной)переменнойхизависимой(результативной)переменнойубудет существовать корреляционная связь. Если предположить, что форма связи между переменными линейная, а значения величин х и у распределены по нормальному закону, то тесноту связи между этими переменными можно оценить по следующей формуле:
(18)
Значение коэффициента корреляции изменяется в интервале от –1 до 1. Чем ближе значение r по абсолютной величине к 1, тем сильнее связь между переменными. Напротив, чем ближе значение r к 0, тем эта связь слабее. При значении |r| = 1 связь между переменной является функциональной.
Если коэффициент корреляции показывает существование достаточно тесной связи между факторной и результативной переменными, имеется возможность использовать для прогнозирования модель линейной регрессии.
Линейная модель влияния фактора x на процесс y запишется в виде:
y (t) = a0 + a1 x (t)
Оценку параметров модели регрессии обычно выполняется методом наименьших квадратов.
19
(19)
Прогнозную оценку переменной у получают подстановкой в модель соответствующего прогнозного значения х для того же момента или периода наблюдений.
Выполним расчеты регрессионной модели и оценим ее точность. Воспользуемся тем же примером, дополненным значениями фактора x:
Пример 2.
Имеются данные относительно среднеквартальной цены корпоративной акции и среднеквартального индекса фондового рынка за девятнадцать периодов наблюдений. Требуется оценить параметры линейной зависимости цены акции у от индекса фондового рынка х. Исходные данные приведены в таблице 4.
Таблица 4
Исходные данные для построения прогностической модели зависимости цены корпоративной акции (у.е.) от индекса фондового рынка (пунктов)
Номер |
Цена акции, yt |
Индекс фондового |
|
наблюдения, t |
рынка, xt |
||
|
|||
1 |
25 |
45 |
|
2 |
25 |
45 |
|
3 |
34 |
47 |
|
4 |
39 |
46 |
|
5 |
42 |
50 |
|
6 |
37 |
48 |
|
7 |
45 |
48 |
|
8 |
44 |
50 |
|
9 |
55 |
54 |
|
10 |
52 |
55 |
|
11 |
67 |
57 |
|
12 |
63 |
56 |
20
