Пример

В таблице (4.5) представлены поквартальные данные о количестве реализованных единиц товара за три года.

Таблица 5. Данные о реализации товара

Год	Квартал	Период, t	Реализация товара	Тренд	Оценка сезонной компоненты
				308,5	0,972
			310,4	1,031
			312,4	1,040
			314,3	0,939
				316,2	0,980
			318,1	1,022
			320,1	1,062
			322,0	0,947
				323,9	0,973
			325,8	1,028
			327,8	1,068
			329,7	0,940

Сначала следует оценить тренд. Для этого к данным табл. 4.5. подберем линейную модель тренда по методу наименьших квадратов.

(4.9.2)

Результаты расчетов трендовых значений для t=1,..,12 представлены в пятом столбце таблицы.

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

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

Рассчитаем средние индексы сезонности. Для первого квартала индекс сезонности равен

Для второго квартала

Для третьего квартала

Для четвертого квартала

Взаимная погашаемость сезонных воздействий в мультипликативной модели выражается в том, что сумма значений сезонной компоненты по отдельным периодам должна быть равна числу периодов в цикле. В нашем примере число периодов в цикле равно 4. Найдем сумму средних оценок сезонной компоненты

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

Найдем прогноз на четвертый год. Для первого квартала (t=13) прогноз размера реализации равен

(306,6+1,9

Для второго квартала

(306,6+1,9

Для третьего квартала

(306,6+1,9

Для четвертого квартала

(306,6+1,9

4.10. ДИАГНОСТИЧЕСКАЯ ПРОВЕРКА АДЕКВАТНОСТИ МОДЕЛЕЙ. КРИТЕРИЙ ДАРБИНА-УОТСОНА.

Формально качество модели определяется ее адекватностью и точностью прогноза.

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

Для проверки свойства нулевого среднего рассчитывается среднее значение остатков

.

Если близко к нулю, то можно считать, что модель не содержит систематической ошибки и адекватна по критерию нулевого среднего. Этот вопрос может быть решен со статистических позиций.

.

Гипотеза принимается, если

,

где

- среднее квадратичное отклонение остатков;

- табличное значение критерия Стьюдента с n-1 степенью свободы уровня значимости .

Далее осуществляется проверка остаточного ряда по критерию случайности.

Проверка независимости последовательных остатков может осуществляться также по критерию Дарбина-Уотсона:

. (4.10.1)

Несложные вычисления позволяют проверить, что статистика Дарбина- Уотсона следующим образом связана с выборочным коэффициентом корреляции между соседними наблюдениями.

(4.10.2)

В самом деле

При большом числе наблюдений n сумма значительно меньше и

,

откуда и следует приближенное равенство (4.10.2), так как в силу условия

Естественно, что в случае отсутствия автокорреляции выборочный коэффициент окажется не сильно отличающимся от нуля, а значение статистики d будет близко к двум.

Для рядов с тесной взаимосвязью между последовательными уровнями значение d близко к нулю; это свидетельствует о том, что регулярная составляющая не полностью отражена в модели тренда, т. е. модель не адекватна реальному процессу. Если последовательные остатки независимы, то d близко к 2, и это свидетельствует о хорошем качестве модели. При отрицательной автокорреляции остатков (строго периодическом чередовании их знаков) d близко к 4.

Заметим, что тест Дарбина-Уотсона, вообще говоря, не представляет собой статистический критерий в том смысле, что нельзя указать критическую область, которая позволяла бы отвергнуть гипотезу об отсутствии корреляции, если бы оказалось что в эту область попало наблюдаемое значение статистики d.

Тест работает следующим образом.

Если d , то значение d сравниваетсяс табличными значениями d_L и d_{U :}

- если , то гипотеза о независимости отвергается;

- если , то гипотеза о независимости принимается;

- если , то нет достаточных оснований для принятия решения.

Если , то с критическими значениями d_L и d_U, взятыми из таблиц, сравнивается не d, а 4- d и решение принимается по тем же правилам.

Изобразим результат Дарбина –Уотсона в виде таблицы

0<
Н0 отвергается (положительная автокорреляция)	Зона неопреде-ленности	Н0 принимается (отсутствие автокорреляции)	Зона неопреде-ленности	Н0 отвергается (отрицательная автокорреляция)

Совокупный критерий согласия Бокса-Пирса рассматривает совокупность первых k автокорреляций ряда остатков - ₁ (e), ₂ (e),…, _k (e):

Статистика Q имеет c² - распределение с k-m-1 степенями свободы. Пороговое значение для величины Q определяется по таблице.

Нормальность ряда остатков проверяется с целью использования этого свойства в дальнейшем при построении доверительных интервалов прогноза.

Ввиду малого числа наблюдений в большинстве временных рядов (меньше 30) это свойство может быть проверено лишь посредством вычисления оценок коэффициентов асимметрии и эксцесса для ряда остатков

; ,

где S – средняя квадратичная ошибка

.

Для нормального распределения и .

Гипотеза Н₀ о нормальном распределении ряда остатков (; ) принимается, если выполняются соотношения:

;

.

В этом случае доверительные интервалы для прогноза будут достаточно достоверными.

4.11. ОСНОВНЫЕ ПРОБЛЕМЫ ИДЕНТИФИКАЦИИ