Построение, анализ и применение для прогнозирования тренд-сезонной модели временного ряда
Предварительный анализ заданного временного ряда на предмет наличия тренда. Обоснование наличия сезонности по графическому представлению одноименных элементов ряда разных лет. Применение модели для прогноза. Выбор типа остатков и корректировка модели.
Рубрика | Программирование, компьютеры и кибернетика |
Вид | контрольная работа |
Язык | русский |
Дата добавления | 12.09.2011 |
Размер файла | 218,8 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ
Факультет экономики и менеджмента
Кафедра Экономической кибернетики
Контрольная работа
По дисциплине
Методы прогнозирования
Тема работы
Построение, анализ и применение для прогнозирования тренд-сезонной модели временного ряда
Симферополь, 2010г.
Задание
1) Предварительный анализ заданного временного ряда на предмет наличия тренда:
a) графическое представление ряда;
b) укрупнение;
c) сглаживание.
2) Обоснование выбора типа модели основного тренда по свойствам кривых роста. Экономический смысл параметров выбранной модели.
3) Обоснование наличия сезонности по графическому представлению одноименных элементов ряда разных лет. Рекомендуется предварительно исключить из исходного ряда основную тенденцию развития (тренд, укрупненные или сглаженные значения) вычитанием или делением. Описать характер сезонных колебаний и по возможности объяснить причины.
4) Построение модели декомпозиции временного ряда. Обоснование формы модели: аддитивная, мультипликативная или смешанная.
5) Анализ точности модели.
6) Проверка адекватности моделей данному временному ряду. Анализ свойств остаточной компоненты.
7) Применение модели для прогноза (точечный и интервальный прогноз на два периода сезонности с заданной доверительной вероятностью Р = 0,95).
8) Верификация прогноза (сопоставление с контрольными данными). Выводы о возможности дальнейшего применения модели.
9) Возможно, изменение модели.
Исходные данные
Житлові будівлі, тис. м2 загальної площі |
|||
I кв |
|
1234,7 |
|
I-II кв |
|
2709,4 |
|
I-III кв |
|
4280,0 |
|
за год |
|
7566,3 |
|
I кв |
|
1453,3 |
|
I-II кв |
|
2893,6 |
|
I-III кв |
|
4524,0 |
|
за год |
|
7816,0 |
|
I кв |
|
1690,9 |
|
I-II кв |
|
3188,7 |
|
I-III кв |
|
5046,6 |
|
за год |
|
8628,4 |
|
I кв |
|
1761,1 |
|
I-II кв |
|
3947,6 |
|
I-III кв |
|
6130,3 |
|
за год |
|
10243,7 |
|
I кв |
|
1897,2 |
|
I-II кв |
|
4044,6 |
|
I-III кв |
|
6483,4 |
|
за год |
|
10495,6 |
Предварительный анализ заданного временного ряда на предмет наличия тренда и сезонности
Тренд - это длительная тенденция изменения экономических показателей, т.е. изменение, определяющее общее направление развития, основную тенденцию временных рядов. Тренд характеризует основные закономерности движения во времени, в некоторой мере свободные от случайных воздействий. При разработке моделей прогнозирования тренд оказывается основной составляющей прогнозируемого временного ряда, на которую накладываются другие составляющие. Результат при этом связывается исключительно с ходом времени. Предполагается, что через время можно выразить влияние всех основных факторов. Если тренд является монотонным (устойчиво возрастает или устойчиво убывает), то анализировать такой ряд обычно нетрудно. Если временные ряды содержат значительную ошибку, то первым шагом выделения тренда является сглаживание.
Представляем исходные значения на точечном графике
Выполним сглаживание
Сглаживание всегда включает некоторый способ локального усреднения данных, при котором несистематические компоненты взаимно погашают друг друга. Самый общий метод сглаживания - скользящее среднее, в котором каждый член ряда заменяется простым или взвешенным средним n соседних членов, где n - ширина окна. Вместо среднего можно использовать медиану значений, попавших в окно.
Для построения модели тренда проведем сглаживание временного ряда с периодом 4. Для этого воспользуемся формулой средней хронологической. В нашем случае данная формула будет выглядеть следующим образом.
Где t = 3,…,n-2; к = 4 - длина цикла сезонности.
Обоснование выбора типа модели основного тренда по свойствам кривых роста
Построив 4 прироста по сглаженному ряду наблюдаем следующую картину
Только четвертые приросты колеблются вокруг постоянного уровня, скорее всего мы можем взять полином третьей степени.
С помощью функции ЛИНЕЙН построим уравнение модели:
ЛИНЕЙН |
-0,533959337 |
19,1 |
-147,7185881 |
2244,95 |
|
0,114276663 |
3,6 |
34,83892337 |
98,01 |
||
0,987638656 |
34,4 |
#Н/Д |
#Н/Д |
||
319,5894121 |
12,0 |
#Н/Д |
#Н/Д |
||
1135607,793 |
14213,3 |
#Н/Д |
#Н/Д |
||
Уравнение модели: |
y = -,53x3+19,1x2-147,72x+2244,95 |
Выбор типа модели: аддитивная, мультипликативная или смешанная
Общая модель. Основная идея сезонной декомпозиции проста. В общем случае временной ряд типа того, который описан выше, можно представить себе состоящим из трех различных компонент: (1) сезонной компоненты (обозначается St, где t обозначает момент времени), (2) тренда (Tt), (3) циклической компоненты (Ct). Разница между циклической и сезонной компонентой состоит в том, что последняя имеет регулярную (сезонную) периодичность, тогда как циклические факторы обычно имеют более длительный эффект, который к тому же меняется от цикла к циклу. Конкретные функциональные взаимосвязи между этими компонентами могут иметь самый разный вид. Однако, можно выделить два основных способа, с помощью которых они могут взаимодействовать: аддитивно и мультипликативно.
Аддитивная модель:
Xt = TCt + St + It
Мультипликативная модель:
Xt = TCt*St*It
Здесь Xt обозначает значение временного ряда в момент времени t. Если имеются какие-то априорные сведения о циклических факторах, влияющих на ряд (например, циклы деловой конъюнктуры), то можно использовать оценки для различных компонент для составления прогноза будущих значений ряда.
тренд сезонность временный ряд
Так как график не имеет ярко выраженной склонности к возрастанию, а колеблется примерно в одном и том же диапазоне, я взял аддитивную модель.
Построим аддитивную модель.
Для того чтобы определить порядок сезонности, выполним ряд вычислений:
Квартал |
Сумма разностей |
Число повторений квартала |
Средняя разность |
Вычетание среднего |
|
I |
-3076,17 |
5 |
-615,23 |
-596,00 |
|
II |
-2462,47 |
5 |
-492,49 |
-473,26 |
|
III |
-1655,05 |
5 |
-331,01 |
-311,77 |
|
IV |
6808,96 |
5 |
1361,79 |
1381,03 |
|
|
|
Среднее |
-19,24 |
0,00 |
Сезонная компонента S Квартал |
S |
|
I |
-596,00 |
|
II |
-473,26 |
|
III |
-311,77 |
|
IV |
1381,03 |
Построим график Y и T+S чтобы посмотреть на качество аддитивной модели
Выбор типа остатков и корректировка модели
Для каждого значения t определим ошибку модели (остаточную компоненту временного ряда) по следующей формуле:
Для того, чтобы выбрать тип остатков: e = y - y или e = (y/y) - 1, построим оба графика
Так как на первом графике димперсия более постоянна, выберем e = y - y.
С помощью процедуры «Поиск решения» определим новые параметры модели
Новые параметры модели: |
|
|||
-0,414011172 |
13,53396 |
-76,8736 |
2000,126 |
y = -0,41*x3+13,53*x2-76,87*x+2000,13
Проверка гипотезы о нормальном законе распределения остатков по ассиметрии и эксцессу:
0,47 |
|||||||
0,22 |
|||||||
0,59 |
|||||||
-0,76 |
0,76 |
||||||
H0: |
остатки распределены в соответствии с нормальным законом |
||||||
0,22 |
< |
0,71 |
0,55 |
< |
1,14 |
||
H1: |
остатки не распределены в соответствии с нормальным законом |
||||||
Следовательно, остатки нормально распределены |
Проверка постоянства дисперсии остатков:
n1 = |
10 |
D1 = |
10892,72291 |
??? |
0,05 |
||
n2 = |
10 |
D2 = |
17834,30861 |
|
|
||
F = |
0,61077349 |
||||||
F1кр = |
4,025994158 |
Fкр = |
3,178893105 |
||||
F2кр = |
0,314574906 |
|
|
||||
H0: |
????????? |
H0: |
????????? |
||||
H1: |
?????????? |
H1: |
????????? |
||||
Следовательно, остатки гомокедастичны. |
Проверка среднего значения остатков:
t |
e |
e - e1ср |
e - e2ср |
||||
1 |
-0,079 |
-0,086 |
|
n1 = |
10 |
||
2 |
0,036 |
0,028 |
|
n2 = |
10 |
||
3 |
0,001 |
-0,006 |
|
||||
4 |
0,007 |
-0,001 |
|
e1ср = |
0,007623145 |
||
5 |
0,112 |
0,105 |
|
e2ср = |
0,006 |
||
6 |
-0,016 |
-0,023 |
|
||||
7 |
-0,025 |
-0,032 |
-0,031 |
??? |
0,05 |
||
8 |
-0,038 |
-0,045 |
-0,044 |
||||
9 |
0,122 |
0,115 |
0,116 |
S = |
0,001102894 |
||
10 |
-0,118 |
-0,125 |
-0,124 |
||||
11 |
-0,037 |
-0,045 |
-0,043 |
H0: |
??????? |
||
12 |
-0,030 |
-0,038 |
-0,036 |
H1: |
???????? |
||
13 |
-0,012 |
-0,020 |
-0,018 |
||||
14 |
0,111 |
0,104 |
0,105 |
t = |
2,33343192 |
||
15 |
0,000 |
-0,008 |
-0,007 |
||||
16 |
0,049 |
0,042 |
0,043 |
t1кр = |
-2,100922037 |
||
17 |
-0,039 |
-0,047 |
-0,046 |
t2кр = |
2,100922037 |
||
18 |
0,016 |
0,008 |
0,010 |
||||
19 |
0,073 |
0,065 |
0,066 |
||||
20 |
0,017 |
0,009 |
0,011 |
||||
Суммкв: |
0,074666557 |
0,054814471 |
|||||
Среднее значение остатков непостоянно |
Проверка гипотезы об отсутствии автокорелляции:
? |
r |
n - ? |
t |
t1кр |
t2кр |
|
|
1 |
-0,27 |
19 |
-1,14 |
-2,11 |
2,11 |
+ |
|
2 |
-0,02 |
18 |
-0,08 |
-2,12 |
2,12 |
+ |
|
3 |
-0,11 |
17 |
-0,43 |
-2,13 |
2,13 |
+ |
|
4 |
-0,09 |
16 |
-0,34 |
-2,14 |
2,14 |
+ |
|
5 |
0,17 |
15 |
0,63 |
-2,16 |
2,16 |
+ |
|
6 |
-0,19 |
14 |
-0,69 |
-2,18 |
2,18 |
+ |
|
7 |
0,17 |
13 |
0,59 |
-2,20 |
2,20 |
+ |
|
8 |
-0,54 |
12 |
-2,01 |
-2,23 |
2,23 |
+ |
Автокорелляция присутствует.
Применение модели для прогноза
Прогноз на основании трендовых моделей (кривых роста) содержит два элемента: точечный и интервальный прогнозы. Точечный прогноз -- это прогноз, которым называется единственное значение прогнозируемого показателя. Это значение определяется подстановкой в уравнение выбранной кривой роста величины времени t, соответствующей периоду упреждения: t=n+1; t=n+2 и т. д. Очевидно, что точное совпадение фактических данных в будущем и прогностических точечных оценок маловероятно. Поэтому точечный прогноз должен сопровождаться двусторонними границами, т.е. интервальным прогнозом.
Интервальный прогноз на базе трендовых моделей осуществляется путем расчета доверительного интервала -- такого интервала, в котором с определенной вероятностью можно ожидать появления фактического значения прогнозируемого экономического показателя. Методы, разработанные для статистических совокупностей, позволяют определить доверительный интервал, зависящий от стандартной ошибки оценки прогнозируемого показателя, от времени упреждения прогноза, от количества уровней во временном ряду и от уровня значимости (ошибки) прогноза.
Точечный прогноз на два периода сезонности:
|
t |
t2 |
t3 |
T |
S |
y |
Y |
|
2004 |
1 |
1 |
1 |
1936,372498 |
-596,00 |
1340,38 |
1234,7 |
|
|
2 |
4 |
8 |
1897,202698 |
-473,26 |
1423,94 |
1474,6 |
|
|
3 |
9 |
27 |
1880,132678 |
-311,77 |
1568,36 |
1570,7 |
|
|
4 |
16 |
64 |
1882,678372 |
1381,03 |
3263,71 |
3286,3 |
|
2005 |
5 |
25 |
125 |
1902,355711 |
-596,00 |
1306,36 |
1453,3 |
|
|
6 |
36 |
216 |
1936,68063 |
-473,26 |
1463,42 |
1440,3 |
|
|
7 |
49 |
343 |
1983,16906 |
-311,77 |
1671,40 |
1630,4 |
|
|
8 |
64 |
512 |
2039,336936 |
1381,03 |
3420,37 |
3292,0 |
|
2006 |
9 |
81 |
729 |
2102,700189 |
-596,00 |
1506,70 |
1690,9 |
|
|
10 |
100 |
1000 |
2170,774753 |
-473,26 |
1697,52 |
1497,8 |
|
|
11 |
121 |
1331 |
2241,076561 |
-311,77 |
1929,30 |
1857,9 |
|
|
12 |
144 |
1728 |
2311,121546 |
1381,03 |
3692,15 |
3581,8 |
|
2007 |
13 |
169 |
2197 |
2378,425641 |
-596,00 |
1782,43 |
1761,1 |
|
|
14 |
196 |
2744 |
2440,504778 |
-473,26 |
1967,25 |
2186,5 |
|
|
15 |
225 |
3375 |
2494,874891 |
-311,77 |
2183,10 |
2182,7 |
|
|
16 |
256 |
4096 |
2539,051913 |
1381,03 |
3920,08 |
4113,4 |
|
2008 |
17 |
289 |
4913 |
2570,551777 |
-596,00 |
1974,56 |
1897,2 |
|
|
18 |
324 |
5832 |
2586,890415 |
-473,26 |
2113,63 |
2147,4 |
|
|
19 |
361 |
6859 |
2585,583761 |
-311,77 |
2273,81 |
2438,8 |
|
|
20 |
400 |
8000 |
2564,147747 |
1381,03 |
3945,18 |
4012,2 |
|
2009 |
21 |
441 |
9261 |
2520,098307 |
-596,00 |
1924,10 |
|
|
|
22 |
484 |
10648 |
2450,951373 |
-473,26 |
1977,69 |
|
|
|
23 |
529 |
12167 |
2354,222879 |
-311,77 |
2042,45 |
|
|
|
24 |
576 |
13824 |
2227,428757 |
1381,03 |
3608,46 |
|
|
2010 |
25 |
625 |
15625 |
2068,08494 |
-596,00 |
1472,09 |
|
|
|
26 |
676 |
17576 |
1873,707362 |
-473,26 |
1400,45 |
|
|
|
27 |
729 |
19683 |
1641,811955 |
-311,77 |
1330,04 |
|
|
|
28 |
784 |
21952 |
1369,914653 |
1381,03 |
2750,94 |
|
Графическое представление:
Интервальный прогноз на два периода сезонности:
L = |
8 |
??? |
0,05 |
Sy = |
0,064905 |
||
n = |
20 |
t? = |
2,100922037 |
|
|
|
t |
t2 |
t3 |
t4 |
y |
Uy |
||
2004 |
1 |
1 |
1 |
1 |
1340,376 |
1340,376049 |
1340,376 |
|
|
2 |
4 |
8 |
16 |
1423,945 |
1423,944764 |
1423,945 |
|
|
3 |
9 |
27 |
81 |
1568,359 |
1568,358863 |
1568,359 |
|
|
4 |
16 |
64 |
256 |
3263,707 |
3263,706569 |
3263,707 |
|
2005 |
5 |
25 |
125 |
625 |
1306,359 |
1306,359263 |
1306,359 |
|
|
6 |
36 |
216 |
1296 |
1463,423 |
1463,422696 |
1463,423 |
|
|
7 |
49 |
343 |
2401 |
1671,395 |
1671,395245 |
1671,395 |
|
|
8 |
64 |
512 |
4096 |
3420,365 |
3420,365133 |
3420,365 |
|
2006 |
9 |
81 |
729 |
6561 |
1506,704 |
1506,70374 |
1506,704 |
|
|
10 |
100 |
1000 |
10000 |
1697,517 |
1697,516819 |
1697,517 |
|
|
11 |
121 |
1331 |
14641 |
1929,303 |
1929,302746 |
1929,303 |
|
|
12 |
144 |
1728 |
20736 |
3692,15 |
3692,149743 |
3692,15 |
|
2007 |
13 |
169 |
2197 |
28561 |
1782,429 |
1782,429192 |
1782,429 |
|
|
14 |
196 |
2744 |
38416 |
1967,247 |
1967,246844 |
1967,247 |
|
|
15 |
225 |
3375 |
50625 |
2183,101 |
2183,101076 |
2183,101 |
|
|
16 |
256 |
4096 |
65536 |
3920,08 |
3920,080111 |
3920,08 |
|
2008 |
17 |
289 |
4913 |
83521 |
1974,555 |
1974,555328 |
1974,555 |
|
|
18 |
324 |
5832 |
104976 |
2113,632 |
2113,632481 |
2113,632 |
|
|
19 |
361 |
6859 |
130321 |
2273,81 |
2273,809946 |
2273,81 |
|
|
20 |
400 |
8000 |
160000 |
3945,176 |
3945,175944 |
3945,176 |
|
2009 |
21 |
441 |
9261 |
194481 |
1924,102 |
1924,27099 |
1923,933 |
|
|
22 |
484 |
10648 |
234256 |
1977,693 |
1977,868154 |
1977,519 |
|
|
23 |
529 |
12167 |
279841 |
2042,449 |
2042,630093 |
2042,268 |
|
|
24 |
576 |
13824 |
331776 |
3608,457 |
3608,645039 |
3608,269 |
|
2010 |
25 |
625 |
15625 |
390625 |
1472,088 |
1472,284381 |
1471,893 |
|
|
26 |
676 |
17576 |
456976 |
1400,449 |
1400,653871 |
1400,245 |
|
|
27 |
729 |
19683 |
531441 |
1330,038 |
1330,251881 |
1329,824 |
|
|
28 |
784 |
21952 |
614656 |
2750,943 |
2751,166625 |
2750,719 |
|
? |
210 |
2870 |
44100 |
722666 |
|
|
|
|
|
406 |
7714 |
164836 |
3756718 |
|
|
|
Графическое представление
Размещено на Allbest.ru
Подобные документы
"Наивная" модель прогнозирования. Прогнозирование методом среднего и скользящего среднего. Метод опорных векторов, деревьев решений, ассоциативных правил, системы рассуждений на основе аналогичных случаев, декомпозиции временного ряда и кластеризации.
курсовая работа [2,6 M], добавлен 02.12.2014Исходный текст программы и ее экранная форма. Программа вычисления и выдачи на печать суммы/произведения элементов бесконечного числового ряда, вычисления числового ряда для известного числа членов ряда. Значение максимального элемента в матрице.
контрольная работа [29,0 K], добавлен 07.12.2010Последовательность разработки чертежа и модели с типоразмерами из параметрического ряда. Построение таблицы переменных в соответствии с исходными данными. Проектирование параметрической модели в системе Компас-3D, внешние переменные для чертежа детали.
практическая работа [5,9 M], добавлен 14.04.2016Проверка наличия линейной связи между соответствующими показателями деятельности коммерческих банков Украины в модуле Multiple Regression ППП Statistica. Расчет теоретических значений зависимой переменной и ошибки модели, вид графика линейной функции.
лабораторная работа [1,5 M], добавлен 19.05.2011Особенности и основные этапы разработка компьютерной модели расчета траектории неуправляемого летательного аппарата, анализ модели динамики. Метод Эйлера как линейное приближение, использующее первые два члена ряда Тейлора, способы решения задач.
курсовая работа [930,6 K], добавлен 14.01.2013Анализ матрицы коэффициентов парной корреляции. Выбор факторных признаков для построения двухфакторной регрессионной модели. Оценка параметров регрессии по методу наименьших квадратов. Нахождение определителей матриц. Применение инструмента Регрессия.
контрольная работа [1,0 M], добавлен 13.01.2013Приведена оптимизация расходов и трудозатрат персонала. Реализация модели ARIMA (модели Бокса-Дженкинса), являющейся интегрированной композицией метода авторегрессии и модели скользящего среднего. Применение средств программного продукта Matlab 2013a.
дипломная работа [876,7 K], добавлен 19.09.2019Построение концептуальной модели и метод имитационного моделирования. Определение переменных уравнений математической модели и построение моделирующего алгоритма. Описание возможных улучшений системы и окончательный вариант модели с результатами.
курсовая работа [79,2 K], добавлен 25.06.2011Применение комплексного математического моделирования в проектировании. Обзор численных методов в моделировании. Решение дифференциальных уравнений в MathCAD. Анализ исходных и результирующих данных. Описание реализации базовой модели в MathCAD.
курсовая работа [240,5 K], добавлен 18.12.2011Общая характеристика ателье "Вита", схема модели рабочего процесса. Исследование заданной системы с помощью моделирования динамических рядов, модели типа "система массового облуживания". Построение имитационной модели деятельности данного ателье.
курсовая работа [1,4 M], добавлен 01.06.2016