Файл: В.М. Волков Эконометрика. Методические указания и пример выполнения контрольной работы для студентов заочной формы обучения экономических специальностей 060500,060400.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 03.06.2024

Просмотров: 221

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

6

Стьюдента приведены в табл. 2.

Таблица 2

k

3

5

7

10

13

16

20

30

tкр

3,18

2,57

2,45

2,23

2,16

2,12

2,09

2,04

1,96

В случае, если t< tкр, гипотеза о случайности временного ряда (равенства средних) принимается. В противном случае – отвергается, что говорит о значимости различия средних первой и второй половины ряда, неслучайном поведении ряда и наличии временного тренда.

4. Автокорреляционный анализ временных рядов

Задачей автокорреляционного анализа временного ряда является установление степени и временного интервала влияния последующих членов ряда от предыдущих. Наличие корреляционной связи между последующими и предыдущими членами ряда также служит информативным признаком временного тренда.

Для этого последовательно рассчитывают коэффициенты автокорреляции rk между первыми и последними (N – k) членами ряда ut

и ut+k (k = 1, 2,...):

r = (

),

ut ut +k

ut

ut +k

(9)

k

s1

s2

где ut , ut +k , ut ut +k – средние значения, а s1 и s2 – средние квадратические отклонения рядов ut и ut+k.

Средние значения величин ut , ut +k , ut ut +k вычисляются как

средние арифметические этих значений по формуле (2), например:

ut ut +k = N∑−k(ut ut +k ) (N −k) . t =1

Среднеквадратическое отклонение s = s 2 , где s 2 – дисперсия

величины иt равна:

s2 2 = N∑−k (ut +k −ut +k )2 (N −k). t =1

После вычисления коэффициентов автокорреляции по формуле (9) проверяется значимость коэффициентов автокорреляции сравнением этих значений с критическими значениями коэффициента корреляции rкр. Если rк < rкр, то корреляция на временном интервале в k


7

единиц отсутствует. При 5-процентном уровне значимости критические значения коэффициента корреляции приведены в табл. 3.

Таблица 3

N–k–2

3

5

7

9

11

13

15

20

30

40

rкр

0,88

0,75

0,67

0,60

0,55

0,51

0,48

0,42

0,35

0,30

Полученные результаты оформляют в виде графической зависимости rк от k, которая носит название коррелограммы. Если для первых k значений выполняется условие rк> rкр, то имеется значимая зависимость между первыми (N–k) и последними (N–k) членами ряда. Таким образом, временная длина зависимости составляет k временных единиц, что говорит о возможности долгосрочного прогноза на k временных шагов вперед.

5. Модели краткосрочного прогноза

Для краткосрочного прогноза на один временной шаг используют 1, 2, 3, 4 или 5 последних значений временного ряда. Прогнозное значение определяют по следующим формулам:

1) прогноз по одному последнему значению

un+1(1) = un ;

(10)

2)

прогноз по двум последним значениям

3)

un+1(2) = 2un–un-1;

(11)

прогноз по трем последним значениям

4)

un+1(3) = (4un+ un-1 –2un-2)/3;

(12)

прогноз по четырем последним значениям

5)

un+1(4) = (2un+ un-1 –un-3)/2;

(13)

прогноз по пяти последним значениям

un+1(5) = (8un+ 5 un-1 + 2un-2–un-3 – 4un-4)/10.

(14)

6. Оценка точности и достоверности краткосрочного прогноза

Оценку точности краткосрочного прогноза проводят на основе сравнения прогнозируемых значений ряда un+1( к ) с известными зна-

чениями un+1 .

Для этого вычисляют абсолютную погрешность прогнозного значения по выражению

=

u

( к ) −и

.

(15)

к

n+1

n+1


8

В первой модели погрешности вычисляют для значений ряда, начиная со второго по формуле (10). Во второй модели для вычисления погрешностей используют прогнозные значения, рассчитанные по формуле (11), и данные значения, начиная с третьего, и т.д.

Затем задают предельное значение погрешности ∆кр, с которым сравнивают рассчитанные абсолютные погрешности. Предельное значение погрешности может выбираться как некоторое абсолютное значение (например, 5 копеек; 1,5 тонны) или как относительное значение – 5 %, 10 % от значений временного ряда, что определяется существом решаемой задачи. Если ∆к< ∆кр, то прогноз считается точным, в противном случае – неточным. Для каждой модели производят подсчет числа точных прогнозов К+. Далее оценивают достоверность каждой модели прогноза, для чего рассчитывают процентное отношение точных прогнозов К+ к общему числу прогнозов К= N–k, где k

– номер модели. То есть достоверность определяют как

δ = (К+/ К) 100%.

(16)

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

7. Определение степени полиномиального тренда методом переменных разностей

Для выделения полиномиального тренда степени p предварительно находят значение этой степени по следующей процедуре, которая аналогична дифференцированию полинома. Так, очевидно, что вторая производная полинома первой степени, третья производная полинома второй степени и т.д. равны нулю.

Для временного ряда операция дифференцирования заменяется вычислением переменных разностей, а условие равенства производной нулю – проверкой гипотезы о равенстве дисперсий предыдущих и последующих разностей.

Таким образом, сначала вычисляют первые разности:

∆1u = u

−u ,

(17)

t

t +1

t

где t = 1, ..., N – 1.

Затем по первым разностям вычисляют вторые разности:

∆2ut = ∆ 1ut +1 −∆ 1ut ,

(18)

где t = 1, ..., N – 2.

И далее последовательно – разности 3-го, ..., n-го порядков:


9

∆ nut = ∆ n-1ut +1 −∆ n-1ut ,

(19)

где t = 1, ..., N – n.

Под разностями нулевого порядка понимается сам временной ряд. На каждом шаге, начиная с n = 0, вычисляют дисперсии разностей

n-го порядка по формуле

N∑−n (∆ n ut −

)2

2

∆ nut

(n!)2

s

=

t =1

(20)

n

(N −n −1) (2n)!

(при n= 0 имеем дисперсию всего временного ряда).

На каждом шаге для каждых двух (предыдущей и последующей) дисперсий проверяют гипотезу о равенстве дисперсий по критерию Фишера в соответствии с (7).

Если Fn < Fкр, то можно принять, что дисперсии отличаются незначимо. В противном случае процедура вычисления разностей и их дисперсий продолжается. Здесь

F

= (s

n−1

)2 s

2

,

(21)

n

n

а Fкр = F(α, k1 , k2), где α – принятый уровень значимости; k1 = N– n, k2 = N – n – 1 – степени свободы.

Для 5 % уровня значимости критические значения распределения Фишера Fкр приведены в табл. 1.

Последовательность дисперсий (20) убывает с ростом n, и при некотором значении p = n – 1 выполняется неравенство Fn < Fкр. Полученное значение p и является степенью полиномиального тренда.

Дисперсия sp2 называется дисперсией случайности, а разности порядка p являются случайной компонентой временного ряда.

8. Выделение полиномиального тренда

После определения степени полиномиального тренда методом наименьших квадратов находят уравнение тренда (оценка его коэффициентов).

Для p = 1 – линейного тренда

yt = at +b ,

(22)

оценки коэффициентов тренда находятся из системы линейных уравнений: