Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

 
Добавить ответ в эту темуОткрыть тему
> Доверительный или прогнозный интервал, для данных из временного ряда
nokh
сообщение 17.06.2022 - 08:19
Сообщение #1





Группа: Пользователи
Сообщений: 1187
Регистрация: 13.01.2008
Из: Челябинск
Пользователь №: 4704



Уважаемые участники форума. Ситуация такая: есть многолетний временной ряд, в котором сильная сезонная составляющая. Её хорошо видно при усреднении данных по месяцам. Хочется к этим средним добавить интервал, возможно прогнозный даже лучше, пусть он и шире. Поскольку данные не независимые, это нужно сделать грамотно, а то, что я смог быстро нагуглить - доверительные и прогнозные интервалы для регрессии и/или ряда целиком. Для этого такую зависимость нужно сначала отмоделировать, а классическими методами для этого я не владею, да и пока задача другая: просто снабдить многолетние средние интервалом. Ситуация осложняется тем, что если смотреть тренды для каждого месяца отдельно, то в августе и сентябре по критерию Манна - Кендалла на тренд таковой статистически значим, хотя снижение совсем небольшое. Как посоветуете представить такие данные? Данные гидрохимические, но обобщается и на медицинские, и на ветеринарные временные серии (недавно смотрел посуточные индивидуальные удои на автоматическом комплексе, где коровы сами подходят доиться когда захотят, так там тоже временной ряд получается...).
Ещё вопрос: насколько корректно смотреть PCA для такого массива данных: 25 строк (годы) х 12 столбцов (месяцы). Можно, что-то почитать умное. Спасибо!
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
100$
сообщение 17.06.2022 - 12:53
Сообщение #2





Группа: Пользователи
Сообщений: 830
Регистрация: 23.08.2010
Пользователь №: 22694



Nokh, для временных рядов нет ДИ, есть ДИ для прогноза. Они зависят от применяемой модели (линейная регрессия, ARIMA, GARCH, etc.). При этом две последние модели применяются к стационарным временным рядам. Для этого необходимо посмотреть на коррелограмму временного ряда и провести формальный тест на стационарность (классика: тесты Дики - Фуллера, Квятковского-Филлипса-Шмидта-Шина (KPSS), Лейбурна-МакКейба, и т.д. Мильёны их).

По поводу теста Манна - Кендалла замечу, что тестировать им - это все равно что в 2022 г. добывать огонь трением.

Если в ряду есть сезонная волна - 12 лаг коррелограммы будет статистически значим.

По этому временному ряду невозбранно проползти "гусеницей" (сингулярным спектральным анализом). У Игоря в мануале к Аттестату все написано, со ссылками на хрестоматийные публикации Голяндиной с соавт.

А вообще, взглянуть бы на этот ряд хоть одним глазком.

Сообщение отредактировал 100$ - 17.06.2022 - 12:57
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
passant
сообщение 17.06.2022 - 13:45
Сообщение #3





Группа: Пользователи
Сообщений: 209
Регистрация: 27.04.2016
Пользователь №: 28223



Небольшое уточнение.
Цитата(100$ @ 17.06.2022 - 12:53) *
Nokh, для временных рядов нет ДИ, есть ДИ для прогноза. Они зависят от применяемой модели (линейная регрессия, ARIMA, GARCH, etc.). При этом две последние модели применяются к стационарным временным рядам.

Вообще-то ARIMA - есть адаптация метода ARMA именно для случая, когда ряд изначально нестационарен, но его удается привести к стационарному виду взятием разностей некоторого порядка. А если речь идет - как у ТС - "ряд, в котором сильная сезонная составляющая. Её хорошо видно при усреднении данных по месяцам". - то в таком случае есть еще одно расширение модели, а именно SARIMA, которая именно для таких задач и ситуаций и предназначена. Причем часто работает вполне удовлетворительно.

Сообщение отредактировал passant - 17.06.2022 - 13:52
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
100$
сообщение 17.06.2022 - 14:19
Сообщение #4





Группа: Пользователи
Сообщений: 830
Регистрация: 23.08.2010
Пользователь №: 22694



Цитата(passant @ 17.06.2022 - 13:45) *
Небольшое уточнение.


Это, конечно, так, только я имел в виду уже непосредственное оценивание модели. Для которого в случае нестационарности ряда остационаривание является обязательной операцией, предшествующей собственно оцениванию.

Про SARIM'у вы хорошо вспомнили. Хотя сезонная волна к стационарности не имеет отношения и ее можно сгладить и по Фурье.

Сообщение отредактировал 100$ - 17.06.2022 - 14:28
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Добавить ответ в эту темуОткрыть тему