Анализ и прогнозирование временных рядов в профессиональном статистическом пакете Statgraphics Centurion

Введение  

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

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

     Основными критериями выбора подходящей модели являются:

- адекватность  построения модели (когда модель  успешно проходит тесты);

- модель  имеет наименьшее количество  ошибок.

     STATGRAPHICS, программный продукт американской  фирмы STSC, Inc., уверенно занимает  одно из лидирующих мест в мире среди аналогичных систем и является наиболее распространенной в России зарубежной статистической и графической системой (в основном, в виде нелегальных копий).

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

     В целом STATGRAPHICS нацелен на использование  в различных областях: в промышленности, образовании, маркетинге, туризме  и

др.

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

1. Сущность и методы  прогнозирования  в ИГиТ. Их реализация  в пакете Statgraphics Plus. 
 

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

     Прогнозы бывают:

- краткосрочными,

- среднесрочными,

- долгосрочными.

На самом  деле не существует единства в понимании временных границ прогнозирования.

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

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

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

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

     Он состоит из 4 компонентов:

- цикла (с периодических колебаний с периодом более одного года),

- тренда (основного или эволюторного направления развития показателя),

- сезонных изменений (периодических колебаний с периодом, равным одному году),

- случайной компоненты или остатков (непрогнозируемой компоненты).

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

- аддитивная (когда найденные индексы сезонности добавляются к компоненту тренд-цикла),

- мультипликативная (индекс сезонности умножается на значение тренд-цикла).

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

     Мультипликативная модель используется, если эффект сезонности усиливается год от года (или уменьшается).

Существует много методов коллективных экспертных оценок. Наиболее популярными являются:

- метод круглого стола;

- метод Дельфи;

- метод эвристического прогнозирования;

- метод коллективной генерации идей;

- метод мозговой атаки. 

2. Анализ временного  ряда продаж шампанского. 

Для изучения продаж шампанского нам понадобятся  данные о динамике: 

Динамика  ежемесячных продаж шампанского.

месяц 1997 1998 1999 2000 2001 2002 2003
январь 2,815 2,541 3,113 5,375 3,633 4,016 2,639
февраль 2,672 2,475 3,006 3,088 4,292 3,957 2,899
март 2,755 3,031 4,047 3,718 4,154 4,51 3,37
апрель 2,721 3,266 3,523 4,514 4,121 4,276 3,74
май 2,946 3,776 3,937 4,52 4,647 4,968 2,927
июнь 3,036 3,23 3,986 4,539 4,753 4,677 3,986
июль 2,282 3,028 3,26 3,663 3,965 3,523 4,217
август 2,212 1,759 1,573 1,643 1,723 1,821 1,738
сентябрь 2,922 3,595 3,528 4,739 5,048 5,222 5,221
октябрь 4,301 4,474 5,211 5,428 6,922 6,872 6,424
ноябрь 5,764 6,838 7,614 8,314 9,858 10,803 9,842
декабрь 7,312 8,357 9,254 10,651 11,331 13,916 13,076
 

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

Временные ряды делятся на:

- дискретные,

- непрерывные.

     Дискретные ряды делятся на ряды с равноотстоящими и произвольными моментами наблюдения.

Рассмотрим  более подробно описательные методы анализа временного ряда на примере  ежемесячных продаж шампанского, представленного  в приложение 1.

Модуль  анализа временных рядов Statgraphics располагает пятью основными процедурами:

Descriptive Methods (описательный/дескриптивный метод);

Smoothing (сглаживание);

Seasonal Decomposition (сезонная декомпозиция);

Forecasting (прогнозирование);

Automatic Forecasting (автоматическое прогнозирование). 

     Рассмотрим дескриптивную модель. Для этого введем имя переменной, установим переключатель Month(s), в поле Starting At установим начальную дату переменной bubbly.

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

 

Рис.1

     На рисунке 1 показан горизонтальный  график ежемесячных продаж шампанского.

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

Предполагаемые  Автокорреляции для  шампанского 

      Lower 95,0% Upper 95,0%
Lag Autocorrelation Stnd. Error Prob. Limit Prob. Limit
1 0,489328 0,109109 -0,21385 0,21385
2 0,0578721 0,132687 -0,260062 0,260062
3 -0,154985 0,132987 -0,26065 0,26065
4 -0,250012 0,13512 -0,264831 0,264831
5 -0,0390596 0,140519 -0,275413 0,275413
6 0,0364701 0,140648 -0,275666 0,275666
7 -0,0377318 0,140761 -0,275887 0,275887
8 -0,246331 0,140881 -0,276123 0,276123
9 -0,181323 0,145919 -0,285996 0,285996
10 -0,0030676 0,148577 -0,291206 0,291206
11 0,373327 0,148578 -0,291207 0,291207
12 0,80455 0,159354 -0,312329 0,312329
13 0,406056 0,202004 -0,395922 0,395922
14 0,0254493 0,211498 -0,414529 0,414529
15 -0,173233 0,211534 -0,414601 0,414601
16 -0,244184 0,213217 -0,417898 0,417898
17 -0,0560871 0,21652 -0,424373 0,424373
18 0,0291966 0,216693 -0,424712 0,424712
19 -0,0333897 0,21674 -0,424803 0,424803
20 -0,206322 0,216801 -0,424923 0,424923
21 -0,146823 0,219126 -0,42948 0,42948
22 -0,0129479 0,220294 -0,43177 0,43177
23 0,278689 0,220303 -0,431787 0,431787
24 0,601805 0,224461 -0,439936 0,439936
 

The StatAdvisor

     Эта таблица показывает предполагаемую автокорреляцию ценности  шампанского. Пределы вероятности составляют 95,0 %. В этом случае 5 из  24 автокоэффициентов корреляции являются статистически существенными на доверительном уровне на 95,0 %. График выходит за пределы, это обозначает, что нужно открыть сезонную декампозицию и разобраться, в чём дело.

Рис. 2 

     Statgraphics построил график интегрированной периодограммы, мы видим это на Рисунке 2. На графике имеются 95% и 99% границы по Колмогорову-Смирнову для равномерного распределения данных. Этот график позволяет проверить, являются ли значения временного ряда случайными.

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

     Проведем  анализ сезонной декомпозиции.

     Для этого необходимо зайти в меню Seasonal Decomposition и выбрать там:

  1. Once every: 1;
  2. Month (s);
  3. Starting At: 1/1990;
  4. Seasonality: 12.
 

Seasonal Decomposition - bubbly

Data variable: bubbly (MONTHLY CHAMPAGNE SALES) 

Number of observations = 84

Start index = 1.1990

Sampling interval = 1,0 month(s)

Length of seasonality = 12 

Seasonal Decomposition

Method: Multiplicative 

The StatAdvisor

This procedure applies a multiplicative seasonal decomposition to bubbly. The purpose of the decomposition is to separate bubbly into trend-cycle, seasonal, and random components. The data cover 84 time periods.

Each of the tables and graphs shows different aspects of the decomposition.  

Рис. 3

Seasonal Indices for bubbly

Seasonal decomposition method: Multiplicative 

Season Index
1 76,4727
2 70,3034
3 81,8533
4 83,3609
5 87,6277
6 87,0381
7 71,4805
8 40,4459
9 89,8232
10 119,213
11 175,451
12 216,93
 

The StatAdvisor

This table shows the seasonal indices for each month, scaled so that an average month equals 100. The indices range from a low of 40,4459 in month 8 to a high of 216,93 in month 12. This indicates that there is a seasonal swing from 40,4459% of average to 216,93% of average throughout the course of one complete cycle. 

Рис. 4 

     В таблице видны 2 критических месяца, они выделены цветом, которые приходятся на август и декабрь. Нам нужно избавиться от них. Для этого в меню  Graphs  надо выбрать horizontal time sequence plot и integrated periodogram.

     Также в Analitic Options выбираем (Math) выбираем None; Seasonal None и удаляем квадратический тренд.

Рис. 5

     График входит в рамки, это хорошо.

     На рисунке 5 мы видим интегрированную периодограмму скорректированной переменной. Её значения не выходят за 95%-ные границы Колмогорова-Смирнова. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

3. Автоматическое и ручное прогнозирование продаж шампанского 

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

     В меню автоматического прогнозирования ставим Seasonality=24, всё остаьное остаётся, как было.

     В Tables выбираем: Analysis Summary, Forecast Table, Model Comparison, Residual Tests for Randomns. В Graphs: Time Sequence Plot, Forecast Plot Residual Integrated Periodogram 

Получается график: 

 

Рис. 6 
 

На рисунке 6 показан прогноз на 2 года. 

 

Рис. 7 

На рисунке 7 видим, что значения периодограммы  не выходят за 95% границы. Это хорошо. 

     Ручное прогнозирование

 

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

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

Рис. 8

Forecast Table for bubbly

Model: Winter's exp. smoothing with alpha = 0,4159, beta = 0,0182, gamma = 0,4166 (см. таблицу 2. в приложениях) 

    Lower 95,0% Upper 95,0%
Period Forecast Limit Limit
1.1997 3,87751 3,08619 4,66884
2.1997 3,91417 2,97639 4,85195
3.1997 4,50158 3,28857 5,71459
4.1997 4,70563 3,31526 6,09599
5.1997 4,77706 3,25512 6,29899
6.1997 5,28965 3,49404 7,08527
7.1997 4,61259 2,95903 6,26616
8.1997 2,23283 1,3933 3,07235
9.1997 6,24723 3,7971 8,69736
10.1997 8,03609 4,76318 11,309
11.1997 12,0051 6,94641 17,0637
12.1997 14,6836 8,30197 21,0653
1.1998 4,50934 2,49334 6,52534
2.1998 4,54343 2,42279 6,66406
3.1998 5,2157 2,72505 7,70635
4.1998 5,44238 2,78775 8,09701
5.1998 5,51536 2,77137 8,25935
6.1998 6,09679 3,00688 9,18669
7.1998 5,30758 2,57056 8,04459
8.1998 2,56508 1,22056 3,9096
9.1998 7,16545 3,35138 10,9795
10.1998 9,20294 4,23268 14,1732
11.1998 13,7274 6,21102 21,2437
12.1998 16,7653 7,4652 26,0655
 

 

Рис. 9 

Model Comparison

Data variable: bubbly

Number of observations = 84

Start index = 1.1990

Sampling interval = 1,0 month(s)

Length of seasonality = 12 

Models

(A) Random walk with drift = 0,0282736

Seasonal adjustment: Multiplicative

(B) Winter's exp. smoothing with alpha = 0,4159, beta = 0,0182, gamma = 0,4166

(C) Simple exponential smoothing with alpha = 0,2032

Seasonal adjustment: Multiplicative

(D) Brown's linear exp. smoothing with alpha = 0,1022

Seasonal adjustment: Multiplicative

(E) Winter's exp. smoothing with alpha = 0,4159, beta = 0,018, gamma = 0,4165 

Estimation Period

Model RMSE MAE MAPE ME MPE
(A) 0,725943 0,464052 12,1618 0,019122 -1,70134
(B) 0,683328 0,490358 12,8366 -0,169971 -7,74267
(C) 0,663364 0,422294 10,0186 0,129023 0,425901
(D) 0,65552 0,414178 10,0297 0,0736636 -0,477326
(E) 0,683429 0,490643 12,8496 -0,171457 -7,77584
 
Model RMSE RUNS RUNM AUTO MEAN VAR
(A) 0,725943 OK OK ** OK OK
(B) 0,683328 OK OK OK OK **
(C) 0,663364 OK OK OK OK OK
(D) 0,65552 OK OK OK OK OK
(E) 0,683429 OK OK OK OK **
 

     Прогноз  на Рисунке 9 нам полностью  подходит, Это значит, что мы справились  с задачей. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Заключение: 

     Statgraphics Centurion обладает рядом свойств,  которые существенно облегчают  статистическую обработку, анализ и прогнозирование данных. Эта программа применима в различных сферах, не исключая и туризм.

     Существует  огромное количество методов  прогнозирования. Прогнозирование  является очень важным элементом  управления на рынке.

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

Анализ и прогнозирование временных рядов в профессиональном статистическом пакете Statgraphics Centurion