Файл: В.М. Волков Эконометрика. Методические указания и пример выполнения контрольной работы для студентов заочной формы обучения экономических специальностей 060500,060400.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 03.06.2024
Просмотров: 56
Скачиваний: 0
10 |
|
|
a∑t2 +b∑t = ∑t |
u |
(23) |
a∑t +bN = ∑ut . |
t |
|
|
|
9. Проверка адекватности трендовой модели
Для получения надежного, долговременного прогноза необходимо проверить трендовую модель на адекватность, то есть выяснить, не являются ли ошибки выбранной аппроксимации также трендовой моделью. А это означает, что случайная составляющая в выбранной модели не была исключена. Поэтому рассматривают ряд остатков –
разность значений ряда и значений тренда: |
|
εt = ut – yt. |
(24) |
Таким образом, проверяют следующие гипотезы:
а) о случайности ряда остатков методом сравнения средних первой и второй половин ряда остатков в соответствии с формулой (8). Если гипотеза о случайности ряда остатков отвергается, то трендовую модель следует считать неадекватной; б) о равенстве математического ожидания ряда остатков нулю по статистике
t =εt n / sε , |
(25) |
где εt – среднее значение ряда остатков; sε – среднее квадратическое
ряда остатков.
На 5 % уровне значимости вычисленное значение t сравнивается с критическим значением, взятым из табл. 2, с n–1 степенями свободы. Если гипотеза отвергается, то модель считается неадекватной на 5 % уровне значимости;
в) отсутствие автокорреляции ряда остатков; для проверки этой гипотезы используется критерий Дарбина-Уотсона со статистикой
n |
n |
2 . |
|
D= ∑(εt −εt −1 )2 |
∑εt |
(26) |
|
t =2 |
t =1 |
|
|
Если D [2;4], следует использовать вспомогательную статистику |
D1=4–D.
Расчетное значение D или D1 сравнивается с верхним d2 и нижним d1 критическими значениями статистики Дарбина-Уотсона, представленными в табл. 4 для различной длины ряда N и числа определяемых параметров модели (k+1) на уровне значимости 5 %.
11
|
|
|
|
|
|
Таблица 4 |
|
N |
k=1 |
|
k=2 |
|
k=3 |
|
|
|
d1 |
d2 |
d1 |
d2 |
d1 |
|
d2 |
10 |
0,88 |
1,32 |
0,70 |
1,64 |
0,53 |
|
2,02 |
15 |
1,08 |
1,36 |
0,95 |
1,54 |
0,82 |
|
1,75 |
20 |
1,20 |
1,41 |
1,10 |
1,54 |
1,00 |
|
1,68 |
30 |
1,35 |
1,49 |
1,28 |
1,57 |
1,21 |
|
1,65 |
Если расчетное значение критерия D больше верхнего табличного значения d2, то гипотеза о независимости уровней остаточной последовательности, т.е. об отсутствии в ней автокорреляции, принимается. Если значение D меньше нижнего табличного значения d1, то эта гипотеза отвергается и модель неадекватна. Если значение D находится между значениями d1 и d2, включая сами эти значения, то считается, что нет достаточных оснований сделать тот или иной вывод и необходимы дальнейшие исследования, например по большему числу наблюдений.
Трендовая модель считается адекватной, если подтверждены все три гипотезы а), б), в).
10. Краткосрочный и долгосрочный прогнозы значений временного ряда и оценка относительной погрешности прогнозов
Кратковременный прогноз значений временного ряда на один шаг для t = N+1 находится с помощью метода скользящих средних по m последним значениям ряда на основе выбора наиболее достоверной модели из п. 5, 6.
Долговременное прогнозирование значений временного ряда на k
шагов вперед осуществляется по уравнению тренда: |
|
~ |
(27) |
uN +k = a(N + k)+b . |
Относительная погрешность прогнозов определяется по следующей формуле
|
|
~ |
|
|
|
|
|
|
|
|
|
||
Пк = |
|
un+k −un+k |
|
100% , (k =1, 2, 3, 4, 5) . |
(28) |
|
un+k |
||||||
|
|
|
|
12
Пример выполнения контрольного задания
Временной ряд представляет среднюю заработную плату работников угольной промышленности Кузбасса за 10 месяцев
1992 года в тыс. рублей: 32, 33, 36, 41, 68, 57, 96, 113, 132, 113. Кон-
трольные значения средней заработной платы за следующие
3месяца 128, 137, 134 тыс. рублей.
1.Строим график временного ряда.
140 |
|
120 |
|
100 |
|
80 |
|
60 |
|
40 |
|
20 |
|
0 |
t |
0 1 2 3 4 5 6 7 8 9 10
Вычисляем среднее значение
1 10
ut = 10 t∑=1ut = 72 ,
дисперсию
s2 = |
1 10∑(u −72)2 |
=1461,9, |
|
9 t =1 t |
|
среднее квадратическое отклонение s= s2 = 1461,9 = 38,2.
2. Проводим линейное сглаживание временного ряда по m = 5 точкам по формулам (5) и (6).
Заданные и сглаженные значения временного ряда заносим в таблицу. Строим графики этих рядов.
t |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
ut |
32 |
33 |
36 |
41 |
68 |
57 |
96 |
113 |
132 |
113 |
~ |
26 |
34 |
42 |
47 |
60 |
75 |
93 |
102 |
117 |
132 |
ut |
|
|
|
|
|
|
|
|
|
13 |
140 |
|
|
|
|
|
|
|
|
|
120 |
|
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
|
|
|
60 |
|
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
t |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 10 |
График сглаженного ряда показывает монотонное возрастание значений ряда во времени.
3. Проверяем гипотезу о случайности ряда на основе сравнения средних первой и второй половин ряда. Предварительно вычисляем величины:
u |
= |
1 |
∑5 |
u |
= 42 , u |
= |
1 |
10∑u |
= 102,2 ; |
|
|
|||||||||
5 |
|
|
|
|||||||||||||||||
|
п |
|
|
1 |
t |
|
в |
|
5 6 |
t |
|
|
|
|
|
|
|
|||
|
2 |
|
1 |
5 |
2 |
|
|
|
|
|
|
2 |
|
1 |
10 |
|
2 |
|||
s |
п |
= |
|
|
|
∑(u − 42) |
= 223,5 , |
s |
в |
= |
|
∑ |
(u |
−102,2) = 800,7 . |
||||||
|
|
|
|
|||||||||||||||||
|
|
4 |
1 |
t |
|
|
|
|
|
|
|
4 |
6 |
t |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Проверяем гипотезу о равенстве дисперсий первой и второй половины ряда по критерию Фишера (9):
F = sв2 = 800,7 = 3,58. sп2 223,5
Из табл. 1 находим Fкр= 5,0. Так как 3,58<5,0, гипотеза о равенстве дисперсий принимается.
Вычисляем величину t по формуле (10):
t = |
42−102,2 |
|
5 5 (5 +5 −2) |
=− 4,21 . |
4 |
223,5 +4 800,7 |
|
5 +5 |
|
По табл. 2 определяем tкр =t8;0,05 =2,38. Так как t > tкр, то различия между средними первой и второй половинами ряда значимы, ряд не является случайным и временной тренд существует на уровне значимости 5 %.
4. Проводим автокорреляционный анализ временного ряда. Вычисляем коэффициенты автокорреляции по формуле (9), пред-
варительно составив следующие таблицы. В последнем столбце каж-
14
дой таблицы вычислено среднее значение. k = 1
ut |
|
|
|
|
32 |
|
|
|
33 |
|
|
36 |
|
|
|
41 |
|
|
|
68 |
|
|
57 |
|
96 |
|
113 |
|
132 |
|
68 |
|||||||||||
ut+1 |
|
|
|
|
33 |
|
|
|
36 |
|
|
41 |
|
|
|
68 |
|
|
|
57 |
|
|
96 |
|
113 |
|
132 |
|
113 |
|
77 |
|||||||||||
ut ut+1 |
|
|
|
1056 |
|
1188 |
|
1476 |
2788 |
3876 |
5472 |
10848 |
14916 |
14916 |
6282 |
|||||||||||||||||||||||||||
Для |
величины ut дисперсия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
s |
2 |
= 1 ∑9 (u |
−68)2 |
=1255,6 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
1 |
|
|
|
9 t =1 |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Для величины ut+1 |
дисперсия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
s |
2 |
= 1 |
∑9 (u |
|
−77)2 |
=1263,6 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
2 |
|
|
|
9 t =1 |
|
t |
+1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Коэффициент корреляции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
r 1 |
= |
|
6282 −68 77 |
= 0,83 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
1255,6 |
1263,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
k = 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ut |
|
|
|
32 |
|
|
33 |
|
|
36 |
|
|
|
41 |
|
|
68 |
|
57 |
|
|
|
96 |
|
|
113 |
|
|
|
60 |
|
|||||||||||
Ut+2 |
|
|
|
36 |
|
|
41 |
|
|
68 |
|
|
|
57 |
|
|
96 |
|
113 |
|
|
132 |
|
113 |
|
|
|
82 |
|
|||||||||||||
ut ut+2 |
|
|
1152 |
|
1353 |
2448 |
2337 |
6528 |
6441 |
12672 |
12769 |
|
|
5713 |
||||||||||||||||||||||||||||
|
|
= 1 |
8 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
8 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||||
s |
|
|
|
|
−60) |
= |
828,3 ; s |
|
|
2 = |
|
|
|
|
−82) |
=1154,5; |
|
|||||||||||||||||||||||||
2 |
∑(u |
|
|
∑(u |
|
|
||||||||||||||||||||||||||||||||||||
1 |
|
|
|
8 t =1 |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
8 t =1 |
|
t +2 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
= |
5713 −60 82 |
= 0,81. |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
28,8 34,0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
k = 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
ut |
|
|
32 |
|
33 |
|
36 |
|
|
41 |
|
|
68 |
|
|
|
57 |
|
96 |
|
|
|
|
|
|
|
|
52 |
|
|||||||||||||
ut+3 |
|
|
41 |
|
68 |
|
57 |
|
|
96 |
|
|
113 |
|
132 |
|
|
113 |
|
|
|
|
|
|
|
|
89 |
|
||||||||||||||
ut ut+3 |
|
1312 |
2244 |
2052 |
3936 |
7684 |
7524 |
10848 |
|
|
|
|
|
|
|
5086 |
||||||||||||||||||||||||||
|
|
= 1 |
7 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
7 |
|
|
|
|
|
2 |
|
|
|
|
|
|
||||||||
s |
|
|
|
|
−52) |
= |
479,2; s |
|
2 = |
|
(u |
|
−89) |
= 973,0; |
|
|
|
|||||||||||||||||||||||||
2 |
∑(u |
|
|
∑ |
|
|
|
|
||||||||||||||||||||||||||||||||||
1 |
|
|
|
7 t =1 |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
7 t =1 |
t +3 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r3 = |
5086 −52 89 |
= 0,67. |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
21,9 31,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|