Автор: Пользователь скрыл имя, 20 Марта 2011 в 16:56, курсовая работа
Целью моей курсовой работы является проведение экономического анализа заданных временных рядов для прогнозирования их значений.
Исходя из поставленной цели, мною были сформулированы следующие задачи:
1. Построить графики временных рядов.
2. Провести первичный статистический анализ временных рядов.
3. Построить модель временного ряда
4. Сделать выводы по всем полученным результатам.
5. Вычислить прогнозные значения по наиболее оптимальной модели.
ВВЕДЕНИЕ…………………………………………………………………………..4
1 ПЕРВИЧНЫЙ СТАТИСТИЧЕСКИЙ АНАЛИЗ……………………………..….6
1.1 Построение графиков временных рядов……………………………………….6
1.2 Вычисление среднего значения, дисперсии, меры разброса………………….8
1.3Вычисление автоковариационной и автокорреляционной функции. Построение коррелограммы………………………………………………………...8
2 ПОСТРОЕНИЕ МОДЕЛИ ВРЕМЕННОГО РЯДА……………………….…….11
2.1 Построение модели для неслучайных компонент……………………...…….11
2.2 Проверка остатков на автокорреляцию с помощью критерия Дарбина – Уотсона…………………………………………………………………………...…23
2.3 Структурные изменения……………………………………………………….25
ЗАКЛЮЧЕНИЕ…………………………………………………………………….32
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ………………………………34
Сформулируем гипотезы:
структурная стабильность
структурное изменение.
1 уравнение: Т(1) =
t:1;2;…;22
T=X-S: 11,548; 9,797;…; 6,941.
Сделаем матричные обозначения и найдём , , .
Оценка
неизвестных параметров
2 уравнение: Т(2) =
t:23;24;…;55
T=X-S:
10,955; 9,589;…;6,186.
Сделаем матричные обозначения и найдём , , .
Оценка неизвестных параметров
В таблице 2.13 приведены расчеты по тесту Чоу.
Таблица 2.11
Тест Чоу
Уравнение | Число наблюдений | ESS | Число параметров | Число степеней свободы |
(1) | 22 | 4,776 | 3 | 19 |
(2) | 33 | 7,409 | 3 | 30 |
(3) | 55 | 38,135 | 3 | 52 |
Остаточная сумма квадратов для кусочной модели [1]:
= + = 12,185
Соответствующее ей число степеней свободы составит:
= 49
Сокращение остаточной дисперсии при переходе от единого уравнения тренда к кусочной модели можно определить следующим образом:
∆ESS = ESS3 – ESS = 38,135-12,185=25,95 0.
Число степеней свободы, соответствующее , будет равно: .
Найдем
F–статистику по следующей формуле:
при k=3, N=55.
Fkp (k, N-2*k)= Fkp(3,49) = 2,8
Так как , следовательно, гипотеза (о структурной стабильности ряда) отвергается при уровне значимости . Можно сделать вывод о том, что влияние структурных изменений на динамику изучаемого показателя признается значимым.
Таким
образом, моделирование тенденции
временного ряда, следует осуществлять
с помощью кусочной модели (Рис.2.7).
Рис.
2.7. Изменение характера тенденции
временного ряда
Проверка ряда на наличие структурных изменений с помощью подхода Гуйарати [2].
Этот
подход основан на использование
фиктивной переменной, которая принимает
два значения: 0 и 1.
Общее уравнение Гуйарати имеет вид:
Yt =
При
, где - параметры 1-го уравнения.
При
, где - параметры 2-го уравнения
- разница между свободными членами уравнений 1) и 2): если θ1 – значимый, то изменение – значимо.
т.е. есть разница между параметрами θ11 и θ12 уравнений (1) и (2): если θ3 – значимый, то изменение – значимо.
т.е. есть разница между параметрами θ21 и θ22 уравнений (1) и (2): если θ5 – значимый, то изменение – значимо.
Оценка статистической значимости различий , а также и эквивалентна оценке статистической значимости параметров уравнения Гуйарати.
Чтобы сделать проверку значимости сначала найдём
t:1;2;…;55
T=X-S: 11,548;9,797 ;…;6,186.
Сделаем
матричные обозначения и найдём
.
Оценка
неизвестных параметров
Проверим значимость параметров , в каждом уравнении регрессии, используя критерий Стьюдента.
Выдвигаем
гипотезы:
=
9,54176 | -9,541758 | -0,00212 | 0,002 | -222,8 | 222,8 |
-9,5418 | 9,7793689 | 0,00212 | -0,003 | 222,8 | -223,2 |
-0,0021 | 0,0021197 | 5E-07 | -5E-07 | 0,0483 | -0,048 |
0,00212 | -0,002806 | -5E-07 | 3E-06 | -0,048 | 0,049 |
-222,8 | 222,7983 | 0,04832 | -0,048 | 5297,4 | -5297 |
222,798 | -223,2377 | -0,04832 | 0,049 | -5297 | 5299 |
S2 = 0,2487
Таблица 2.12
Проверка значимости параметров уравнения
Коэффициенты | Стандартная ошибка | t-статистика | |
8,772 | 1,559 | 5,625 | |
-0,0022 | 9E-04 | -2,446 | |
-248,33 | 36,3 | -6,841 |
Таким образом, параметры , , признаются значимыми.
В
целом по подходу Гуйарати можно
сделать вывод о влияние
Прогнозирование
значений месячного уровня осадков
будем осуществлять по следующей
модели:
X1(t)=30,714+Si+E
Результаты расчетов представлены в таблице 2.15.
Таблица 2.13
Расчет прогнозных значений Х1
t | Si | а | а+S |
56 | 1,931 | 30,714 | 32,645 |
57 | 1,606 | 30,714 | 32,32 |
58 | 8,563 | 30,714 | 39,277 |
59 | 4,8981 | 30,714 | 35,6121 |
60 | -4,145 | 30,714 | 26,569 |
График результатов прогнозирования на 5 месяцев представлен на рисунке 2.8.
Рис.2.8.
Результаты прогнозирования для
месячного уровня осадков
Количество осадков зависит только от сезона года и график это полностью подтверждает.
Наилучшая модель, описывающая среднемесячные удои молока, является кусочная модель:
1)Х(t) = ,
2)Х(t) = ,
Таблица 2.14
Расчет прогнозных значений для Х2
t | Si | T | T+S |
56 | 0,02731 | 6,58715 | 6,6145 |
57 | 0,16377 | 6,64968 | 6,8135 |
58 | 1,55856 | 6,71746 | 8,2760 |
59 | 1,94502 | 6,79035 | 8,7354 |
60 | -0,3893 | 6,86822 | 6,4789 |
На
рисунке 2.9 представлен результат прогноза
на 5 месяцев.
Рис.2.9. Результаты прогнозирования для среднемесячных удоев молока
Таким образом, первичный статистический анализ показал, что в среднем на необитаемом острове выпадает 30,48 миллиметров осадков, среднемесячные удои молока за весь изучаемый период составляют 7,22 галлонов.
Проверка гипотезы о наличие неслучайных компонент в разложение временного ряда месячного уровня осадков подтвердила предположение, выдвинутое при анализе графика, об отсутствие трендовой составляющий. Также было выявлено, что отсутствует и сезонная компонента, но дальнейший анализ показал, что эта компонента все-таки присутствует в разложение данного временного ряда. Присутствие сезонной компоненты объясняется тем, что среднемесячное количество осадков зависит от времени года: в летний период на острове выпадает больше количество осадков по сравнению с зимним.