Построение моделей временных рядов

Автор работы: Пользователь скрыл имя, 15 Сентября 2011 в 13:15, курсовая работа

Описание работы

Цель: проведение эконометрического анализа заданных временных рядов для прогнозирования их значений.


1.Построить графики временных рядов. Для каждого временного ряда провести первичный статистический анализ, включая:
◦Вычисление среднего значения, дисперсии, меры разброса;
◦Вычисление автоковариационной и автокорреляционной функций;
◦Построение коррелограммы.

Файлы: 1 файл

Курсач Яша.docx

— 174.04 Кб (Скачать файл)

(N) =     (55)=6 

V(N) > 30,306;    τ(N) < 6 => гипотеза не отвергается , следовательно неслучайные компоненты могут и не присутствовать. 

Критерий  Фостера-Стьюарта 

mt =  

lt =  

dt = mt – lt , rt = mt + lt ; 

D = = -2 

R = = 10 

tD = , tR, 

= = 2,6775

= = 2,1423

= 2 = 3.6227

tD = -0.747;  tR = 1.32146 

| tD | < 2,0057

| tR| < 2.0057, следовательно, гипотеза не отвергается, а значит в структуре временного ряда отсутствуют трендовые компоненты. 
 
 

     Временной ряд X2 (t) 

Критерий  серий. 

Xmed =   
 
 

Хmed = 6,8; V(N) = 20; τ(N) = 7

Условия: τ(N) ≥ 1,43Ln(N+1); V(N) ≤ 0,5(N+2-1,96), 

Если выполняется  хотя бы одно условие, то гипотеза отвергается.

τ(N) > 5,756; V(N) < 21,298 , следовательно гипотеза отвергается, а значит неслучайные компоненты присутствуют. 

Критерий  «восходящих» и «нисходящих» серий 

“+”: x(t+1) – x(t) >0  
“-”: x(t+1) – x(t) < 0
 

V(N) = 35; τ(N) = 3

Условия: V(N) < ; τ(N) ≥ τ0 (N),

Если выполняется  хотя бы одно условие, то гипотеза  отвергается. 

(N) =     (55)=6 

V(N) > 30,306; τ(N) < 6 => гипотеза не отвергается , следовательно неслучайные компоненты могут и не присутствовать. 

Критерий  Фостера-Стьюарта 

mt =  

lt =  

dt = mt – lt , rt = mt + lt ; 

D = = -8

R = = 12 

tD = , tR, 

= = 2,6775

= = 2,1423 

= 2 = 7,169

tD = -2,98785;  tR = 2,25504 

| tD | > 2,0057

| tR| > 2.0057, следовательно, гипотеза отвергается, а значит в структуре временного ряда присутствуют трендовые компоненты.

В силу того, что критерий «восходящих» и «нисходящих» серий говорит о возможном отсутствии неслучайных компонент во временном ряде Х2, проведем еще одну проверку ряда на наличие неслучайной компоненты при помощи критерия Аббе.

Критерий  Аббе 

Данный  критерий служит для выявления систематического смещения значений временного ряда. Для  проверки гипотезы Н0 необходимо вычислить статистику вида

,

где

Гипотеза  Н0 отвергается с вероятностью ошибки 0,05, если выполняется условие

,

Где - квантиль стандартного нормального распределения для доверительной вероятности . 

Для модели времянного ряда X2: 

= 2,32944; 

,

,

 – неравенство выполняется, следовательно гипотеза Н0 отвергается, а значит делаем вывод о наличии в структуре временного ряда Х2 неслучайных, зависящих от времени компонент. 
 

     Проверив  данные временные ряды на наличие  неслучайной компоненты можно сделать  следующие выводы:

     При проверке гипотезы о наличии неслучайных  компонент во временном ряду месячного  уровня осадков получились, что  по критерию серий и критерию восходящих и нисходящих серий, критерию Фостера-Cтьюарта неслучайная компонента в модели отсутствует.

     При проверке гипотезы о наличии неслучайных  компонент во временном ряду среднемесячного  удоя молока получились, что  по критерию серий и по критерию Аббе неслучайная компонента в модели присутствует. По критерию восходящих и нисходящих серий неслучайная компонента может и не присутствовать. Критерий Фостера-Стьюарта показывает наличие трендовой компоненты. В целом мы можем сделать вывод о присутствии неслучайной компоненты во временном ряду среднемесячного удоя молока.

 

Временной ряд X1(t) 

Шаг 1. Произведем выравнивание исходного временного ряда методом скользящей средней. Найдем годовую добычу дичи Робинзоном, после чего найдем скользящие средние. Для приведения в соответствие с фактическими моментами времени, найдем центрированные скользящие средние. 

месяц кол-во скользящие средние за год центрированные  скользящие средние оценка сезонной компоненты
1 27,50      
2 31,90      
3 29,80      
4 31,70      
5 29,60      
6 25,00 30,45    
7 24,10 30,46667 30,45833 -6,36
,,, ,,, ,,, ,,, ,,,
49 27,30 30,61667 30,68333 -3,38333
50 33,20      
51 30,30      
52 32,10      
53 30,80      
54 24,60      
55 22,50      
 

Шаг 2. Оценим сезонные компоненты как разность между фактическими элементами ряда и центрированными скользящими средними.

Скорректированные оценки сезонной компоненты определяются путем вычитания из средней оценки сезонной компоненты для месяца корректирующего  коэффициента. 

k =  0,010475 

год\мес. 1 2 3 4 5 6 7 8 9 10 11 12
1             -6,35833 2,583333 1,2 8,070833 4,258333 -3,32917
2 -2,87083 2,5625 0,066667 1,0625 -0,57083 -6,29583 -7,2375 2,333333 1,129167 8,1125 4,904167 -3,5625
3 -1,98333 1,9625 0,379167 0,804167 -0,34583 -4,925 -7,77917 2,375 1,4 8,016667 4,579167 -5,33333
4 -2,87917 2,191667 -0,0875 1,091667 0,345833 -5,53333 -6,66667 3,520833 1,208333 7,7 4,991667 -4,775
5 -3,38333                      
Siсредн -2,77917 2,238889 0,119444 0,986111 -0,19028 -5,58472 -7,01042 2,703125 1,234375 7,975 4,683333 -4,25
Si -2,78964 2,228414 0,10897 0,975637 -0,20075 -5,5952 -7,02089 2,69265 1,2239 7,964525 4,672859 -4,26047
 
 

Шаг 3. Устраним сезонную компоненту S из исходных уравнений ряда и получим выровненные данные T+ =X(t)-S

Месяц Х(t) Si X(t) - Si
1 27,50 -2,78964 30,29
2 31,90 2,228414 29,67159
3 29,80 0,10897 29,69103
53 30,80 -0,20075 31,00075
54 24,60 -5,5952 30,1952
55 22,50 -7,02089 29,52089

T=X(t) -S

 

     После исключения сезонной компоненты из временного ряда можно говорить об отсутствии трендовой компоненты в разложении временного ряда. 

Временной ряд X2(t) 

Шаг 1. Произведем выравнивание исходного временного ряда методом скользящей средней. Найдем годовую добычу дичи Робинзоном, после чего найдем скользящие средние. Для приведения в соответствие с фактическими моментами времени, найдем центрированные скользящие средние. 

Месяц Кол-во Скользящие средние за год Центрированные  скользящие средние Оценка сезонной компоненты
1 11,00      
2 9,70      
3 8,50      
4 8,70      
5 7,60      
6 6,70 8,066667    
7 6,10 7,65 7,858333 -1,76
,,, ,,, ,,, ,,, ,,,
49 6,50 7,166667 7,25 -0,75
50 7,50      
51 6,90      
52 7,60      
53 6,60      
54 5,30      
55 4,80      

Информация о работе Построение моделей временных рядов