16.02.15 Виктор Гаврилов
44859 0
Временным рядом называется последовательность значений, изменяемых во времени. О некоторых простых, но эффективных подходах к работе с подобными последовательностями я попробую рассказать в данной статье. Примеров таких данных можно встретить очень много – котировки валют, объемы продаж, обращения клиентов, данные в различных прикладных науках (социология, метеорология, геология, наблюдения в физике) и многое другое.
Ряды являются распространенной и важной формой описания данных, так как позволяют наблюдать всю историю изменения интересующего нас значения. Это даёт нам возможность судить о «типичном» поведении величины и об отклонениях от такого поведения.
Передо мной встала задача выбрать набор данных, на котором можно было бы наглядно продемонстрировать особенности временных рядов. Я решил воспользоваться статистикой пассажиропотока на международных авиалиниях, поскольку этот набор данных весьма нагляден и стал своего рода стандартным (http://robjhyndman.com/tsdldata/data/airpass.dat , источник Time Series Data Library, R. J. Hyndman). Ряд описывает количество пассажиров международных авиалиний в месяц (в тысячах) за период с 1949 по 1960 года.
Поскольку у меня всегда под рукой , в которой есть интересный инструмент « » для работы с рядами, я воспользуюсь именно им. Перед импортом данных в файл нужно добавить столбец с датой, чтобы была привязка значений ко времени, и столбец с именем ряда для каждого наблюдения. Ниже видно, как выглядит мой исходный файл, который я импортировал в Prognoz Platform с помощью мастера импорта непосредственно из инструмента анализа временных рядов.
Первое, что мы обычно делаем с временным рядом, это отображаем его на графике. Prognoz Platform позволяет построить график, просто «перетащив» ряд в рабочую книгу.
Временной ряд на графике
Символ ‘M’ в конце имени ряда означает, что ряд имеет месячную динамику (интервал между наблюдениями равен одному месяцу).
Уже из графика мы видим, что ряд демонстрирует две особенности:
- тренд – на нашем графике это долгосрочный рост наблюдаемых значений. Видно, что тренд практически линейный.
- сезонность – на графике это периодические колебания величины. В следующей статье на тему временных рядов мы узнаем, как можно вычислить период.
Наш ряд достаточно «аккуратный», однако часто встречаются ряды, которые помимо двух описанных выше характеристик демонстрируют ещё одну – наличие «шума», т.е. случайных вариаций в той или иной форме. Пример такого ряда можно увидеть на графике ниже. Это синусоидальный сигнал, смешанный со случайной величиной.
При анализе рядов нас интересует выявление их структуры и оценка всех основных компонентов – тренда, сезонности, шума и других особенностей, а также возможность строить прогнозы изменения величины в будущих периодах.
При работе с рядами наличие шума часто затрудняет анализ структуры ряда. Чтобы исключить его влияние и лучше увидеть структуру ряда, можно использовать методы сглаживания рядов.
Самый простой метод сглаживания рядов – скользящее среднее. Идея заключается в том, что для любого нечётного количества точек последовательности ряда заменять центральную точку на среднее арифметическое остальных точек:
где x i – исходный ряд, s i – сглаженный ряд.
Ниже можно увидеть результат применения данного алгоритма к двум нашим рядам. Prognoz Platform по умолчанию предлагает использовать сглаживание с размером окна в 5 точек (k в нашей формуле выше будет равно 2). Обратите внимание, что сглаженный сигнал уже не так подвержен влиянию шума, однако вместе с шумом, естественно, пропадает и часть полезной информации о динамике ряда. Также видно, что у сглаженного ряда отсутствуют первые (и также последние) k точек. Это связано с тем, что сглаживание выполняется для центральной точки окна (в нашем случае для третьей точки), после чего окно сдвигается на одну точку, и вычисления повторяются. Для второго, случайного ряда, я использовал сглаживание с окном равным 30, чтобы лучше выявить структуру ряда, так как ряд «высокочастотный», точек очень много.
Метод скользящего среднего имеет определённые недостатки:
- Скользящее среднее неэффективно в вычислении. Для каждой точки среднее необходимо перевычислять по новой. Мы не можем переиспользовать результат, вычисленный для предыдущей точки.
- Скользящее среднее нельзя продлить на первые и последние точки ряда. Это может вызвать проблему, если нас интересуют именно эти точки.
- Скользящее среднее не определено за пределами ряда, и как следствие, не может использоваться для прогнозирования.
Экспоненциальное сглаживание
Более продвинутый метод сглаживания, который также можно использовать для прогнозирования – экспоненциальное сглаживание, также иногда называемое методом Хольта-Уинтерса (Holt-Winters) в честь имён его создателей.
Существует насколько вариантов данного метода:
- одинарное сглаживание для рядов, у которых нет тренда и сезонности;
- двойное сглаживание для рядов, у которых есть тренд, но нет сезонности;
- тройное сглаживание для рядов, у которых есть и тренд, и сезонность.
Метод экспоненциального сглаживания вычисляет значения сглаженного ряда путём обновления значений, рассчитанных на предыдущем шаге, используя информацию с текущего шага. Информация с предыдущего и текущего шагов берётся с разными весами, которыми можно управлять.
В простейшем варианте одинарного сглаживания соотношение такое:
Параметр α определяет соотношение между несглаженным значением на текущем шаге и сглаженным значением с предыдущего шага. При α =1 мы будем брать только точки исходного ряда, т.е. никакого сглаживания не будет. При α =0 ряд мы будем брать только сглаженные значения с предыдущих шагов, т.е. ряд превратится в константу.
Чтобы понять, почему сглаживание называется экспоненциальным, нам нужно раскрыть соотношение рекурсивно:
Из соотношения видно, что все предыдущие значения ряда вносят вклад в текущее сглаженное значение, однако их вклад угасает экспоненциально за счёт роста степени параметра α .
Однако, если в данных есть тренд, простое сглаживание будет «отставать» от него (либо придётся брать значения α близкими к 1, но тогда сглаживание будет недостаточным). Нужно использовать двойное экспоненциальное сглаживание.
Двойное сглаживание использует уже два уравнения – одно уравнение оценивает тренд как разницу между текущим и предыдущим сглаженным значениями, потом сглаживает тренд простым сглаживанием. Второе уравнение выполняет сглаживание как в случае простого варианта, но во втором слагаемом используется сумма предыдущего сглаженного значения и тренда.
Тройное сглаживание включает ещё один компонент – сезонность, и использует ещё одно уравнение. При этом различаются два варианта сезонного компонента – аддитивный и мультипликативный. В первом случае амплитуда сезонного компонента постоянна и со временем не зависит от базовой амплитуды ряда. Во втором случае амплитуда меняется вместе с изменением базовой амплитуды ряда. Это как раз наш случай, как видно из графика. С ростом ряда амплитуда сезонных колебаний увеличивается.
Так как наш первый ряд имеет и тренд, и сезонность, я решил подобрать параметры тройного сглаживания для него. В Prognoz Platform это довольно просто сделать, потому что при обновлении значения параметра платформа сразу же перерисовывает график сглаженного ряда, и визуально можно сразу увидеть, насколько хорошо он описывает наш исходный ряд. Я остановился на следующих значениях:
Как я вычислил период, мы рассмотрим в следующей статье о временных рядах.
Обычно в качестве первых приближений можно рассматривать значения между 0,2 и 0,4. Prognoz Platform также использует модель с дополнительным параметром ɸ , который дэмпфирует тренд так, что он приближается к константе в будущем. Для ɸ я взял значение 1, что соответствует обычной модели.
Также я сделал прогноз значений ряда данным методом на последние 2 года. На рисунке ниже я пометил точку начала прогноза, проведя через неё черту. Как видно, исходный ряд и сглаженный весьма неплохо совпадают, в том числе и на периоде прогнозирования – неплохо для такого простого метода!
Prognoz Platform также позволяет автоматически подобрать оптимальные значения параметров, используя систематический поиск в пространстве значений параметров и минимизируя сумму квадратов отклонений сглаженного ряда от исходного.
Описанные методы весьма просты, их легко применять, и они являются хорошей отправной точкой для анализа структуры и прогнозирования временных рядов.
Еще больше о временных рядах читайте в следующей статье.
Эконометрические модели начали использоваться для экономического прогнозирования в 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-модель используется для аппроксимации данных, имеющих вид различий первого периода.
Ценность новаторских идей Р. Ингла и К. Грэнджера заключается не только в том, что они предложили новые методы моделирования экономических зависимостей, но и в том, что разработанные ими модели открыли новые сферы исследований. При этом нобелианты фундаментально обосновали использование таких моделей, доказали корректность эконометрической оценки их параметров в случае нарушения ряда классических прогнозов. Важно и то, что каждый из предложенных методов подтвердил теоретические результаты.
3.3.1. Методы анализа и прогнозирования временных рядов
Модели стационарных и нестационарных временных рядов. Пусть Рассмотрим временной ряд X (t ). Пусть сначала временной ряд принимает числовые значения. Это могут быть, например, цены на батон хлеба в соседнем магазине или курс обмена доллара на рубли в ближайшем обменном пункте. Обычно в поведении временного ряда выявляют две основные тенденции - тренд и периодические колебания.
При этом под трендом понимают зависимость от времени линейного, квадратичного или иного типа, которую выявляют тем или иным способом сглаживания (например, экспоненциального сглаживания) либо расчетным путем, в частности, с помощью метода наименьших квадратов. Другими словами, тренд - это очищенная от случайностей основная тенденция временного ряда.
Временной ряд обычно колеблется вокруг тренда, причем отклонения от тренда часто обнаруживают правильность. Часто это связано с естественной или назначенной периодичностью, например, сезонной или недельной, месячной или квартальной (например, в соответствии с графиками выплаты заплаты и уплаты налогов). Иногда наличие периодичности и тем более ее причины неясны, и задача статистика - выяснить, действительно ли имеется периодичность.
Элементарные методы оценки характеристик временных рядов обычно достаточно подробно рассматриваются в курсах "Общей теории статистики" (см., например, учебники ), поэтому нет необходимости подробно разбирать их здесь. О некоторых современных методах оценивания длины периода и самой периодической составляющей речь пойдет ниже в подразделе 3.3.2.
Характеристики временных рядов. Для более подробного изучения временных рядов используются вероятностно-статистические модели. При этом временной ряд X (t ) рассматривается как случайный процесс (с дискретным временем). Основными характеристиками X (t ) являются математическое ожидание X (t ), т.е.
дисперсия X (t ), т.е.
и автокорреляционная функция временного ряда X (t )
т.е. функция двух переменных, равная коэффициенту корреляции между двумя значениями временного ряда X (t ) и X (s ).
В теоретических и прикладных исследованиях рассматривают широкий спектр моделей временных рядов. Выделим сначала стационарные модели. В них совместные функции распределения для любого числа моментов времени k , а потому и все перечисленные выше характеристики временного ряда не меняются со временем . В частности, математическое ожидание и дисперсия являются постоянными величинами, автокорреляционная функция зависит только от разности t - s. Временные ряды, не являющиеся стационарными, называются нестационарными.
Линейные регрессионные модели с гомоскедастичными и гетероскедастичными, независимыми и автокоррелированными остатками. Как видно из сказанного выше, основное - это "очистка" временного ряда от случайных отклонений, т.е. оценивание математического ожидания. В отличие от простейших моделей регрессионного анализа, рассмотренных в главе 3.2, здесь естественным образом появляются более сложные модели. Например, дисперсия может зависеть от времени. Такие модели называют гетероскедастичными, а те, в которых нет зависимости от времени - гомоскедастичными. (Точнее говоря, эти термины могут относиться не только к переменной "время", но и к другим переменным.)
Далее, в главе 3.2 предполагалось, что погрешности независимы между собой. В терминах настоящей главы это означало бы, что автокорреляционная функция должна быть вырожденной - равняться 1 при равенстве аргументов и 0 при их неравенстве. Ясно, что для реальных временных рядов так бывает отнюдь не всегда. Если естественный ход изменений наблюдаемого процесса является достаточно быстрым по сравнению с интервалом между последовательными наблюдениями, то можно ожидать "затухания" автокорреляции" и получения практически независимых остатков, в противном случае остатки будут автокоррелированы.
Идентификация моделей. Под идентификацией моделей обычно понимают выявление их структуры и оценивание параметров. Поскольку структура - это тоже параметр, хотя и нечисловой, то речь идет об одной из типовых задач прикладной статистики - оценивании параметров.
Проще всего задача оценивания решается для линейных (по параметрам) моделей с гомоскедастичными независимыми остатками. Восстановление зависимостей во временных рядах может быть проведено на основе методов наименьших квадратов и наименьших модулей оценивания параметров в моделях линейной (по параметрам) регрессии. На случай временных рядов переносятся результаты, связанные с оцениванием необходимого набора регрессоров, в частности, легко получить предельное геометрическое распределение оценки степени тригонометрического полинома.
Однако на более общую ситуацию такого простого переноса сделать нельзя. Так, например, в случае временного ряда с гетероскедастичными и автокоррелированными остатками снова можно воспользоваться общим подходом метода наименьших квадратов, однако система уравнений метода наименьших квадратов и, естественно, ее решение будут иными. Формулы в терминах матричной алгебры, о которых упоминалось в главе 3.2, будут отличаться. Поэтому рассматриваемый метод называется "обобщенный метод наименьших квадратов (ОМНК)".
Замечание. Как уже отмечалось в главе 3.2, простейшая модель метода наименьших квадратов допускает весьма далекие обобщения, особенно в области системам одновременных эконометрических уравнений для временных рядов. Для понимания соответствующей теории и алгоритмов необходимо владение методами матричной алгебры. Поэтому мы отсылаем тех, кому это интересно, к литературе по системам эконометрических уравнений и непосредственно по временным рядам , в которой особенно много интересуются спектральной теорией, т.е. выделением сигнала из шума и разложением его на гармоники. Подчеркнем еще раз, что за каждой главой настоящей книги стоит большая область научных и прикладных исследований, вполне достойная того, чтобы посвятить ей много усилий. Однако из-за ограниченности объема книги мы вынуждены изложение сделать конспективным.
Системы эконометрических уравнений. В качестве первоначального примера рассмотрим эконометрическую модель временного ряда, описывающего рост индекса потребительских цен (индекса инфляции). Пусть I (t ) - рост цен в месяц t (подробнее об этой проблематике см. главу 7 в ). По мнению некоторых экономистов естественно предположить, что
I (t ) = с I (t - 1) + a + bS (t - 4) + e , (1)
где I (t -1) - рост цен в предыдущий месяц (а с - некоторый коэффициент затухания, предполагающий, что при отсутствии внешний воздействий рост цен прекратится), a - константа (она соответствует линейному изменению величины I (t ) со временем), bS (t- 4) - слагаемое, соответствующее влиянию эмиссии денег (т.е. увеличения объема денег в экономике страны, осуществленному Центральным Банком) в размере S (t- 4) и пропорциональное эмиссии с коэффициентом b , причем это влияние проявляется не сразу, а через 4 месяца; наконец, e - это неизбежная погрешность.
Модель (1), несмотря на свою простоту, демонстрирует многие характерные черты гораздо более сложных эконометрических моделей. Во-первых, обратим внимание на то, что некоторые переменные определяются (рассчитываются) внутри модели, такие, как I (t ). Их называют эндогенными (внутренними). Другие задаются извне (это экзогенные переменные). Иногда, как в теории управления, среди экзогенных переменных, выделяют управляемые переменные - те, с помощью выбора значений которых можно привести систему в нужное состояние.
Во-вторых, в соотношении (1) появляются переменные новых типов - с лагами, т.е. аргументы в переменных относятся не к текущему моменту времени, а к некоторым прошлым моментам.
В-третьих, составление эконометрической модели типа (1) - это отнюдь не рутинная операция. Например, запаздывание именно на 4 месяца в связанном с эмиссией денег слагаемом bS (t- 4) - это результат достаточно изощренной предварительной статистической обработки. Далее, требует изучения вопрос зависимости или независимости величин S (t- 4) и I(t ) в различные моменты времени t . От решения этого вопроса зависит, как выше уже отмечалось, конкретная реализация процедуры метода наименьших квадратов.
С другой стороны, в модели (1) всего 3 неизвестных параметра, и постановку метода наименьших квадратов выписать нетрудно:
Проблема идентифицируемости. Представим теперь модель тапа (1) с большим числом эндогенных и экзогенных переменных, с лагами и сложной внутренней структурой. Вообще говоря, ниоткуда не следует, что существует хотя бы одно решение у такой системы. Поэтому возникает не одна, а две проблемы. Есть ли хоть одно решение (проблема идентифицируемости)? Если да, то как найти наилучшее решение из возможных? (Это - проблема статистической оценки параметров.)
И первая, и вторая задача достаточно сложны. Для решения обеих задач разработано множество методов, обычно достаточно сложных, лишь часть из которых имеет научное обоснование. В частности, достаточно часто пользуются статистическими оценками, не являющимися состоятельными (строго говоря, их даже нельзя назвать оценками).
Коротко опишем некоторые распространенные приемы при работе с системами линейных эконометрических уравнений.
Система линейных одновременных эконометрических уравнений. Чисто формально можно все переменные выразить через переменные, зависящие только от текущего момента времени. Например, в случае уравнения (1) достаточно положить
H (t ) = I (t- 1), G (t) = S (t- 4).
Тогда уравнение примет вид
I (t ) = с H (t ) + a + bG (t ) + e . (2)
Отметим здесь же возможность использования регрессионных моделей с переменной структурой путем введения фиктивных переменных. Эти переменные при одних значениях времени (скажем, начальных) принимают заметные значения, а при других - сходят на нет (становятся фактически равными 0). В результате формально (математически) одна и та же модель описывает совсем разные зависимости.
Косвенный, двухшаговый и трехшаговый методы наименьших квадратов. Как уже отмечалось, разработана масса методов эвристического анализа систем эконометрических уравнений. Они предназначены для решения тех или иных проблем, возникающих при попытках найти численные решения систем уравнений.
Одна из проблем связана с наличием априорных ограничений на оцениваемые параметры. Например, доход домохозяйства может быть потрачен либо на потребление, либо на сбережение. Значит, сумма долей этих двух видов трат априори равна 1. А в системе эконометрических уравнений эти доли могут участвовать независимо. Возникает мысль оценить их методом наименьших квадратов, не обращая внимания на априорное ограничение, а потом подкорректировать. Такой подход называют косвенным методом наименьших квадратов.
Двухшаговый метод наименьших квадратов состоит в том, что оценивают параметры отдельного уравнения системы, а не рассматривают систему в целом. В то же время трехшаговый метод наименьших квадратов применяется для оценки параметров системы одновременных уравнений в целом. Сначала к каждому уравнению применяется двухшаговый метод с целью оценить коэффициенты и погрешности каждого уравнения, а затем построить оценку для ковариационной матрицы погрешностей. После этого для оценивания коэффициентов всей системы применяется обобщенный метод наименьших квадратов.
Менеджеру и экономисту не следует становиться специалистом по составлению и решению систем эконометрических уравнений, даже с помощью тех или иных программных систем, но он должен быть осведомлен о возможностях этого направления эконометрики, чтобы в случае производственной необходимости квалифицированно сформулировать задание для специалистов по прикладной статистике.
От оценивания тренда (основной тенденции) перейдем ко второй основной задаче эконометрики временных рядов - оцениванию периода (цикла).
| Предыдущая |
АНАЛИЗ ВРЕМЕННЫХ РЯДОВ
ВВЕДЕНИЕ
ГЛАВА 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).
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/
Загрузка...
