Построили график каждого временного ряда.
а) График временного ряда для количества дичи:

В среднем количество убитых уток увеличивается (т.к. левая часть графика ниже правой, что говорит о наличии тренда);
Ярко прослеживаются сезонные колебания, так как функция в 4х точках достигает своего максимума в среднем через 12 месяцев и своего минимума в среднем через 12 месяцев, что говорит о сезонности.
б) график временного ряда для веса Робинзона:

При динамике веса выражены сезонные колебания, а роста среднего значения веса в течении 5 лет не прослеживается.
Провели первичный статистический анализ временных рядов, включая вычисление среднего значения, меры разброса. Сделали соответствующие выводы.
1) Х1ср = 75,69091
Х2ср = 75,74545
Значит в среднем Робинзон добывал ~76 дичи, а его вес в среднем составлял ~76 кг.
2)
- дисперсия временного ряда


3. Проверим гипотезу Н0 и ее альтернативу Н1:
H0(0): 
H0(1): 
H1(0): 
H1(1): 
Проверка гипотез с помощью критерия серий:
| Критерий серий | |||||||||
| t | X1 | t | X1 | t | X2 | t | X2 | ||
| - | - | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| - | |||||||||
| - | - | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| - | + | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| - | - | ||||||||
| - | + | ||||||||
| - | - | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| - | - | ||||||||
| - | - | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| - | + | ||||||||
| + | - | ||||||||
| + | |||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| + | |||||||||
| + | - | ||||||||
| + | + | ||||||||
| - | + | ||||||||
| + | |||||||||
| + | + | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | - | ||||||||
| + | + | ||||||||
| + | + |
Найдем статистики
и
и проверим неравенства:
1) 
2)
.
Если хоть одно из них не выполняются, то гипотеза
отвергается с вероятностью ошибки a, такой, что 0,05 < a < 0,0975, что подтверждается наличие зависящей от времени неслучайной составляющей в
.
-количество серий;
- максимальное количество элементов с серии;
·

·

1)
→ не выполняется (8<21,3);
2)
→ не выполняется (12>5,8);
1)
→ выполняется (30>21,3)
2)
→ выполняется (6>5,8)
Следовательно, гипотезы
и
отвергаются.
Проверка гипотезы с помощью критерия "восходящих" и "нисходящих" серий:
| Критерий восходящих и нисходящих серий | |||||
| Х1 | Х2 | ||||
| + | + | ||||
| - | + | ||||
| - | - | ||||
| - | + | ||||
| + | - | ||||
| - | + | ||||
| - | + | ||||
| - | + | ||||
| + | - | ||||
| - | + | ||||
| + | - | ||||
| - | + | ||||
| + | - | ||||
| + | + | ||||
| - | + | ||||
| - | - | ||||
| - | + | ||||
| - | + | ||||
| - | - | ||||
| - | - | ||||
| + | |||||
| + | + | ||||
| + | - | ||||
| + | + | ||||
| + | |||||
| - | + | ||||
| - | - | ||||
| - | + | ||||
| - | - | ||||
| + | - | ||||
| - | + | ||||
| + | + | ||||
| + | - | ||||
| + | - | ||||
| + | - | ||||
| + | + | ||||
| + | - | ||||
| + | - | ||||
| - | + | ||||
| - | - | ||||
| - | - | ||||
| - | + | ||||
| - | + | ||||
| + | - | ||||
| + | + | ||||
| + | - | ||||
| + | + | ||||
| + | - | ||||
| + | + | ||||
| + | - | ||||
| - | + | ||||
| - | + | ||||
| - | - | ||||
| - | + | ||||
| - | + | ||||
| v1(N)кр= | 30,30636 | v1(N)кр= | 30,30636 | ||
| t1(N)кр= | t1(N)кр= | ||||
| v1(N)= | отв | v1(N)= | не отв | ||
| t1(N)= | отв | t1(N)= | не отв | ||
| отвергается | не отвергается |
Найдем статистики
и
. Если не выполняется хоть одно из неравенств:
1)
;
2)
,
то гипотеза
отвергается.
Найдем:


1)
=6;
2)
не выполняется (18<30,3);
1)
не выполняется (7>6);
2)
выполняется (37>30,3);
выполняется (3<6)
Следовательно, гипотеза
отвергается, а гипотеза
не отвергается.
Задание 4. Построить уравнение для неслучайных компонент, присутствие которых в модели было доказано. Провести сравнительный анализ моделей А, В и С, выбирая типы моделей по таблице 2. Обосновать по результатам эконометрического анализа выбор наилучшей модели.
Построим уравнения для неслучайных компонент временного ряда
:
Шаг 1: Используя метод скользящего среднего, выровняем исходный ряд
. Для этого суммируем элементы ряда последовательно за каждые 12 месяцев со сдвигом на один момент времени и, разделив полученные суммы на 12, найдем скользящие средние. Для приведения в соответствие с фактическими моментами времени найдем центрированные скользящие средние. Получим
| скользящие средние за год | центророванные скользящие средние |
| 39,33333 | |
| 42,16667 | 40,75 |
| 45,58333 | 43,875 |
| 48,08333 | 46,83333 |
| 51,33333 | 49,70833 |
| 53,66667 | 52,5 |
| 55,25 | 54,45833 |
| 56,125 | |
| 58,91667 | 57,95833 |
| 59,45833 | |
| 62,08333 | 61,04167 |
| 62,75 | 62,41667 |
| 65,16667 | 63,95833 |
| 67,08333 | 66,125 |
| 67,58333 | 67,33333 |
| 69,66667 | 68,625 |
| 70,75 | 70,20833 |
| 70,875 | |
| 72,25 | 71,625 |
| 73,125 | |
| 77,41667 | 76,70833 |
| 79,08333 | 78,25 |
| 81,41667 | 80,25 |
| 82,83333 | 82,125 |
| 82,91667 | |
| 85,08333 | 84,04167 |
| 85,83333 | 85,45833 |
| 87,41667 | 86,625 |
| 89,75 | 88,58333 |
| 91,58333 | 90,66667 |
| 92,25 | 91,91667 |
| 94,16667 | 93,20833 |
| 95,5 | 94,83333 |
| 96,91667 | 96,20833 |
| 96,66667 | 96,79167 |
| 97,5 | 97,08333 |
| 98,25 | 97,875 |
| 98,66667 | 98,45833 |
| 100,0833 | 99,375 |
| 100,8333 | 100,4583 |
| 102,0833 | 101,4583 |
| 103,25 | 102,6667 |
| 105,0833 | 104,1667 |
Шаг 2: Найдем оценки сезонной компоненты S, которые определяются как разность между фактическими элементами ряда и центрированными скользящими средними.
Получим
| Шаг 2 | ||||||||||||
| k= | 0,10 | |||||||||||
| год\мес. | ||||||||||||
| -20,75 | -27,88 | -2,83 | -15,71 | 12,50 | 4,54 | |||||||
| 28,88 | 31,04 | 10,54 | 7,96 | -1,42 | -12,96 | -25,13 | -28,33 | -11,63 | -11,21 | 2,13 | 16,38 | |
| 34,88 | 20,00 | 18,29 | 3,75 | -16,25 | -16,13 | -20,92 | -21,04 | -11,46 | -7,63 | 12,42 | 14,33 | |
| 18,08 | 26,79 | 9,17 | 4,79 | -4,79 | -9,08 | -27,88 | -12,46 | -9,38 | -4,46 | -3,46 | 12,33 | |
| 14,83 | ||||||||||||
| Siсредн | 24,17 | 25,94 | 12,67 | 5,50 | -7,49 | -12,72 | -23,67 | -22,43 | -8,82 | -9,75 | 5,90 | 11,90 |
| Si | 24,07 | 25,84 | 12,57 | 5,40 | -7,59 | -12,82 | -23,77 | -22,53 | -8,92 | -9,85 | 5,80 | 11,80 |
, где 

График сезонности
Выводы по графику:
На графике видно, что Робинзону лучше всего ходить на охоту в начале года (с 1го по 5й месяцы), а также в конце года (с середины 10го месяца по 12й).
В эти периоды ему следует запасаться на середину года (с 5го по 10е месяцы).
Шаг 3: Исключим влияние сезонной компоненты, вычитая ее значения из каждого элемента временного ряда. Получим
.
Шаг 4: Используя регрессионный анализ, проведем аналитическое выравнивание ряда X-S. Построим три модели:
1)
, 2)
, 3)
. Найдем оценки параметров
и
для каждого уравнения с помощью МНК.
СНУ:


Делаем замену 
Получим
СНУ:


Делаем замену 
Получим
СНУ:

Получим три уравнения
,
и 
Построим соответствующие графики

Исходя из анализа графика, видим, что наиболее точно зависимость количества убитой дичи от времени описывают модели
и

Проверим значимость найденных моделей
Если
, то гипотеза
отвергается с вероятностью 0,95.
Если
, то гипотеза
не отвергается с вероятностью 0,95.

Для 


R2= 0,928
F= 685,31>Fкрит =4,02, следовательно,
отвергается, что говорит о значимости нашей модели, и ее можно использовать для прогнозирования
Для 


R2= 0,952
F= 1106,08>Fкрит =4,02, следовательно,
отвергается, что говорит о значимости нашей модели, и ее можно использовать для прогнозирования.
Для 


R2= 0,4
F= 35,72>Fкрит =4,02, следовательно,
отвергается, что говорит о значимости нашей модели, и ее можно использовать для прогнозирования.
Сравнивая эти полученные значения R2, делаем вывод, что лучшая модель
, так как R2= 0,952 в данном случае наибольшее, следовательно, эта модель лучше всего описывает зависимость количества убитой дичи от времени.
5. Выводы:
После построения графика сделали следующие выводы:
· в среднем количество убитых уток увеличивается; прослеживаются сезонные колебания;
· при динамике веса отсутствуют сезонные колебания, и рост среднего значения веса не прослеживается.
· Провели первичный статистический анализ временных рядов, включая вычисление среднего значения, меры разброса. Таким образом, среднее отклонение количества добытых уток от среднего значения равно 28,6, а среднее отклонение веса 3,3. Следовательно колебания количества добытой дичи в течение времени весьма значительно, а вот колебание веса Робинзона невелико.
· Построили график сезонности на год проживания Робинзона на острове, чтобы дать рекомендации Робинзону в какой период ему лучше охотиться.
· Построили 3 модели, и в дальнейшем вывели, что модели А и В наиболее точно описывает исходные данные (количество убитой дичи).