Эконометрическое прогнозирование по временным рядам
 
Предыдущий раздел Меню навигации Следующий раздел

Экспоненциальная модель

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

где b0 – множитель, а b1 – оценка годовой нормы прироста показателя. Эта модель существенно нелинейна, но может быть сведена к линейной посредством логарифмирования.

Теперь зависимость логарифма значения временного ряда от момента наблюдения становится линейной, так что снова можно применить простейший вариант регрессионного анализа. При этом предварительно следует прологарифмировать значения Y (будет ли это десятичный, или натуральный логарифм – безразлично; мы будем использовать натуральные логарифмы) и учесть, что после построения регрессии константа уравнения в действительности будет равна логарифму множителя, коэффициент пропорциональности соответственно логарифму нормы прироста показателя. Прогноз по модели можно будет выполнить аналогично линейному случаю, но важно, что после получения численного значения придется взять его экспоненту, чтобы перейти к тем же натуральным показателям, что и в значениях исходного ряда. Соответственно ошибки прогноза будут рассчитаны не для самого показателя, а для его логарифма, и для получения интервальной оценки также потребуется потенцирование. Тем не менее все эти операции не представляют особой сложности в электронных таблицах.

Вы можете сразу посмотреть итоговые результаты экспоненциального моделирования объема продаж фирмы Kodak и выводы, либо загрузить таблицу в формате Quattro или в формате Excel и выполнить упражнение, а затем сравнить свои результаты с представленными ниже.

  • Расположите окна с табличным процессором и браузером так, чтобы Вам было удобно с ними работать. Если разрешение монитора невелико, просто переключайтесь между окнами.
  • В колонке C пронумеруйте все имеющиеся значения отклика – от 1 до 27 (последние 4 значения будут использоваться для прогнозирования).
  • Заполните значениями колонку натуральных логарифмов (колонка D). Для этого введите в D7 формулу @ln(B7) при работе в Quattro или =ln(B7) при работе в Excel и прокопируйте вниз до 29-й строки таблицы включительно.
  • Вызовите функцию построения регрессии.
  • В поле независимых переменных введите блок ячеек с номерами значений от C7 до C29.
  • В поле зависимых переменных введите блок ячеек с логарифмами значений от D7 до D29.
  • Определите блок вывода результатов B36..D38 (в загруженной таблице этот блок выделен линией сверху и слева).
  • Нажмите клавишу Ok для проведения расчетов. Результаты будут выведены в указанный Вами блок.
  • В ячейку F7 введите формулу для расчета по полученным коэффициентам регрессии и скопируйте ее вниз по колонке вплоть до ячейки F33, при этом четыре последних значения этой колонки дают логарифм прогноза на 4 периода вперед.
  • В ячейку G7 введите формулу @exp(F7) для Quattro Pro или =exp(F7) для Excel и скопируйте ее вниз по колонке вплоть до ячейки G33.
  • Вызовите функцию построения графика.
  • В качестве оси X укажите года, выделив все ячейки от A7 до A33.
  • В качестве первой серии укажите все ячейки от B7 до B29.
  • В качестве второй серии укажите все ячейки от G7 до G33.
  • При необходимости измените цвета кривых и отмените вывод маркеров на графике.
  • Сравните вид Вашего графика и Ваши данные с тем, что представлено на странице пособия. Можно при этом округлить рассчитанные значения так же, как это сделано в пособии.
  • Сохраните таблицу с расчетами в нужном Вам каталоге или закройте систему без сохранения результатов.

Результаты расчета коэффициентов регрессии

Результат регрессионного анализа
Константа 1.061
Оценка стандартной ошибки Y 0.109
Коэффициент вариации 0.970
Число наблюдений 23
Степени свободы 21
Коэффициент(ы) X 0.090
Стандартная ошибка коэффициента 0.003
Критерий Стьюдента 26.048

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

Наблюдаемые и рассчитанные значения отклика

Год Измеренный объем выпуска Логарифм объема выпуска Значение X Логарифм регрессионной оценки Рассчитанное значение
1970 2.8 1.030 1 1.150 3.159
1971 3.0 1.099 2 1.240 3.455
1972 3.5 1.253 3 1.329 3.778
1973 4.0 1.386 4 7.419 4.132
1974 4.6 1.526 5 1.508 4.520
1975 5.0 1.609 6 1.598 4.943
1976 5.4 1.686 7 1.688 5.406
1977 6.0 1.792 8 1.777 5.913
1978 7.0 1.946 9 1.867 6.467
1979 8.0 2.079 10 1.956 7.073
1980 9.7 2.272 11 2.046 7.736
1981 10.3 2.332 12 2.135 8.461
1982 10.8 2.380 13 2.225 9.254
1983 10.2 2.322 14 2.315 10.121
1984 10.6 2.361 15 2.404 11.069
1985 10.6 2.361 16 2.494 12.106
1986 11.5 2.442 17 2.583 13.241
1987 13.3 2.588 18 2.673 14.482
1988 17.0 2.833 19 2.762 15.839
1989 18.4 2.912 20 2.852 17.323
1990 18.9 2.939 21 2.942 18.946
1991 19.4 2.965 22 3.031 20.722
1992 20.1 3.001 23 3.121 22.663
Прогнозные значения
1993 24 3.210 24.787
1994 25 3.300 27.110
1995 26 3.389 29.650
1996 27 3.479 32.429

 

В данном случае модель предсказывает объем чистых продаж на 1993 год на уровне 24.8 миллиарда долларов, что заметно выше прогноза по прочим моделям. Качество экспоненциальной модели можно грубо оценить по коэффициенту корреляции (0.97 - очень высокое значение, хотя и сравнимое со случаем квадратичной модели). Но прямое сравнение этих двух моделей по критерию Фишера непосредственно из результатов регрессионного анализа невозможно, поскольку в квадратичной модели отклик задан в натуральных величинах, а в экспоненциальной - в их логарифмах, и стандартная ошибка модели рассчитана именно для логарифма отклика. И тем не менее дисперсионное сопоставление выполнить можно: следует только рассчитать дефект модели в каждой точке ряда, возвести разности в квадрат, просуммировать по всему ряду и поделить на число степеней свободы – получится значение дисперсии адекватности модели именно для натуральных показателей ряда. В нашем случае дисперсия адекватности экспоненциальной модели выше дисперсии для квадратичной, но ниже, чем для линейной, хотя все эти модели неразличимы с точки зрения критерия Фишера при доверительной вероятности 95%. Что касается применения критерия Стьюдента к коэффициенту экспоненциальной модели, то он показывает отчетливую значимость при всех разумных доверительных вероятностях.

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

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


Предыдущий раздел Меню навигации Следующий раздел


Created 17.08.2007, Revised 13.04.2008
Используются технологии uCoz