Анализ динамики производства древесноволокнистых плит в России за 2000-2009 годы
Автокорреляционная функция временного ряда темпов роста производства древесноволокнистых плит в Российской Федерации. Расчет значений сезонной компоненты в аддитивной модели и коэффициента автокорреляции третьего порядка по логарифмам уровней ряда.
Рубрика | Экономико-математическое моделирование |
Вид | контрольная работа |
Язык | русский |
Дата добавления | 15.11.2014 |
Размер файла | 300,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Размещено на http://www.allbest.ru/
КАЗАНСКИЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. А.Н. ТУПОЛЕВА - КАИ
ФИЛИАЛ «ВОСТОК»
КОНТРОЛЬНАЯ РАБОТА
ПО ДИСЦИПЛИНЕ «Эконометрика»
на тему «Анализ динамики производства древесноволокнистых плит в РФ за 2000-2009 гг.»
Выполнил:
Стаценко Н.А.
Чистополь
Имеем исходные данные о производстве древесноволокнистых плит в РФ за период 2000-2009 гг. (табл. 1). Требуется выбрать наилучший вид тренда и определить его параметры.
Таблица 1. Производство древесноволокнистых плит в РФ (миллионов м2)
2000 |
2005 |
2006 |
2007 |
2008 |
2009 |
||
Производство древесноволокнистых плит в РФ |
278 |
413 |
439 |
481 |
479 |
373 |
Задание 1
Проведем анализ динамики производства древесноволокнистых плит в РФ за период 2000-2009 гг.
1. Построим график данного временного ряда (рис. 1).
Рис. 1. Динамика темпов роста производства древесноволокнистых плит в РФ
На графике рис. 1 наглядно видно наличие возрастающей тенденции и небольшие колебания. 2. Для дальнейшего анализа определим коэффициенты автокорреляции по уровням ряда и их логарифмам (табл. 2).
Таблица 2. Автокорреляционная функция временного ряда темпов роста производства древесноволокнистых плит в РФ за 2000-2009 гг.
Лаг |
Автокорреляционная функция |
||
по уровням ряда |
по логарифмам уровней ряда |
||
1 |
0,198 |
0,201 |
|
2 |
-0,246 |
-0,218 |
|
3 |
-0,637 |
-0,629 |
Для нахождения коэффициента корреляции были проведены следующие расчеты: 1. Рассчитываем коэффициент автокорреляции первого порядка по уровням ряда (табл. 3).
Таблица 3. Коэффициент автокорреляции первого порядка
t |
yt |
yt-1 |
yt-y1c |
yt-1 - y2c |
(yt-y1c)*(yt-1-y2c) |
(yt-y1c)^2 |
(yt-1 - y2c)^2 |
|
1 |
- |
- |
- |
- |
- |
- |
- |
|
2 |
413 |
278 |
-24 |
-140 |
3360 |
576 |
19600 |
|
3 |
439 |
413 |
2 |
-5 |
-10 |
4 |
25 |
|
4 |
481 |
439 |
44 |
21 |
924 |
1936 |
441 |
|
5 |
479 |
481 |
42 |
63 |
2646 |
1764 |
3969 |
|
6 |
373 |
479 |
-64 |
61 |
-3904 |
4096 |
3721 |
|
сумма |
2185 |
2090 |
3016 |
8376 |
27756 |
|||
среднее |
437 |
418 |
R1 = 0,198
2. Рассчитываем коэффициент автокорреляции второго порядка по уровням ряда (табл. 4).
Таблица 4. Коэффициент автокорреляции второго порядка
t |
yt |
yt-2 |
yt-y3c |
yt-2 - y4c |
(yt-y3c)*(yt-2-y4c) |
(yt-y3c)^2 |
(yt-2 - y4c)^2 |
|
1 |
- |
- |
- |
- |
- |
- |
- |
|
2 |
- |
- |
- |
- |
- |
- |
- |
|
3 |
439 |
278 |
-4 |
-124,75 |
499 |
16 |
15562,563 |
|
4 |
481 |
413 |
38 |
10,25 |
389,5 |
1444 |
105,0625 |
|
5 |
479 |
439 |
36 |
36,25 |
1305 |
1296 |
1314,0625 |
|
6 |
373 |
481 |
-70 |
78,25 |
-5477,5 |
4900 |
6123,0625 |
|
сумма |
1772 |
1611 |
-3284 |
7656 |
23104,75 |
|||
среднее |
443 |
402,75 |
R2 = -0,246
3. Рассчитываем коэффициент автокорреляции третьего порядка по уровням ряда (табл. 5).
Таблица 5. Коэффициент автокорреляции третьего порядка
t |
yt |
yt-3 |
yt-y5c |
yt-3 - y6c |
(yt-y5c)*(yt-3-y6c) |
(yt-y3c)^2 |
(yt-3 - y6c)^2 |
|
1 |
- |
- |
- |
- |
- |
- |
- |
|
2 |
- |
- |
- |
- |
- |
- |
- |
|
3 |
- |
- |
- |
- |
- |
- |
- |
|
4 |
481 |
278 |
36,67 |
-98,67 |
-3617,78 |
1344,44 |
9735,11 |
|
5 |
479 |
413 |
34,67 |
36,33 |
1259,56 |
1201,78 |
1320,11 |
|
6 |
373 |
439 |
-71,33 |
62,33 |
-4446,44 |
5088,44 |
3885,44 |
|
сумма |
1333 |
1130 |
-6804,67 |
7634,67 |
14940,67 |
|||
среднее |
444,33 |
376,67 |
R3 = -0,637
4. Рассчитываем коэффициент автокорреляции по логарифмам уровней ряда (табл. 6).
Исходные данные по логарифмам уровней ряда
yt |
Lnyt |
|
278 |
5,63 |
|
413 |
6,02 |
|
439 |
6,08 |
|
481 |
6,18 |
|
479 |
6,17 |
|
373 |
5,92 |
Таблица 6. Коэффициент автокорреляции первого порядка по логарифмам уровней ряда
Lnyt |
lnyt-1 |
Lnyt - lny1c |
lnyt-1 - lny2 |
(lnyt-lny1) *(lnyt-1-lny2c) |
(lnyt-lny1c)^2 |
lnyt-1 - lny2c |
||
1 |
- |
- |
- |
- |
- |
- |
- |
|
2 |
6,02 |
5,63 |
-0,054 |
-0,386 |
0,020844 |
0,002916 |
0,148996 |
|
3 |
6,08 |
6,02 |
0,006 |
0,004 |
0,000024 |
3,6E-05 |
1,6E-05 |
|
4 |
6,18 |
6,08 |
0,106 |
0,064 |
0,006784 |
0,011236 |
0,004096 |
|
5 |
6,17 |
6,18 |
0,096 |
0,164 |
0,015744 |
0,009216 |
0,026896 |
|
6 |
5,92 |
6,17 |
-0,154 |
0,154 |
-0,02372 |
0,023716 |
0,023716 |
|
сумма |
30,37 |
30,08 |
0,01968 |
0,04712 |
0,20372 |
|||
среднее |
6,074 |
6,016 |
R1 = 0,201
5. Рассчитываем коэффициент второго порядка по логарифмам уровней ряда (табл. 7).
Таблица 7. Коэффициент автокорреляции первого порядка по логарифмам уровней ряда
Lnyt |
lnyt-2 |
Lnyt - lny3c |
lnyt-2 - lny4с |
(lnyt-lny3с) *(lnyt-1-lny4c) |
(lnyt-lny3c)^2 |
lnyt-2 - lny4c |
||
1 |
- |
- |
- |
|||||
2 |
- |
- |
- |
- |
- |
- |
- |
|
3 |
6,08 |
5,63 |
-0,0075 |
-0,3475 |
0,002606 |
5,625E-05 |
0,120756 |
|
4 |
6,18 |
6,02 |
0,0925 |
0,0425 |
0,003931 |
0,0085562 |
0,001806 |
|
5 |
6,17 |
6,08 |
0,0825 |
0,1025 |
0,008456 |
0,0068062 |
0,010506 |
|
6 |
5,92 |
6,18 |
-0,1675 |
0,2025 |
-0,03392 |
0,0280563 |
0,041006 |
|
сумма |
24,35 |
23,91 |
-0,01893 |
0,043475 |
0,174075 |
|||
среднее |
6,0875 |
5,9775 |
R2 = -0,218
6. Рассчитываем коэффициент автокорреляции третьего порядка по логарифмам уровней ряда (табл. 8).
автокорреляция древесноволокнистый плита логарифм
Таблица 8. Коэффициент автокорреляции третьего порядка по логарифмам уровней ряда
Lnyt |
lnyt-3 |
Lnyt - lny5c |
lnyt-3 - lny6с |
(lnyt-lny5с) *(lnyt-3-lny6c) |
(lnyt-lny6c)^2 |
lnyt-3 - lny6с |
||
1 |
- |
- |
- |
|||||
2 |
- |
- |
- |
- |
- |
- |
- |
|
3 |
- |
- |
- |
- |
- |
- |
- |
|
4 |
6,18 |
5,63 |
0,09 |
-0,28 |
-0,0252 |
0,0081 |
0,0784 |
|
5 |
6,17 |
6,02 |
0,08 |
0,11 |
0,0088 |
0,0064 |
0,0121 |
|
6 |
5,92 |
6,08 |
-0,17 |
0,17 |
-0,0289 |
0,0289 |
0,0289 |
|
сумма |
18,27 |
17,73 |
-0,0453 |
0,0434 |
0,1194 |
|||
среднее |
6,09 |
5,91 |
R3 = -0,629
Строим автокоррелограммы по уровням ряда (рис. 2) и по логарифмам уровней ряда (рис. 3).
Рис. 2. Автокоррелограмма по уровням ряда
Рис. 3. Автокоррелограмма по логарифмам уровней ряда
Низкие значения коэффициентов автокорреляции первого, второго и третьего порядков свидетельствуют о том, что ряд содержит тенденцию. Приблизительно равные значения коэффициентов автокорреляции по уровням ряда и логарифмам уровней позволяют сделать вывод: если ряд содержит нелинейную тенденцию, то она выражена в неявной форме. Поэтому для моделирования его тенденции целесообразно использовать линейный тренд. 3. Построим тренд на основе линейной функции yt = a+bt. Для этого строим вспомогательную таблицу (табл. 9).
Таблица 9. Вспомогательная таблица
Годы |
y |
t |
t^2 |
t*y |
yt |
|
2000 |
278 |
-3 |
9 |
-834 |
361,3214 |
|
2005 |
413 |
-2 |
4 |
-826 |
377,7143 |
|
2006 |
439 |
-1 |
1 |
-439 |
394,1071 |
|
2007 |
481 |
1 |
1 |
481 |
426,8929 |
|
2008 |
479 |
2 |
4 |
958 |
443,2857 |
|
2009 |
373 |
3 |
9 |
1119 |
459,6786 |
|
сумма |
2463 |
0 |
28 |
459 |
2463 |
Параметры линейного уравнения находим из системы уравнений.
a = b = = .
Записываем уравнение: yt = 410,5 + 16,4*t. Темпы роста производства древесноволокнистой плиты за 2000 - 2009 гг. изменялись от уровня 410,5 м2 со средним за год абсолютным приростом, равным 16,4 м2.
Задание 2
Построить аддитивную модель. Шаг 1. Проведем выравнивание исходных уровней ряда методом скользящей средней. Для этого: а) просуммируем уровни ряда последовательно за каждые два года со сдвигом на один момент времени и определим условные объемы производства древесноволокнистых плит (гр. 3 табл. 10); b) разделив полученные суммы на 2, найдем скользящие средние (гр. 4 табл. 10). Выровненные значения уже не содержат сезонной компоненты;c) приведем эти значения в соответствие с фактическими моментами времени, для чего найдем средние значения из двух последовательных скользящих средних - центрированные скользящие (гр. 5 табл. 10)
Таблица 10. Расчет оценок сезонной компоненты по аддитивной модели
Год t |
Производство древесн. плит |
Итого за 2 года |
Скользящая средняя за 2 года |
Центрированная скользящая средняя |
Оценка сезонной компоненты |
|
1 |
278 |
- |
- |
- |
- |
|
2 |
413 |
691 |
345,5 |
- |
- |
|
3 |
439 |
852 |
426 |
385,75 |
53,25 |
|
4 |
481 |
920 |
460 |
443 |
38 |
|
5 |
7831 |
15153 |
7576,5 |
6655,75 |
1175,25 |
|
6 |
5787 |
13618 |
6809 |
7192,75 |
-1405,75 |
Шаг 2. Найдем оценки сезонной компоненты как разность между фактическими уровнями ряда и центрированными скользящими средними (гр. 6 табл. 10). Используем эти оценки для расчета значений сезонной компоненты S (табл. 11). Для этого найдем средние за каждый год оценки сезонной компоненты S. В моделях с сезонной компонентной обычно предполагается, что сезонные воздействия за период взаимопогашаются. В аддитивной модели это выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна нулю.
Таблица 11. Расчет значений сезонной компоненты в аддитивной модели
показатели |
год |
|||
1 |
2 |
|||
1 |
- |
53,25 |
||
2 |
38 |
9 |
||
3 |
-80 |
- |
||
итого |
-42 |
62,25 |
||
средняя оценка |
-14 |
20,75 |
||
скорректированная |
-17,375 |
17,375 |
Для данной модели имеем: -14+20,75 = 6,75. Определим корректирующий коэффициент: k = 6.75/2 = 3.375. Рассчитаем скорректированные значения сезонной компоненты как разность между ее средней оценкой и корректирующим коэффициент k. Проверим условие равенства нулю суммы значений сезонной компоненты: -17,375 + 17,375 = 0. Таким образом, получены следующие значения сезонной компоненты:S1 = - 17.375; S2 = 17.375. Занесем значения в табл. 11.
Шаг 3. Элиминируем влияние сезонной компоненты, вычитая ее значение из каждого уровня исходного временного ряда. Получим величины T+E=Y-S (гр. 4 табл. 12). Эти значения рассчитываются за каждый момент времени и содержат только тенденцию и случайную компоненту.
Таблица 12. Расчет выровненных значений T и ошибок E в аддитивной модели
t |
yt |
Si |
T+E-yt-Si |
T |
T+S |
E |
E^2 |
|
1 |
278 |
-17,375 |
295,375 |
4967,714 |
4950,339286 |
-4672,34 |
21830754 |
|
2 |
413 |
17,375 |
395,625 |
7180,75 |
7198,125 |
-6785,13 |
46037921 |
|
3 |
439 |
-17,375 |
456,375 |
7606,964 |
7589,589286 |
-7150,59 |
51130927 |
|
4 |
481 |
17,375 |
463,625 |
8295,464 |
8312,839286 |
-7831,84 |
61337707 |
|
5 |
479 |
-17,375 |
496,375 |
8262,679 |
8245,303571 |
-7766,3 |
60315471 |
|
6 |
373 |
17,375 |
355,625 |
6525,036 |
6542,410714 |
-6169,41 |
38061629 |
|
сумма |
2463 |
278714409 |
||||||
среднее |
410,5 |
Шаг 4. Определяем компоненту T данной модели. Для этого проводили аналитическое выравнивание ряда (T+S) с помощью линейного тренда: T = 410,5 + 16,4*t. График уравнения тренда приведен на рис. 4
Рис. 4. Производство древесноволокнистых плит в РФ
Шаг 5. Найдем значения уровней ряда, полученные по аддитивной модели. Для этого прибавим к уровням T значения сезонной компоненты для соответствующих годов. Графические значения (T+S) представлены на рис. 4. Шаг 6. В соответствии с методикой построения аддитивной модели расчет ошибки производится по формуле E=Y - (T+S). Это абсолютная ошибка. Численные значения абсолютных ошибок приведены в гр. 7 табл. 12. По аналогии с моделью регрессии для оценки качества построения модели или для выбора наилучшей модели можно применять сумму квадратов полученных абсолютных ошибок. Для данной аддитивной модели сумма квадратов абсолютных ошибок равна 278714409. По отношению к общей сумме квадратов отклонений уровней ряда от его среднего уровня, равной (2463-410,5)2=4212756, эта величина составляет -6515,964: (1-278714409/4212756)*100 = - 6515,964.
Прогнозирование по аддитивной модели.
По данным требуется дать прогноз производства древесноволокнистых плит в РФ на ближайшие 3 года. Прогнозное значение Ft уровня временного ряда в аддитивной модели в соответствии с соотношением Y =T +S + E есть сумма трендовой и сезонной компонент. Для определения трендовой компоненты воспользуемся уравнением тренда T = 410,5 + 16,4*t. Получим: Т4 = 410,5 + 16,4 *4 = 476,1Т5 = 410,5 + 16,4 *5 = 492,5Т6 = 410,5 + 16,4 *6 = 508,9 Значения сезонной компоненты равны:S1 = - 17.375; S2 = 17.375.
Таким образом, F4 = T4+S1 = 476,1+(- 17.375) = 458,73, F5 = T5+S2 = 492,5+17,375= 509,88, F6 = T6+S1 = 508,9 + (-17.375) = 491,53
Прогноз объема производства древесноволокнистых плит в РФ на ближайшие 3 года составит: 458,73; 509,88; 491,53.
Размещено на Allbest.ru
Подобные документы
Анализ автокорреляции уровней временного ряда, характеристика его структуры; построение аддитивной и мультипликативной модели, отражающую зависимость уровней ряда от времени; прогноз объема выпуска товаров на два квартала с учетом выявленной сезонности.
лабораторная работа [215,7 K], добавлен 23.01.2011Построение графика временного ряда. Тренд - устойчивое систематическое изменение процесса в течение продолжительного времени. Динамика продаж бензина на АЗС. Выявление сезонной составляющей и тренда. Коррелограмма, построенная в программе Statistica.
курсовая работа [1,2 M], добавлен 15.11.2013Изучение понятия имитационного моделирования. Имитационная модель временного ряда. Анализ показателей динамики развития экономических процессов. Аномальные уровни ряда. Автокорреляция и временной лаг. Оценка адекватности и точности трендовых моделей.
курсовая работа [148,3 K], добавлен 26.12.2014Сглаживание с помощью метода скользящей средней. Анализ исходного ряда на наличие стационарности. Тест Дики-Фуллера. Выделение сезонной компоненты в аддитивной и мультипликативной модели. Составление уравнения тренда в виде полинома пятой степени.
лабораторная работа [2,6 M], добавлен 17.02.2014Аддитивная модель временного ряда. Мультипликативная модель временного ряда. Одномерный анализ Фурье. Регрессионная модель с переменной структурой. Сущность адаптивной сезонной модели Тейла – Вейджа. Прогнозирование естественного прироста населения.
курсовая работа [333,1 K], добавлен 19.07.2010Теория и анализ временных рядов. Построение линии тренда и прогнозирование развития случайного процесса на основе временного ряда. Сглаживание временного ряда, задача выделения тренда, определение вида тенденции. Выделение тригонометрической составляющей.
курсовая работа [722,6 K], добавлен 09.07.2019Этапы и проблемы эконометрических исследований. Параметры парной линейной регрессии. Оценка тесноты связи с помощью показателей корреляции и детерминации. Расчет коэффициентов автокорреляции второго порядка для временного ряда расходов на потребление.
контрольная работа [60,3 K], добавлен 05.01.2011Прогнозирование, его основные подходы и виды. Текущее состояние российского кинематографа, его проблемы и тенденции. Прогнозирование числа выходящих кинофильмов в Российской Федерации методом экстраполяции временного ряда и методом наименьших квадратов.
курсовая работа [280,0 K], добавлен 20.06.2014Эконометрическое моделирование стоимости квартир в московской области. Матрица парных коэффициентов корреляции. Расчет параметров линейной парной регрессии. Исследование динамики экономического показателя на основе анализа одномерного временного ряда.
контрольная работа [298,2 K], добавлен 19.01.2011Эффективная оценка по методу наименьших квадратов. Корелляционно-регрессионный анализ в эконометрическом моделировании. Временные ряды в эконометрических исследованиях. Моделирование тенденции временного ряда. Расчет коэффициента автокорреляции.
контрольная работа [163,7 K], добавлен 19.06.2015