Временной ряд интегральной характеристики с пропуском...
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Временной ряд интегральной характеристики с пропуском...
Здравствуйте. Есть временной ряд некой инертной расчётной среднемноголетней интегральной характеристики (скажем, теплозапас большой толщи океанских вод) с пропуском в несколько значений подряд идущих лет... Скажем так:
1-ый интервал) за 1948...1995 гг. обеспечен ряд таковыми значениями;
2-ой интервал) за 1996...2005 гг. ПРОПУСК;
3-ий интервал) за 2006...2015 гг. обеспечен ряд таковыми значениями...
С помощью каких R-package'-ей попробовать "поиграться" с восстановлением значений во 2-ом интервале таковой интегральной расчётной характеристики?
1-ый интервал) за 1948...1995 гг. обеспечен ряд таковыми значениями;
2-ой интервал) за 1996...2005 гг. ПРОПУСК;
3-ий интервал) за 2006...2015 гг. обеспечен ряд таковыми значениями...
С помощью каких R-package'-ей попробовать "поиграться" с восстановлением значений во 2-ом интервале таковой интегральной расчётной характеристики?
-
- Активный участник
- Сообщения: 205
- Зарегистрирован: 20 фев 2013, 21:48
- Репутация: 30
Re: Временной ряд интегральной характеристики с пропуском...
Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.
Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
А потом по найденной зависимости восстановить пропуски.
Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Временной ряд интегральной характеристики с пропуском...
sergsh, спасибо! - да, конечно, - это в планах - хоть и мало измеренного, но пытаемся... а пока вот так, - как я спросилsergsh писал(а):Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.
Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
-
- Гуру
- Сообщения: 4056
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1054
- Ваше звание: программист
- Откуда: Казань
Re: Временной ряд интегральной характеристики с пропуском...
у вас там наверняка дикая инерционность, так что МСМС ваше все
Из внешних переменных я бы в первую очередь посмотрел на SST, его, насколько я понимаю, считали еще до Терры/Аквы по метеоспутникам, а разрешения в несколько км вам хватит. Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты, наверняка и другие пакеты есть, тот же mice (я с временными рядами не очень дружу).
Т.е. вы строите случайные процесс, например типа RW1/2, с внешними переменными (возможно с лагами - нужно смотреть реальную задачу), и получаете распределение для пропущенных данных с доверительным интервалом. Можно для начала прописать такой процесс в среде WinBUS/JAGS и посмотреть, там вообще все просто, особенно если моделировать только отклонение от регрессии.
Из внешних переменных я бы в первую очередь посмотрел на SST, его, насколько я понимаю, считали еще до Терры/Аквы по метеоспутникам, а разрешения в несколько км вам хватит. Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты, наверняка и другие пакеты есть, тот же mice (я с временными рядами не очень дружу).
Т.е. вы строите случайные процесс, например типа RW1/2, с внешними переменными (возможно с лагами - нужно смотреть реальную задачу), и получаете распределение для пропущенных данных с доверительным интервалом. Можно для начала прописать такой процесс в среде WinBUS/JAGS и посмотреть, там вообще все просто, особенно если моделировать только отклонение от регрессии.
- Игорь Черниенко
- Активный участник
- Сообщения: 137
- Зарегистрирован: 28 мар 2009, 01:05
- Репутация: 11
- Откуда: Хабаровск, Южно-Сахалинск
-
- Гуру
- Сообщения: 4056
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1054
- Ваше звание: программист
- Откуда: Казань
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Временной ряд интегральной характеристики с пропуском...
gamm, Игорь, спасибо!
-
- Активный участник
- Сообщения: 117
- Зарегистрирован: 31 окт 2011, 00:18
- Репутация: 14
- Откуда: Кривий Ріг
Re: Временной ряд интегральной характеристики с пропуском...
Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Временной ряд интегральной характеристики с пропуском...
ymr3R9Jge, обязательно гляну... спасибо!ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.
-
- Активный участник
- Сообщения: 117
- Зарегистрирован: 31 окт 2011, 00:18
- Репутация: 14
- Откуда: Кривий Ріг
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Временной ряд интегральной характеристики с пропуском...
Замечательно! Спасибо! Теперь понятно связь этой статьи с вопросом Вашим: viewtopic.php?f=31&t=19602 Что-то так и не собрался в своё время ответить - извините - у самого сразу не получилось.. Ещё раз спасибо!
-
- Активный участник
- Сообщения: 117
- Зарегистрирован: 31 окт 2011, 00:18
- Репутация: 14
- Откуда: Кривий Ріг
Re: Временной ряд интегральной характеристики с пропуском...
Да, вопрос задавал, когда только приступал к "разбору полётов". Но на функцию na.approx сам вышел чисто случайно! До этого долго не удавалось найти удобный инструмент для подобных задач. Первые наши расчёты загрязнений были не очень хорошими, поскольку пришлось выкинуть треть значений. Но теперь с интерполяцией - можно два временных ряда соединять "ноздря в ноздрю".
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Временной ряд интегральной характеристики с пропуском...
gamm писал(а): ... Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты...
- Вложения
-
- v63i20.pdf
- Статья об R-INLA
- (1.5 МБ) 471 скачивание
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 12 гостей