Автор работы: Пользователь скрыл имя, 08 Декабря 2013 в 14:17, контрольная работа
В ходе работы был проведен эконометрический анализ двух временных рядов - Х1 – месячного уровня осадков и Х2 – среднемесячных удоев молока. В конце нашей работы был выполнен прогноз значений временных рядов на ближайшие пять месяцев. Можно предположить, что сделанный прогноз достаточно точен, так как ошибки в данной модели малы.
В ходе своих наблюдений Робинзон заметил, что удои резко сократились в некоторый момент времени. Он пришел к выводу, что необходимо построить новое пастбище для своих коз, и поэтому огородил новую местность. Это изменение привело к резкому скачку в удоях, что отразилось на временном ряде и вызвало структурное изменение – соответственно мы перешли к разделенной модели, которая и является оптимальной для составления прогнозов.
Вычислим характеристики данной последовательности: количество серий – ν, длину максимальной серии – τ. Для данных рядов они примут следующие значения:
ν (Х1) = 25 ν (Х2) = 18
τ (Х1) = 5 τ (Х2) = 8
Проверим удовлетворяют ли эти значения неравенствам:
ν (55) > 0,5(55 + 2 - 1,96 ) = 21,232
τ (55) < 1,43ln(55 + 1) = 5,76
В первом случае оба неравенства выполняются, то есть не отвергается, что говорит об отсутствии трендовой компоненты в разложении ряда. Во втором оба неравенства не выполняются, это означает, что во временном ряде удоев присутствуют неслучайные компоненты. Речь идет в первую очередь о трендовой компоненте.
Проверим ту же гипотезу, используя критерий «восходящих» и «нисходящих» серий, чтобы «уловить» постепенное смещение среднего значения не только монотонного, но и периодического характера.
Вместо исходных элементов временного ряда Х(t) сформируем последовательность знаков: +, если Х(t+1) > Х(t),
−, если Х(t+1) < Х(t).
Полученные результаты для временного ряда осадков оформим в виде таблицы:
t |
X1 |
t |
X1 |
t |
X1 |
t |
X1 |
||||
1 |
29,2 |
+ |
15 |
30,6 |
+ |
29 |
30 |
- |
43 |
24,6 |
+ |
2 |
32,6 |
- |
16 |
31,2 |
- |
30 |
24,7 |
- |
44 |
32,7 |
- |
3 |
29,9 |
+ |
17 |
31,2 |
- |
31 |
22,1 |
+ |
45 |
31,9 |
+ |
4 |
31,3 |
- |
18 |
25 |
- |
32 |
33 |
- |
46 |
38,8 |
- |
5 |
31,3 |
- |
19 |
23,4 |
+ |
33 |
31,8 |
+ |
47 |
35,3 |
- |
6 |
24,8 |
- |
20 |
32,4 |
- |
34 |
39,3 |
- |
48 |
27 |
+ |
7 |
23,3 |
+ |
21 |
32,1 |
+ |
35 |
35,7 |
- |
49 |
28,2 |
+ |
8 |
32,2 |
- |
22 |
38,6 |
- |
36 |
27 |
+ |
50 |
32,5 |
- |
9 |
31,8 |
+ |
23 |
35,1 |
- |
37 |
27,8 |
+ |
51 |
30,8 |
+ |
10 |
38,7 |
- |
24 |
26,7 |
+ |
38 |
32,5 |
- |
52 |
32 |
- |
11 |
34,8 |
- |
25 |
28 |
+ |
39 |
30,3 |
+ |
53 |
30,7 |
- |
12 |
26,7 |
+ |
26 |
33,5 |
- |
40 |
31,8 |
- |
54 |
24,9 |
- |
13 |
27,7 |
+ |
27 |
31,6 |
+ |
41 |
30,7 |
- |
55 |
24,1 |
- |
14 |
33,1 |
- |
28 |
32,6 |
- |
42 |
24,9 |
- |
56 |
Для временного ряда удоев результаты имеют вид:
t |
X2 |
t |
X2 |
t |
X2 |
t |
X2 |
||||
1 |
11,2 |
- |
15 |
6,6 |
+ |
29 |
12 |
- |
43 |
5,6 |
+ |
2 |
9,8 |
- |
16 |
7,2 |
- |
30 |
8,5 |
- |
44 |
7,1 |
- |
3 |
8,6 |
- |
17 |
6,7 |
- |
31 |
7,3 |
+ |
45 |
7,1 |
+ |
4 |
8,4 |
- |
18 |
5,5 |
- |
32 |
9,2 |
- |
46 |
8,3 |
- |
5 |
8,1 |
- |
19 |
5,3 |
+ |
33 |
8,1 |
+ |
47 |
7,6 |
- |
6 |
6,8 |
- |
20 |
7 |
- |
34 |
9,7 |
- |
48 |
5,8 |
+ |
7 |
5,9 |
+ |
21 |
6,8 |
+ |
35 |
8,7 |
- |
49 |
6 |
+ |
8 |
7,7 |
- |
22 |
8,1 |
- |
36 |
6,5 |
+ |
50 |
6,8 |
- |
9 |
7,6 |
+ |
23 |
7,6 |
- |
37 |
6,8 |
+ |
51 |
6,5 |
+ |
10 |
8,7 |
- |
24 |
5,3 |
+ |
38 |
7,6 |
- |
52 |
7 |
- |
11 |
7,7 |
- |
25 |
5,7 |
+ |
39 |
7,2 |
- |
53 |
6,3 |
- |
12 |
5,8 |
+ |
26 |
7,1 |
- |
40 |
6,9 |
- |
54 |
5,1 |
- |
13 |
6,6 |
+ |
27 |
6,6 |
- |
41 |
6,9 |
- |
55 |
4,9 |
|
14 |
7,5 |
- |
28 |
6,4 |
+ |
42 |
5,6 |
- |
56 |
Вычислим характеристики данной последовательности: количество серий – ν, длину максимальной серии – τ. Для данных рядов они примут следующие значения:
ν (Х1) = 36 ν (Х2) = 31
τ (Х1) = 3 τ (Х2) = 8
Проверим, удовлетворяют ли эти значения неравенствам:
ν (55) > 1/3(55×2 – 1) - 1,96 = 30,31
τ (55) < 6
В первом случае выполняются оба неравенства, это означает, что не отвергается, то есть не отвергается предположении об отсутствии неслучайных компонент в разложении временного ряда. Однако графический и первичный статистический анализ отвергают данное предположение. Поэтому предположим, что в разложении временного ряда присутствует сезонность с периодом 12 месяцев.
Во втором одно из условий не выполняется, поэтому гипотеза отвергается, т.е. в модели присутствуют неслучайные компоненты. Данный критерий говорит о присутствии сезонной или циклической компоненты.
3.2. ПОСТРОЕНИЕ УРАВНЕНИЙ
В ходе анализа было выявлено, что в разложении временного ряда осадков присутствует только сезонная компонента. Значит, аддитивная модель данного временного ряда будет иметь следующий вид:
Как было определено выше, в уравнении временного ряда удоев присутствуют как трендовая, так и сезонная компоненты. Тогда аддитивная модель данного временного ряда выглядит следующим образом:
Далее оценим сезонные и трендовые компоненты временных рядов.
3.3. ОЦЕНКА СЕЗОННОЙ КОМПОНЕНТЫ
С помощью метода скользящих средних оценим сезонную компоненту во временном ряду осадков.
Учитывая, что при графическом и статистическом анализе был определен период колебаний, равный году, а также исходя из условия задачи, установим период скольжения в 12 месяцев. Рассчитаем скользящие средние и центрируем их. Результаты расчетов представлены в таблице:
X1 |
Xск ср |
Xцент |
X1 |
Xск ср |
Xцент |
X1 |
Xск ср |
Xцент |
X1 |
Xск ср |
Xцент |
29,2 |
31,2 |
30,575 |
30,571 |
22,1 |
30,775 |
30,767 |
38,8 |
30,767 |
30,775 | ||
32,6 |
31,2 |
30,567 |
30,579 |
33 |
30,758 |
30,717 |
35,3 |
30,783 |
30,783 | ||
29,9 |
25 |
30,592 |
30,592 |
31,8 |
30,675 |
30,621 |
27 |
30,783 |
30,783 | ||
31,3 |
23,4 |
30,592 |
30,604 |
39,3 |
30,567 |
30,533 |
28,2 |
30,783 |
30,763 | ||
31,3 |
32,4 |
30,617 |
30,633 |
35,7 |
30,500 |
30,529 |
32,5 |
30,742 |
|||
24,8 |
32,1 |
30,650 |
30,692 |
27 |
30,558 |
30,567 |
30,8 |
||||
23,3 |
30,550 |
30,488 |
38,6 |
30,733 |
30,792 |
27,8 |
30,575 |
30,679 |
32 |
||
32,2 |
30,425 |
30,446 |
35,1 |
30,850 |
30,800 |
32,5 |
30,783 |
30,771 |
30,7 |
||
31,8 |
30,467 |
30,496 |
26,7 |
30,750 |
30,738 |
30,3 |
30,758 |
30,763 |
24,9 |
||
38,7 |
30,525 |
30,521 |
28 |
30,725 |
30,671 |
31,8 |
30,767 |
30,746 |
24,1 |
||
34,8 |
30,517 |
30,513 |
33,5 |
30,617 |
30,642 |
30,7 |
30,725 |
30,708 |
|||
26,7 |
30,508 |
30,517 |
31,6 |
30,667 |
30,654 |
24,9 |
30,692 |
30,692 |
|||
27,7 |
30,525 |
30,529 |
32,6 |
30,642 |
30,671 |
24,6 |
30,692 |
30,708 |
|||
33,1 |
30,533 |
30,542 |
30 |
30,700 |
30,725 |
32,7 |
30,725 |
30,725 |
|||
30,6 |
30,550 |
30,563 |
24,7 |
30,750 |
30,763 |
31,9 |
30,725 |
30,746 |
|||
33 |
30,583 |
30,563 |
30,1 |
30,617 |
30,633 |
33,2 |
30,775 |
30,779 |
|||
31,1 |
30,542 |
30,567 |
25,2 |
30,650 |
30,650 |
31,9 |
30,783 |
30,758 |
Найдем разность между элементами ряда и центрированными средними. С помощью осреднения определим оценки сезонной компоненты. Рассчитанные разности, а также значения сезонной компоненты приведены в таблице:
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
- |
- |
- |
- |
- |
- |
-7,188 |
1,754 |
1,304 |
8,179 |
4,287 |
-3,817 |
-2,829 |
2,558 |
0,038 |
0,629 |
0,621 |
-5,592 |
-7,204 |
1,767 |
1,408 |
7,808 |
4,300 |
-4,038 |
-2,671 |
2,858 |
0,946 |
1,929 |
-0,725 |
-6,063 |
-8,667 |
2,283 |
1,179 |
8,767 |
5,171 |
-3,567 |
-2,879 |
1,729 |
-0,463 |
1,054 |
-0,008 |
-5,792 |
-6,108 |
1,975 |
1,154 |
8,025 |
4,517 |
-3,783 |
-2,563 |
- |
- |
- |
- |
- |
- |
|||||
-2,735 |
2,382 |
0,174 |
1,204 |
-0,038 |
-5,815 |
-7,292 |
1,945 |
1,261 |
8,195 |
4,569 |
-3,801 |
Для выполнения условия взаимопогашаемости сезонных воздействий приведем сумму значений сезонной компоненты к нулю, отниманием от каждого значения величины:
K = 0,048611/12 = 0,004051
Получим следующие значения сезонной компоненты:
-2,739 |
2,378 |
0,170 |
1,200 |
-0,042 |
-5,819 |
-7,296 |
1,941 |
1,257 |
8,191 |
4,565 |
-3,805 |
Таким образом, получена оценка сезонной компоненты временного ряда осадков.
Для определения постоянной компоненты данного ряда исключим из него сезонность:
t |
X1(t)-S(t) |
t |
X1(t)-S(t) |
t |
X1(t)-S(t) |
t |
X1(t)-S(t) |
t |
X1(t)-S(t) |
t |
X1(t)-S(t) |
1 |
31,939 |
11 |
30,235 |
21 |
30,843 |
31 |
29,396 |
41 |
30,742 |
51 |
30,630 |
2 |
30,222 |
12 |
30,505 |
22 |
30,409 |
32 |
31,059 |
42 |
30,719 |
52 |
30,800 |
3 |
29,730 |
13 |
30,439 |
23 |
30,535 |
33 |
30,543 |
43 |
31,896 |
53 |
30,742 |
4 |
30,100 |
14 |
30,722 |
24 |
30,505 |
34 |
31,109 |
44 |
30,759 |
54 |
30,719 |
5 |
31,342 |
15 |
30,430 |
25 |
30,739 |
35 |
31,135 |
45 |
30,643 |
55 |
31,396 |
6 |
30,619 |
16 |
30,000 |
26 |
31,122 |
36 |
30,805 |
46 |
30,609 |
56 |
|
7 |
30,596 |
17 |
31,242 |
27 |
31,430 |
37 |
30,539 |
47 |
30,735 |
57 |
|
8 |
30,259 |
18 |
30,819 |
28 |
31,400 |
38 |
30,122 |
48 |
30,805 |
58 |
|
9 |
30,543 |
19 |
30,696 |
29 |
30,042 |
39 |
30,130 |
49 |
30,939 |
59 |
|
10 |
30,509 |
20 |
30,459 |
30 |
30,519 |
40 |
30,600 |
50 |
30,122 |
60 |