Страница 1 из 1
Временной ряд интегральной характеристики с пропуском...
Добавлено: 21 ноя 2015, 13:14
nickleb
Здравствуйте. Есть временной ряд некой инертной расчётной среднемноголетней интегральной характеристики (скажем, теплозапас большой толщи океанских вод) с пропуском в несколько значений подряд идущих лет... Скажем так:
1-ый интервал) за 1948...1995 гг. обеспечен ряд таковыми значениями;
2-ой интервал) за 1996...2005 гг. ПРОПУСК;
3-ий интервал) за 2006...2015 гг. обеспечен ряд таковыми значениями...
С помощью каких R-package'-ей попробовать "поиграться" с восстановлением значений во 2-ом интервале таковой интегральной расчётной характеристики?
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 21 ноя 2015, 16:39
sergsh
Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.
Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 21 ноя 2015, 16:48
nickleb
sergsh писал(а):Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.
Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
sergsh, спасибо! - да, конечно, - это в планах - хоть и мало измеренного, но пытаемся... а пока вот так, - как я спросил
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 21 ноя 2015, 21:49
gamm
у вас там наверняка дикая инерционность, так что МСМС ваше все
Из внешних переменных я бы в первую очередь посмотрел на SST, его, насколько я понимаю, считали еще до Терры/Аквы по метеоспутникам, а разрешения в несколько км вам хватит. Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты, наверняка и другие пакеты есть, тот же mice (я с временными рядами не очень дружу).
Т.е. вы строите случайные процесс, например типа RW1/2, с внешними переменными (возможно с лагами - нужно смотреть реальную задачу), и получаете распределение для пропущенных данных с доверительным интервалом. Можно для начала прописать такой процесс в среде WinBUS/JAGS и посмотреть, там вообще все просто, особенно если моделировать только отклонение от регрессии.
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 23 ноя 2015, 08:48
Игорь Черниенко
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 23 ноя 2015, 11:31
gamm
это и есть наследие WinBUGS (который нынче OpenBUGS)
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 23 ноя 2015, 12:35
nickleb
gamm, Игорь, спасибо!
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 23 ноя 2015, 15:52
ymr3R9Jge
Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете
zoo есть функция
na.approx для восстановления пропущенных значений во временных рядах.
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 23 ноя 2015, 18:37
nickleb
ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете
zoo есть функция
na.approx для восстановления пропущенных значений во временных рядах.
ymr3R9Jge, обязательно гляну... спасибо!
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 02 дек 2015, 13:59
ymr3R9Jge
ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете
zoo есть функция
na.approx для восстановления пропущенных значений во временных рядах.
В
статье на сайте описан пример интерполяции временного ряда.
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 02 дек 2015, 19:57
nickleb
ymr3R9Jge писал(а):ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете
zoo есть функция
na.approx для восстановления пропущенных значений во временных рядах.
В
статье на сайте описан пример интерполяции временного ряда.
Замечательно! Спасибо! Теперь понятно связь этой статьи с вопросом Вашим:
viewtopic.php?f=31&t=19602 Что-то так и не собрался в своё время ответить - извините - у самого сразу не получилось.. Ещё раз спасибо!
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 02 дек 2015, 22:10
ymr3R9Jge
Да, вопрос задавал, когда только приступал к "разбору полётов". Но на функцию na.approx сам вышел чисто случайно! До этого долго не удавалось найти удобный инструмент для подобных задач. Первые наши расчёты загрязнений были не очень хорошими, поскольку пришлось выкинуть треть значений. Но теперь с интерполяцией - можно два временных ряда соединять "ноздря в ноздрю".
Re: Временной ряд интегральной характеристики с пропуском...
Добавлено: 20 сен 2016, 13:41
nickleb
gamm писал(а):
... Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты...