Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
уч_комплекс_эконометр.doc
Скачиваний:
33
Добавлен:
04.04.2013
Размер:
2.81 Mб
Скачать

Анализ сезонных колебаний

При анализе рядов динамики важное значение имеет выявление сезонных колебаний. Этим колебаниям свойственны более или менее устойчивые изменения уровней ряда по временным периодам: годам, кварталам, месяцам, неделям. Гармонический тренд – один из способов представления сезонных явлений.

Формального математического определения сезонности нет, но зато есть огромное количество всевозможных практических примеров проявления сезонности (например, спрос на меховые шапки подвержен сезонным колебаниям, количество пересылаемых по почте писем также связано с сезонностью – в конце каждого квартала организации пересылают по почте свою отчетность в налоговые инспекции). В сезонности выделяют 2 параметра, которые необходимо оценить – период и форму.

Замечание: для выявления сезонности не обязательно производить наблюдения строго периодически.

Для выявления сезонных колебаний обычно анализируются месячные, годовые и квартальные уровни ряда динамики за год или за несколько лет. При изучении сезонных колебаний используются модели сезонных волн. Например, в следующем виде:

где k – определяет гармонику тригонометрического многочлена,

t – время.

T – объем выборки.

Целью гармонического анализа является выявление и измерение периодических колебаний во временных ряда. Это уравнение представляет собой ряд Фурье. Обычно высчитывают не более четырех гармоник.

Введем понятия частота ипериод k-й гармоники

Будем строить регрессоры из функций вида ив соответствие с вышеуказанной идеей – т.е. необходимо построитьT линейно независимых регрессоров.

На практике описанное выше представление очень громоздко, пользоваться им нельзя. Оно нужно для того, чтобы построить т.н. спектрограмму. Сначала ищутся оценки всех коэффициентов в сумме. Тогда значения спектрограммы Спектрограмма – это графикrk относительно частот.

Пики спектрограммы соответствуют частотам, вносящим наиболее яркий и значимый вклад в прогноз. Наибольшая частота на графике спектрограмме всегда = ½, это максимальная частота, различимая при дискретных по времени наблюдениях.

Спектрограмма позволяет связать теорию рядов Фурье с практикой прогнозирования. Она дает возможность оценить те частоты, которые вносят наибольший вклад, т.е. аппроксимировать громоздкую сумму разложения Фурье. Иногда ее удается сократить до 2-3 слагаемых.

Замечание: на графике спектрограммы следует выделить не только частоты, соответствующие пикам, но и близкие к ним частоты. В соответствии с общей методологией анализа лучше сначала взять больше частот, чем нужно, а затем исключить избыточные.

Модель ряда с гармоническим трендом:

,

- наиболее значимые частоты.

Как и полиномиальные тренды, гармонические сводятся к модели линейной регрессии

Время (t) выражается в радиальной мере или градусах:

Месяцы (t)

1

2

3

4

5

6

7

8

9

10

11

12

Радиальная мера

0

π/6

π/3

π/2

2π/3

5π/6

π

7π/6

4π/3

3π/2

5π/3

11π/6

Градусы

0

30

60

90

120

150

180

210

240

270

300

320

Уровни ()

Эта процедура перехода от номеров времени к градусам называется преобразованием Фурье.

Задача 3. По данным о продаже картофеля на рынках сельхозпродуктов города проанализировать сезонные колебания продаж.

Данные и промежуточные результаты представлены в таблице:

Месяц, t

Продано,

y

ycost

ysint

1

30

30,0

0

35,60

2

40

36,64

20

40,51

3

43

21,5

37,24

45,11

4

54

0

54,00

48,03

5

67

-33,5

58,02

48,55

6

29

-25,11

14,5

46,53

7

35

-35,00

0

42,52

8

34

-29,44

-17,00

38,25

9

45

-22,50

-38,97

33,05

10

35

0

-35,00

30,13

11

29

14,50

-25,11

29,61

12

28

24,25

-14,00

31,63

Итого

469

-20,66

53,68

469,52

На основе полученных итогов данных таблицы находим:

Отсюда Подставляя в это уравнение значенияcost и sint, получаем теоретические значения количества проданного картофеля .

Параметры гармоник второго и высшего порядка рассчитываются аналогично, и их значения последовательно присоединяются к значениям первой гармоники.

Например, опуская расчеты, запишем уравнение для нашего ряда с учетом второй гармоники:

Подставив в это уравнение значения cos2t и sin2t, получим значенияпо месяцам. Сравнив эти две модели с помощью, например,R2, можно сделать вывод, какая гармоника ряда Фурье близка к фактическим уровням ряда.

Для уменьшения времени вычисления преобразования Фурье разработан алгоритм, получивший название быстрого преобразования Фурье. Это преобразование реализовано в Microsoft Excel. Однако у этой процедуры есть один существенный недостаток: входной интервал должен содержать число значений, равных степени 2. Один из способов преодоления данного недостатка – добавление в ряд констант (например, нулей) до тех пор, пока длина ряда не станет равной степени 2.

Задача 4. Даны данные о динамике урожайности зерновых культур в одном из хозяйств области.

Год

1983

1984

1985

1986

1987

1988

1989

1990

Урожайность

17,6

18,1

17,4

16,8

16

15,4

14

16,6

Год

1991

1992

1993

1994

1995

1996

1997

1998

Урожайность

14,4

14,2

14,6

13,8

13,4

14,2

13,2

13,2

Для представленного временного ряда требуется провести гармонический анализ динамики отклонений от основной тенденции.

Решение задачи начнем с построения трендовой модели ряда. В Microsoft Excel данную операцию удобнее всего проводить с помощью инструмента "Подбор линии тренда" из мастера диаграмм. Анализ трендовых моделей показывает, что в качестве рабочей модели модно выбрать линейную модель: . Такой выбор обусловлен тем, что, во-первых, коэффициент детерминацииR2=0,82 имеет достаточно высокое значение ( лишь очень незначительно уступает коэффициенту детерминации R2=0,85 для полинома 2-го порядка); во-вторых, все коэффициенты модели значимы; в-третьих, при прочих равных условиях данная модель наиболее проста для вычислений и наиболее "прозрачна" для последующей экономической интерпретации.

График уравнения тренда показан на рисунке:

Для более детального анализа построенной модели можно использовать режим "Регрессия". Показатели, рассчитанные в данном режиме, представлены в таблице:

ВЫВОД ИТОГОВ

Регрессионная статистика

Множественный R

0,907548626

R-квадрат

0,823644509

Нормированный R-квадрат

0,811047688

Стандартная ошибка

0,726020725

Наблюдения

16

Дисперсионный анализ

df

SS

MS

F

Значимость F

Регрессия

1

34,4648897

34,46489

65,3851

1,21E-06

Остаток

14

7,37948529

0,527106

Итого

15

41,844375

Коэффициенты

Стандартная ошибка

t-статистика

P-Значение

Нижние 95%

Верхние 95%

Y-пересечение

17,8875

0,38072848

46,9823

8,3E-17

17,0709

18,7040

t

-0,318382353

0,03937400

-8,08611

1,21E-06

-0,4028

-0,2339

ВЫВОД ОСТАТКА

Наблюдение

Предсказан-

ное Урожайность

Остатки

1

17,56911765

0,03088235

2

17,25073529

0,84926470

3

16,93235294

0,46764705

4

16,61397059

0,18602941

5

16,29558824

-0,2955882

6

15,97720588

-0,5772058

7

15,65882353

-1,6588235

8

15,34044118

1,25955882

9

15,02205882

-0,6220588

10

14,70367647

-0,5036764

11

14,38529412

0,21470588

12

14,06691176

-0,2669117

13

13,74852941

-0,3485294

14

13,43014706

0,76985294

15

13,11176471

0,08823529

16

12,79338235

0,40661764

Столбец Остатки содержит значения отклонений от основной тенденции (разность между эмпирическими и теоретическими значениями). Гармонический анализ проведем с помощью режима "Анализ Фурье" для остатков. Тогда мы построим обобщенный тренд, описывающий в наших данных и основную тенденцию, и сезонности.

Результат применения режима работы указан в таблице.

Остатки

Y

Действительная часть, Yd

Мнимая часть, Ym

R2

0,0308823

0

0

a0

0,00

0,0000

0,8492647

3,2179246833565+0,985461907212836i

3,217925

0,985

a1

0,20

b1

0,06

0,0442

0,4676470

1,39644406072503-1,12572977822617i

1,396444

-1,126

a2

0,09

b2

-0,07

0,0126

0,1860294

-2,2329138282402-1,2706916705285i

-2,23291

-1,271

a3

-0,14

b3

-0,08

0,0258

-0,2955882

-0,3470588235293

+1,0470588235294i

-0,34706

1,047

a4

-0,02

b4

0,07

0,0048

-0,5772058

0,71036905643527-3,48957546372124i

0,710369

-3,490

a5

0,04

b5

-0,22

0,0495

-1,6588235

-1,2905617077838

+3,3801525747150i

-1,29056

3,380

a6

-0,08

b6

0,21

0,0511

1,2595588

0,916384794330886-1,02165718009751i

0,916385

-1,022

a7

0,06

b7

-0,06

0,0074

-0,6220588

-4,2470588235294

-4,24706

0,000

a8

-0,27

b8

0,00

0,0705

-0,5036764

0,916384794330889+1,0216571800975I

0,916385

1,022

a9

0,06

b9

0,06

0,0074

0,2147058

-1,2905617077838-3,38015257471502i

-1,29056

-3,380

a10

-0,08

b10

-0,21

0,0511

-0,2669117

0,710369056435279+3,4895754637212I

0,710369

3,490

a11

0,04

b11

0,22

0,0495

-0,3485294

-0,34705882352939

-1,04705882352941i

-0,34706

-1,047

a12

-0,02

b12

-0,07

0,0048

0,7698529

-2,2329138282402

+1,2706916705285i

-2,23291

1,271

a13

-0,14

b13

0,08

0,0258

0,0882352

1,39644406072503+1,12572977822617i

1,396444

1,126

a14

0,09

b14

0,07

0,0126

0,4066176

3,21792468335649-0,98546190721284I

3,217925

-0,985

a15

0,20

b15

-0,06

0,0442

Столбец Действительная часть содержаться данные, полученные с помощью математической функции МНИМ.ВЕЩ., примененной к столбцу Y, а столбец Мнимая часть содержит значения, полученные с помощью математической функции МНИМ.ЧАСТЬ., примененной к столбцу Y.

Действительные и мнимые части связаны с гармоническими коэффициентами следующим образом:

. R2 = . ГрафикR2 представляет собой спектрограмму. Посмотрев на график, можно выделить 3 основные гармоники (они соответствуют 3-м основным пикам).

Их и следует включить в модель.

Тогда для наших данных мы получим следующую модель: y = -0,3184 t + 17,888 + 0,4 cos t –0,12 sin t + 0,04 cos 5t + 0,44 sin 5t – 0,27 cos 8t + 0 sin 8t + 0,04 cos 11t – 0,44 sin 11t .

Проверить значимость гармоник можно, построив линейную регрессию остатков на эти 4 гармоники (8 переменных) и посмотрев значимость коэффициентов. Если хотя бы одна переменная из двух выбранных пар sin и cos значима, оставляем всю гармонику.

Таким образом, получаем гармоническую модель, описывающую наши данные. Сложив трендовую и гармоническую модели, получаем общую модель.

На рисунке представлены графики реальных и подогнанных значений (полученных по обобщенной модели) нашего ряда для остатков после удаления тренда. Причем на рисунке 2 графика подогнанных, оценочных значений – в первом случае мы просто взяли первые четыре гармоники, а для второго графика мы включили 5, 8 и 11 гармоники. Ниже в таблице приведены расчеты.

t

u1

u2

u3

u4

y прогнозное 1

Остат

ки

u5

u8

u11

y прогноз

ное 2

0,000

0,40224

0,08727

-0,13956

-0,02169

0,328

0,03088

0,04439

-0,26544

0,04439

-0,17664

0,393

0,32453

0,16120

0,09322

-0,1309

0,448

0,84926

0,38621

0,26544

0,38554

1,03719

0,785

0,19747

0,14078

0,21101

0,02148

0,571

0,46764

-0,33928

-0,26544

-0,34104

-0,94576

1,178

0,04037

0,03797

0,06850

0,13093

0,278

0,18602

-0,12716

0,26543

-0,12315

0,01513

1,570

-0,1228

-0,0870

-0,158

-0,0212

-0,39

-0,2955

0,0445

-0,2654

0,4357

0,2149

1,963

-0,2674

-0,1611

-0,19

-0,1309

-0,75

-0,5772

-0,0578

0,2654

-0,2121

-0,0046

2,355

-0,3712

-0,1409

0,0128

0,0210

-0,478

-1,6588

-0,2790

-0,2654

-0,2725

-0,8170

2,748

-0,4187

-0,0382

0,1998

0,131

-0,126

1,2595

0,4190

0,2654

0,4218

1,1063

3,140

-0,4024

0,0868

0,1403

-0,0208

-0,196

-0,6220

-0,0409

-0,2654

-0,0520

-0,3583

3,533

-0,3249

0,1610

-0,0923

-0,1310

-0,387

-0,5036

-0,3878

0,2654

-0,3818

-0,5042

3,925

-0,1980

0,1410

-0,2110

0,0206

-0,247

0,2147

0,3370

-0,2654

0,3458

0,4174

4,318

-0,0410

0,0384

-0,0694

0,1310

0,059

-0,2669

0,1304

0,2654

0,1157

0,5116

4,710

0,1222

-0,086

0,1578

-0,0204

0,173

-0,3485

-0,436

-0,2653

-0,4348

-1,1369

5,103

0,2668

-0,1610

0,1904

-0,131

0,165

0,7698

0,2029

0,2653

0,2188

0,6871

5,495

0,3709

-0,141

-0,0118

0,0202

0,238

0,0882

0,2817

-0,2653

0,2664

0,2828

5,888

0,4186

-0,0387

-0,1995

0,1311

0,311

0,4066

-0,4180

0,2653

-0,4238

-0,5765

Где

u1=

u2=…,

u11=.

{y прогнозное 1} = a0 + u1 +u2 +u3 +u4

{y прогнозное 2} = a0 + u5 +u8 +u11.

Мы подбираем гармоническую модель для остатков. Потом, выбрав из двух моделей одну, наиболее точную, мы складываем ее с трендовой моделью, и получаем обобщенную модель.

График ее выглядит следующим образом:

Из анализа графика можно увидеть, что получившаяся модель достаточно хорошо описывает данные, достаточно точная.