Введение
В данной главе рассматриваются задачи описания упорядоченных данных, полученных последовательно (во времени). Вообще говоря, упорядоченность может иметь место не только во времени, но и в пространстве, например, диаметр нити как функция её длины (одномерный случай), значение температуры воздуха как функция пространственных координат (трёхмерный случай).
В отличие от регрессионного анализа, где порядок строк в матрице наблюдений может быть произвольным, во временных рядах важна упорядоченность, а следовательно, интерес представляет взаимосвязь значений, относящихся к разным моментам времени.
Если значения ряда известны в отдельные моменты времени, то такой ряд называют дискретным , в отличие от непрерывного , значения которого известны в любой момент времени. Интервал между двумя последовательными моментами времени назовём тактом (шагом) . Здесь будут рассматриваться в основном дискретные временные ряды с фиксированной протяжённостью такта, принимаемой за единицу счёта. Заметим, что временные ряды экономических показателей, как правило, дискретны.
Значения ряда могут быть измеряемыми непосредственно (цена, доходность, температура), либо агрегированными (кумулятивными) , например, объём выпуска; расстояние, пройдённое грузоперевозчиками за временной такт.
Если значения ряда определяются детерминированной математической функцией, то ряд называют детерминированным . Если эти значения могут быть описаны лишь с привлечением вероятностных моделей, то временной ряд называют случайным .
Явление, протекающее во времени, называют процессом , поэтому можно говорить о детерминированном или случайном процессах. В последнем случае используют часто термин “стохастический процесс” . Анализируемый отрезок временного ряда может рассматриваться как частная реализация (выборка) изучаемого стохастического процесса, генерируемого скрытым вероятностным механизмом.
Временные ряды возникают во многих предметных областях и имеют различную природу. Для их изучения предложены различные методы, что делает теорию временных рядов весьма разветвленной дисциплиной. Так, в зависимости от вида временных рядов можно выделить такие разделы теории анализа временных рядов:
– стационарные случайные процессы, описывающие последовательности случайных величин, вероятностные свойства которых не изменяются во времени. Подобные процессы широко распространены в радиотехнике, метереологии, сейсмологии и т. д.
– диффузионные процессы, имеющие место при взаимопроникновении жидкостей и газов.
– точечные процессы, описывающие последовательности событий, таких как поступление заявок на обслуживание, стихийных и техногенных катастроф. Подобные процессы изучаются в теории массового обслуживания.
Мы ограничимся рассмотрением прикладных аспектов анализа временных рядов, которые полезны при решении практических задач в экономике, финансах. Основной упор будет сделан на методы подбора математической модели для описания временного ряда и прогнозирования его поведения.
1.Цели, методы и этапы анализа временных рядов
Практическое изучение временного ряда предполагает выявление свойств ряда и получение выводов о вероятностном механизме, порождающем этот ряд. Основные цели при изучении временного ряда следующие:
– описание характерных особенностей ряда в сжатой форме;
– построение модели временного ряда;
– предсказание будущих значений на основе прошлых наблюдений;
– управление процессом, порождающим временной ряд, путем выборки сигналов, предупреждающих о грядущих неблагоприятных событиях.
Достижение поставленных целей возможно далеко не всегда как из-за недостатка исходных данных (недостаточная длительность наблюдения), так из-за изменчивости со временем статистической структуры ряда.
Перечисленные цели диктуют в значительной мере, последовательность этапов анализа временных рядов:
1) графическое представление и описание поведения ряда;
2) выделение и исключение закономерных, неслучайных составляющих ряда, зависящих от времени;
3) исследование случайной составляющей временного ряда, оставшейся после удаления закономерной составляющей;
4) построение (подбор) математической модели для описания случайной составляющей и проверка ее адекватности;
5) прогнозирование будущих значений ряда.
При анализе временных рядов используются различные методы, наиболее распространенными из которых являются:
1) корреляционный анализ, используемый для выявления характерных особенностей ряда (периодичностей, тенденций и т. д.);
2) спектральный анализ, позволяющий находить периодические составляющие временного ряда;
3) методы сглаживания и фильтрации, предназначенные для преобразования временных рядов с целью удаления высокочастотных и сезонных колебаний;
5) методы прогнозирования.
2.Структурные компоненты временного ряда
Как уже отмечалось, в модели временного ряда принято выделять две основные составляющие: детерминированную и случайную (рис.). Под детерминированной составляющей временного ряда
понимают числовую последовательность , элементы которой вычисляются по определенному правилу как функция времени t . Исключив детерминированную составляющую из данных, мы получим колеблющийся вокруг нуля ряд, который может в одном предельном случае представлять чисто случайные скачки, а в другом – плавное колебательное движение. В большинстве случаев будет нечто среднее: некоторая иррегулярность и определенный систематический эффект, обусловленный зависимостью последовательных членов ряда.В свою очередь, детерминированная составляющая может содержать следующие структурные компоненты:
1) тренд g, представляющий собой плавное изменение процесса во времени и обусловленный действием долговременных факторов. В качестве примера таких факторов в экономике можно назвать: а) изменение демографических характеристик популяции (численности, возрастной структуры); б) технологическое и экономическое развитие; в) рост потребления.
2) сезонный эффект s , связанный с наличием факторов, действующих циклически с заранее известной периодичностью. Ряд в этом случае имеет иерархическую шкалу времени (например, внутри года есть сезоны, связанные с временами года, кварталы, месяцы) и в одноименных точках ряда имеют место сходные эффекты.
Рис. Структурные компоненты временного ряда.
Типичные примеры сезонного эффекта: изменение загруженности автотрассы в течение суток, по дням недели, временам года, пик продаж товаров для школьников в конце августа - начале сентября. Сезонная компонента со временем может меняться, либо носить плавающий характер. Так на графике объема перевозок авиалайнерами (см рис.) видно, что локальные пики, приходящиеся на праздник Пасхи «плавают» из-за изменчивости ее сроков.
Циклическая компонента c , описывающая длительные периоды относительного подъема и спада и состоящая из циклов переменной длительности и амплитуды. Подобная компонента весьма характерна для рядов макроэкономических показателей. Циклические изменения обусловлены здесь взаимодействием спроса и предложения, а также наложением таких факторов, как истощение ресурсов, погодные условия, изменения в налоговой политике и т. п. Отметим, что циклическую компоненту крайне трудно идентифицировать формальными методами, исходя только из данных изучаемого ряда.
«Взрывная» компонента i , иначе интервенция, под которой понимают существенное кратковременное воздействие на временной ряд. Примером интервенции могут служить события «черного вторника» 1994г., когда курс доллара за день вырос на несколько десятков процентов.
Случайная составляющая ряда отражает воздействие многочисленных факторов случайного характера и может иметь разнообразную структуру, начиная от простейшей в виде «белого шума» до весьма сложных, описываемых моделями авторегрессии-скользящего среднего (подробнее дальше).
После выделения структурных компонент необходимо специфицировать форму их вхождения во временной ряд. На верхнем уровне представления с выделением лишь детерминированной и случайной составляющих обычно используют аддитивную либо мультипликативную модели.
Аддитивная модель имеет вид
;мультипликативная –
Эконометрические модели начали использоваться для экономического прогнозирования в 60-е годы ХХ в. С этого времени структура экономики развитых стран и методы эконометрического анализа претерпели кардинальные изменения. В то же время проблема прогнозирования будущего состояния экономики остается нерешенной, что требует усовершенствования эконометрических моделей.
Специалисты сосредоточены на исследованиях, связанных с коинтеграцией (метод определения долговременной взаимосвязи в группе переменных динамичных рядов); на прогнозировании и оценке параметров, меняющихся во времени. В частности, разработка американским экономистом Р. Инглом проблемы коинтеграции меняет подход экономистов-практиков к изучению временных рядов.
Временные ряды - последовательность наблюдений за экономическими изменениями за одинаковые временные интервалы.
Анализ временных рядов - основной инструмент экономической науки и одна из самых плодотворных сфер анализа для экономистов. Временные ряды необходимы для анализа эволюции во времени экономических и социальных связей между переменными (например, эконометрическая модель поведения совокупной безработицы, которая базируется на временных рядах, может дать ценную информацию об ее эволюции во времени, хотя не дает сведений о структуре или продолжительности безработицы). Большая часть использующихся данных имеет вид временных рядов, массив которых постоянно расширяется.
Одним из известнейших исследователей в этой области является К. Грэнджер.
Грэнджер (Granger) Клив (также Клайв) (род в 1934) - американский экономист, лауреат Нобелевской премии (2003). Родился в г. Суонс (Уэльс, Великобритания). Учился в Ноттингемском университете, где в 1955 г. защитил бакалаврскую работу по математике, а в 1959 г. - докторскую диссертацию по статистике. Работал профессором Калифорнийского университета (г. Сан-Диего).
Он автор больше десяти книг, свыше двухсот научных статей.
К. Грэнджер - член Британской национальной академии наук, Американского эконометрического общества, Американской и Финской академий искусств и наук; заслуженный член Американской экономической ассоциации, почетный доктор Ноттингемского, Мадридского, Лафборского университетов и Стокгольмской школы экономики, заслуженный профессор Калифорнийского университета.
Ингл (Engle) Роберт (род. в 1942) - американский экономист, лауреат Нобелевской премии (2003). Родился в г. Сиракьюс (штат Нью-Йорк, США). Учился в Корнелльском университете. В 1969 г. защитил докторскую диссертацию по экономике. В течение 1969-1974 гг. работал ассистентом профессора Массачусетского технологического института; в 1975 г. - адъюнкт-профессор Калифорнийского университета г. Сан-Диего. Через два года занял должность профессора. На протяжении 1990-1994 гг. был деканом экономического факультета этого же университета, позже - профессор менеджмента финансового факультета Нью-Йоркского университета.
Р. Ингл - известный эксперт по анализу временных рядов в течение долгосрочных периодов на финансовых рынках. Его исследования посвящены таким инновационным статистическим методам, как ARCH-моделирование, коинтеграция, взаимосвязанные спектральные регрессии. В своих исследованиях использует методы финансовой эконометрии для проведения операций с акциями, валютными и процентными ставками, опционами.
Он член Американского эконометрического общества и Американской академии искусств и наук.
Разработка анализа временных рядов (и на его основе - прогнозирование и контроль) основала новое направление в методах прогнозирования, стала теоретической основой ARIMA-анализа, по которому определенный временной ряд моделируют лишь с помощью его прошлых значений и экзогенной случайной величины, и методологии, необходимым условием которой является стационарность рассматриваемого временного ряда. Такая методология является сравнительно новым поколением средств прогнозирования, основанных на анализе вероятностных (стохастических) особенностей временных рядов. При этом определенный временной ряд моделируется лишь с помощью его прошлых значений (лагов) и экзогенной случайной величины. Необходимым условием внедрения ARIMA-методологии является стационарность временного ряда - математического ожидания (среднее), дисперсия и автоковариация (в разных промежутках) которого не зависят от времени. Если он стационарный, то его можно смоделировать разными способами, в частности с помощью двух составляющих - авторегрессийной (AR) и скользящего среднего (MA). Соответственно сама модель является комбинацией этих двух составляющих.
Поскольку ARIMA-методология используется только для стационарных рядов, то первым шагом в идентификации процесса становится проверка временного ряда на стационарность. Необходимость того, чтобы временные ряды были стационарными при ARIMA-моделировании, обусловлена тем, что эти модели используются для прогноза, а прогнозировать можно поведение только тех процессов, основные характеристики которых (средняя, дисперсия и коэффициенты автоковариации) не зависят от времени. Невозможно предусмотреть поведение того процесса, в основе которого нестационарный временной ряд (математическое ожидание, дисперсия и автоковариация его меняются в зависимости от времени). В таком случае сложно найти постоянные средней и дисперсии, поэтому следует искать возможные преобразования ряда, которые могут свести его к стационарному. Такими преобразованиями и является операция различий.
Моделирование экономических процессов с помощью ARIMA-моделей дает возможность выявить динамическую связь между поточными и лаговыми значениями исследуемого показателя. Эти модели являются удобным инструментом кратко- и среднесрочного прогнозирования отдельных временных рядов. Однако современные исследования сосредоточены на разработке аппарата одновременного моделирования нескольких временных рядов с помощью системы динамических уравнений ARIMA-процессов, что дает возможность включать и исследовать взаимообратные связи между показателями и их лаговыми значениями.
Таким образом, VAR-модели (векторная авторегрессионная модель) является расширением концепции ARIMA-моделирования отдельного временного ряда. Термин «вектор» в этом случае указывает, что моделируются одновременно два или более временных ряда. Термин «авторегрессионная» означает включение лаговых значений зависимых переменных в правую часть каждого отдельного уравнения системы. Стабильность VAR-моделей является необходимым условием их практического использования. Она предусматривает, что последовательность внешних шоков для VAR-системы имеет конечный падающий эффект, то есть если шоки затухают со временем, то VAR-модель является стационарной.
В 90-е годы ХХ в. активно развивается новое направление моделирования с помощью моделей корректирования ошибки (error correction model - ECM). Эти модели являются структурной формой VAR-моделей, которая включает нестационарные переменные. Для оценки таких систем необходимы дополнительные знания, в частности коинтеграции временных рядов. Коинтеграция переменных дает возможность строить корректные модели даже в случае их нестационарности, не преобразуя временные ряды оператором различий в стационарные. Это важно для прикладных исследований, так как, используя оператор различий, утрачивается ценная «долгосрочная» информация о динамике поведения временного ряда. Поэтому преобразовывать ряды целесообразно только при необходимости.
Построение и корректное внедрение ЕСМ предусматривает определенную последовательность.
1. Проверка рядов на стационарность. Если они не стационарны, то необходимо определить порядок интеграции. При одинаковом порядке интеграции можно переходить к проверке рядов на коинтеграцию.
И только тогда, когда ряды коинтегрируют, можно строить ЕСМ (она является не чем иным, как VAR в структурной форме), и оценивать ее неизвестные параметры.
Именно Р. Ингл и К. Грэнджер предложили собственное понимание коинтеграции: если между рассматриваемыми переменными существует долгосрочная связь, то очевидно долгосрочное равновесие достигается, когда:
γ1Υ1t + γ2Υ2t + ... + γkΥkt = 0,
или в матричном виде:
γΥt = 0, где γt = {γ1, γ2, ..., γk}, Υt = {Υ1t, Υ2t, ..., Υkt}.
Отклонения от долгосрочного равновесия называют «ошибкой равновесия», что, соответственно, равняется et = γΥt.
Если равновесие есть, то необходимо, чтобы ошибка равновесия была стационарным процессом.
Исходя из приведенных формул, Р. Ингл и К. Грэнджер утверждают: компоненты вектора Υt = {Υ1t, Υ2t, ..., Υkt} являются коинтегрированными порядка d,b: ~ CI (d,b), если:
- все компоненты Υt имеют одинаковый порядок интеграции d;
- существует вектор коэффициентов γt = {γ1, γ2, ..., γk} такой, что линейная комбинация γ1Υ1t + γ2Υ2t + ... + γkΥkt является интегрированной величиной порядка (d - b), b > 0.
Вектор γ = {γ1, γ2, ..., γk} называют «коинтеграционным вектором». Очевидно, что если γt = {γ1, γ2, ..., γk} является коинтеграционным вектором, то для любого дополнительного Ø значения Øγ = {Øγ1, Øγ2, ..., Øγk} является также коинтеграционным вектором. Поэтому на практике одна из переменных используется для нормализации коинтеграционного вектора, то есть соответствующий коэффициент должен быть равен единице.
2. Проверка временных рядов на коинтеграцию. Есть несколько принципиально разных подходов к тестированию временных рядов на коинтеграцию. Р. Ингл и К. Грэнджер предложили сначала оценить уравнение долгосрочного равновесия, рассчитать избытки, то есть получить соответствующий временной ряд избытков, а потом, если излишки окажутся стационарным рядом, можно сделать вывод о коинтеграции.
На рубеже ХХ и ХХI в. в рамках теории валютного курса основные усилия ученых были направлены на изучение долгосрочных взаимосвязей относительных цен и валютных курсов. Они использовали современные методы эконометрического анализа: определение порядка интегрированности временных рядов и их тестирование на коинтеграцию.
Если два временных ряда коинтегрированы, то это означает, что их индивидуальные тренды взаимозависимы и не могут значительно отклоняться один от другого. Согласно разработке Р. Ингла и К. Грэнджера для коинтегрированных переменных существует ЕС-распределение (error correction representation). Этот механизм улавливает краткосрочное корригирование переменных относительно долгосрочного равновесия. То есть если номинальные валютные курсы и относительные цены коинтегрированы, то паритет покупательской способности является условием подобного равновесия валютных курсов, и они в своей долгосрочной динамике приближаются к нему.
В большинстве исследований, проводимых с использованием этой методологии, было установлено коинтеграцию валютных курсов и относительных цен. Однако в отношении разных групп стран эти результаты не были однозначными. Так, американские ученые К. Хабермайер и М. Месквита нашли подтверждение теории паритета покупательской способности для развитых стран, но не смогли доказать возможность ее использования развивающимися странами.
Значимыми являются также исследования Р. Ингла и К. Грэнджера и в сфере эконометрического прогнозирования. Известны разнообразные конкурирующие прогнозы с разными информационными множественностями и разными стратегиями моделирования. Эти прогнозы можно сравнивать с «прогностической способностью», то есть сопоставлять суммы квадратов ошибок прогнозирования. С помощью комбинирования разнообразных прогнозов также получают хорошие прогнозы. Такое комбинирование можно осуществить, рассчитывая регрессии фактических значений ряда от разнообразных прогнозов, константы и лаговых значений того же ряда. Прогноз, который не «вписывается» в такую регрессию, можно отбросить, поскольку над ним доминируют другие прогнозы.
Точечное прогнозирование имеет небольшую ценность для принятия решений без каких-либо указаний на неопределенность. Относительно большинства традиционных экономических прогнозов 95 %-ные интервалы вокруг прогнозированной точки необычайно велики, поэтому иногда рекомендуются 50 %-е интервалы. Еще одна проблема заключается в том, что дисперсии ошибок прогнозирования могут меняться во времени. Как и условная средняя fn,h, условная дисперсия может быть функцией использованной информационной множественности In:
h2n = E[(xn+h - fn,h)2|In].
Методы моделирования h2n менее разработаны, чем методы моделирования fn,h ошибки прогнозирования еп,1 = xn+1 - fn,1 часто являются белым шумом, но квадратичная ошибка может выявиться не такой, которая указывает на то, что условные дисперсии могут быть прогнозированными.
Обозначив через εt = xt - ft-1 одношаговые ошибки прогнозирования, Р. Ингл рассмотрел спецификацию:
и выявленный процесс назвал «авторегрессионным условным гетероскедастичным процессом» (такой, который предусматривает переменный разброс). Если дисперсия меняется во времени прогнозируемо, то преимущество ее моделирования заключается в том, что при учете гетероскедастичности удается достичь более точных оценок параметров в ft, а также получить более точные оценки интервалов вокруг прогноза средней.
Р. Ингл рассмотрел разные формы для ht, сделав вывод об их особенностях и методике оценки, а также использовав метод множителей Лагранжа для проверки авторегрессивной условной гетероскедастичности (построенные на основе этого подхода модели названы «ARCH-моделями»). Он использовал этот метод для анализа данных об инфляции в Великобритании и выявил четкие признаки прогнозируемости дисперсий: стандартное отклонение инфляции выросло за несколько лет с 0,6 до 1,5 % в меру движения экономики из предусматриваемых 60-х в хаотические 70-е годы.
Приведенное выше выражение для ht можно использовать для включения наблюдаемых управляющих переменных. Как пример К. Грэнджер исследовал связь между розничными и оптовыми ценами, причем в каждом уравнении дисперсии были специфицированы вышеописанным порядком, но с добавлением квадратичных лаговых значений моделированных и других цен, а также квадратичных ошибок прогнозирования других показателей. Обогащение спецификации ARCH обусловило появление лучших (по коэффициентам правдоподобности) моделей, а также более интересных интерпретаций моделей. Было выявлено, что и средние значения, и дисперсии оптовых цен влияют соответственно на средние значения и дисперсии потребительских цен. А квадраты потребительских цен не влияют на дисперсию оптовых цен. Если бы эти модели были построены без учета ARCH, то создавалась бы видимость влияния потребительских цен на оптовые цены. Однако с учетом ARCH эта причинная связь стала слабой.
Поскольку на практике дисперсии меняются во времени предсказуемо, то использование моделей ARCH можно рекомендовать для случаев, когда доверительным интервалам прогноза уделяется значительное внимание. Другие сферы анализа сосредотачиваются на тех отраслях экономической теории, где дисперсию используют как показатель риска (например, финансовая теория).
В последнее время инструментарий анализа временных рядов стремительно развивался. Но если брать для проверки на коинтеграцию две переменных, то лучше и далее использовать тест Ингла-Грэнджера (если проверять больше двух, то можно использовать технику Йохансена).
Исследование методов анализа экономических временных рядов в условиях изменчивости временной зависимости (ARCH) Р. Ингл и К. Грэнджер проводили на основе математической модели, которая дает возможность прогнозировать тенденции изменений ВВП, потребительских цен, процентных ставок, биржевого курса не только на следующий день, а даже на год вперед. Дело в том, что на финансовых рынках случайные отклонения показателей от постоянного значения (волатильность) являются необычайно важными, поскольку стоимость акций, опционов и других финансовых инструментов зависит от рисков. Отклонения могут значительно изменяться во времени: после периодов значительных перемен наступают периоды незначительных. Помимо того, что реальная волатильность изменчива, экономисты долгое время внедряли статистические методы, которые предусматривают ее постоянность.
И только выявленная в 1982 г. Р. Инглом авторегрессивная гетероскедастическая модель точно описывает множество временных рядов, которые встречаются в экономике.
Результаты исследования волатильности широко используют на практике, в частности:
а) с 1996 г. международные соглашения (так называемые Базельские правила) обязывают использовать показатели стоимости, поддающейся риску, при контроле необходимого капитала банков. Использование метода ARCH в этих и других ситуациях сделало его необходимым инструментом для оценки риска в финансовой сфере;
б) ими воспользовались эксперты для введения евро. Так, проект экономического и валютного союза, касающийся интересов ряда государств, был детально проанализирован академическими экономистами США и Великобритании.
Их интересовали вопросы, вырастут или уменьшатся флуктуации (случайные отклонения величины) параметров системы, то есть обменного курса, вследствие введения евро, вырастут или уменьшатся при переходе к единой валюте флуктуации платежного баланса, чего можно ожидать от курса доллар США/евро.
С помощью волатильности обменного курса было доказано, что флуктуаций станет меньше. Между странами-участницами Еврозоны они исчезнут вообще. А поскольку зона евро рассматривается как неизменное во времени творение, то будут равны нулю все форвардные премии и исчезнет разница в процентных ставках; останутся только ножницы в налоговых ставках и рисках дефолта. Содружество государств валютного союза станет великой зоной валютной стабильности.
Специалисты тоже пришли к выводу, что колебания платежных балансов при общей валюте станут меньшими, чем те, которые наблюдались при плавающих курсах. Исчезнут два источника нестабильности:
1) не будет колебаться обменный курс, движение которого стимулируют потоки капитала (спекулятивные потоки капитала исчезнут или существенно ослабнут);
2) в монетарной политике профициты платежных балансов, которые будут меньше или больше от желаемого уровня, автоматически будут корректироваться механизмом перелива резервов.
От платежных балансов внутри стран зоны евро не откажутся, но их корректировка будет программироваться ранее и окажется внешне не наблюдаемой за исключением экстраординарных случаев.
В отношении курса доллар США/евро отмечается, что он станет важнейшим ценовым фактором в мире. Некоторые считали, что этот курс должен иметь большие колебания, чем курс доллар США/немецкая марка, поскольку экономика Евросоюза более замкнута, чем объединенные в союз национальные экономики. Однако специалисты отклонили такое мнение. Если ориентироваться не на отношения импорта или экспорта к ВВП, а на общий баланс платежей, и прежде всего на движение капитала, то с устранением спекулятивных мотивов в зоне евро исчезнут и дестабилизирующие сдвиги от «более слабых» валют к «более сильным».
В целом сегодня уже невозможно изучать ключевые моменты в стабильности мировой денежной системы, не используя волатильности обменного курса. Кроме того, модель Ингла является незаменимой не только для ученых, но и для финансовых и рыночных аналитиков, которые используют ее при оценке собственности и рисков портфельных инвестиций.
Специалисты считают, что во многих аспектах экономические преобразования 90-х годов подобны преобразованиям первого десятилетия ХХ века. Эффект от осторожной финансовой политики одинаковый.
И все же, по мнению Р.-А. Манделла, мироустройство изменилось в худшую сторону: из-за постоянной изменчивости (волатильность) обменных курсов при отсутствии мировой валюты. От волатильности обменных курсов особенно страдают страны, которые стремятся поодиночке путем введения собственных масштабов и индексов достичь стабильности цен. Поэтому волатильность является мерилом тех изменений, которые претерпевают реальные обменные курсы, и отражает дисфункциональные перекосы внутреннего и международного развития отраслей, что еще больше усиливает свойственную финансовым рынкам нестабильность.
Последние разработки в области анализа нестационарных временных рядов уже влияют на методы прогнозирования. Р. Ингл и К. Грэнджер рассматривают свойства двух и большего количества объединенных переменных, каждая из которых является интегрированной первого порядка, в то время как их комбинация является стационарной (то есть интегрированной нулевого порядка). Такие переменные называются «коинтегрированными».
Коинтеграция играет важную роль в экономическом моделировании и прогнозировании. Во-первых, если переменные уравнения не коинтегрированы, то, поскольку ошибки не стационарны, связь между переменными может быть неправильно специфицирована (или получить в значительной мере достоверную оценку параметров будет сложно). Во-вторых, Р. Ингл и К. Грэнджер доказали, что если х и у являются интегрированными первого порядка, имеют постоянные средние и коинтегрированы, то существует механизм, который корректирует ошибки генерирования данных (модель корректировки ошибок), выражаемый аналитически следующим образом:
Δyt = -α1ut-1 + лаговые значения (Δy, Δx) + d(L)ε1t,
Δxt = - α2ut-1 + лаговые значения (Δy, Δx) + d(L)ε2t, (6.1)
где ut = yt - βxt, (6.2)
а Δ - оператор первых разниц. Здесь d(L) является конечным полиномом лагового оператора L, а εi - случайный процесс, причем
│α1│+│α2│ ≠ 0. (6.3)
Интерпретация (6.1) облегчается рассмотрением равновесной ситуации, при которой разницы в формуле (6.1) нулевые, и выражение (6.1) преобразовывается в (6.2) при иt = 0, то есть в равновесии в пропорциональный х. Отсюда, согласно выражению (6.2), и - это отклонение от равновесного значения, и поскольку и является стационарным с нулевой средней, то отклонение от равновесия в период t - 1 частично корректируется в период t. Значит, механизм корректировки ошибок в экономической интерпретации обеспечивает связь между структурными моделями и моделями временных рядов. Такой механизм корректировки ошибок является важнейшим для прогнозирования, поскольку он означает, что модель, включающая только различия переменных первого порядка, будет неправильно специфицирована по коинтегрированным переменным. Это может произойти, если, например, VAR-модель используется для аппроксимации данных, имеющих вид различий первого периода.
Ценность новаторских идей Р. Ингла и К. Грэнджера заключается не только в том, что они предложили новые методы моделирования экономических зависимостей, но и в том, что разработанные ими модели открыли новые сферы исследований. При этом нобелианты фундаментально обосновали использование таких моделей, доказали корректность эконометрической оценки их параметров в случае нарушения ряда классических прогнозов. Важно и то, что каждый из предложенных методов подтвердил теоретические результаты.
1 Виды и методы анализа временных рядов
Временным рядом называется ряд наблюдений за значениями некоторого показателя (признака), упорядоченный в хронологической последовательности, т.е. в порядке возрастания переменной t- временного параметра. Отдельные наблюдения временного ряда называются уровнями этого ряда.
1.1 Виды временных рядов
Временные ряды делятся на моментные и интервальные. В моментных временных рядах уровни характеризуют значения показателя по состоянию на определенные моменты времени. Например, моментными являются временные ряды цен на определенные виды товаров, временные ряды курсов акций, уровни которых фиксируются для конкретных чисел. Примерами моментных временных рядов могут служить также ряды численности населения или стоимости основных фондов, т.к. значения уровней этих рядов определяются ежегодно на одно и то же число.
В интервальных рядах уровни характеризуют значение показателя за определенные интервалы (периоды) времени. Примерами рядов этого типа могут служить временные ряды производства продукции в натуральном или стоимостном выражении за месяц, квартал, год и т.д.
Иногда уровни ряда представляют собой не непосредственно наблюдаемые значения, а производные величины: средние или относительные. Такие ряды называются производными. Уровни таких временных рядов получаются с помощью некоторых вычислений на основе непосредственно наблюдаемых показателей. Примерами таких рядов могут служить ряды среднесуточного производства основных видов промышленной продукции или ряды индексов цен.
Уровни ряда могут принимать детерминированные или случайные значения. Примером ряда с детерминированными значениями уровней служит ряд последовательных данных о количестве дней в месяцах. Естественно, анализу, а в дальнейшем и прогнозированию, подвергаются ряды со случайными значениями уровней. В таких рядах каждый уровень может рассматриваться как реализация случайной величины - дискретной или непрерывной.
1.2 Методы анализа временных рядов
Методы анализа временных рядов. Для решения этих задач существует большое количество различных методов. Из них наиболее распространенными являются следующие:
1. Корреляционный анализ, позволяющий выявить существенные периодические зависимости и их лаги (задержки) внутри одного процесса (автокорреляция) или между несколькими процессами (кросскорреляция);
2. Спектральный анализ, позволяющий находить периодические и квазипериодические составляющие временного ряда;
3. Сглаживание и фильтрация, предназначенные для преобразования временных рядов с целью удаления из них высокочастотных или сезонных колебаний;
5. Прогнозирование, позволяющее на основе подобранной модели поведения временного рада предсказывать его значения в будущем.
2 Основы прогнозирования развития перерабатывающих отраслей и торговых организаций
2.1 Прогнозирование развития перерабатывающих предприятий
Сельскохозяйственная продукция производится на предприятиях различных организационных форм. Здесь она может храниться, сортироваться и готовиться к переработке, вместе с тем могут быть и специализированные предприятия хранения. Дальше продукция транспортируется на перерабатывающие предприятия, где производится разгрузка, хранение, сортировка, переработка, фасовка; отсюда осуществляется транспортировка в торговые предприятия. На самих же предприятиях торговли производится реализация послепродажная упаковка и доставка.
Все виды перечисленных технологически и организационных операций должны прогнозироваться и планироваться. При этом используются различные приемы и методы.
Но надо отметить, что пищевые перерабатывающие предприятия имеют некоторую специфику планирования.
Пищевая перерабатывающая промышленность занимает важное место в системе АПК. Сельскохозяйственное производство обеспечивает эту промышленность сырьевыми ресурсами, то есть по существу, имеется жесткая технологическая связь между сферами 2 и 3 АПК.
В зависимости от вида используемого сырья и особенностей реализации конечной продукции сложились три группы отраслей пищевой и перерабатывающей промышленности: первичной и вторичной переработки сельскохозяйственных ресурсов и добывающей пищевой промышленности. В первую группу входят отрасли, которые перерабатывают малотранспортабельную сельскохозяйственную продукцию (крахмалопаточная, плодоовощеконсервная, спиртовая и др.), во вторую – отрасли использующие сельскохозяйственное сырьё, которое прошло первичную переработку (хлебопекарная, кондитерская, пищеконцетратная, производство сахара рафинада и др.). К третьей группе относятся соленая и рыбная отрасли.
Предприятия первой группы располагаются ближе к районам производства сельскохозяйственное продукции, здесь производство носит сезонный характер. Предприятия второй группы тяготеют, как правила, к районам потребления этой продукции; они работают ритмично на протяжении всего года.
Наряду с общими особенностями предприятия всех трех групп имеют свои внутренние, обусловленные номенклатурой выпускаемой продукции, в используемых технических средствах, технологиях, организации труда и производства и др.
Важным исходным началом прогнозирования этих отраслей является учет внешних и внутренних особенностей, специфики каждой отрасли промышленности.
В состав пищевых и перерабатывающих отраслей АПК входят зерноперерабатывающая, хлебопекарная и макаронная, сахарная, маложирная, кондитерская, плодоовощная, пищеконцетратная и др.
2.2 Прогнозирование развития торговых организаций
В торговле при ее прогнозировании используются те же методы, что и в других отраслях народного хозяйства. Перспективными являются создание рыночных структур в виде сети оптовых продовольственных рынков, совершенствование фирменной торговли, а также создание широкой информационной сети. Оптовая торговля позволяет сократить количество посредников при доведении продукции от товаропроизводителя до потребителя, создать альтернативные каналы реализации, точнее прогнозировать потребительский спрос и предложение.
В большинстве случаев план экономического и социального развития торгового предприятия состоит в основном из пяти разделов: розничный и оптовый товарооборот и товарное обеспечение; финансовый план; развитие материально-технической базы; социально развитие коллективов; план по труду.
Планы могут разрабатываться в виде долгосрочных – до 10 лет, среднесрочных – от трех до пяти лет, текущих – до одного месяца.
В основе планирования – товарооборот по каждой ассортиментной группе товаров.
Оптовый и розничный товарооборот может прогнозироваться в следующей последовательности:
1. оценивают ожидаемое выполнение плана за текущий год;
2. исчисляют среднегодовые темпы товарооборота за два-три года, предшествовавших периоду прогноза;
3. на основании анализа первых двух позиций экспертным методом устанавливают в процентах темпы роста (снижения) продажи отдельных товаров (товарных групп на прогнозируемый период).
Умножением объема ожидаемого товарооборота за текущий год на прогнозируемый темп роста продажи рассчитывают возможный товарооборот в прогнозируемом периоде.
Необходимые товарные ресурсы состоят из ожидаемого товарооборота и товарных запасов. Товарные запасы могут измеряться в натуральном и денежном выражении или в днях оборота. Товарные запасы обычно планируют на основе экстраполяции данных по четвертому кварталу за ряд лет.
Товарное обеспечение определяют путем сравнения потребности в необходимых товарных ресурсах и их источников. Необходимые товарные ресурсы рассчитывают как сумму товарооборота, вероятного прироста товарных запасов за минусом естественной убыли товаров и их уценки.
Финансовый план торгового предприятия включает кассовый план, кредитный план и сметы доходов и расходов. Кассовый план составляю по квартально, в кредитном плане определяют потребность в различных видах кредита, в смете доходов и расходов – по статьям доходы и поступления денежных средств, расходы и отчисления средств.
Объектами планирования материально-технической базы является торговая сеть, техническое оснащение, складское хозяйство, то есть планируются общая потребность в торговой площади, торговых предприятиях, их размещение и специализация, потребность в механизмах и оборудовании, необходимые складские емкости.
Показатели социального развития коллектива включают разработку планов повышения квалификации, улучшения условий труда и охраны здоровья работников, жилищных и культурно-бытовых условий, развития общественной активности.
Достаточно сложным разделом является план по труду. Необходимо подчеркнуть, что в торговле результатом труда выступает не продукт, а услуга, здесь преобладают затраты живого труда в связи с затруднением механизации большинства трудоемких процессов.
Производительность труда в торговле измеряется показателями среднего товарооборота, приходящегося на одного работника за определенный период времени, то есть сумма товарооборота делится на среднесписочную численность работников. В связи с тем что реализация различных товаров по своей трудоемкости не одинакова, при планировании следует учитывать изменения в товарообороте, индексы цен, ассортиментный состав товаров.
Развитие товарооборота требует увеличения количества предприятий торговли, общественного питания. При расчете количества на плановый период исходя из нормативов обеспеченности населения торговыми предприятиями для городской и сельской местности.
В качестве примера приведем содержание плана экономического и социального развития плодоовощного торгового предприятия. Он включает следующие разделы: исходные данные; основные экономические показатели работы предприятия; техническое и организационное развитие предприятия; план закладки продукции на длительное хранение; план реализации продукции; план розничного товарооборота; распределение издержек по завозу, хранению и оптовой реализации по группам товаров; издержки обращения розничной реализации продукции; затраты на производства продукции, ее переработку и реализацию; численность работников и фонд заработной планы; прибыль от оптовой реализации продукции; план прибыли от всех видов деятельности; распределение дохода; распределение прибыли; социальное развитие коллектива; финансовый план. Методика составления этого плана такая же, как и в других отраслях АПК.
3 Расчет прогноза экономического временного ряда
Имеются данные об экспорте железобетонной продукции товаров (в страны вне СНГ), млрд. долларов США.
Таблица 1
Экспорт товаров за 2002, 2003, 2004, 2005 годы (млрд.долларов США)
Прежде, чем приступить к анализу, обратимся к графическому изображению исходных данных (рис. 1).
Рис. 1. Экспорт товаров
Как видно из построенного графика, четко прослеживается тенденция к увеличению объемов импорта. Проанализировав полученный график можно сделать вывод о нелинейном развитии процесса, предположив об экспоненциальном или параболическом развитии.
Теперь сделаем графический анализ квартальных данных за четыре года:
Таблица 2
Экспорт товаров за кварталы 2002,2003, 2004 и 2005 годов
Рис. 2. Экспорт товаров
Как видно из графика яркое выражение имеет сезонность колебаний. Амплитуда колебания довольно не фиксированная, что указыает на наличие мультипликативной модели.
В исходных данных нам представлен интервальный ряд с равноотстоящими уровнями во времени. Поэтому для определения среднего уровня ряда воспользуемся следующей формулой:
Млрд.долл.
Для количественной оценки динамики явлений применяются следующие основные аналитические показатели:
· абсолютный прирост;
· темпы роста;
· темпы прироста.
Рассчитаем каждый из этих показателей для интервального ряда с равноотстоящими уровнями во времени.
Представим статистические показатели динамики в виде таблицы 3.
Таблица 3
Статистические показатели динамики
| t | y t | Абсолютный прирост, млрд.долларов США | Темп роста, % | Темп прироста, % | |||
| Цепной | Базисный | Цепной | Базисный | Цепной | Базисный | ||
| 1 | 48,8 | - | - | - | - | - | - |
| 2 | 61,0 | 12,2 | 12,2 | 125 | 125 | 25 | 25 |
| 3 | 77,5 | 16,5 | 28,7 | 127,05 | 158,81 | 27,05 | 58,81 |
| 4 | 103,5 | 26 | 54,7 | 133,55 | 212,09 | 33,55 | 112,09 |
Темпы роста были примерно одинаковые. Это говорит о том, что для определения прогнозного значения можно использовать средний темп роста:
Проверим гипотезу о наличии тренда с помощью критерия Фостера-Стюарта . Для этого заполним вспомогательную таблицу 4:
Таблица 4
Вспомогательная таблица
| t | yt | mt | lt | d | t | yt | mt | lt | d |
| 1 | 9,8 | - | - | - | 9 | 16,0 | 0 | 0 | 0 |
| 2 | 11,8 | 1 | 0 | 1 | 10 | 18,0 | 1 | 0 | 1 |
| 3 | 12,6 | 1 | 0 | 1 | 11 | 19,8 | 1 | 0 | 1 |
| 4 | 14,6 | 1 | 0 | 1 | 12 | 23,7 | 1 | 0 | 1 |
| 5 | 12,9 | 0 | 0 | 0 | 13 | 21,0 | 0 | 0 | 0 |
| 6 | 14,7 | 1 | 0 | 1 | 14 | 23,9 | 1 | 0 | 1 |
| 7 | 15,5 | 1 | 0 | 1 | 15 | 26,9 | 1 | 0 | 1 |
| 8 | 17,8 | 1 | 0 | 1 | 16 | 31,7 | 1 | 0 | 1 |
Применим критерий Стьюдента:
Получаем, , то есть 
, следовательно гипотеза Н
0
отвергается, тренд есть.
Проанализируем структуру временного ряда с использованием коэффициента автокорреляции.
Найдем последовательно коэффициенты автокорреляции:
–
коэффициент автокорреляции первого порядка, так как сдвиг во времени равен единице (-лаг).
Аналогично находим остальные коэффициенты.
– коэффициент автокорреляции второго порядка.
– коэффициент автокорреляции третьего порядка.
– коэффициент автокорреляции четвертого порядка.
Таким образом, мы видим, что самым высоким является коэффициент автокорреляции четвертого порядка. Это говорит о том, что во временном ряде присутствуют сезонные колебания с периодичностью в четыре квартала.
Проверим значимость коэффициента автокорреляции. Для этого введем две гипотезы: Н 0: , Н 1: .
Находится по таблице критических значений отдельно для >0 и <0. Причем, если ||>||, то принимается гипотеза Н 1 ,то есть коэффициент значим. Если ||<||, то принимается гипотеза Н 0 и коэффициент автокорреляции незначим. В нашем случае коэффициент автокорреляции достаточно велик, и проверять его значимость необязательно.
Требуется провести сглаживание временного ряда и восстановить потерянные уровни.
Проведем сглаживание временного ряда с помощью простой скользящей средней. Результаты расчетов представим в виде следующей таблицы 13.
Таблица 5
Сглаживание исходного ряда с помощью скользящей средней
| № года | № квартала | t | Импорт товаров, млрд.долларов США, yt | Скользящая средняя, | |
| 1 | I | 1 | 9,8 | - | - |
| II | 2 | 11,8 | - | - | |
| III | 3 | 12,6 | 12 , 59 | 1,001 | |
| IV | 4 | 14,6 | 13,34 | 1,094 | |
| 2 | I | 5 | 12,9 | 14,06 | 0,917 |
| II | 6 | 14,7 | 14,83 | 0,991 | |
| III | 7 | 15,5 | 15,61 | 0,993 | |
| IV | 8 | 17,8 | 16,41 | 1,085 | |
| 3 | I | 9 | 16 | 17,36 | 0,922 |
| II | 10 | 18 | 18,64 | 0,966 | |
| III | 11 | 19,8 | 20,0 | 0,990 | |
| IV | 12 | 23,7 | 21,36 | 1,110 | |
| 4 | I | 13 | 21 | 22,99 | 0,913 |
| II | 14 | 23,9 | 24,88 | 0,961 | |
| III | 15 | 26,9 | - | - | |
| IV | 16 | 31,7 | - | - |
Теперь рассчитаем отношение фактических значений к уровням сглаженного ряда. В результате получим временной ряд, уровни которого отражают влияние случайных факторов и сезонности.
Предварительные оценки сезонной составляющей получим усреднением уровней временного ряда для одноименных кварталов:
Для I квартала:
Для II квартала:
Для II квартала:
Для IV квартала: 
Взаимопогашаемость сезонных воздействий в мультипликативной форме выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна числу фаз в цикле. В нашем случае число фаз равно четырем. Просуммировав средние значения по кварталам, получаем:
Поскольку сумма получилась неравной четырем, необходимо произвести корректировку значений сезонной составляющей. Найдем поправку, на которую надо изменить предварительные оценки сезонности:
Определяем скорректированные значения сезонной, результаты сведем в таблицу 6.
Таблица 6
Оценивание сезонной компоненты в мультипликативной модели .
| № квартала | i | Предварительная оценка сезонной компоненты, | Скорректированное значение сезонной компоненты, |
| I | 1 | 0,917 | 0,921 |
| II | 2 | 0,973 | 0,978 |
| III | 3 | 0,995 | 1,000 |
| IV | 4 | 1,096 | 1,101 |
| 3,981 | 4 |
Проводим сезонную корректировку исходных данных, то есть, удаляем сезонную составляющую.
Таблица 7
Построение мультипликативной тренд сезонной модели.
| t | Импорт товаров , млрд.долларов США | Сезонная компонента, | Десезонализированный импорт товаров, | Расчетное значение, | Расчетное значение импорта товаров,  |
| 1 | 9,8 | 0,921 | 10,6406 | 11,48 | 10,57308 |
| 2 | 11,8 | 0,978 | 12,0654 | 11,85 | 11,5893 |
| 3 | 12,6 | 1 | 12,6 | 12,32 | 12,32 |
| 4 | 14,6 | 1,101 | 13,2607 | 12,89 | 14,19189 |
| 5 | 12,9 | 0,921 | 14,0065 | 13,56 | 12,48876 |
| 6 | 14,7 | 0,978 | 15,0307 | 14,33 | 14,01474 |
| 7 | 15,5 | 1 | 15,5 | 15,2 | 15,2 |
| 8 | 17,8 | 1,101 | 16,1671 | 16,17 | 17,80317 |
| 9 | 16 | 0,921 | 17,3724 | 17,24 | 15,87804 |
| 10 | 18 | 0,978 | 18,4049 | 18,41 | 18,00498 |
| 11 | 19,8 | 1 | 19,8 | 19,68 | 19,68 |
| 12 | 23,7 | 1,101 | 21,5259 | 21,05 | 23,17605 |
| 13 | 21 | 0,921 | 22,8013 | 22,52 | 20,74092 |
| 14 | 23,9 | 0,978 | 24,4376 | 24,09 | 23,56002 |
| 15 | 26,9 | 1 | 26,9 | 25,76 | 25,76 |
| 16 | 31,7 | 1,101 | 28,792 | 27,53 | 30,31053 |
По МНК получаем следующее уравнение тренда:3
Изобразим графически ряд остатков:
Рис. 3. График остатков
Проанализировав полученный график можно сделать вывод о случайности колебаний этого ряда.
Так же качество модели можно проверить с помощью показателей асимметрии и эксцесса остатков. В нашем случае получаем:
,
то гипотеза о нормальном распределении остатков отвергается.
Поскольку одно из неравенств выполняется, то уместен вывод о том, что гипотеза о нормальном характере распределения остатков отвергается.
Заключительным этапом применения кривых роста является расчет прогнозов на базе выбранного уравнения.
Для прогнозирования импорта товаров в следующем году оценим значения тренда при t =17, t =18, t =19 и t =20:
4. Личко Н.М. Планирование на предприятиях АПК. – М., 1996.
5. Финам. События и рынки, – http://www.finam.ru/
Цели анализа временных рядов. При практическом изучении временных радов на основании экономических данных на определенном промежутке времени эконометрист должен сделать выводы о свойствах этого ряда и о вероятностном механизме, порождающем этот ряд. Чаще всего при изучении временных рядов ставятся следующие цели:
1. Краткое (сжатое) описание характерных особенностей ряда.
2. Подбор статистической модели, описывающей временной ряд.
3. Предсказание будущих значений на основе прошлых наблюдений.
4. Управление процессом, порождающим временной ряд.
На практике эти и подобные цели достижимы далеко не всегда и далеко не в полной мере. Часто этому препятствует недостаточный объем наблюдений из-за ограниченного времени наблюдений. Еще чаще – изменяющаяся с течением времени статистическая структура временного ряда.
Стадии анализа временных рядов. Обычно при практическом анализе временных рядов последовательно проходят следующие этапы:
1. Графическое представление и описание поведения временного рада.
2. Выделение и удаление закономерных составляющих временного рада, зависящих от времени: тренда, сезонных и циклических составляющих.
3. Выделение и удаление низко- или высокочастотных составляющих процесса (фильтрация).
4. Исследование случайной составляющей временного ряда, оставшейся после удаления перечисленных выше составляющих.
5. Построение (подбор) математической модели для описания случайной составляющей и проверка ее адекватности.
6. Прогнозирование будущего развития процесса, представленного временным рядом.
7. Исследование взаимодействий между различными временными радами.
Для решения этих задач существует большое количество различных методов. Из них наиболее распространенными являются следующие:
8. Корреляционный анализ, позволяющий выявить существенные периодические зависимости и их лаги (задержки) внутри одного процесса (автокорреляция) или между несколькими процессами (кросскорреляция).
9. Спектральный анализ, позволяющий находить периодические и квазипериодические составляющие временного ряда.
10. Сглаживание и фильтрация, предназначенные для преобразования временных рядов с целью удаления из них высокочастотных или сезонных колебаний.
12. Прогнозирование, позволяющее на основе подобранной модели поведения временного рада предсказывать его значения в будущем.
Модели тренда
простейшие модели тренда. Приведем модели трендов, наиболее часто используемые при анализе экономических временных рядов, а также во многих других областях. Во-первых, это простая линейная модель
где а 0 , а 1 – коэффициенты модели тренда;
t – время.
В качестве единицы времени, может быть, час, день (сутки), неделя, месяц, квартал или год. Модель 269, несмотря на свою простоту, оказывается полезной во многих реальных задачах. Если нелинейный характер тренда очевиден, то может подойти одна из следующих моделей:
1. Полиномиальная:
(270)
где значение степени полинома п в практических задачах редко превышает 5;
2. Логарифмическая:
Эта модель чаще всего применяется для данных, имеющих тенденцию сохранять постоянные темпы прироста;
3. Логистическая:
(272)
4. Гомперца
(273), где
Две последние модели задают кривые тренда S-образной формы. Они соответствуют процессам с постепенно возрастающими темпами роста в начальной стадии и постепенно затухающимитемпами роста в конце. Необходимость подобных моделей обусловлена невозможностью многих экономических процессов продолжительное время развиваться с постоянными темпами роста или по полиномиальным моделям, в связи с их довольно быстрым ростом (или уменьшением).
При прогнозировании тренд используют в первую очередь для долговременных прогнозов. Точность краткосрочных прогнозов, основанных только на подобранной кривой тренда, как правило, недостаточна.
Для оценки и удаления трендов из временных рядов чаще всего используется метод наименьших квадратов. Этот метод достаточно подробно рассматривался во втором разделе пособия в задачах линейного регрессионного анализа. Значения временного ряда рассматриваюткак отклик (зависимую переменную), а время t – какфактор, влияющий на отклик (независимую переменную).
Для временных рядов характерна взаимная зависимость его членов (по крайней мере, не далеко отстоящих по времени) и это является существенным отличием от обычного регрессионного анализа, для которого все наблюдения предполагаются независимыми. Тем не менее, оценки тренда и в этих условиях обычно оказываются разумными, если выбрана адекватная модель тренда и если среди наблюдений нет больших выбросов. Упомянутые выше нарушения ограничений регрессионного анализа сказываются не столько на значениях оценок, сколько наих статистических свойствах. Так, при наличии заметной зависимости между членами временного ряда оценки дисперсии, основанные на остаточнойсумме квадратов, дают неправильные результаты. Неправильными оказываются и доверительные интервалы для коэффициентов модели, и т.д. В лучшем случае их можно рассматривать как очень приближенные.
АНАЛИЗ ВРЕМЕННЫХ РЯДОВ
ВВЕДЕНИЕ
ГЛАВА 1. АНАЛИЗ ВРЕМЕННЫХ РЯДОВ
1.1 ВРЕМЕННОЙ РЯД И ЕГО ОСНОВНЫЕ ЭЛЕМЕНТЫ
1.2 АВТОКОРРЕЛЯЦИЯ УРОВНЕЙ ВРЕМЕННОГО РЯДА И ВЫЯВЛЕНИЕ ЕГО СТРУКТУРЫ
1.3 МОДЕЛИРОВАНИЕ ТЕНДЕНЦИИ ВРЕМЕННОГО РЯДА
1.4 МЕТОД НАИМЕНЬШИХ КВАДРАТОВ
1.5 ПРИВЕДЕНИЕ УРАВНЕНИЯ ТРЕНДА К ЛИНЕЙНОМУ ВИДУ
1.6 ОЦЕНКА ПАРАМЕТРОВ УРАВНЕНИЯ РЕГРЕССИИ
1.7 АДДИТИВНАЯ И МУЛЬТИПЛИКАТИВНАЯ МОДЕЛИ ВРЕМЕННОГО РЯДА
1.8 СТАЦИОНАРНЫЕ ВРЕМЕННЫЕ РЯДЫ
1.9 ПРИМЕНЕНИЕ БЫСТРОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ К СТАЦИОНАРНОМУ ВРЕМЕННОМУ РЯДУ
1.10 АВТОКОРРЕЛЯЦИЯ ОСТАТКОВ. КРИТЕРИЙ ДАРБИНА- УОТСОНА
Введение
Почти в каждой области встречаются явления, которые интересно и важно изучать в их развитии и изменении во времени. В повседневной жизни могут представлять интерес, например, метеорологические условия, цены на тот или иной товар, те или иные характеристики состояния здоровья индивидуума и т. д. Все они изменяются во времени. С течением времени изменяются деловая активность, режим протекания того или иного производственного процесса, глубина сна человека, восприятие телевизионной программы. Совокупность измерений какой-либо одной характеристики подобного рода в течение некоторого периода времени представляют собой временной ряд.
Совокупность существующих методов анализа таких рядов наблюдений называется анализом временных рядов.
Основной чертой, выделяющей анализ временных рядов среди других видов статистического анализа, является существенность порядка, в котором производятся наблюдения. Если во многих задачах наблюдения статистически независимы, то во временных рядах они, как правило, зависимы, и характер этой зависимости может определяться положением наблюдений в последовательности. Природа ряда и структура порождающего ряд процесса могут предопределять порядок образования последовательности.
Цель работы состоит в получении модели для дискретного временного ряда во временной области, обладающей максимальной простотой и минимальным числом параметров и при этом адекватно описывающей наблюдения.
Получение такой модели важно по следующим причинам:
1) она может помочь понять природу системы, генерирующей временные ряды;
2) управлять процессом, порождающим ряд;
3) ее можно использовать для оптимального прогнозирования будущих значений временных рядов;
Временные ряды лучше всего описываются нестационарными моделями, в которых тренды и другие псевдоустойчивые характеристики, возможно меняющиеся во времени, рассматриваются скорее как статистические, а не детерминированные явления. Кроме того, временные ряды, связанные с экономикой, часто обладают заметными сезонными , или периодическими, компонентами; эти компоненты могут меняться во времени и должны описываться циклическими статистическими (возможно, нестационарными) моделями.
Пусть наблюдаемым временным рядом является y 1 , y 2 , . . ., y n . Мы будем понимать эту запись следующим образом. Имеется Т чисел, представляющих собой наблюдение некоторой переменной в Т равноотстоящих моментов времени. Эти моменты для удобства пронумерованы целыми числами 1, 2, . . .,Т. Достаточно общей математической (статистической или вероятностной) моделью служит модель вида:
y t = f(t) + u t , t = 1, 2, . . ., T.
В этой модели наблюдаемый ряд рассматривается как сумма некоторой полностью детерминированной последовательности {f(t)}, которую можно назвать математической составляющей, и случайной последовательности {u t }, подчиняющейся некоторому вероятностному закону. (И иногда для этих двух составляющих используются соответственно термины сигнал и шум). Эти компоненты наблюдаемого ряда ненаблюдаемы; они являются теоретическими величинами. Точный смысл указанного разложения зависит не только от самих данных, но частично и оттого, что понимается под повторением эксперимента, результатом которого являются эти данные. Здесь используется так называемая «частотная» интерпретация. Полагается, что, по крайней мере, принципиально можно повторять всю ситуацию целиком, получая новые совокупности наблюдений. Случайные составляющие, кроме всего прочего, могут включать в себя ошибки наблюдений.
В данной работе рассмотрена модель временного ряда, в которой на тренд накладывается случайная составляющая, образующая случайный стационарный процесс. В такой модели предполагается, что течение времени никак не отражается на случайной составляющей. Точнее говоря, предполагается, что математическое ожидание (то есть среднее значение) случайной составляющей тождественно равно нулю, дисперсия равна некоторой постоянной и что значения u t в различные моменты времени некоррелированны. Таким образом, всякая зависимость от времени включается в систематическую составляющую f(t). Последовательность f(t) может зависеть от некоторых неизвестных коэффициентов и от известных величин, меняющихся со временем. В этом случае её называют «функцией регрессии». Методы статистических выводов для коэффициентов функции регрессии оказываются полезными во многих областях статистики. Своеобразие же методов, относящихся именно к временным рядам, состоит в том, что здесь исследуются те модели, в которых упомянутые выше величины, меняющиеся со временем, являются известными функциями t.
Глава 1. Анализ временных рядов
1.1 Временной ряд и его основные элементы
Временной ряд –это совокупность значений какого-либо показателя за несколько последовательных моментов или периодов времени. Каждый уровень временного ряда формируется под воздействием большого числа факторов, которые условно можно подразделить на три группы:
· факторы, формирующие тенденцию ряда;
· факторы, формирующие циклические колебания ряда;
· случайные факторы.
При различных сочетаниях в изучаемом процессе или явлении этих факторов зависимость уровней ряда от времени может принимать различные формы. Во-первых , большинство временных рядов экономических показателей имеют тенденцию, характеризующую долговременное совокупное воздействие множества факторов на динамику изучаемого показателя. Очевидно, что эти факторы, взятые в отдельности, могут оказывать разнонаправленное влияние на исследуемый показатель. Однако в совокупности они формируют его возрастающую или убывающую тенденцию.
Во-вторых, изучаемый показатель может быть подвержен циклическим колебаниям. Эти колебания могут носить сезонный характер, поскольку деятельность ряда отраслей экономики и сельского хозяйства зависит от времени года. При наличии больших массивов данных за длительные промежутки времени можно выявить циклические колебания, связанные с общей динамикой временного ряда.
Некоторые временные ряды не содержат тенденции и циклической компоненты, а каждый следующий их уровень образуется как сумма среднего уровня ряда и некоторой(положительной или отрицательной) случайной компоненты.
В большинстве случаев фактический уровень временного ряда можно представить как сумму или произведение трендовой, циклической и случайной компонент. Модель, в которой временной ряд представлен как сумма перечисленных компонент, называется аддитивной моделью временного ряда. Модель, в которой временной ряд представлен как произведение перечисленных компонент, называется мультипликативной моделью временного ряда. Основная задача статистического исследования отдельного временного ряда – выявление и придание количественного выражения каждой из перечисленных выше компонент с тем чтобы использовать полученную информацию для прогнозирования будущих значений ряда.
1.2 Автокорреляция уровней временного ряда и выявление его структуры
При наличии во временном ряде тенденции и циклических колебаний значения каждого последующего уровня ряда зависят от предыдущих. Корреляционную зависимость между последовательными уровнями временного ряда называют автокорреляцией уровней ряда .
Количественно её можно измерить с помощью линейного коэффициента корреляции между уровнями исходного временного ряда и уровнями этого ряда, сдвинутыми на несколько шагов во времени.
Одна из рабочих формул для расчёта коэффициента автокорреляции имеет вид:
(1.2.1)В качестве переменной х мы рассмотрим ряд y 2 , y 3 , … , y n ; в качестве переменной у – ряд y 1 , y 2 , . . . ,y n – 1 . Тогда приведённая выше формула примет вид:
(1.2.2)Аналогично можно определить коэффициенты автокорреляции второго и более высоких порядков. Так, коэффициент автокорреляции второго порядка характеризует тесноту связи между уровнями у t и y t – 1 и определяется по формуле
(1.2.3)Число периодов, по которым рассчитывается коэффициент автокорреляции, называют лагом . С увеличением лага число пар значений, по которым рассчитывается коэффициент автокорреляции, уменьшается. Некоторые авторы считают целесообразным для обеспечения статистической достоверности коэффициентов автокорреляции использовать правило – максимальный лаг должен быть не больше (n/4).