Прогнозирование с помощью временных рядов

Все методы прогнозирования базируются на информации об объекте прогнозирования и его прошлом развитии. Прогноз, получающийся в результате применения методов прогнозирования, определяет ожидаемые варианты экономического развития. При этом предполагается, что основные факторы и тенденции прошлого периода сохранятся на период прогноза или что можно обосновать и учесть направление их изменений в рассматриваемой перспективе. Такую гипотезу выдвигают исходя из инерционности развития социально-экономических явлений и процессов. Инерционность проявляется во взаимосвязях, т. e. сохраняются зависимости, корреляции прогнозируемой переменной от совокупности факторных признаков, темпы и направление развития, вариация показателей на протяжении длительного периода времени. Инерционность развития экономики связана с длительно действующими факторами (структура основных фондов, их возраст и эффективность, степень устойчивости технологических взаимосвязей отраслей производства и др.).

Для того чтобы в прогнозе содержалось не только правильное качественное предсказание, но и наиболее вероятное количественное значение прогнозируемого признака, необходимо, чтобы прогностическая модель допускала малую ошибку прогноза. Ошибка прогноза будет тем меньше, чем меньше срок упреждения и чем длиннее прошлый период, на информации из которого построена прогностическая модель, т. e. чем длиннее база прогноза. Нет общих правил определения допустимого срока упреждения при заданной точности прогноза, и наоборот: нельзя указать точность прогноза в зависимости от срока упреждения. В большинстве случаев срок упреждения не должен превышать третьей части длины базы прогноза. Например, для прогноза на 5 лет желательно, чтобы база прогноза (ряд динамики) содержала не менее 15 уровней. В каждом конкретном случае соотношение длины базы прогноза и срока упреждения необходимо обосновывать, используя имеющуюся информацию.

При прогнозировании нужно взвешивать все существующие методы, чтобы воспользоваться тем из них, который наиболее полно отвечает данным обстоятельствам. Прежде всего, следует рассмотреть метод, в котором исследуемый динамический ряд экстраполируется. При этом тренд, краткосрочную осцилляцию, сезонный эффект объединяют сложением или умножением, в зависимости от обстоятельств, с тем чтобы сформировать прогноз. Затем исследуют ошибки прогноза, т. e. вычисляют стандартную ошибку оценки или доверительный интервал оценки, выражая на языке вероятностей степень уверенности в том, что оценка лежит в заданной области. Все эти действия основываются на том, что исследуемая выборка извлечена случайным образом из генеральной совокупности. При прогнозировании, осуществляя разложение (на тренд, краткосрочную осцилляцию, сезонную и случайную компоненты), строят модель. Ошибки прогноза проявляются и вследствие ошибок спецификации этой модели.

Наиболее точный способ оценивания надежности метода прогноза состоит в исследовании его “работы” за какой-либо период. По ряду рассчитанных ошибок можно сформировать хорошую эмпирическую оценку ошибки, которая, вероятно, встретится в будущем. Однако такой метод оценки надежности требует большего труда. Разумным компромиссом оказывается вычисление ошибок для прошлых моментов времени на основе текущих значений. Это может привести к недооценке истинных ошибок, но, по крайней мере, будет получено некоторое представление об ошибках в будущем.

Число единиц времени, на которое делается прогноз, называется горизонтом прогнозирования.

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

9.1. Прогнозирование с использованием показателей средних характеристик ряда динамики. Одним из наиболее распространенных методов краткосрочного прогнозирования социально-экономических явлений и процессов является экстраполяция, т. e. распространение прошлых и настоящих закономерностей, связей, соотношений на будущее. Наиболее простым методом экстраполяции одномерных рядов динамики является использование средних характеристик: среднего уровня, среднего абсолютного прироста и среднего темпа роста.

При использовании среднего уровня ряда динамики в прогнозировании социально-экономических явлений прогнозируемый уровень принимается равным среднему значению уровней ряда в прошлом:

 

,

 

Прогноз вычисляется на моментов времени вперед (период упреждения), т. e. до момента (горизонт прогнозирования). Получается прогностическая точечная оценка, которая, вообще говоря, не совпадает с фактическими данными. Поэтому для средней указывается доверительный интервал прогноза

,

где табличное значение -критерия Стьюдента с = n - 1 степенями свободы и уровнем доверия ; средняя квадратичная ошибка средней:

.

 

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

 

(9.1)

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

 

.

 

Если общая тенденция развития динамического ряда является линейной или выполняется неравенство: где - остаточная дисперсия, не объясненная экстраполяцией по среднему абсолютному приросту; - общий прирост показателя от начального уровня до конечного, то выполняется экстраполяция по среднему абсолютному приросту. Прогнозное значение уровня определяют по формуле:

 

где - уровень ряда динамики, принятый за базу экстраполяции; - средний абсолютный прирост; - период упреждения.

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

 

где - средний темп роста; - уровень ряда динамики, принятый за базу экстраполяции.

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

Отметим, что если уровни ряда динамики непрерывно возрастают за рассматриваемый период, то средний темп роста вычисляют по формуле

 

 

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

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

 

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

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

Пример 9.1. Выпуск цемента за период с 1975 по 1990 г. характеризуется динамическим рядом, представленным в табл. 9.1.

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

При экстраполяции на основе среднего уровня ряда используется принцип, при котором прогнозируемый уровень принимается равным среднему значению уровней в прошлом:

 

Таблица 9.1

Год, t	Производство цемента, млн.т..	Год, t	Производство цемента, млн.т

 

Доверительный интервал прогноза для средней вычислим по формуле (9.1):

 

.

 

Табличное значение t-статистики Стьюдента с = n - 1 = 15 степенями свободы при уровне доверия = 0,05 равно = 2,13. Среднее квадратичное отклонение, связанное с выборочной средней и варьированием уровней ряда вокруг средней, равно:

 

 

Подставив найденные значения в формулу (9.1), полущим доверительный интервал (116,1639; 143,9561), который с доверительной вероятностью 0,95 включает прогнозируемое значение производства цемента равно: млн.т.

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

 

Средний абсолютный прирост

 

 

Тогда прогнозное значение уровня на 1991г.

 

(млн.т.)

 

Экстраполяция по среднему темпу роста осуществляется по формуле

где

За базу экстраполяции примем среднее арифметическое трех последних уровней, т. e. 140,3. В этом случае прогнозируемый уровень ряда равен:

(млн.т.)

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

9.2. Прогнозирование динамики социально-экономических явлений по трендовым моделям. Прогнозирование с помощью трендов - также один из простейших и распространенных методов статистического прогнозирования. Суть этого метода заключается во временной экстраполяции. При этом предполагается, что:

- период, для которого построен тренд, достаточен для выявления тенденции;

- анализируемый процесс устойчив и обладает инерционностью;

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

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

, (9.2)

где - точечный прогноз на момент + ; - табличное значение -критерия Стьюдента с степенями свободы при уровне доверия ; - число параметров тренда; - средняя квадратичная ошибка тренда:

 

.

 

В основу расчета доверительного интервала прогноза положен показатель, определяющий колеблемость ряда заданных значений признака. Чем выше эта колеблемость тем менее определено положение тренда и тем уже должен быть интервал для вариантов прогноза при одном и том же уровне доверия. В качестве такого показателя ряда наблюдаемых значений признака обычно рассматривается среднее квадратичное отклонение фактических наблюдений от расчетных, полученных при выравнивании динамического ряда, т. e. средняя квадратичная ошибка тренда.

Доверительный интервал (9.1) учитывает неопределенность, связанную с положением тренда. Но он должен учитывать также и возможность отклонения от тренда, т. e. среднюю квадратичную ошибку прогноза . Тогда доверительный интервал прогноза будет иметь вид

 

. (9.3)

 

Стандартная ошибка прогноза, когда тренд описывается прямой , вычисляется по формуле:

 

 

и доверительный интервал (9.2) примет вид

 

(9.4)

 

где - среднее квадратичное отклонение фактических уровней динамического ряда от расчетных, называемая стандартной ошибкой тренда; К - величина, зависящая только от длины динамического ряда и периода упреждения ; табличное значение - критерия Стьюдента с = n-2 степенями свободы при уровне доверия .

С увеличением значения К уменьшаются, а с увеличением - увеличиваются. Поэтому достаточно надежный прогноз получается при относительно большом числе наблюдений (для линейного тренда n = 6, для параболического второй степени n = 13, для кубического n = 23), когда период упреждения не очень большой. При одном и том же с ростом доверительный интервал прогноза увеличивается. Кроме того, доверительный интервал прогноза при одной и той же величине средней квадратичной ошибки будет тем шире, чем выше степень полинома, характеризующего тренд.

Доверительные интервалы для линейного тренда, изображены на рис. 9.1.

Рис. 9.1

Проиллюстрируем нахождение прогноза по уравнению тренда и построение доверительного интервала на примере.

Пример 9.2. Рассмотрим динамический ряд, характеризующий производительность труда с февраля 1988г. по апрель 1989г. (см. пример 3.7). Для данного ряда наилучшей функцией, характеризующей тренд, была признана прямая , поэтому прогностическая модель имеет вид: .

Прогнозирование с помощью этой модели осуществляется весьма просто: необходимо вместо в уравнение подставить нужное значение и найти прогноз. Так, для прогнозирования производительности труда в апреле 1989г. нужно подставить = 15, вследствие чего

Если прогноз необходимо сделать в году , a период упреждения равен , то в прогностическую модель подставляется значение , где n = 14 соответствует марту 1989г. Доверительные интервалы найдем, используя выражение (9.3). Необходимые для этого значения можно определить из соответствующей таблицы (см. [19,прил. 7]). Результаты вычислений сведены в табл. 9.2.

 

 

Таблица 9.2

Месяц и год, t	Факти- ческий уровень	Выравнен- ное значе- ние и про- гноз			90% доверительный интервал
нижняя граница	верхняя граница
02.1988		22,72	2,72	-	-	-
03.1988		24,63	0,63	-	-	-
04.1988		26,54	1,46	-	-	-
05.1988		28,45	1,55	-	-	-
06.1988		30,36	0,64	-	-	-
07.1988		32,27	0,73	-	-	-
08.1988		34,18	0,18	-	-	-
09.1988		36,09	0,91	-	-	-
10.1988		38,00	0,0	-	-	-
11.1988		39,91	0,09	-	-	-
12.1988		41,82	0,82	-	-	-
01.1989		43,73	0,73	-	-	-
02.1989		45,64	0,64	-	-	-
03.1989		47,55	0,45	-	-	-
04.1989	-	49,46	1,15	2,0462	47,11	51,81
05.1989	-	51,37	1,15	2,1000	48,95	53,78
06.1989	-	53,28	1,15	2,1590	50,80	55,76

 

 

Среднее квадратичное отклонение равно 1,15. Из таблицы видим, что доверительные интервалы оказались достаточно широкими - свыше 9% прогнозируемого уровня.

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

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

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

 

- для линейного тренда

;

 

- для тренда, определяемого параболой второго порядка,

;

 

- для тренда, определяемого параболой третьего порядка,

 

 

.

 

Сопоставляя подкоренные выражения стандартных ошибок, видим, что при одном и том же значении доверительный интервал прогноза тем шире, чем выше степень полинома, характеризующего тренд. Это объясняется тем, что дисперсия уравнения тренда определяется как взвешенная сумма дисперсий соответствующих параметров уравнений. Хотя средняя квадратичная ошибка тренда является не единственной характеристикой, определяющей ширину доверительного интервала, однако она оказывает преобладающее влияние на эту величину.

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

В практике криволинейного выравнивания широко распространены два вида преобразований: логарифмирование и обратное преобразование . При этом возможно преобразование как зависимой переменной , так и независимой или одновременно той и другой.

Рассмотрим процедуру построения доверительного интервала прогноза для модифицированной экспоненты Прологарифмировав уравнение, получим линейную функцию от :

.

Пусть асимптота задана, т.е. , и не содержит ошибки. Обозначим . Тогда доверительный интервал прогноза для модифицированной экспоненты будет определяться как доверительный интервал (9.3) для прямой, т.е. ,

где среднее квадратичное отклонение от тренда . Зная границы доверительного интервала для , легко определить доверительные границы прогноза для :

. (9.5)

Так как экспоненциальная кривая логарифмированием преобразуется к виду , то доверительный интервал прогноза имеет вид

. (9.6)

 

В доверительных интервалах (9.4) и (9.5)

.

 

Пример 9.3. В табл. 9.3 дан динамический ряд, характеризующий объем продаж.

Таблица 9.3

Год	Объем продаж , млн.р.	Год	Объем продаж , млн.р.

 

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

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

Используя формулу (9.5), строим доверительные интервалы прогноза для периода упреждения: Значения коэффициента будут равны соответственно:

 

, ,

 

 

Значение квантиля , взятое из таблицы распределения Стьюдента при заданном уровне значимости и степенях свободны, равно . Вычислим для значения :

 

Тогда доверительные интервалы, для , будут иметь вид:

Прогноз для составит соответственно:

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

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

9.5. Прогнозирование методом экспоненциального сглаживания. Вначале введем понятие экспоненциальной средней. При вычислении скользящих средних – простой и взвешенной – всем уровням динамического ряда присваивались одинаковые веса. Вес отдельного наблюдения указывает на часть вклада его значения в значение средней. В случае простой скользящей средней эта часть равна для наблюдений, входящих в среднюю, и нулю для наблюдений, отсутствующих в ней. При этом недавние данные имеют тот же вес, что и данные, относящиеся к далекому прошлому (старые). Однако понятно, что недавние данные имеют более важное значение и должны иметь больший вес. Поэтому предлагается процедура усреднения с разными весами. При этом система весов образует ряд, в котором веса убывают во времени по экспоненциальному закону:

, .

Сумма этого ряда стремится к единице при неограниченном увеличении числа слагаемых.

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

 

, (9.7)

 

которую можно преобразовать к виду:

 

(9.8)

 

На основе этого уравнения строятся другие модели экспоненциального сглаживания.

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

 

(9.9)

 

Из этой формулы легко получаются выражения:

 

 

Экспоненциально взвешенная средняя имеет ряд преимуществ перед традиционной скользящей средней.

1. Для вычисления экспоненциально взвешенной средней используется предыдущая экспоненциально взвешенная средняя и последнее значение уровня .

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

3. В экспоненциально взвешенной средней значения весов убывают со временем, т.е. нет такой точки, на которой веса обрываются.

Метод экспоненциально взвешенной средней разработан для анализа динамических рядов, состоящих из большего числа наблюдений. Поэтому, если динамические ряды слишком короткие (уровней) и в случае, когда темпы роста и прироста велики, метод не “успевает” отразить все изменения. Метод тем точнее, чем больше число наблюдений (уровней динамического ряда).

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

Предположим, что динамический ряд представлен в виде

,

где - тренд; - случайная компонента. Если на изучаемом интервале времени коэффициенты уравнения, описывающего тренд, остаются неизменными, то для построения модели прогноза можно использовать метод наименьших квадратов. Однако в течение анализируемого периода коэффициенты уравнения тренда изменяются во времени. И так как динамические ряды, характеризующие экономические процессы, содержат небольшое число уровней, применение метода наименьших квадратов для оценки параметров модели прогноза может привести к существенным ошибкам. Поэтому применяется метод экспоненциально взвешенной средней, в котором новым данным придаются большие веса, чем старым.

Пусть тренд определяется линейной функцией . Как показал Р.Г. Браун, оценки коэффициентов и выражаются через экспоненциально взвешенные средние по формулам:

.

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

. (9.10)

 

Чтобы воспользоваться формулой (9.9) для прогнозирования, нужно определить значения параметров и , которые выражаются через экспоненциально взвешенные средние. А из формул (9.7) и (9.8) следует, что для вычисления и необходимо задать начальные значения и или в общем случае , которые будем называть в дальнейшем начальными условиями.

Начальные условия либо задают исходя из экономических соображений (например, из величины лага), либо вычисляют по формулам:

 

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

 

Затем вычисляются экспоненциально взвешенные средние первого и второго порядков:

 

, , …

 

, , …

 

Ошибка прогноза при использовании доверительного интервала (9.3) определяется по формуле

,

 

где – средняя квадратичная ошибка, характеризующая отклонения от линейного тренда: .

 

При использовании прогностической модели (9.10) одной из основных проблем является выбор оптимального значения параметра сглаживания , где . От численного значения зависит, насколько быстро будет уменьшаться вес предшествующих наблюдений, т.е. насколько быстро будет уменьшаться степень их влияния на сглаженный уровень. Это значит, что чувствительность экспоненциально взвешенной средней в целях повышения адекватности прогностической модели может быть в любой момент изменена путем изменения значений . Чем больше , тем выше чувствительность средней. Чем меньше значение , тем устойчивее становится экспоненциально взвешенная средняя. Если подходящими оказываются более высокие значения , это указывает на нарушение условий стационарности и означает, что экспоненциально взвешенная средняя становится неприемлемой для прогнозирования. Значения при условии равенства среднего значения степени старения данных можно выбирать, используя формулу или . Значения , используемые в области экономического прогнозирования, находятся в пределах от 0,05 до 0,3. Длина усреднения в скользящем среднем с точки зрения чувствительности прогноза может быть найдена в соответствии с из таблицы

	0,05	0,1	0,2	0,3

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

Пример 9.4. Рассмотрим динамический ряд, характеризующий производство цемента (таблица 9.4).

.

Таблица 9.4

Год	Производство цемента , млн. т.	Год	Производство цемента , млн. т.	Год	Производство цемента , млн. т
				2081

 

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

 

 

Решая ее, находим: . Уравнение тренда имеет вид:

.

Для прогноза выпуска цемента на 1991 г. воспользуемся формулой (9.10). Оценки коэффициентов и найдем из выражений:

которые содержат экспоненциально взвешенные средние и и параметр . Параметр сглаживания положим равным 0,15, так как для рекомендуется брать в нашем примере . Вычисление и осуществим по рекуррентной формуле (9.9), предварительно определив начальные условия и

,

,

где, и - коэффициенты уравнения тренда. Тогда:

,

 

Затем вычисляем и

,

и осуществляем прогноз на 1976 г. Далее по рекуррентной формуле (9.9) вычисляем новые и :

, ,

и по ним находим и которые используем для прогноза производства цемента на 1977 г., и т.д. В таблице 9.5 приведены:

- экспоненциально взвешенные средние, вычисленные по формуле (9.9);

-соответствующие коэффициенты и ;

- результаты прогноза и отклонения фактических уровней от прогнозируемых в случае ретроспективного прогноза;

- указан прогноз производства цемента на 1991 г.

Таблица 9.5

Год	Производст-во цемента, млн. т					Прогноз	Отклонение
		125,82	123,36	128,28	0,43	128,71	-4,71
		125,60	123,70	127,50	0,34	127,84	-0,84
		125,81	124,02	127,60	0,32	127,95	-0,92
		125,39	124,23	126,55	0,20	126,75	-3,75
		125,33	124,39	126,27	0,17	126,44	-1,44
		125,58	124,57	126,59	0,18	126,77	0,23
		125,39	124,69	125,99	0,11	126,10	-2,10
		125,74	124,85	126,63	0,16	126,69	1,31
		126,38	125,08	127,68	0,23	127,91	2,09
		127,07	125,38	128,76	0,30	129,06	1,94
		128,26	125,81	130,71	0,43	131,14	3,86
		129,57	126,37	132,77	0,56	133,33	3,67
		130,98	127,06	134,90	0,69	135,59	3,41
		132,33	127,85	136,81	0,79	137,60	2,40
		133,78	128,74	138,82	0,89	139,71	2,99
	-	135,01	129,68	140,34	0,94	141,28	-

 

Для прогноза производства цемента на 1991 г. использовались следующие значения экспоненциально взвешенных средних: , и оценки коэффициентов модели , . Ошибку прогноза вычислим по формуле

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

 

Тогда доверительный интервал прогноза определяется в виде (9.3), где - квантиль распределения Стьюдента при уровне доверия и числе степеней свободы , равный . Подставив значения и , получим доверительный интервал прогноза

.

 

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

Рассмотрим временной ряд , математическая модель которого имеет вид ,

где - неслучайная функция времени (тренд); - стационарная случайная компонента.

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

Для применения метода гармонических весов ряд динамики разбивается на интервалы, каждый из которых содержит 3 – 5 уровней. Число интервалов меньше . Для каждого интервала определяется линейный тренд

Причем для для

для .

Оценивание параметров скользящего тренда осуществляется методом наименьших квадратов. Вычислив оценки параметров и , получим уравнение. Вычислим далее значение в точках , где . Для каждого уравнения получим число значений функции , равное числу уровней, содержащихся в интервале скольжения. Образуем множества , значений функции - , для которых . Эти функции обозначим , а число таких функций - . Вычислим средние функций, содержащихся в построенных множествах:

(9.11)

 

Соединив точки отрезками прямой, получим тренд исследуемого динамического ряда в виде ломаной линии.

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

 

(9.12)

где Величину называют сдвигом. Сдвиг, которому соответствует наибольший коэффициент автокорреляции, называют временным лагом. График нормированной автокорреляционной функции называют коррелограммой.

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

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

Значения автокорреляционной функции образуют ряд , , …, (верхний индекс означает число наблюдений, для которого вычисляется автокорреляционная функция). Затем исключают из динамического ряда первый или последний уровень и вычисляют значения автокорреляционной функции , , …, . Продолжая указанный процесс, исключают уровней динамического ряда и вычисляют значения автокорреляционных функций. Таким образом, получают групп коэффициентов автокорреляции, в каждой из которых будет коэффициентов. Отклонения от скользящего тренда образуют стационарный в широком смысле процесс, если коэффициенты автокорреляции, входящие в одну и ту же группу, однородны.

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

, .

Затем для этой группы находят среднюю и вычисляют величину

, (9.13)

которая распределена по закону - квадрат с степенями свободы. Тогда, сравнивая вычисленное значение величины (9.13) с табличным, при с вероятностью принимаем гипотезу об однородности рассматриваемой группы коэффициентов автокорреляции. Аналогичную проверку однородности проводим для всех групп коэффициентов автокорреляции. Если гипотеза об однородности принимается для всех групп, то делаем вывод о том, что отклонения от скользящего тренда образуют стационарный в широком смысле случайный процесс. Кроме того, если значения автокорреляционной функции, вычисленные для ряда отклонений от скользящего тренда, уменьшаются, это значит, что более поздняя информация сильнее отражается на прогнозируемой величине, чем более ранняя.

Установив, что отклонения образуют стационарный процесс, вычисляем приросты, равные разностям средних значений скользящих трендов:

.

 

Средняя приростов вычисляется по формуле

 

,

где гармонические коэффициенты, удовлетворяющие следующим условиям:

, . (9.14)

Гармонические коэффициенты определяем так, чтобы более поздним наблюдениям придавались большие веса. Для этого полагаем

,

где , т.е.

. Следовательно, и гармонические коэффициенты удовлетворяют условиям (9.14).

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

 

. Применив неравенство Чебышева, можно записать: ,

где заданное положительное число; . Так как значения коррелированны между собой, то в неравенстве Чебышева является величиной переменной, вычисляемой по формуле .

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

.

 

Пример 9.5. Рассмотрим динамический ряд, характеризующий производство цемента (см. таблицу 9.5). Этот динамический ряд не имеет скачкообразных изменений и достаточно хорошо описывается линейным трендом. Поэтому для дисконтирования уровней ряда динамики с целью определения прогноза производства цемента в 1991 г. применим метод гармонических весов. Построим вначале скользящий тренд. Для этого разобьем исходный ряд динамики на интервалы, каждый из которых содержит 5 уровней. Для каждого интервала скольжения строим методом наименьших квадратов линейный тренд . Так, первый интервал скольжения состоит из уровней 122, 124, 127, 127, 123, то система нормальных уравнений для оценки параметров, имеет вид

откуда находим , . Следовательно, линейный тренд для первого интервала скольжения выражается уравнением

Аналогично определяем параметры уравнений для всех интервалов скольжения:

 

 

 

 

 

С помощью построенных уравнений определим значения скользящего тренда по формуле (9.11).

При имеем

 

При имеем два значения функций, для которых:

, ,

 

откуда

 

Аналогично находим все остальные значен