Моделирование и прогнозирование естественного прироста населения в РФ
Моделирование и прогнозирование естественного прироста населения в РФ
В данной курсовой работе методы и модели эконометрического анализа используется с целью моделирования и прогнозирования естественного прироста населения в Российской Федерации. Значения естественного прироста населения представлены в виде одномерного временного ряда ежемесячных данных в период с января 2006г. по декабрь 2008 г.
Естественный прирост
- разность между числами родившихся и умерших.
Актуальность исследования тесно связана со сложившейся тяжелой демографической ситуацией в России. Такая ситуация начала складываться в начале 90-х годов XX века и совпала с экономическим кризисом в стране.
Сокращение населения затронуло практически все территории Российской Федерации и почти все этнические группы. Депопуляция в Российской Федерации обусловливается не только низкой рождаемостью. Крайне острой проблемой является высокая смертность населения. Некоторое улучшение ситуации со смертностью в 1995-1998 годах оказалось непродолжительным. С 1999 года смертность населения страны вновь начала расти.
Данное исследование может послужить основой для формирования демографической политики в субъектах Российской Федерации и разработки региональных планов по улучшению демографической ситуации.
Предметом исследования являются ежемесячные данные об уровне естественного прироста населения в абсолютном выражении, т.к. по абсолютным данным модели и прогнозы являются более достоверными.
Целью данной курсовой работы является описание предмета исследования и эконометрических методов его анализа, моделирование и выявление тенденций изменения изучаемого показателя, прогнозирование и анализ полученных прогнозных значений.
В теоретической части курсовой работы подробно раскрывается содержание вопросов темы - предпрогнозная ориентация (программа исследования), включающая в себя анализ внешней обстановки; анализ внутренней обстановки; подробную характеристику объекта и предмета исследования; определение времени основания и учреждения в прогнозе; обоснование и описание методов, используемых в ходе исследования; организацию проведения исследования; формулирование рабочих гипотез.
В практической части проводится исследование временных рядов естественного прироста населения с помощью мультипликативной и аддитивной моделей, рядов Фурье, адаптивной сезонной модели, фиктивных переменных. По полученным данным определяется наиболее адекватная и точная модель, далее строится итоговый прогноз на период с января 2009г. по декабрь 2009г.
Для нахождения параметров, оценок моделей, проведения тестов на адекватность и значимость, а также для прогнозирования использовались пакеты MicrosoftExcel и STATISTICA.
Программа исследования естественного прироста населения в РФ в период с января 2006г. по декабрь 2008г. предполагает содержание подробной характеристики проблемы состоящую из следующих пунктов:
· анализ внешней обстановки (характеристика экономической, политической, социальной ситуации в государстве);
· анализ внутренней обстановки (тенденции внутри объекта исследования);
· подробная характеристика объекта и предмета исследования;
· определение времени основания и упреждения в прогнозе: текущий, краткосрочный;
· обоснование и описание методов, используемых в ходе исследования;
· формулирование рабочих гипотез.
Остановимся подробнее на каждом пункте исследования.
В России в основном завершен переход к рыночной экономической системе. Создана система базовых правовых норм и других институтов, обеспечивающих развитие рыночных отношений.
Достигнута высокая степень открытости российской экономики. Внешнеторговый оборот в 2007 году составил 45 процентов валового внутреннего продукта, что является одним из наиболее высоких показателей для стран с развитой экономикой.
В целом обеспечена макроэкономическая стабильность. Экономика защищена от внешних шоковых воздействий международными резервными активами Российской Федерации.
Сформировался мощный слой развивающихся компаний, успешно конкурирующих на внутреннем и внешнем рынках и активно привлекающих капитал для своего развития. Российский фондовый рынок стал важным фактором привлечения инвестиций и обеспечения экономического роста страны. В условиях развивающегося мирового финансового кризиса российская финансовая система (при активной поддержке государства) показала свою устойчивость.
Преодолены тенденции социальной конфронтации в обществе, наблюдавшиеся в 90-е годы. Развиваются институты гражданского общества.
Показатели социально-экономического развития Российской Федерации, достигнутая макроэкономическая стабильность и финансовая устойчивость свидетельствуют о результативности проведенных системных преобразований. Последние 5 лет темпы роста российской экономики значительно превышают среднемировой уровень (около 107 процентов и 104,6 процента соответственно).
Основанием для прогноза служит упорядоченная ежемесячная выборка, предоставленная Федеральной службой государственной статистики, в период с января 2006г. по декабрь 2008г. и содержащая 36 наблюдений.
Временем упреждения прогноза является период с января по декабрь 2009 г. – прогноз является краткосрочным.
Для данного исследования выбран достаточно широкий диапазон моделей исследования:
· Аддитивная модель;
· Мультипликативная модель;
· Одномерный анализ Фурье;
· Регрессионная модель с переменной структурой (фиктивные переменные);
· Адаптивная сезонная модель.
Выбор производился исходя из особенностей выборки – в исходных данных присутствует тенденция временного ряда, а также наблюдается наличие периодических колебаний. Выбранные модели достаточно хорошо описывают такие процессы.
Для аддитивной и мультипликативной моделей временного ряда необходимо произвести сглаживание методом простой скользящей средней.
Оценка параметров каждой из модели производится методом наименьших квадратов (МНК).
В период с января по декабрь 2009г. динамика изменения естественного прироста населения будет иметь возрастающую тенденцию и, следовательно, значения показателя будут увеличиваться.
Рассмотрим график временного ряда исходных данных естественного прироста населения РФ в период с января 2006 по декабрь 2008 года (Приложение 1). Проанализировав график, делаем вывод о наличии сезонных колебаний с периодичностью 12 месяцев и возрастающей тенденцией, что наглядно отражено в построенном графике сезонной волны (Приложение 2). Подтверждение данному факту отражено в АКФ и ЧАКФ (Таблица 1).
Таблица 1 - Значения АКФ и ЧАКФ
Лаг
АКФ
ЧАКФ
1
0,664
0,664
2
0,537
0,173
3
0,337
-0,135
4
0,242
0,011
5
0,065
-0,164
6
-0,058
-0,120
7
-0,051
0,155
8
0,044
0,225
9
0,104
0,067
10
0,216
0,152
11
0,247
-0,031
12
0,369
0,132
13
0,208
-0,300
14
0,162
-0,025
15
0,024
-0,036
Наибольшее значение достигается на 1 лаге, следовательно, присутствует тенденция временного ряда. Выбросы по АКФ – 1 и 12 лаг, по ЧАКФ – 1 и 13 лаг – гипотеза о сезонных колебаниях с периодичностью 12 месяцев подтверждается. Качество каждой модели будем оценивать по показателям среднеквадратической ошибки и средней ошибки аппроксимации. После построения всех моделей сделаем по каждой из них прогноз и проанализируем полученные результаты.
По графику временного ряда можно установить наличие приблизительно равной амплитуды колебаний. Это свидетельствует о соответствии этого ряда аддитивной модели. Рассчитаем ее компоненты.
Расчетная таблица модели приведена в Приложении 3.
Шаг 1.
Проведем выравнивание исходных уровней ряда методом простой скользящей средней. Для этого:
· Просуммируем уровни ряда последовательно за каждые 12 месяцев со сдвигом на один момент времени и определим условные годовые объемы показателя;
· Разделив полученные суммы на 12, найдем скользящие средние. Отметим, что полученные таким образом выравненные значения уже не содержат сезонной компоненты;
· Приведем эти значения в соответствие с фактическими моментами времени, для чего найдем средние значения из двух последовательных скользящих средних – центрированные скользящие средние.
Шаг 2.
Найдем оценки сезонной компоненты как разность между фактическими уровнями временного ряда и центрированными скользящими средними. Используем эти оценки для расчета значений сезонной компоненты S. Для этого найдем средние за каждый месяц (по всем годам) оценки сезонной компоненты Si. В аддитивной модели это выражается в том, что сумма значений сезонной компоненты по всем месяцам должна быть равна 0.
Таким образом получены следующие значения скорректированной сезонной компоненты (Таблица 2):
Таблица 2 - Значения скорректированной сезонной компоненты
Январь
S1
-20480,000
Июль
S7
13869,083
Февраль
S2
92,000
Август
S8
17015,208
Март
S3
-10291,958
Сентябрь
S9
14070,708
Апрель
S4
-6639,813
Октябрь
S10
5001,646
Май
S5
-11262,333
Ноябрь
S11
-142,500
Июнь
S6
-355,333
Декабрь
S12
-876,708
Занесем полученные значения для соответствующих месяцев каждого года.
Шаг 3.
Элиминируем влияние скорректированной сезонной компоненты, вычитая ее значение из каждого уровня исходного временного ряда. Получим: T + E = Y – S. Эти значения рассчитываются для каждого момента времени и содержат только тенденцию и случайную компоненту.
Шаг 4.
Определим компоненту T данной модели. Для этого проведем аналитическое выравнивание ряда (T + E) с помощью линейного тренда. Результаты выравнивания следующие:
T = -63923,013 + 1156,975 ·t; R2 = 0,889
Таблица 3 - Статистика уравнения тренда
Коэффициенты
Стандартная ошибка
t-статистика
P-Значение
Y-пересечение
-63923,013
1487,618
-42,970
3,23251E-31
t
1156,975
70,114
16,501
8,46925E-18
Уравнение описывает на 88,9% вариацию исходного показателя естественного прироста, при этом уравнение является статистически значимым при уровне надежности 95%.
Подставляя в это уравнение значения t = 1,…, 36, найдем уровни T для каждого момента времени.
Шаг 5.
Найдем значения уровней ряда, полученные по аддитивной модели. Для этого прибавим к уровням T значения сезонной компоненты для соответствующих месяцев.
Шаг 6.
В соответствие с методикой построения аддитивной модели расчет абсолютной ошибки производится по формуле:
E = Y – (T + S).
Исходя из значений выше приведенных показателей качества, можно сделать вывод о том, что модель обладает высокой точностью и пригодна для прогнозирования.
Шаг 1.
Проведем выравнивание исходных уровней временного ряда методом простой скользящей средней. Методика, применяемая для мультипликативной модели, полностью совпадает с методикой аддитивной модели. Расчетная таблица модели приведена в Приложении 4.
Шаг 2.
Найдем оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние. Используем эти оценки для расчета значений скорректированной сезонной компоненты. Для этого найдем средние за каждый месяц оценки сезонной компоненты. Взаимопогашаемость сезонных воздействий в мультипликативной модели выражается в том, что сумма значений скорректированной сезонной компоненты по всем кварталам должна быть равно числу периодов в цикле, т.е. двенадцати, так как в нашем случае число периодов одного цикла равно 12 месяцам. В результате имеем следующие сезонные компоненты: 1,526 + 1,010 + 1,280 + 1,183 + 1,329 + 1,032 + 0,712 + 0,631 + 0,680 + 0,874 + 1,008 + 1,021 = 12,287.
Определим скорректированные значения сезонной компоненты, умножив ее на средние оценки на корректирующий коэффициент k:
Получим следующие значения скорректированной сезонной компоненты (Таблица 4):
Таблица 4 - Значения скорректированной сезонной компоненты
Январь
S1
1,491
Июль
S7
0,695
Февраль
S2
0,987
Август
S8
0,616
Март
S3
1,250
Сентябрь
S9
0,664
Апрель
S4
1,155
Октябрь
S10
0,853
Май
S5
1,298
Ноябрь
S11
0,985
Июнь
S6
1,008
Декабрь
S12
0,997
Занесем полученные значения для соответствующих месяцев каждого года.
Шаг 3.
Разделим каждый уровень исходного ряда на соответствующие значения скорректированной сезонной компоненты. Получим: T·E = Y / S, которые содержат только тенденцию и случайную компоненту.
Шаг 4.
Определим компоненту T в мультипликативной модели. Для этого рассчитаем параметры линейного тренда, используя уровни (T·E). Результаты аналитического выравнивания:
Т = -64018,071 + 1201,064 · t; R2 = 0,863
Таблица 5 - Статистика уравнения тренда
Коэффициенты
Стандартная ошибка
t-статистика
P-Значение
Y-пересечение
-64018,071
1744,938
-36,688
6,24E-29
t
1201,064
82,242
14,604
3,27E-16
Уравнение описывает на 86,3% вариацию исходного показателя естественного прироста, при этом уравнение является статистически значимым при уровне надежности 95%.
Подставляя в это уравнение значения t = 1,…,36, найдем уровни Т для каждого момента времени.
Шаг 5.
Найдем уровни ряда по мультипликативной модели, умножив уровни Т на значения скорректированной сезонной компоненты для соответствующих месяцев.
Шаг 6.
Расчет ошибки в мультипликативной модели производится по формуле:
E = Y / (T·S)
Исходя из значений выше приведенных показателей качества, можно сделать вывод о том, что модель обладает высокой точностью и пригодна для прогнозирования.
2.4 Одномерный анализ Фурье
Выполним одномерный анализ Фурье для показателя естественного прироста населения РФ. Расчетная таблица ряда Фурье представлена в Приложении 5.
Переменными для составления модели будут следующие: t, cos(2·Π·t/12); sin(2·Π·t/12); cos(4·Π·t/12); sin(4·Π·t/12); cos(6·Π·t/12); sin(6·Π·t/12).
Значение знаменателя каждой дроби обусловлено периодичностью сезонных колебаний.
Построив модель с включением данных переменных, получаем следующее уравнение:
Таблица 6 - Статистика уравнения для модели ряда Фурье
Коэффициенты
Стандартная ошибка
t-статистика
P-Значение
Y-пересечение
-64314,412
2390,342
-26,906
1,49932E-21
t
1178,132
113,994
10,335
4,62652E-11
cos(2Pi*t/12)
-5360,004
1595,524
-3,359
0,002267582
sin(2Pi*t/12)
-12253,175
1647,329
-7,438
4,22248E-08
cos(4Pi*t/12)
-4098,437
1595,524
-2,569
0,015828755
sin(4Pi*t/12)
1894,178
1603,648
1,181
0,247471866
cos(6Pi*t/12)
933,424
1595,524
0,585
0,56321529
sin(6Pi*t/12)
-5109,257
1595,524
-3,202
0,003385848
Уравнение описывает на 89,7% вариацию исходного показателя естественного прироста, при этом уравнение является статистически значимым при уровне надежности 95%.
Но коэффициенты перед переменными sin(4·Π·t/12) и cos(6·Π·t/12) не удовлетворяют данному уровню надежности.
Исключим их из модели и перестроим уравнение регрессии.
Таблица 7 - Статистика уравнения для модели ряда Фурье
Коэффициенты
Стандартная ошибка
t-статистика
P-Значение
Y-пересечение
-64096,083
2361,646
-27,140
1,138E-22
t
1166,330
112,370
10,379
1,90859E-11
cos(2Pi*t/12)
-5348,202
1588,773
-3,366
0,002101931
sin(2Pi*t/12)
-12297,219
1639,342
-7,501
2,31486E-08
cos(4Pi*t/12)
-4086,636
1588,773
-2,572
0,015299504
sin(6Pi*t/12)
-5121,059
1588,773
-3,223
0,003049779
Уравнение описывает на 89,1% вариацию исходного показателя естественного прироста, уравнение статистически значимо при уровне надежности 95%. Все коэффициенты уравнения статистически значимы при аналогичном уровне надежности.
Модель имеет высокие показатели среднеквадратической ошибки и средней ошибки аппроксимации, но может быть использована для прогнозирования.
Рассмотрим еще один метод моделирования временного ряда, содержащего сезонные колебания, - построение модели регрессии с включением фактора времени и фиктивных переменных. Количество фиктивных переменных в такой модели должно быть на единицу меньше числа моментов (периодов) времени внутри одного цикла колебаний. В данном случае при моделировании ежемесячных данных модель должна включать двенадцать независимых переменных – фактор времени и одиннадцать фиктивных переменных. Каждая фиктивная переменная отражает сезонную (циклическую) компоненту временного ряда для какого-либо одного периода. Она равна единице для данного периода и нулю для всех остальных периодов.
Построим модель регрессии с включением фактора времени и фиктивных переменных для данных о естественном приросте населения в РФ. В данной модели двенадцать независимых переменных: t, D2, D3, D4, D5, D6, D7, D8, D9, D10, D11, D12 и результативная переменная Y. Составим матрицу исходных данных (Приложение 6).
Уравнение описывает на 96,0% вариацию исходного показателя естественного прироста, уравнение статистически значимо при уровне надежности 95%.
Оценим параметры уравнения регрессии обычным МНК. Результаты оценки приведены в Таблице 8.
Таблица 8 - Статистика уравнения для модели с фиктивными переменными
Коэффициенты
Стандартная ошибка
t-статистика
P-Значение
Y-пересечение
-89444,083
2879,238
-31,065
2,76051E-20
t
1132,083
79,218
14,291
6,2844E-13
D2
24047,583
3803,309
6,323
1,8825E-06
D3
17218,167
3805,783
4,524
0,000152385
D4
21431,750
3809,903
5,625
1,00291E-05
D5
15077,333
3815,664
3,951
0,000634609
D6
26904,583
3823,058
7,037
3,59459E-07
D7
40734,833
3832,075
10,630
2,38467E-10
D8
43809,083
3842,705
11,401
6,10242E-11
D9
38606,667
3854,934
10,015
7,43321E-10
D10
32848,917
3868,747
8,491
1,52119E-08
D11
26662,833
3884,126
6,865
5,33422E-07
D12
24437,083
3901,054
6,264
2,16207E-06
Проанализируем эти результаты. Все коэффициенты уравнения и само уравнение статистически значимы при уровне надежности 95%.
Исходя из значений выше приведенных показателей качества, можно сделать вывод о том, что модель обладает высокой точностью и пригодна для прогнозирования.
Рассмотрим аддитивную модель сезонных явлений с линейным ростом, предложенную Г. Тейлом и С. Вейджем. Параметры адаптации определим методом последовательных итераций, исходя из принципа минимизации средней ошибки аппроксимации модели. В результате получим следующие значения: α1
= 0,9; α2
= 0,1; α3
= 0,1.
Тренд – линейный, уравнение тренда выглядит следующим образом:
T = -67660,089 + 1358,979 ·t; R2 = 0,579
Начальные условия для нулевого цикла представлены в таблице 9:
Таблица 9 -Начальные условия
i
ĝi0
i
ĝi0
1
-24733,642
7
14639,816
2
-912,954
8
17487,170
3
-7969,267
9
12057,857
4
-3982,580
10
6073,211
5
-10563,892
11
-339,768
6
1036,462
12
-2792,414
Исходя из значений выше приведенных показателей качества, можно сделать вывод о том, что модель обладает высокой точностью и пригодна для прогнозирования.
Рассмотрим прогнозные значения естественного прироста населения в РФ по вышеописанным моделям, сравним полученные значения с фактическими, и выберем наиболее адекватную и точную модель для целей прогнозирования (Таблица 10).
Для оценки точности каждого прогноза рассчитаем среднюю относительную ошибку прогноза по формуле:
Таблица 10 - Прогнозные значения
Адд. модель
Мультипл. модель
Ряд Фурье
Модель с фикт. переменными
Адапт. модель Тейла-Вейджа
Фактические значения
Январь 2009
-41595
-29184
-38887
-47557
-46805
-47976
Февраль 2009
-19866
-18134
-31056
-22377
-22944
-24401
Март 2009
-29093
-21475
-21699
-28075
-29994
-32121
Апрель 2009
-24284
-18459
-23375
-22729
-26006
-27017
Май 2009
-27749
-19177
-24958
-27951
-32588
-28463
Июнь 2009
-15685
-13678
-13849
-14992
-20989
-19821
Июль 2009
-304
-8598
-86
-30
-7384
-4237
Август 2009
3999
-6886
2590
4177
-4539
1050
Сентябрь 2009
2212
-6624
-348
106
-9968
-3263
Октябрь 2009
-5701
-7484
-426
-4519
-15951
-12170
Ноябрь 2009
-9688
-7452
-4684
-9573
-22368
-25891
Декабрь 2009
-9265
-6350
-17547
-10667
-24818
-25116
Средняя относит. ошибка прогноза (%)
66,260
111,627
56,422
62,296
74,758
-
Исходя из показателя средней относительной ошибки прогноза, можно сделать вывод о том, что показатель естественного прироста населения наиболее точно прогнозируется рядом Фурье.
В ходе работы было проведено моделирование и прогнозирование естественного прироста населения в РФ. Исследование было проведено с помощью следующих моделей:
· Аддитивная модель;
· Мультипликативная модель;
· Одномерный анализ Фурье;
· Регрессионная модель с переменной структурой (фиктивные переменные);
· Адаптивная сезонная модель.
Выдвинутая гипотеза о возрастающей тенденции динамики изменения естественного прироста населения в РФ в 2009 году подтверждается.
По каждой модели сделан прогноз на 2009 год, при этом следует отметить, что наиболее точный прогноз дает модель с использованием ряда Фурье, в тоже время вариацию исходного показателя наиболее точно описывает адаптивная модель Тейла – Вейджа. Также можно сделать вывод о том, что для получения наиболее достоверного прогноза показателя необходимо комбинировать прогнозные значения нескольких наиболее точных моделей.
3. Концепция демографической политики Российской Федерации на период до 2015 года // http://www.akdi.ru/econom/program/demogr.htm
4. Концепция долгосрочного социально-экономического развития Российской Федерации на период до 2020 года // http://www.youngscience.ru/753/820/978/index.shtml
5. МеньшоваИ.В. Моделирование экономических процессов методами регрессионного анализа.- Воскресенск: Издательский дом «Лира», 2009. - 113 с.
6. Президент России молодым ученым и специалистам // http://www.youngscience.ru/
7. Федеральная служба государственной статистики // http://www.gks.ru/
8. Эконометрика: учебник / И.И. Елисеева, С.В. Курышев, Т.В. Костеева и др.; под ред. И.И. Елисеевой. – 2-е изд. - М.: Финансы и статистика, 2007. – 576 с.