
Metodichka_lab_Ekonometria
.pdfЕсли наиболее высоким окажется коэффициент автокорреляции порядка k , то ряд содержит циклические колебания с периодичностью в k моментов времени.
Если ни один из коэффициентов автокорреляции не является значимым, то можно сделать одно из двух предположений относительно структуры этого ряда: либо ряд не содержит тенденции и сезонных колебаний, либо содержит сильную нелинейную тенденцию, для выявления которой нужно провести дополнительные исследования.
Поэтому коэффициент автокорреляции уровней и автокорреляционную функцию целесообразно использовать для выявления во временном ряде наличия или отсутствия трендовой компоненты Q(t ) и
сезонной компоненты S (t ).
Замечание. Временной ряд нельзя исследовать как регрессионную модель методом наименьших квадратов. В регрессионной модели считается, что значение объясняемой переменой получается, как мгновенная реакция на изменение объясняющей переменной. Во временных рядах значение объясняемой переменной получается, как правило, с опозданием на некоторое время, которое называется лагом. Существенным является и то, что кроме тенденции временные ряды имеют и сезонные колебания. Если же описывать временной ряд регрессионной моделью, то вместе с ошибками мы погасим и колебания.
Рассмотрим пример выполнения заданий лабораторной работы 8.
9.2 Методы сглаживания временных рядов
Собраны сведения о потреблении электроэнергии жильцами 107-го дома по ул. Партизанской г. Донецк (кВт/час), yt за три года (январь 2009
г. – декабрь 2011 г.). Данные внесены в табл. 9.1.
Задание 1
1.1.Сгладить ряд с помощью скользящей средней.
1.2.Сгладить ряд с помощью экспоненциального сглаживания.
1.3.Сделать прогноз потреблении электроэнергии на январь 2012 г.
Таблица 9.1
Потребление электроэнергии
t |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
y |
153 |
136 |
153 |
158 |
151 |
178 |
174 |
178 |
211 |
233 |
234 |
202 |
t |
||||||||||||
t |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
y |
174 |
153 |
168 |
170 |
159 |
181 |
174 |
182 |
218 |
246 |
253 |
202 |
t |
||||||||||||
t |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
yt |
180 |
155 |
169 |
180 |
171 |
203 |
198 |
210 |
236 |
274 |
280 |
232 |
119

Решение. Скользящие средние определим по формуле (9.4), в которой выбираем интервал (лаг) k =12 .
1.Вносим данные в Excel в виде столбца Дата, Порядковый номер t и столбца yt (табл. 9.2).
2.Строим график yt с помощью Excel, Диаграмма (рис. 9.1).
электроэнергии |
300 |
|
|
|
|
250 |
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
|
|
Потребление |
150 |
|
|
|
|
100 |
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
|
|
|
Месяцы |
|
|
|
|
|
|
|
|
|
|
|
|
Исходные данные |
|
|
Рис. 9.1 – График исходных данных
Таблица 9.2
Исходные данные и скользящие средние
Дата |
t |
yt |
Прогноз, Ft |
Остатки, et |
1 |
2 |
3 |
4 |
5 |
01.01.09 |
1 |
153 |
#Н/Д |
#Н/Д |
01.02.09 |
2 |
136 |
#Н/Д |
#Н/Д |
01.03.09 |
3 |
153 |
#Н/Д |
#Н/Д |
01.04.09 |
4 |
158 |
#Н/Д |
#Н/Д |
01.05.09 |
5 |
151 |
#Н/Д |
#Н/Д |
01.06.09 |
6 |
178 |
#Н/Д |
#Н/Д |
01.07.09 |
7 |
174 |
#Н/Д |
#Н/Д |
01.08.09 |
8 |
178 |
#Н/Д |
#Н/Д |
01.09.09 |
9 |
211 |
#Н/Д |
#Н/Д |
01.10.09 |
10 |
233 |
#Н/Д |
#Н/Д |
01.11.09 |
11 |
234 |
#Н/Д |
#Н/Д |
01.12.09 |
12 |
202 |
179,8 |
#Н/Д |
01.01.10 |
13 |
174 |
181,5 |
#Н/Д |
01.02.10 |
14 |
153 |
183,0 |
#Н/Д |
01.03.10 |
15 |
168 |
184,2 |
#Н/Д |
01.04.10 |
16 |
170 |
185,3 |
#Н/Д |
120

|
|
|
|
Продолжение таблицы 9.2 |
||
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
|
5 |
|
01.05.10 |
17 |
159 |
186,0 |
|
#Н/Д |
|
01.06.10 |
18 |
181 |
186,2 |
|
#Н/Д |
|
01.07.10 |
19 |
174 |
186,2 |
|
#Н/Д |
|
01.08.10 |
20 |
182 |
186,5 |
|
#Н/Д |
|
01.09.10 |
21 |
218 |
187,1 |
|
#Н/Д |
|
01.10.10 |
22 |
246 |
188,2 |
|
#Н/Д |
|
01.11.10 |
23 |
253 |
189,8 |
|
30,4 |
|
01.12.10 |
24 |
202 |
189,8 |
|
29,9 |
|
01.01.11 |
25 |
180 |
190,3 |
|
30,0 |
|
01.02.11 |
26 |
155 |
190,5 |
|
30,4 |
|
01.03.11 |
27 |
169 |
190,5 |
|
30,7 |
|
01.04.11 |
28 |
180 |
191,3 |
|
30,6 |
|
01.05.11 |
29 |
171 |
192,3 |
|
30,2 |
|
01.06.11 |
30 |
203 |
194,2 |
|
30,3 |
|
01.07.11 |
31 |
198 |
196,2 |
|
30,1 |
|
01.08.11 |
32 |
210 |
198,5 |
|
30,2 |
|
01.09.11 |
33 |
236 |
200,1 |
|
30,7 |
|
01.10.11 |
34 |
274 |
202,5 |
|
33,1 |
|
01.11.11 |
35 |
280 |
204,7 |
|
35,1 |
|
01.12.11 |
36 |
232 |
207,2 |
|
35,6 |
|
3. Выбираем Анализ данных – Скользящие средние, ОК. Вносим требуемые данные в диалоговое окно, рис. 9.2.
Рис. 9.2 – Диалоговое окно ввода данных скользящих средних
121

Получили столбцы Прогноз, Остатки (табл. 9.2), графики |
||||||
фактических и прогнозных значений (рис. 9.3). |
|
|
||||
электроэнергии |
300 |
|
|
|
|
|
250 |
|
|
|
|
||
200 |
|
|
|
|
||
Потребление |
150 |
|
|
|
|
|
100 |
|
|
|
|
||
0 |
10 |
20 |
30 |
40 |
||
|
||||||
|
|
|
Месяцы |
|
|
|
|
|
|
Фактический |
Прогноз |
|
Рис. 9.3 – Скользящие средние из Анализа данных
Остатки вычисляются по формуле
|
|
1 |
t |
|
|
et = |
∑ (yi − Fi )2 , 23 ≤t ≤ 36 . |
(9.10) |
|||
12 |
|||||
|
i=t−11 |
|
|||
|
|
|
|
Остатки определяются как среднеквадратическое отклонение уровней ряда от сглаженных.
4. Скользящие средние используем для прогноза. Для этого в столбце Прогноз выделяем 12 последних значений и определяем среднее их значение, получаем y37прогноз = y01.01.12прогноз =196,722 (СРЗНАЧ).
5. Ряд, тренд, прогнозные значения и прогноз на будущий месяц можно изобразить графически с помощью Excel, Диаграмма (рис. 9.4).
С помощью остатков можно оценивать качество прогноза и находить доверительные интервалы прогноза, но при этом будем иметь очень большую погрешность, что видно из рис. 9.4.
122

электроэнергии |
300 |
|
|
|
|
|
250 |
|
|
|
|
||
200 |
|
|
|
|
||
|
|
|
|
|
||
Потребление |
150 |
|
|
|
|
|
100 |
|
|
|
|
||
0 |
10 |
20 |
30 |
40 |
||
|
||||||
|
|
|
Месяцы |
|
|
|
|
|
Исходный |
Прогноз |
Прогноз на будущее |
Рис. 9.4 – Скользящие средние
Экспоненциальное сглаживание осуществляется по формуле (9.5).
1.Вносим данные в Excel в виде столбца Дата, Порядкового номера t
истолбца yt .
2.Строим график yt с помощью Диаграммы (рис. 9.1).
3.Выбираем Анализ данных, Экспоненциальное сглаживание,
ОК. Вносим требуемые данные в диалоговое окно (рис. 9.5).
Рис. 9.5 – Диалоговое окно ввода данных экспоненциального сглаживания
Получили столбцы Прогноз, Остатки, графики фактических и прогнозных значений, табл. 9.3.
123

Таблица 9.3
Исходные данные и экспоненциально слаженные уровни
Дата |
t |
y |
Прогноз, F |
Остатки, e |
|
|
t |
t |
t |
|
|
|
|
#Н/Д |
01.01.09 |
1 |
153 |
#Н/Д |
|
01.02.09 |
2 |
136 |
153 |
#Н/Д |
01.03.09 |
3 |
153 |
147,9 |
#Н/Д |
01.04.09 |
4 |
158 |
149,43 |
#Н/Д |
01.05.09 |
5 |
151 |
152,001 |
11,37915 |
01.06.09 |
6 |
178 |
151,7007 |
5,786677 |
01.07.09 |
7 |
174 |
159,5905 |
15,98020 |
01.08.09 |
8 |
178 |
163,9133 |
17,32329 |
01.09.09 |
9 |
211 |
168,1393 |
19,12870 |
01.10.09 |
10 |
233 |
180,9975 |
27,34413 |
01.11.09 |
11 |
234 |
196,5983 |
39,74807 |
01.12.09 |
12 |
202 |
207,8188 |
44,49787 |
01.01.10 |
13 |
174 |
206,0732 |
37,13490 |
01.02.10 |
14 |
153 |
196,4512 |
28,64399 |
01.03.10 |
15 |
168 |
183,4158 |
31,36109 |
01.04.10 |
16 |
170 |
178,7911 |
32,42603 |
01.05.10 |
17 |
159 |
176,1538 |
27,09821 |
01.06.10 |
18 |
181 |
171,0076 |
14,24995 |
01.07.10 |
19 |
174 |
174,0053 |
12,53505 |
01.08.10 |
20 |
182 |
174,0037 |
11,46152 |
01.09.10 |
21 |
218 |
176,4026 |
7,38890 |
01.10.10 |
22 |
246 |
188,8818 |
24,45596 |
01.11.10 |
23 |
253 |
206,0173 |
41,05593 |
01.12.10 |
24 |
202 |
220,1121 |
48,99048 |
01.01.11 |
25 |
180 |
214,6785 |
43,96176 |
01.02.11 |
26 |
155 |
204,2749 |
35,29882 |
01.03.11 |
27 |
169 |
189,4924 |
36,32567 |
01.04.11 |
28 |
180 |
183,3447 |
36,74486 |
01.05.11 |
29 |
171 |
182,3413 |
30,87149 |
01.06.11 |
30 |
203 |
178,9389 |
13,65958 |
01.07.11 |
31 |
198 |
186,1572 |
15,47846 |
01.08.11 |
32 |
210 |
189,7101 |
16,81083 |
01.09.11 |
33 |
236 |
195,7970 |
19,41536 |
01.10.11 |
34 |
274 |
207,8579 |
26,88376 |
01.11.11 |
35 |
280 |
227,7006 |
46,19788 |
01.12.11 |
36 |
232 |
243,3904 |
53,93293 |
124

Здесь |
F2 = y1 =153 , |
|
|
(9.11) |
|||
|
|
|
|||||
|
Ft |
= 0,3yt−1 + 0,7Ft−1, 3 ≤t ≤36 . |
(9.12) |
||||
Остатки вычисляются по формуле |
|
|
|||||
|
|
|
1 |
t |
− Fi )2 , 5 |
|
|
|
et |
= |
∑ (yi |
≤ t ≤ 36 . |
(9.13) |
||
|
|
|
4 |
i=t−3 |
|
|
|
|
|
|
|
|
|
|
|
электроэнергии |
300 |
|
|
|
|
|
|
250 |
|
|
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Потребление |
150 |
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
0 |
|
10 |
20 |
30 |
40 |
||
|
|
|
|
|
Месяцы |
|
|
|
|
|
|
Фактический |
Прогноз |
|
Рис. 9.6 – Экспоненциальное сглаживание из Анализа данных
4.Для получения прогноза продолжаем столбец прогнозных значений ещё на одну клетку и получаем y37прогноз = 239,97 .
5.Ряд, тренд, прогнозные значения и прогноз на будущий месяц можно изобразить графически с помощью Excel, Диаграмма (рис. 9.7).
Можно, с помощью остатков, оценивать качество прогноза и находить доверительные интервалы прогноза, но при этом будем иметь очень большую погрешность, что видно и из рис 9.7.
9.3Выявление тренда во временном ряде и построение кривых
роста Задание 2
2.1Определить наличие тренда в исследуемом ряде.
125

2.2 Найти кривую роста. |
|
|
|
||
2.3 Сделать прогноз на следующий квартал. |
|
|
|||
электроэнергии |
300 |
|
|
|
|
250 |
|
|
|
|
|
200 |
|
|
|
|
|
|
|
|
|
|
|
Потребление |
150 |
|
|
|
|
100 |
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
|
|
|
|
|
|
|
|
|
|
Месяцы |
|
|
|
Исходный |
Тренд |
Прогноз |
Прогноз на будущее |
Рис. 9.7 – Скользящие средние
2.4 Найти доверительные интервалы прогноза.
Решение. 1. Ряд разбиваем на две приблизительно равные части по 18 наблюдений. Проверяем гипотезу о равенстве дисперсий обоих частей ряда с помощью F-критерия Фишера. Для этого бóльшую дисперсию делим на меньшую и сравниваем с табличным значением при уровне значимости 5 %:
Fрасч. = 1389,794768,8824 =1,807 < Fтабл. = Fα; n1 −1; n2 −1=F0,05; 17; 17 = 2,28. (9.14)
Дисперсии в выбранных частях ряда различаются незначимо. С вероятностью 0,95 расхождение между ними есть величина случайная.
2. Наличие тренда в исследуемом ряде определяем с помощью критерия
|
|
|
|
1 − |
|
2 |
|
|
|
n n (n + n − 2) |
|
|
|
||||
tнабл. = |
|
|
y |
y |
= 3,03 |
, |
(9.15) |
||||||||||
|
|
|
|
|
|
|
|
1 |
2 |
1 |
2 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
(n1 |
−1)Sy2 |
+ (n2 |
−1)Sy2 |
|
n1 + n2 |
|
|
|
||||||||
|
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
tтабл. = tα,n−2 = t0,05; 34 = 2,03 ,
126

tнабл. > tтабл. − принимаем гипотезу о наличии тренда. 3. Находим линию тренда (рис. 9.8)
yˆt = 1,9371t + 156,64 . |
(9.16) |
4. По лини тренда находим прогнозные значения yˆt , подставляя
вместо t : 37, 38,39,40, табл. 9.4. 5. По формуле
U |
|
= yˆ |
|
±U |
|
= yˆ |
|
± S t 1+ |
1 |
+ |
3(n + 2k −1)2 |
(9.17) |
y |
n+k |
k |
n+k |
|
n(n2 −1) |
|||||||
|
|
|
|
e α |
n |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
находим доверительные интервалы прогноза.
Стандартную ошибку Se находим с помощью: надстройки «Анализ данных», функция «Регрессия». Se =30,734 , tα, n−2 = t0,05, 34 =2,03 . Uk −
полуширина доверительных интервалов.
Результаты вычислений вносим в табл. 9.4, графическая иллюстрация, рис. 9.8.
Таблица 9.4
Прогноз и доверительный интервал по кривой роста
Прогноз
t |
Uk |
yˆt –Uk |
yˆt |
yˆt +Uk |
37 |
65,906 |
162,407 |
228,313 |
294,218 |
38 |
66,194 |
164,056 |
230,250 |
296,444 |
39 |
66,496 |
165,691 |
232,187 |
298,683 |
40 |
66,811 |
167,313 |
234,124 |
300,936 |
9.4 Выявление сезонных колебаний во временном ряде. Построение аддитивной и мультипликативной моделей временного ряда с сезонными колебаниями.
Задание 3
3.1Проверить наличие сезонной составляющей во временном ряде.
3.2Выделить сезонные составляющие.
3.3Найти линию тренда после удаления сезонной составляющей.
3.4Найти прогноз на следующий квартал с учётом сезонной составлящей.
127

электроэнергии |
300 |
|
|
|
|
|
250 |
|
|
|
|
|
|
200 |
|
|
|
|
|
|
Потребление |
150 |
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
|
|
|
|
Месяцы |
|
|
|
|
Исходный |
|
Тренд |
|
|
|
|
Прогноз |
|
Доверительныйинтервал |
|
||
|
Доверительныйинтервал |
|
|
|
|
|
Рис. 9.8 – Тренд, прогнозные значения по кривой роста и доверительный |
||||||
|
|
интервал |
|
|
|
|
Решение. Наличие сезонной составляющей определяем с помощью |
||||||
коррелограммы. |
|
|
|
|
|
|
Вычисляем |
коэффициенты |
автокорреляции |
порядков |
1–33. |
||
Коэффициенты корреляции находим с помощью Функция, Корреляция, |
||||||
табл. 9.5. Коэффициенты автокорреляции |
вычисляются по формуле (9.8), |
|||||
схема (9.9). |
|
|
|
|
|
Таблица 9.5
Значения коэффициента автокорреляции
τ |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
r |
–0,5 |
–0,5 |
–0,5 |
–0,1 |
–0,4 |
–0,4 |
–0,4 |
–0,4 |
–0,1 |
0,3 |
0,7 |
τ |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
r |
1,0 |
0,8 |
0,4 |
0,1 |
–0,2 |
–0,3 |
–0,4 |
–0,4 |
–0,4 |
–0,2 |
0,3 |
τ |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
r |
0,8 |
1,0 |
0,8 |
0,7 |
0,6 |
0,7 |
0,6 |
0,2 |
0,4 |
–0,6 |
–0,6 |
Строим коррелограмму, рис. 9.9.
128