Построение временного ряда уровня безработицы за 10 месяцев 2010 года в России
Построение временной ряда величины по данным об уровне безработицы в России за 10 месяцев 2010 г., вычисление ее числовых характеристик. Регрессионная модель временного тренда. Краткосрочный и долгосрочный прогнозы изменения рассматриваемой величины.
| Рубрика | Экономико-математическое моделирование |
| Вид | контрольная работа |
| Язык | русский |
| Дата добавления | 26.02.2012 |
| Размер файла | 118,1 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Содержание
- Начальные данные
- 1. Построить временной ряд величины и вычислить ее основные числовые характеристики
- 2. Сгладить ряд методом скользящих средний
- 3. Проверить гипотезы о случайности временного ряда
- 4. Автокорреляционный анализ временного ряда
- 5. Оценить модели краткосрочного прогноза
- 6. Определить степень полиномиального тренда методом переменных разностей
- 7. Построить регрессионную модель временного тренда
- 8. Оценить адекватность трендовой модели
- 9. Дать краткосрочный и долгосрочный прогнозы изменения рассматриваемой величины
- Список литературы
Начальные данные
Временной ряд представляет уровень безработицы за 10 месяцев 2010 года в России. Данные взяты с сайта "Статистическая база данных по российской экономике" (http://stat. hse.ru/hse/indexn.html)
|
месяц |
Уровень безработицы, % |
|
|
1 |
9,2 |
|
|
2 |
8,6 |
|
|
3 |
8,6 |
|
|
4 |
8,2 |
|
|
5 |
7,3 |
|
|
6 |
6,8 |
|
|
7 |
7,0 |
|
|
8 |
6,9 |
|
|
9 |
6,6 |
|
|
10 |
6,8 |
1. Построить временной ряд величины и вычислить ее основные числовые характеристики
Строим график временного ряда.
Вычисляем среднее значение
Дисперсию
среднее квадратическое отклонение
2. Сгладить ряд методом скользящих средний
Проводим линейное сглаживание временного ряда по точкам по формулам
Заданные и сглаженные значения временного ряда заносим в таблицу.
|
t |
ut |
||
|
1 |
9,2 |
||
|
2 |
8,6 |
||
|
3 |
8,6 |
||
|
4 |
8,2 |
||
|
5 |
7,3 |
||
|
6 |
6,8 |
||
|
7 |
7,0 |
||
|
8 |
6,9 |
||
|
9 |
6,6 |
||
|
10 |
6,8 |
Строим графики этих рядов.
График сглаженного ряда показывает монотонное убывание значений ряда во времени.
3. Проверить гипотезы о случайности временного ряда
Проверяем гипотезу о случайности ряда на основе сравнения средних значений первой и второй половины ряда. Предварительно вычисляем величины:
Проверяем гипотезу о равенстве дисперсий первой и второй половины ряда по критерию Фишера (9):
Из табл.1 находим Fкр= 5,0. Так как 21,89 > 5,0, гипотеза о равенстве дисперсий не принимается.
Так как дисперсии не равны, проверку гипотезы о случайности ряда выполняем методом поворотных точек. Число поворотных точек (т.е. точек, в которых значение больше или меньше значения в соседних точках) равно 3.
Для случайного ряда среднее число поворотных точек и их дисперсия равны
Вычисляем статистику
Следовательно, гипотеза о случайности ряда принимается, и временного тренда нет.
4. Автокорреляционный анализ временного ряда
Проводим автокорреляционный анализ временного ряда.
Вычисляем коэффициенты автокорреляции по формуле
,
Где , , - средние значения,
а s1 и s2 - средние квадратические отклонения рядов ut и ut+k..
Предварительно составим вспомогательные таблицы.
В последнем столбце каждой таблицы вычислено среднее значение.
k = 1
|
ut |
ut+1 |
ut·ut+1 |
|
|
9,2 |
8,6 |
79,12 |
|
|
8,6 |
8,6 |
73,96 |
|
|
8,6 |
8,2 |
70,52 |
|
|
8,2 |
7,3 |
59,86 |
|
|
7,3 |
6,8 |
49,64 |
|
|
6,8 |
7,0 |
47,6 |
|
|
7,0 |
6,9 |
48,3 |
|
|
6,9 |
6,6 |
45,54 |
|
|
6,6 |
6,8 |
44,88 |
|
|
7,689 |
7,422 |
57,713 |
Для величины ut дисперсия
Для величины ut+1 дисперсия
.
Коэффициент корреляции
k = 2
|
ut |
ut+2 |
ut·ut+2 |
|
|
9,2 |
8,6 |
79,12 |
|
|
8,6 |
8,2 |
70,52 |
|
|
8,6 |
7,3 |
62,78 |
|
|
8,2 |
6,8 |
55,76 |
|
|
7,3 |
7,0 |
51,1 |
|
|
6,8 |
6,9 |
46,92 |
|
|
7,0 |
6,6 |
46,2 |
|
|
6,9 |
6,8 |
46,92 |
|
|
7,825 |
7,275 |
57,415 |
Для величины ut дисперсия
Для величины ut+2 дисперсия
.
Коэффициент корреляции
временной ряд величина безработица
k = 3
|
ut |
ut+3 |
ut·ut+3 |
|
|
9,2 |
8,2 |
75,44 |
|
|
8,6 |
7,3 |
62,78 |
|
|
8,6 |
6,8 |
58,48 |
|
|
8,2 |
7,0 |
57,4 |
|
|
7,3 |
6,9 |
50,37 |
|
|
6,8 |
6,6 |
44,88 |
|
|
7,0 |
6,8 |
47,6 |
|
|
7,957 |
7,086 |
56,707 |
Для величины ut дисперсия
Для величины ut+3 дисперсия
.
Коэффициент корреляции
k = 4
|
ut |
ut+4 |
ut·ut+4 |
|
|
9,2 |
7,3 |
67,16 |
|
|
8,6 |
6,8 |
58,48 |
|
|
8,6 |
7,0 |
60,2 |
|
|
8,2 |
6,9 |
56,58 |
|
|
7,3 |
6,6 |
48,18 |
|
|
6,8 |
6,8 |
46,24 |
|
|
8,117 |
6,900 |
56,140 |
Для величины ut дисперсия
Для величины ut+4 дисперсия
.
Коэффициент корреляции
k = 5
|
ut |
ut+5 |
ut·ut+5 |
|
|
9,2 |
6,8 |
62,56 |
|
|
8,6 |
7,0 |
60,2 |
|
|
8,6 |
6,9 |
59,34 |
|
|
8,2 |
6,6 |
54,12 |
|
|
7,3 |
6,8 |
49,64 |
|
|
8,380 |
6,820 |
57,172 |
Для величины ut дисперсия
Для величины ut+5 дисперсия
.
Коэффициент корреляции
Первые пять значений коэффициентов автокорреляции имеют следующие значения: r1 = 1,33; r2 = 1,37; r3 = 1,79; r4 = 4,18; r5 = 2,88.
По табл.3 находим критические значения для этих коэффициентов при 5-% уровне значимости: r1кр = 0,67; r2кр = 0,71; r3кр = 0,75; r4кр = 0,81; r5кр = 0,87. Коэффициент автокорреляции значим, если его значение больше соответствующего критического.
Построим коррелограмму:
Здесь сплошной линией обозначена автокорреляционная функция, а пунктирной - критический уровень коэффициентов автокорреляции.
Так как для первых 5 значений выполняется условие rк rкр, то имеется значимая зависимость между первыми 5 и последними 5 членами ряда. Таким образом, временная длина зависимости составляет 5 временных единиц, что говорит о возможности долгосрочного прогноза на 5 временных шагов вперед.
5. Оценить модели краткосрочного прогноза
Даем оценку моделей краткосрочного прогноза по следующим формулам:
1) прогноз по одному последнему значению un+1 (1) = un
2) прогноз по двум последним значениям un+1 (2) = 2un-un-1
3) прогноз по трем последним значениям un+1 (3) = (4un+ un-1 - 2un-2) /3
4) прогноз по четырем последним значениям un+1 (4) = (2un+ un-1 - un-3) /2
5) прогноз по пяти последним значениям
un+1 (5) = (8un+ 5 un-1 + 2un-2-un-3 - 4un-4) /10
|
t |
ut |
ut (1) |
ut (2) |
ut (3) |
ut (4) |
ut (5) |
кр |
|
|
1 |
9,2 |
0,92 |
||||||
|
2 |
8,6 |
9,2 (+) |
0,86 |
|||||
|
3 |
8,6 |
8,6 (+) |
8 (+) |
0,86 |
||||
|
4 |
8,2 |
8,6 (+) |
8,6 (+) |
8, 20 (+) |
0,82 |
|||
|
5 |
7,3 |
8,2 |
7,8 (+) |
8,07 |
7,9 (+) |
0,73 |
||
|
6 |
6,8 |
7,3 (+) |
6,4 (+) |
6,73 (+) |
7,1 (+) |
7,12 (+) |
0,68 |
|
|
7 |
7,0 |
6,8 (+) |
6,3 (+) |
6,03 |
6,15 |
6,43 (+) |
0,7 |
|
|
8 |
6,9 |
7 (+) |
7,2 (+) |
6,73 (+) |
6,3 (+) |
6,2 |
0,69 |
|
|
9 |
6,6 |
6,9 (+) |
6,8 (+) |
7,00 (+) |
6,75 (+) |
6,37 (+) |
0,66 |
|
|
10 |
6,8 |
6,6 (+) |
6,3 (+) |
6,43 (+) |
6,65 (+) |
6,53 (+) |
0,68 |
|
|
89% |
100% |
71% |
83% |
80% |
В качестве критических величин погрешностей прогнозных значенийкр выбираем 10% от заданных значений ряда.
Числа со знаками "+" показывают достоверные прогнозы. Процент достоверных прогнозов приведен в последней колонке.
Из таблицы следует, что наиболее достоверной, в данном случае, является модель краткосрочного прогноза по двум последним точкам.
6. Определить степень полиномиального тренда методом переменных разностей
Определяем степень полиномиального тренда методом переменных разностей по формулам
Д1ut = ut+1, - ut, где t = 1, …, N - 1.
Д2ut = Д1ut+1, - Д1ut, где t = 1, …, N - 2.
……………………………….
Дnut = Дn-1ut+1, - Дn-1ut, где t = 1, …, N - n.
Под разностями нулевого порядка понимается сам временной ряд.
На каждом шаге, начиная с n = 0, вычисляют дисперсии разностей n - го порядка по формуле
На каждом шаге для каждых двух (предыдущей и последующей) дисперсий проверяют гипотезу о равенстве дисперсий по критерию Фишера
Если Fкр = F (, N - n, N - n - 1), то можно принять, что дисперсии отличаются незначимо. В противном случае процедура вычисления разностей и их дисперсий продолжается.
|
t |
0ut |
1ut |
2ut |
|
|
1 |
9,2 |
-0,6 |
0,6 |
|
|
2 |
8,6 |
0 |
-0,4 |
|
|
3 |
8,6 |
-0,4 |
-0,5 |
|
|
4 |
8,2 |
-0,9 |
0,4 |
|
|
5 |
7,3 |
-0,5 |
0,7 |
|
|
6 |
6,8 |
0,2 |
-0,3 |
|
|
7 |
7,0 |
-0,1 |
-0,2 |
|
|
8 |
6,9 |
-0,3 |
0,5 |
|
|
9 |
6,6 |
0,2 |
||
|
10 |
6,8 |
|||
|
7,6 |
-0,267 |
0,100 |
Дисперсия разностей нулевого и первого порядков:
Вычислим значение параметра распределения Фишера и сравним его с критическим значением из табл.1.
Так как F1 F1кр, то значимо отличается от .
Проводим далее аналогично сравнение дисперсий и .
Вычислим значение параметра распределения Фишера и сравним его с критическим значением из табл.1.
Так как F2 < F2кр, то незначимо отличается от , и степень полиномиального тренда р = 1.
7. Построить регрессионную модель временного тренда
Строим полиномиальный тренд временного ряда степени p = 1:
оценки коэффициентов а и b находятся из системы линейных уравнений:
В нашем случае:
откуда находим а = - 0,292 и b = 9, 207.
Тогда уравнение временного тренда:
Строим графики временного ряда и тренда.
8. Оценить адекватность трендовой модели
Проверяем адекватность трендовой модели. Сначала вычисляем трендовые значения и значения остатков ряда по формуле t = ut - yt. Результаты заносим в таблицу.
|
t |
ut |
yt |
t = ut - yt |
|
|
1 |
9,2 |
8,91 |
0,29 |
|
|
2 |
8,6 |
8,62 |
-0,02 |
|
|
3 |
8,6 |
8,33 |
0,27 |
|
|
4 |
8,2 |
8,04 |
0,16 |
|
|
5 |
7,3 |
7,75 |
-0,45 |
|
|
6 |
6,8 |
7,45 |
-0,65 |
|
|
7 |
7,0 |
7,16 |
-0,16 |
|
|
8 |
6,9 |
6,87 |
0,03 |
|
|
9 |
6,6 |
6,58 |
0,02 |
|
|
10 |
6,8 |
6,29 |
0,51 |
а) Проверяем гипотезу о случайности ряда остатков методом поворотных точек. Здесь
Число поворотных точек заданного ряда находим из графика ряда, d = 5 Вычисляем статистику . Гипотеза о случайности ряда остатков принимается.
б) Далее проверяем гипотезу о равенстве нулю математического ожидания остатков ряда по статистике
где - среднее значение ряда остатков, - среднее квадратическое ряда остатков
Так как вычисленное значение статистики t меньше критического tкр = 2,23 из табл.2, то гипотеза принимается.
Трендовая модель адекватна, потому что обе гипотезы приняты.
9. Дать краткосрочный и долгосрочный прогнозы изменения рассматриваемой величины
Кратковременный прогноз временного ряда на один шаг (один месяц) вперед выполняем по формуле прогноза по двум последним значениям
un+1 (2) = 2un-un-1, u11 (2) = 2·u10-u9 = 2·6,8 - 6,6 = 7,0
Долговременный прогноз на три шага вперед временного ряда производим по найденному уравнению линейного тренда
Таким образом, средняя прогнозируемый уровень безработицы за 11 месяц 2010 года составит 7%, а за 1 месяц 2011 года 5,41%.
Список литературы
1. Методические указания к работе
2. Степанов В.Г. Эконометрика. М., 2010.
Размещено на Allbest.ru
Подобные документы
Теория и анализ временных рядов. Построение линии тренда и прогнозирование развития случайного процесса на основе временного ряда. Сглаживание временного ряда, задача выделения тренда, определение вида тенденции. Выделение тригонометрической составляющей.
курсовая работа [722,6 K], добавлен 09.07.2019Изучение понятия имитационного моделирования. Имитационная модель временного ряда. Анализ показателей динамики развития экономических процессов. Аномальные уровни ряда. Автокорреляция и временной лаг. Оценка адекватности и точности трендовых моделей.
курсовая работа [148,3 K], добавлен 26.12.2014Показатели статистики занятости и безработицы, а также баланс трудовых ресурсов. Изучение межрегиональной вариации уровня безработицы. Построение уравнения регрессии. Регрессионная модель зависимости уровня безработицы и внутреннего валового продукта.
курсовая работа [604,2 K], добавлен 16.09.2014Построение графика временного ряда. Тренд - устойчивое систематическое изменение процесса в течение продолжительного времени. Динамика продаж бензина на АЗС. Выявление сезонной составляющей и тренда. Коррелограмма, построенная в программе Statistica.
курсовая работа [1,2 M], добавлен 15.11.2013Аддитивная модель временного ряда. Мультипликативная модель временного ряда. Одномерный анализ Фурье. Регрессионная модель с переменной структурой. Сущность адаптивной сезонной модели Тейла – Вейджа. Прогнозирование естественного прироста населения.
курсовая работа [333,1 K], добавлен 19.07.2010Выборка и генеральная совокупность. Модель множественной регрессии. Нестационарные временные ряды. Параметры линейного уравнения парной регрессии. Нахождение медианы, ранжирование временного ряда. Гипотеза о неизменности среднего значения временного ряда.
задача [62,0 K], добавлен 08.08.2010Анализ автокорреляции уровней временного ряда, характеристика его структуры; построение аддитивной и мультипликативной модели, отражающую зависимость уровней ряда от времени; прогноз объема выпуска товаров на два квартала с учетом выявленной сезонности.
лабораторная работа [215,7 K], добавлен 23.01.2011Двойственные оценки как мера влияния ограничений на функционал. Построение экономико-математической модели задачи. Выявление аномальных уровней временного ряда с использованием метода Ирвина. Построение графика общих годовых затрат по выгодному способу.
контрольная работа [282,7 K], добавлен 16.01.2012Построение вариационного (статистического) ряда, гистограммы и эмпирической функции распределения. Определение выборочных оценок числовых характеристик случайной величины. Расчет матрицы парных коэффициентов корреляции и создание модели парной регрессии.
контрольная работа [2,0 M], добавлен 05.04.2014Анализ временных рядов с помощью статистического пакета "Minitab". Механизм изменения уровней ряда. Trend Analysis – анализ линии тренда с аппроксимирующими кривыми (линейная, квадратическая, экспоненциальная, логистическая). Декомпозиция временного ряда.
методичка [1,2 M], добавлен 21.01.2011
