Методические указания по лабораторным работам По дисциплине (стр. 4 из 6)
Рисунок 6 – Входящая панель процедуры сглаживания
Система выведет в рабочей области сводку простого пятиточного сглаживания, установленного по умолчанию. В табличных опциях устанавливаем флажок в поле Date Table (таблица данных), а в графических – Time Sequence Plot (график временной последовательности) (рис. 7, 8, 9, 10). Получим отчет (рис. 11).
Рис. 7
Рис. 8
Рис. 9
Рис. 10
Щелкнем правой кнопкой мыши на второй табличной панели – появится всплывающее меню, в котором выберем пункт Pane Options (опции панели). Система предоставит возможность изменить метод сглаживания. Установим переключатель в поле EWMA (рис. 11). Система рассчитает значения и построит график для взвешенного экспоненциального сглаживания с параметром 0,1.
Рис. 11 – Сглаживание уровней рядов динамики методом временной последовательности
Проведем аналогичные расчеты с помощью устойчивого нелинейного сглаживания.
Анализ графиков позволяет сделать вывод, что изменение темпов прироста выработки лучше всего аппроксимирует нелинейное устойчивое сглаживание. Кривые простого скользящего и взвешенного экспоненциального сглаживания менее информативны.
Рис. 12 – Панель установки опций сглаживания
Определив общие закономерности изменения выработки, можно приступить к подбору модели и расчету прогнозных значений моделируемого показателя. С этой целью воспользуемся процедурой Forecasting.
Свернем в пинктограмму окно анализа с результатами сглаживания и выберем указанную процедуру. Заполним верхнюю панель, введя имя переменной и количество лет прогноза, равное пяти (рис. 3). С целью проверки адекватности модели используем три последних наблюдения. Поэтому в поле Withhold for Validation (число точек для проверки правильности модели) введем цифру три.
Система поместит в рабочую область сводку прогноза по модели случайной выборки.
Учитывая, что STATGRAPHICS может сравнивать одновременно пять типов моделей, оптимизируя их параметры, выберем для анализа линейную модель (Liner trend), параболу (Quadratic trend), линейное экспоненциальное сглаживание Брауна (Brown/s linear exp. Smoothing), линейное экспоненциальное сглаживание Хольта (Holt/s linear exp. Smoothing), квадратическое экспоненциальное сглаживание Брауна (Brown/s quadratic exp. Smoothing). Напомним, что в STATGRAPHICS реализовано три типа сглаживания Брауна. Простое сглаживание основано на предположении стационарности изучаемого процесса, линейное предполагает линейный тренд в данных. Квадратическое сглаживание базируется на том, что моделируемый показатель может быть основан на том, что моделируемый показатель может быть описан полиномом второго порядка, т.е. параболой.
Указанные модели выбираются с помощью панели Model Specification Options (опции спецификации модели) (рис. 4). В области Model щелкнем на пункт А, а в области Type установим флажок Linear trend. Затем выберем пункт В и установим флажок Quadratic trend. Для модели С выберем Brown/s linear exp. Smoothing; для моделей D и E установим флажки на полях Holt/s linear exp. Smoothing и Brown/s quadratic exp. Smoothing. Остальные поля оставим со значениями по умолчанию, в т.ч. с флажком Optimize, активным для всех экспоненциальных моделей. Щелкаем ОК с установленным для модели С переключателем. Это означает, что все расчеты система выполнит для этой модели (рис. 13).
Рис. 13 – Анализ линейной модели Брауна
Результаты сравнительного анализа можно вывести, установив флажок Model Comparisons (сравнение моделей) табличных опций.
Получим листинг сравнения моделей (рис. 13). Кратко опишем его. В верхней части листинга приводится информация о данных и уравнениях или коэффициентах построенных моделей.
Наибольший интерес представляют таблицы со статиками прогнозирования, позволяющими оценить адекватность полученных зависимостей. К таким характеристикам относится средняя арифметическая ошибка (МЕ), описывающая отклонения фактических значений от выровненных. Чем ближе она к нулю, тем точнее осуществлена аппроксимация. Средняя квадратическая ошибка (MSE) и средняя абсолютная ошибка (МАЕ) используются для сравнения разных процедур прогнозирования.
Среднепроцентная ошибка (МРЕ) и среднеабсолютная процентная ошибка (МАРЕ) рассчитываются по остаткам одношагового выравнивания, которые делятся на фактическое значение выработки.
Листинг содержит также стандартную ошибку остатков (RMSE) и пять тестов RUNS, RUNM, AUTO, MEAN, VAR:
RUNS – тест на чрезмерное количество пиков и впадин (Test for excessive runs up and down) – рассчитывает число повышений или падений в последовательности анализируемых данных. Тест чувствителен к долгосрочным циклам;
RUNM – тест на чрезмерное количество отклонений от медианы (Test for excessive runs above and below median) – рассчитывает число наблюдений, значение которых выше или ниже медианы, значение которых выше или ниже медианы, и игнорируют значения, которые являются равными медиане. Тест чувствителен к наличию тренда в данных;
AUTO – тест на чрезмерную автокорреляцию (Box-Pierce test for excessive autocorrelation) – рассчитывает коэффициент сериальной корреляции Бокса-Пирса;
MEAN – тест на существенность разности средних (test for difference in mean 1st half to 2nd half) – служит для определения тенденции среднего значения;
VAR – тест на существенность разности дисперсий (test for difference in variance 1st half to 2nd half) – позволяет установить тенденцию вариабельности.
Модель прогнозирования будет адекватной, если все тесты будут иметь значение OK. Т.о., тесты остатков несущественны. Знак звездочки означает, что тест статистически существенен. Количество звездочек определяет уровень существенности критерия. Три звездочки означают, что тест значим с вероятностью, превышающей 99%.
Данные листинга, приведенного на рисунке 14, показывают, что линейная экспоненциальная модель Брауна наиболее удачно аппроксимирует фактические данные. Поэтому для расчетов используется эта модель. Текстовые и графические результаты прогнозирования можно вывести на экран, установив, например, флажки на полях Forecast Table (таблица прогнозов) табличных опций и Time Sequence Plot (график временной последовательности).
Рис. 14 – Сравнительный анализ моделей
На рисунке 15 показаны фактические и прогнозные значения темпов прироста выработки, полученные по линейной экспоненциальной модели Брауна (a=0,2145). Рассчитаем прогнозы по другим моделям и результаты сведем в таблицу 3.
Сравнивая различные варианты прогнозов на 5 лет, следует отметить неоднозначность полученных результатов (таблица 7).
Рис. 15 – Прогноз по линейной экспоненциальной модели Брауна
Таблица 7 – Прогноз темпов прироста производительности труда по различным моделям и методам
| Тип модели | Порядковый номер прогноза | ||||
| 28 | 29 | 30 | 31 | 32 | |
| 7,86123-0,112565*t | 4,71 | 4,6 | 4,48 | 4,37 | 4,25 |
| 9.67752-0.531709*t+0.016766*t2 | 7,93 | 8,36 | 8,81 | 9,31 | 9,83 |
| Brown/s linear wish a=0,2159 | 4,8 | 4,68 | 4,56 | 4,45 | 4,33 |
| Holt/s linear wish a=0,3144, b=0,076 | 4,57 | 4,44 | 4,32 | 4,19 | 4,07 |
| Brown/s quadratic wish a=0,1652 | 4,69 | 4,53 | 4,35 | 4,17 | 3,99 |
В качестве пессимистического варианта прогноза можно рассматривать результаты по квадратической экспоненциальной модели. При этом надо иметь в виду, что остатки у этой модели по тесту MEAN незначительно существенны.