Владимир Брюков - Как предсказать курс доллара. Эффективные методы прогнозирования с использованием Excel и EViews
Однако если посмотреть, как изменялась точность обеих статистических моделей в различные периоды времени, то начиная с 1999 г. стационарная модель дает более точные прогнозы. В частности, в период с января 1999 г. по июнь 2010 г. средняя точность стационарной модели оказалась выше точности нестационарной модели на 0,2 коп. по модулю (см. цифры, выделенные жирным шрифтом в табл. 6.14). А с января 2009 г. по июнь 2010 г. эта разница составила уже 8,7 коп.
Естественно, что и по величине средней точности прогноза (в % по модулю) стационарная модель с января 1999 г. также дает более точные прогнозы. В частности, в период с января 1999 г. по июнь 2010 г. средняя точность стационарной модели (в % по модулю) оказалась выше точности нестационарной модели на 0,08 процентного пункта (см. цифры, выделенные жирным шрифтом в табл. 6.15). В свою очередь с января 2009 г. по июнь 2010 г. эта разница составила уже более 0,26 процентного пункта. С учетом этого можно сделать вывод, что точность стационарной статистической модели за последние 11,5 лет оказалась выше, чем у нестационарной модели.
Воспользовавшись диалоговым мини-окном FORECAST, мы получили не только оценку точности прогноза для стационарной статистической модели log(USDollar) = с + а × log(USDollar(-1)) + МА(1), но и файл с точечными прогнозами USDOLLARF за период с июля 1992 г. по июль 2010 г. Открыв этот файл, мы выяснили, что точечный прогноз на июль 2010 г. оказался равен 31 руб. 19 коп., однако фактический курс доллара в июле 2010 г. был равен 30 руб. 19 коп. Следовательно, разница составила 1 руб. Посмотрим, попал ли фактический курс доллара в диапазон интервального прогноза?
Однако, прежде чем это сделать, проверим остатки, полученные по модели log(USDollar) = с + а × log(USDollar(-l)) + МА(1), на нормальное распределение и на стационарность (см. алгоритм действий № 9).
В первом случае откроем файл RESID и выберем опции VIEW (смотреть)/DESCRIPTIVE STATISTICS (описательная статистика)/ STATS TABLE (таблица со статистикой). При этом следует иметь в виду, что проверку на нормальное распределение остатков целесообразно проводить относительно логарифмических остатков, поскольку наша статистическая модель построена на логарифмическом временном ряде. Логарифмические остатки нетрудно найти, если при составлении прогнозов в диалоговом мини-окне FORECAST (прогноз) поставим «галочку» у файла LOG(USDOLLAR) (см. рис. 6.6). В результате мы получили табл. 6.16.
Судя по тому, что коэффициент асимметрии (Skewness) в табл. 6.16 положителен, можно прийти к выводу, что в распределении остатков, полученных по стационарной модели, наблюдается положительная асимметрия. Отсюда можно сделать вывод, что в динамике курса доллара к рублю чаще наблюдались резкие (вполне очевидно, что незначительные плавные колебания курса легко поддаются прогнозированию) подъемы, чем аналогичные падения. В свою очередь величина коэффициент эксцесса (Kurtosis) существенно выше 3, что свидетельствует об «островершинном» распределении остатков. По сути, это означает, что в этом распределении имеется ярко выраженное ядро плотности распределения, внутри которого диапазон колебаний величины остатков незначителен, и рассеянное «гало», где разброс колебаний величины остатков весьма значителен. Поскольку величина тестовой статистики Жарка — Бера составила 11990,08, а уровень ее значимости (Probability) оказался равен нулю, то, следовательно, мы вынуждены отвергнуть гипотезу о нормальном распределении остатков. Поскольку, как мы уже говорили нашим читателям, при уровне значимости критерия Жарка — Бера (Probability) меньше 0,05 нулевая гипотеза о нормальном распределении отклоняется.
В EViews есть возможность посмотреть в графическом виде оценку ядра плотности распределения с помощью опций DISTRIBUTION/ KERNEL DENSITY GRAPHS… (распределение/графики ядра плотности распределения). В появившемся мини-окне KERNEL DENSITY (ядро плотности распределения) установим опцию EPANECHNICOV. В результате получим рис. 6.11, наглядно показывающий «островершинный» характер распределения остатков с правосторонней асимметрией.
Для проверки остатков на стационарность воспользуемся расширенным тестом Дикки — Фуллера. Следует заметить, что проверку остатков на стационарность также целесообразно проводить относительно логарифмических остатков. После проведения тестирования мы убедились, что получили стационарные остатки (табл. 6.17). Поскольку статистика теста Дикки — Фуллера в этом случае составила -15,61466, а ее значимость (Probability) равна 0,0000, то нулевая гипотеза о том, что D(RESID) имеет единичный корень, отвергается. Следовательно, мы можем принять альтернативную гипотезу о стационарности полученных остатков.
Несмотря на то что остатки, полученные по стационарной модели, нельзя считать нормально распределенными, мы тем не менее уже знаем, что при больших выборках можно строить интервальные прогнозы исходя из их нормального распределения. Поэтому нашим следующим шагом будет расчет интервальных прогнозов не только на июль 2010 г. (курс доллара по этому месяцу не включен в базу данных), но и для всех наблюдений, на основе которых составлена статистическая модель log(USDollar) = с + а × log(USDollar(-1)) + МА(1). Это поможет нам проверить соответствие составленных интервальных прогнозов нормальному распределению, поскольку уровень надежности для интервальных прогнозов рассчитывался исходя из предположения о нормальном распределении остатков. Попутно заметим, что интервальные прогнозы будут построены начиная с июля 1992 г., поскольку первое наблюдение во временном ряде нам потребовалось для создания факторной переменной log(USDollar(-l).
После того как на основе алгоритма действий № 12 будут составлены соответствующие интервальные прогнозы, у нас появится возможность сопоставить заданные уровни надежности с фактической долей точных интервальных прогнозов. Судя по табл. 6.18, доля точных прогнозов оказалась незначительно ниже заданного уровня надежности при 99,9 %-ном и 99 %-ном уровнях надежности и практически ему равна при 95 %-ном (отклонение на -0,1 процентного пункта). Однако при 90 %-ном уровне надежности и более низких уровнях доля фактических прогнозов становится выше заданного уровня. Причем эта положительная разница растет при снижении уровня надежности, достигая своего максимума при 40 %-ном уровне надежности, когда она равна 36 процентным пунктам.
Однако если сравнить табл. 6.18, характеризующую точность интервальных прогнозов для стационарной модели log(USDollar) = с + а × log(USDollar(-l)) + МА(1), с табл. 4.8, которая характеризует их точность для нестационарной модели USDOLLAR= а х USDOLLAR(-l) + а × USDOLLAR(-2), то выяснится, что последняя модель с этой точки зрения несколько точнее.