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

Metodichka_lab_Ekonometria

.pdf
Скачиваний:
85
Добавлен:
13.04.2015
Размер:
29.38 Mб
Скачать

Если наиболее высоким окажется коэффициент автокорреляции порядка 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=t11

 

 

 

 

 

Остатки определяются как среднеквадратическое отклонение уровней ряда от сглаженных.

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,3yt1 + 0,7Ft1, 3 t 36 .

(9.12)

Остатки вычисляются по формуле

 

 

 

 

 

1

t

Fi )2 , 5

 

 

 

et

=

(yi

t 36 .

(9.13)

 

 

 

4

i=t3

 

 

 

 

 

 

 

 

 

 

электроэнергии

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α,n2 = 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α, n2 = 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]