Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Анализ временных рядов и прогнозирование. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
2.
9663
,
24072
6784
,
105
2878
,
23967
=
+
()()(
)
−+=+−⋅+−+⋅+=−=∆
εεε
ε
3.
4.
p
xx
tt
p
y
t
yy
tt
13
p
2
=++=
=+=
;
2522,1431318,2491080,374458,21308 =+=
;
2878,23967119178,17211109,1160457,1767
;
6784,1052522,143738,010*322,5
5. Точечный прогноз фактического значения:
д.е.
7.1.2. Метод последовательных разностей
В ряде случаев вместо аналитического выравнивания вре­менного ряда с целью устранения тенденции можно применить более простой метод – метод последовательных разностей.
Если временной ряд содержит ярко выраженную линей­ную тенденцию, ее можно устранить путем замены исходных уровней ряда цепными абсолютными приростами (первыми раз­ностями). Пусть:
p
(7.3)
ε
+=
t
tt
(7.4)
tbay
+=
где
ε
случайная ошибка;
t
yy
p
t
Тогда
1
btbatbayy
ttttttt
(7.5)
111
Коэффициент b – константа, которая не зависит от време­ни. При наличии сильной линейной тенденции остатки
ε
доста-
t
точно малы и в соответствии с предпосылками МНК носят слу­чайный характер. Поэтому первые разности уровней ряда
не
t
зависят от переменной времени, их можно использовать для дальнейшего анализа.
Если временной ряд содержит тенденцию в форме парабо­лы второго порядка, то для ее устранения можно заменить ис­ходные уровня ряда на вторые разности.
Пусть имеет место соотношение (7.4), однако
p
t
2
. (7.6)
tbtbay
++=
21
– – 111
Тогда
()(
)
(
)
(
)
−−+⋅⋅+−=∆−∆=∆
ε
ε
2
211
2
( )
tbbb
++=
εε
1221
tt
21
Как показывает это соотношение, первые разности
2
11
++++++==
tbtbatbtbayy
εε
1
ttttt
(7.7)
не-
t
посредственно зависят от фактора времени t и, следовательно, содержат тенденцию.
Определим вторые разности:
22
tbbb
( ) ( )( ) ( )
12211
ttttt
2212
btbbb
++=++
(7.8)
εεεεε
21221221
ttttt
Очевидно, что вторые разности
2 не содержат тенденции,
t
поэтому при наличии в исходных уровнях тренда в форме пара­болы второго порядка их можно использовать для дальнейшего анализа. Если тенденции временного ряда соответствует экспо­ненциальный или степенной тренд, метод последовательных разностей следует применять не к исходным уровням ряда, а к их логарифмам.
Пример 7.2. Изучение зависимости расходов на конечное потребление от совокупного дохода с помощью метода последо­вательных разностей.
Обратимся вновь к данным о расходах на конечное по­требление y
и совокупном доходе xt (табл.7.1).Проанализируем
t
зависимость между этими рядами, используя вторые разности, такой выбор объясняется наличием в исходных уровнях тренда в форме параболы второго порядка (табл.7.4).
Проверка временных рядов вторых разностей на автокор­реляцию дает следующие результаты: r
= –0,333; r
ty
= –0,006.
tx
Поскольку полученные ряды не содержат автокорреляции, мы можем использовать их вместо исходных данных для измерения зависимости между расходами на конечное потребление и сово­купным доходом. Коэффициент корреляции рядов вторых раз­ностей составляет r
= 0,827. В данном случае можно под-
tytx
твердить вывод о наличии тесной прямой связи между расхода­ми на конечное потребление и совокупным доходом, приведен-
– – 112
ный в примере 7.1.
Таблица 7.4. Вторые разности временных рядов
t yt хt
1 1746,66 1776,02 – 2 2853,53 2908,00 – 3 3870,68 3983,76 –89,72 –56,23 4 5214,38 5324,52 326,56 265,01 5 6710,06 6829,30 151,97 164,02 6 8644,39 8885,61 438,65 551,53 7 10780,18 10976,29 201,46 34,37 8 13614,23 13818,98 698,26 752,00 9 16706,07 17289,94 257,79 628,28
10 20489,59 21308,58 691,68 547,68
∆′′ty ∆′′tx
Построение уравнения регрессии зависимости расходов на конечное потребление от совокупных доходов по вторым разно­стям привело к следующим результатам:
Показатель Значение
a
0
a
1 2
R
r
tytx
65,5594
0,7456 0,6835 0,8267
Уравнение регрессии приняло вид:
p
t
.
xy
+= 746,05594,65)(
t
В отличие от уравнения регрессии по отклонениям от тренда, параметрам данного уравнения легко дать интерпрета­цию. При изменении прироста дохода на 1 д.е. прирост потреб­ления изменится в среднем на 0,746 д.е. в том же направлении.
При всей своей простоте метод последовательных разно­стей имеет два существенных недостатка. Во-первых, его при­менение сильно связано с сокращением числа пар наблюдений, по которым строится уравнение регрессии, а, следовательно, с потерей числа степеней свободы. Во-вторых, использование
– – 113
вместо исходных уровней временных рядов их приростов, или ускорений, приводит к потере информации, содержащейся в ис­ходных данных. Следовательно, данное расчетное уравнение не позволяет нам сделать прогноз на дальнейшие периоды.
7.2. Методы исключения фактора времени
7.2.1.Представление исходного ряда в виде тренда
Рассмотрим два временных ряда y(t), x(t) и будем изучать зависимость между ними. Особенность этой задачи состоит в том, что обе переменные являются случайными функциями вре­мени t. Как функции от t они могут изменяться в очень согласо­ванной форме. Поэтому их совместное рассмотрение и попытка построить зависимость вида
= b0 + b1 x(t) + ε(t) (7.9)
y(t)
могут привести к неверным результатам и выводам.
Отсюда следует, что зависимость (7.9) или ее аналоги на­прямую применять нельзя. Один из способов проверки зависи­мости между y(t) и x(t) состоит в том, чтобы устранить влияние основной переменной t. Для этого поступают следующим обра­зом.
Будем считать, что y(t) и x(t) описываются трендовыми мо­делями, а именно:
y(t) = g1(t) + x(t) = g2(t) +
ε
(t); (7.10)
1
ε
(t). (7.11)
2
Найдем тренды g1(t) и g2(t) и введем новые переменные (устранение трендов):
x1(t) = x(t) – g2(t); y1 = y(t) – g1(t). (7.12)
Теперь влияние переменной t в значительной степени уст­ранено, и можно изучать линейную зависимость
(t) = b0 + b1 x1 + ε (7.13)
y
1
Для анализа зависимости (7.13), ее значимости или незна­чимости применяются стандартные методики, рассмотренные
– – 114
ранее.
В приложениях используются и другие подходы, позво­ляющие изучать зависимость между временными рядами y(t) и x(t). Например, в зависимости (7.9) явным образом учитывают время. При таком подходе вместо формулы (7.1) рассматривают ее обобщение:
y(t) = b0 + b1 x(t) + b2 t + ε(t), (7.14)
где время t выступает как еще одна объясняющая переменная.
Пример 7.3. Имеется набор данных, отражающих измене- ние во времени следующих переменных: y = y(t) – объем реали­зации предприятия оптовой торговли (у.е.); x = x(t) – размер тор­говой площади (м2); t – время (квартал). Требуется найти линей­ные тренды в динамике x(t)
и y(t), установить или опровергнуть
зависимость изменения объема реализации от размера торговой площади. Исходные данные представлены в табл.7.5.
Таблица 7.5. Зависимость объема реализации и торговой площа-
ди от времени
t xt yt
1 127,2 15,12 2 129 16,62 3 131,2 17,16 4 134,8 17,63 5 136,6 18,14 6 138,6 16,92 7 143,2 17,66 8 144,3 18,4
9 145,6 19,97 10 148,7 21,35 11 150,4 20,16 12 154,7 21,04 13 158,2 20,27 14 159,3 23,33 15 159,6 21,67 16 163,6 23,05 17 164,9 23,9 18 169,4 23,85
– – 115
Будем считать, что выполнены все предположения МНК
относительно случайных составляющих
ε
(t),
ε
1
(t) и ε, входящих
2
в отношения (7.10)-(7.11). Уровень значимости выберем 0,05.
Решение.
1) Предположим, что уравнения (7.10) и (7.11) содержат
линейные тренды:
g1(t) = a0 + a1t, g2(t) = c0 + c1t.
Вычисляя оценки параметров этих трендов, получаем, что
a0 = 15,16; a1 = 0,49; c0 = 124,5; c1 = 2,45.
Сформируем новые переменные x1 и y1 по формуле (7.12) и
сведем все результаты в табл.7.6.
Таблица 7.6. Исходные и расчетные значения переменных
t xt yt g1 g2 x1 y1
1 127,2 15,12 15,64778 126,9468 0,253216 –0,52778 2 129 16,62 16,13523 129,3929 –0,39291 0,484771 3 131,2 17,16 16,62268 131,839 –0,63904 0,53732 4 134,8 17,63 17,11013 134,2852 0,514826 0,519869 5 136,6 18,14 17,59758 136,7313 –0,1313 0,542418 6 138,6 16,92 18,08503 139,1774 –0,57743 –1,16503 7 143,2 17,66 18,57248 141,6236 1,576436 –0,91248 8 144,3 18,4 19,05993 144,0697 0,230306 –0,65993
9 145,6 19,97 19,54739 146,5158 –0,91582 0,422614 10 148,7 21,35 20,03484 148,962 –0,26195 1,315163 11 150,4 20,16 20,52229 151,4081 –1,00808 –0,36229 12 154,7 21,04 21,00974 153,8542 0,845786 0,030261 13 158,2 20,27 21,49719 156,3003 1,899656 –1,22719 14 159,3 23,33 21,98464 158,7465 0,553526 1,345359 15 159,6 21,67 22,47209 161,1926 –1,5926 –0,80209 16 163,6 23,05 22,95954 163,6387 –0,03873 0,090458 17 164,9 23,9 23,44699 166,0849 –1,18486 0,453007 18 169,4 23,85 23,93444 168,531 0,869006 –0,08444
По данным табл.7.6 получим зависимость (7.13): b0 = –
2,81015; b1 = –0,176. Величина F-критерия равна F
= 0,768.
расч
Табличное значение F-критерия с 1 и 16 степенями свободы и α
=
0,05 равно F
= 4,49. Так как F
табл
расч
< F
, зависимость (7.13)
табл
является незначимой. Следовательно, между переменными x1 и
– – 116
y1 нет линейной зависимости, и отсутствует линейная связь ме­жду объемом реализации и размером торговой площади. Обе переменные растут относительно согласованно, но на y(t) влия­ют какие-то другие, не учтенные факторы.
2) Для получения модели (7.14), рассматривающей время как дополнительную влияющую переменную, достаточно при­менить к исходным данным традиционный МНК (или процеду­ру Регрессия из Анализа данных). Результаты расчетов пред­ставлены в табл.7.7.
Таблица 7.7. Результаты расчетов по МНК
t xt yt y
1 127,2 15,12 15,6031 -0,4831 2 129 16,62 16,20456 0,415443 3 131,2 17,16 16,73544 0,424563 4 134,8 17,63 17,01929 0,610709 5 136,6 18,14 17,62075 0,51925 6 138,6 16,92 18,18692 -1,26692 7 143,2 17,66 18,29433 -0,63433 8 144,3 18,4 19,0193 -0,6193
9 145,6 19,97 19,70898 0,26102 10 148,7 21,35 20,08106 1,268942 11 150,4 20,16 20,70016 -0,54016 12 154,7 21,04 20,8605 0,179498 13 158,2 20,27 21,162 -0,892 14 159,3 23,33 21,88697 1,443027 15 159,6 21,67 22,7531 -1,0831 16 163,6 23,05 22,96638 0,083623 17 164,9 23,9 23,65606 0,243941 18 169,4 23,85 23,78111 0,068889
19 170,98 24,4219 20 173,42 24,90935
r
t
yt – y
r
t
Значения коэффициентов регрессии представлены в сле-
дующей таблице:
Коэффициенты
b
0
b1 –0,17645 b2 0,919063
37,12806
– – 117
Статистические параметры полученной модели – в сле-
дующей таблице:
df SS MS F
Регрессия 2 115,5893 57,79464 88,96513 Остаток 15 9,744489 0,649633 Итого 17 125,3338 Коэффициент детерминации R
2
0,9223
Полученные результаты говорят о том, что модель (7.14) адекватно описывает зависимость объема реализации от торго­вой площади и времени (высокие коэффициент детерминации и значение F-критерия, явно превышающее F табличное). Следо­вательно, уравнение
= 37,128 – 0,176 x(t) + 0,919 t
y(t)
можно использовать для прогнозирования объема реализации.
Прогноз на предстоящие два периода представлен в табл.7.7. При этом значения x(t) получены по уравнению (7.11).
7.2.2. Включение в модель регрессии фактора времени
В корреляционно-регрессионном анализе можно устранить воздействие какого-либо фактора, если зафиксировать воздейст­вие этого фактора на результат и другие включенные в модель факторы. Данный прием широко применяется в анализе времен­ных рядов, когда тенденция фиксируется через включение фак­тора времени в модель в качестве независимой переменной.
Модель вида
yt = a + b1 xt + b2 t +
ε
(7.15)
t
относится к группе моделей, включающих фактор времени. Очевидно, что число независимых переменных в такой модели может быть больше единицы. Кроме того, это могут быть не только текущие, но и лаговые значения независимой и результа­тивной переменных.
Преимущество данной модели перед методами отклонений от трендов и последовательных разностей состоит в том, что она позволяет учесть всю информацию, содержащуюся в исходных
– – 118
данных, поскольку значения yt и xt – это уровни исходных вре-
∑∑∑
менных рядов. Кроме того, модель строится по всей совокупно­сти данных за рассматриваемый период в отличие от метода по­следовательных разностей, который приводит к потере числа наблюдений. Параметры модели с включением фактора времени определяются обычным МНК. Покажем это на примере 7.4.
Пример 7.4. Построение модели регрессии с включением фактора времени.
Вернемся к данным табл.7.1. Построим уравнение регрес­сии, описывающее зависимость расходов на конечное потребле­ние yt от совокупного дохода xt и фактора времени t. Для расчета параметров уравнения регрессии (7.15) воспользуемся обычным МНК. Система нормальных уравнений имеет вид:
ytbxban
 
1
2
21
=++
tt
(7.16)
yxxtbxbxa
=++
2
2
21
ttttt
ytxbxtbta
=++
ttt
Решив эту систему относительно а, b1 и b2, находим:
а = –36,242; b1 = 0,919; b2 = 99,571.
Следовательно, уравнение регрессии имеет вид:
p
t
txy
ε
+++= 571,99919,0242,36
tt
Упростить расчет коэффициентов а, b1 и b2 можно с помо­щью процедуры Регрессия (в надстройке Анализ данных). Расчетные значения и их отклонения от исходных приведены в следующей таблице.
Интерпретация параметров этого уравнения следующая. Параметр b1 = 0,919 характеризует, что при увеличении сово­купного дохода на 1 д.е. расходы на конечное потребление воз­растут в среднем на 0,919 д.е. в условиях существования неиз­менной тенденции. Параметр b
= 99,571 означает, что воздейст-
2
вие всех факторов, кроме совокупного дохода, на расходы на конечное потребление приведет к его среднегодовому абсолют-
– – 119
ному приросту на 99,571 д.е.
хt t yt y
1776,02 1 1746,66 1694,65 52,01 2908,00 2 2853,53 2833,98 19,55 3983,76 3 3870,68 3921,66 –50,99 5324,52 4 5214,38 5252,76 –38,37 6829,30 5 6710,06 6734,51 –24,44
8885,61 6 8644,39 8722,86 –78,47 10976,29 7 10780,18 10742,78 37,40 13818,98 8 13614,23 13453,43 160,80 17289,94 9 16706,07 16741,17 –35,11 21308,58 10 20489,59 20531,98 –42,39
r
yt – yr
Применение полученного уравнения для прогноза на пред-
стоящие периоды предполагает знание о значениях зависимой переменной x
в эти периоды. Для этого можно воспользоваться
t
уравнением линейного тренда из табл.7.2.
Например, в нашем случае сделаем прогноз расходов на
конечное потребление для t = 11. Сначала определяем трендовое значение факторного признака:
p
t
2
=++= ttx
,
1348,249108889,1860364,331755,1933
затем прогнозное значение по уравнению:
p
y
t
23939,674211571,991348,24910919,0242,36 =++=
д.е.,
что практически не отличается от прогноза по отклонениям от тренда.
7.3. Автокорреляция в остатках. Критерий Дарбина-Уотсона
7.3.1. Определение критерия Дарбина-Уотсона
Рассмотрим уравнение регрессии вида
k
, (7.17)
++=
xbay
j
=
1
ε
tjtjt
где kчисло независимых переменных модели.
– – 120
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]