
эконометрика / эконометрика / 4.0 Одномерные временные ряды
.pdf
Лабораторная работа № 4. МОДЕЛИРОВАНИЕ ОДНОМЕРНЫХ ВРЕМЕННЫХ РЯДОВ
Постановка задачи: известны статистические данные наблюдений за некоторое количество моментов или периодов времени.
Требуется:
1.Построить график динамики уровней ряда.
2.Рассчитать значения сезонных компонент методом скользящей средней.
3.Устранить сезонную компоненту из исходных уровней ряда. Построить уравнение, моделирующее динамику трендовой компоненты.
4.Найти прогноз фактора y.
Пример 4. Поквартальные данные об объемах потребления электроэнергии жителями региона за 2001 – 2004 гг представлены в табл. 4.7. Построить аддитивную модель временного ряда, рассчитать прогноз объема потребления электроэнергии в 1 квартале 2005 года.
Таблица 4.7
№ квартала |
t |
yt |
1 кв. 2001 года |
1 |
6 |
2 кв. 2001 года |
2 |
4,4 |
3 кв. 2001 года |
3 |
5 |
4 кв. 2001 года |
4 |
9 |
1 кв. 2002 года |
5 |
7,2 |
2 кв. 2002 года |
6 |
4,8 |
3 кв. 2002 года |
7 |
6 |
4 кв. 2002 года |
8 |
10 |
1 кв. 2003 года |
9 |
8 |
2 кв. 2003 года |
10 |
5,6 |
3 кв. 2003 года |
11 |
6,4 |
4 кв. 2003 года |
12 |
11 |
1 кв. 2004 года |
13 |
9 |
2 кв. 2004 года |
14 |
6,6 |
3 кв. 2004 года |
15 |
7 |
4 кв. 2004 года |
16 |
10,8 |
Порядок выполнения работы.
0. Подготовительная работа.
Занесем имеющиеся статистические данные в электронную таблицу (рис. 4.1).
|
|
Потребление электроэнергии |
||||||||||
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 12 13 |
14 15 16 |
|
|
|
|
|
|
|
кварталы |
|
|
|||
Рис. 4.1. Исходные данные |
Рис. 4.2. График динамики потребления электроэнергии |
1

1. Построить график динамики уровней ряда.
Построение графика динамики уровней ряда произведем при помощи функции МАСТЕР ДИАГРАММ. Выделим диапазон ячеек D4:D20 (рис. 4.1) и вызовем команду ВСТАВКА – ДИАГРАММЫ – ГРАФИК.
Построенную автоматически диаграмму (рис. 4.2) можно видоизменять при помощи закладок КОНСТРУКТОР и ФОРМАТ при выделенной диаграмме.
Изучая график потребления электроэнергии (рис. 4.2), заметим наличие сезонных колебаний с периодом 4 квартала и возрастающую тенденцию в уровнях ряда.
2. Рассчитать значения сезонных компонент методом скользящей средней.
Рис. 4.3. Расчет сезонных компонент |
Рис. 4.4. Расчет сезонных компонент |
методом скользящей средней |
|
Составим рабочую таблицу (рис. 4.3). Поскольку в нашем случае известны поквартальные данные (длина цикла равна 4), то найдем скользящие средние для каждой группы из четырех соседних значений уровней ряда. Для этого используем функцию СРЗНАЧ (ФОРМУЛЫ – СТАТИСТИЧЕСКИЕ – СРЗНАЧ).
Центрированные скользящие средние найдем, чтобы привести данные в соответствие с фактическими моментами времени. Они находятся как средние значения между парами соседних скользящих средних. Например, в ячейке F29 рассчитано среднее между значениями ячеек E28 и E29. Затем функцию можно скопировать вниз по столбцу.
Оценки сезонной компоненты для аддитивной модели рассчитаем, как разность между фактическими уровнями ряда и центрированными скользящими средними. В нашем примере G29 = D29 – F29 и т.д. Если нужно построить мультипликативную модель, то фактические уровни ряда делятся на центрированную скользящую среднюю.
В следующей таблице (рис. 4.4) рассчитаем оценки сезонных компонент, распределенных по кварталам. Первые две оценки сезонных компонент в ячейках G29:G30 скопируем (рис. 4.3) и вставим значения в таблицу (рис. 4.4) в ячейках D48:D49. Чтобы в таблицу попали именно значения стоит использовать СПЕЦИАЛЬНУЮ ВСТАВКУ: ГЛАВНАЯ – ВСТАВИТЬ – ЗНАЧЕНИЯ (либо воспользоваться правой клавишей мыши и в контекстном меню выбрать необходимую команду). Аналогично скопируем остальные данные. Таким образом, мы установили соответствие между оценками сезонных компонент, соответствующих определенному кварталу.
Найдем средние оценки сезонных компонент для каждого квартала (по годам) S i при помощи встроенной статистической функции СРЗНАЧ. Например, H46 = СРЗНАЧ (D46:G46) и т.д.
Корректирующий коэффициент в ячейке H50 считается как среднее арифметическое (СРЗНАЧ)
значение по всем S i для аддитивной модели, и как среднее геометрическое (СРГЕОМ) – для мультипликативной. Для аддитивной модели скорректированные значения сезонных компонент Si
рассчитаем как разность между средней оценкой S i и корректирующим коэффициентом. В нашем примере, I46 = H46 – H50 и т.д. Для мультипликативной модели находим частное от деления средней оценки на корректирующий коэффициент.
2

Следует помнить, что в аддитивной модели сумма значений сезонных компонент Si должна быть равна нулю, в мультипликативной – произведение сезонных компонент равно единице.
3. Устранить сезонную компоненту из исходных уровней ряда. Построить уравнение, моделирующее динамику трендовой компоненты.
Составим рабочую таблицу (рис. 4.5). Скопируем статистические данные в ячейки D54:D69. Столбец сезонных компонент St сформируем при помощи команды СПЕЦИАЛЬНАЯ ВСТАВКА. В следующем столбце таблицы рассчитаем выровненные значения ряда: для аддитивной модели yt St ,
для мультипликативной модели yt / St . Мы получили ряд, в котором отсутствует сезонная компонента.
Рис. 4.5. Устранение сезонных компонент из |
Рис. 4.6. Использование команды РЕГРЕССИЯ |
уровней ряда |
для нахождения трендовой компоненты уровней |
|
ряда |
Построим уравнение линейной регрессии, в котором в качестве независимого фактора выступает время t. Получим линейное уравнение вида: T b0 b1t , где T – значение трендовой компоненты, t –
номер наблюдения. Используем команду ДАННЫЕ – АНАЛИЗ ДАННЫХ – РЕГРЕССИЯ (рис 4.6). В меню РЕГРЕССИЯ обязательно поставим флажок Остатки, чтобы оценить наличие автокорреляции в остатках.
Для рассматриваемого примера получим числовые данные, представленные на рис. 4.7.
3

Рис. 4.7. Протокол выполнения регрессионного анализа
Оценим качество полученного уравнения по коэффициенту корреляции, детерминации, t- статистикам. Обратим внимание, что в ячейках D93:D108 находятся расчетные значения трендовой компоненты T для каждого наблюдения, а в ячейках E93:E108 – остатки (рис. 4.7).
Рассчитаем значение коэффициента Дарбина-Уотсона (рис. 4.8). Статистику Дарбина – Уотсона найдем по формуле:
n
(et et 1 )2
DW |
t 2 |
|
. |
(4.1) |
|
n |
|||
|
|
|
|
|
|
|
et2 |
|
|
t 1
Вячейке F109 найдем значение числителя, а в ячейке E109 – значение знаменателя, здесь
использована математическая функция СУММКВ, которая возвращает значение суммы квадратов заданного массива (рис. 4.8).
4

Рис. 4.8. Расчет коэффициента Дарбина-Уотсона
Для рассматриваемого примера все статистики оказались удовлетворительными, кроме коэффициента Дарбина-Уотсона. DW = 1,286, причем dL = 1,1 < 1,286 < 1,37 = dU для 16 наблюдений и 1 независимой переменной t. Таким образом, невозможно определить, присутствует ли автокорреляция первого порядка в остатках построенной модели. На практике в этом случае предполагают наличие автокорреляции.
4.Рассчитать прогноз самостоятельно.
5.Оформить отчет о проделанной работе.
План отчета:
1.Укажите фамилию, имя, название группы, номер варианта.
2.Запишите полученные значения сезонных компонент. Сделайте выводы по ним о поведении ряда.
3.Запишите уравнение линейного тренда в виде
T ... ... t, R2 ..., DW ...
tcm (...) (...)
4.Дайте экономическую интерпретацию параметра модели линейного тренда.
5.Дайте интерпретацию коэффициента детерминации и сделайте вывод о качестве модели.
6.Оцените статистическую значимость параметров модели.
7.Сделайте выводы об автокорреляции остатков модели.
8.Рассчитайте прогноз экономического показателя на ближайший временной промежуток, используя построенную модель.
9.Дополнительное задание: оцените прогнозные свойства модели с помощью средней ошибки аппроксимации.
© Денисейко И.В.
5