Критерий KPSS
Материал из MachineLearning.
(Различия между версиями)
м (ссылки) |
м (оформление) |
||
Строка 30: | Строка 30: | ||
== Реализации == | == Реализации == | ||
- | |||
* MATLAB: В версии 2013b и выше встроен пакет методов [http://www.mathworks.com/help/econ/index.html Econometrics Toolbox], в котором реализована функция [h,pValue] = kpsstest(___) <ref name="kpsstestmatlab"> [http://www.mathworks.com/help/econ/kpsstest.html KPSS test for MATLAB]</ref>. | * MATLAB: В версии 2013b и выше встроен пакет методов [http://www.mathworks.com/help/econ/index.html Econometrics Toolbox], в котором реализована функция [h,pValue] = kpsstest(___) <ref name="kpsstestmatlab"> [http://www.mathworks.com/help/econ/kpsstest.html KPSS test for MATLAB]</ref>. | ||
- | |||
* R: в пакете [http://cran.r-project.org/web/packages/tseries/index.html tseries] реализован метод для вычисления критерия KPSS kpss.test(x) <ref name="kpsstestR"> [http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/library/tseries/html/kpss.test.html KPSS test for R]</ref>. | * R: в пакете [http://cran.r-project.org/web/packages/tseries/index.html tseries] реализован метод для вычисления критерия KPSS kpss.test(x) <ref name="kpsstestR"> [http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/library/tseries/html/kpss.test.html KPSS test for R]</ref>. | ||
- | |||
== Пример использования == | == Пример использования == |
Текущая версия
Критерий KPSS (KPSS test) — критерий, используемый для проверки на стационарность наблюдаемого временного ряда.
Критерий назван по первым буквам ученых Квятковский-Филлипс-Шмидт-Шин (Kwiatkowski–Phillips–Schmidt–Shin), которые ввели его в 1992 году. [1]
Содержание |
Определение
Если рассматриваемый ряд имеет вид:
где
- — коэффициент тренда,
- — некоторый стационарный процесс,
- — некоторый независимый и одинаково распределенный с процесс с математическим ожиданием 0 и дисперсией .
Выдвигаются две конкурирующие гипотезы:
- : временной ряд являются стационарным (или, аналогично ),
- : временной ряд не являются стационарным ().
Вычисляем статистику:
- ,
где
- — размер выборки,
- — стандартная ошибка в форме Ньюи-Уеста (Newey–West estimate) [1]
Реализации
- MATLAB: В версии 2013b и выше встроен пакет методов Econometrics Toolbox, в котором реализована функция [h,pValue] = kpsstest(___) [1].
- R: в пакете tseries реализован метод для вычисления критерия KPSS kpss.test(x) [1].
Пример использования
- a = 1:100;
- b = normrnd(50, 20, 100, 1);
- [~,pValuea] = kpsstest(a);
- [~,pValueb] = kpsstest(b);
Полученные значения p-value 0.1 и 0.001 соответственно, то есть гипотеза о стационарности в первом случае отклоняется, во втором - нет.
Ссылки
- Hamilton, J. D. Time Series Analysis. Princeton, NJ: Princeton University Press, 1994.