Автор работы: Пользователь скрыл имя, 15 Сентября 2011 в 13:15, курсовая работа
Цель: проведение эконометрического анализа заданных временных рядов для прогнозирования их значений.
1.Построить графики временных рядов. Для каждого временного ряда провести первичный статистический анализ, включая:
◦Вычисление среднего значения, дисперсии, меры разброса;
◦Вычисление автоковариационной и автокорреляционной функций;
◦Построение коррелограммы.
(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 |