Выполняя исследование временных рядов в R, я обнаружил, что он arima
предоставляет только значения коэффициентов и их стандартные ошибки для подобранной модели. Тем не менее, я также хочу получить р-значение коэффициентов.
Я не нашел ни одной функции, обеспечивающей значение coef.
Поэтому я хочу рассчитать его самостоятельно, но я не знаю степень свободы в распределении коэффициентов t или chisq. Поэтому мой вопрос заключается в том, как получить p-значения для коэффициентов модели подогнанной аримы в R?
Ответы:
«Значение t» - это отношение коэффициента к стандартной ошибке. Степени свободы (ndf) - это число наблюдений минус максимальный порядок различий в модели минус количество оценочных коэффициентов. Значение «F» будет квадратом значения «t». Для точного вычисления вероятности вам нужно будет вызвать нецентральную функцию хи-квадрат и передать значение F и степени свободы (1, ndf). или, возможно, просто вызвать поиск F-функции.
источник
Поскольку
arima
для оценки используется максимальная вероятность, коэффициенты асимптотически нормальны. Следовательно, разделите коэффициенты на их стандартные ошибки, чтобы получить z-статистику, а затем рассчитать p-значения. Вот пример с в R с первым примером соarima
страницы справки:Последняя строка дает p-значения.
источник
Вы также можете использовать
coeftest
изlmtest
пакета:источник