Временной ряд интегральной характеристики с пропуском...

Вопросы по статистическому пакету R. Не обязательно гео.
Ответить
nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 21 ноя 2015, 13:14

Здравствуйте. Есть временной ряд некой инертной расчётной среднемноголетней интегральной характеристики (скажем, теплозапас большой толщи океанских вод) с пропуском в несколько значений подряд идущих лет... Скажем так:
1-ый интервал) за 1948...1995 гг. обеспечен ряд таковыми значениями;
2-ой интервал) за 1996...2005 гг. ПРОПУСК;
3-ий интервал) за 2006...2015 гг. обеспечен ряд таковыми значениями...
С помощью каких R-package'-ей попробовать "поиграться" с восстановлением значений во 2-ом интервале таковой интегральной расчётной характеристики?

sergsh
Активный участник
Сообщения: 205
Зарегистрирован: 20 фев 2013, 21:48
Репутация: 30

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение sergsh » 21 ноя 2015, 16:39

Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.

Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 21 ноя 2015, 16:48

sergsh писал(а):Если это реально теплозапас, то сначала нужно бы найти зависимость теплозапаса от температуры и количества солнечных дней, и других факторов, влияющих на теплозапас.
А потом по найденной зависимости восстановить пропуски.

Иначе, опираясь только на сам ряд значений, можно сделать грубую ошибку
sergsh, спасибо! - да, конечно, - это в планах - хоть и мало измеренного, но пытаемся... а пока вот так, - как я спросил

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение gamm » 21 ноя 2015, 21:49

у вас там наверняка дикая инерционность, так что МСМС ваше все :-)

Из внешних переменных я бы в первую очередь посмотрел на SST, его, насколько я понимаю, считали еще до Терры/Аквы по метеоспутникам, а разрешения в несколько км вам хватит. Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты, наверняка и другие пакеты есть, тот же mice (я с временными рядами не очень дружу).

Т.е. вы строите случайные процесс, например типа RW1/2, с внешними переменными (возможно с лагами - нужно смотреть реальную задачу), и получаете распределение для пропущенных данных с доверительным интервалом. Можно для начала прописать такой процесс в среде WinBUS/JAGS и посмотреть, там вообще все просто, особенно если моделировать только отклонение от регрессии.

Аватара пользователя
Игорь Черниенко
Активный участник
Сообщения: 137
Зарегистрирован: 28 мар 2009, 01:05
Репутация: 11
Откуда: Хабаровск, Южно-Сахалинск

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение Игорь Черниенко » 23 ноя 2015, 08:48

JAGS+rjags еще

gamm
Гуру
Сообщения: 4046
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1050
Ваше звание: программист
Откуда: Казань

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение gamm » 23 ноя 2015, 11:31

Игорь Черниенко писал(а):JAGS+rjags еще
это и есть наследие WinBUGS (который нынче OpenBUGS) :-)

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 23 ноя 2015, 12:35

gamm, Игорь, спасибо!

ymr3R9Jge
Активный участник
Сообщения: 117
Зарегистрирован: 31 окт 2011, 00:18
Репутация: 14
Откуда: Кривий Ріг

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение ymr3R9Jge » 23 ноя 2015, 15:52

Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 23 ноя 2015, 18:37

ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.
ymr3R9Jge, обязательно гляну... спасибо!

ymr3R9Jge
Активный участник
Сообщения: 117
Зарегистрирован: 31 окт 2011, 00:18
Репутация: 14
Откуда: Кривий Ріг

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение ymr3R9Jge » 02 дек 2015, 13:59

ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.
В статье на сайте описан пример интерполяции временного ряда.

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 02 дек 2015, 19:57

ymr3R9Jge писал(а):
ymr3R9Jge писал(а):Если Вас устроит линейная или сплайновая интерполяция между известными значениями - в пакете zoo есть функция na.approx для восстановления пропущенных значений во временных рядах.
В статье на сайте описан пример интерполяции временного ряда.
Замечательно! Спасибо! Теперь понятно связь этой статьи с вопросом Вашим: viewtopic.php?f=31&t=19602 Что-то так и не собрался в своё время ответить - извините - у самого сразу не получилось.. Ещё раз спасибо!

ymr3R9Jge
Активный участник
Сообщения: 117
Зарегистрирован: 31 окт 2011, 00:18
Репутация: 14
Откуда: Кривий Ріг

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение ymr3R9Jge » 02 дек 2015, 22:10

Да, вопрос задавал, когда только приступал к "разбору полётов". Но на функцию na.approx сам вышел чисто случайно! До этого долго не удавалось найти удобный инструмент для подобных задач. Первые наши расчёты загрязнений были не очень хорошими, поскольку пришлось выкинуть треть значений. Но теперь с интерполяцией - можно два временных ряда соединять "ноздря в ноздрю".

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: Временной ряд интегральной характеристики с пропуском...

Сообщение nickleb » 20 сен 2016, 13:41

gamm писал(а): ... Из пакетов я бы посмотрел в сторону R-INLA, там есть подобные инструменты...
Вложения
v63i20.pdf
Статья об R-INLA
(1.5 МБ) 471 скачивание

Ответить

Вернуться в «R»

Кто сейчас на конференции

Сейчас этот форум просматривают: Bing [Bot] и 11 гостей