Сравнение неких модельных полей с фактическими данными

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

Сравнение неких модельных полей с фактическими данными

Сообщение nickleb » 23 апр 2015, 09:42

Здравствуйте. Есть за некий период осреднения по времени два регулярных растра в одном и том же домейне
в одной и той же равнопромежуточной проекции:
1) поле температуры - некая МОДЕЛЬ 1
2) поле температуры - некая МОДЕЛЬ 2

Нашёл реальные данные температуры за тот же период -
превратил их в
3) растр той же проекции - РЕАЛЬНЫЕ данные - иррегулярная сетка.

Надо сравнить какая из двух МОДЕЛЕЙ годится больше.

В "Сёрфере" Residuals'ы построить - такой вариант?

Или кто посоветует, м.б., какой-то инструмент под эти цели годящийся дающий какие то сводные характеристики и оценки для принятия решения... Желательно из свободного soft'а...
Последний раз редактировалось nickleb 02 июн 2015, 09:41, всего редактировалось 1 раз.


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

Re: Сравнение неких модельных полей с фактическими данными

Сообщение nickleb » 23 апр 2015, 11:18


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

Re: Сравнение неких модельных полей с фактическими данными

Сообщение gamm » 24 апр 2015, 16:50

не мудрите - в том же R постройте МВА-интерполяцию с данных на свои точки, и вычтите ваши значения в них, получите аналог residuals. А дальше обычный анализ выборки - mean(), rmse, sd(), quantiles(), hist(), curtosis(), ...

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

Re: Сравнение неких модельных полей с фактическими данными

Сообщение nickleb » 24 апр 2015, 17:42

gamm писал(а):не мудрите - в том же R постройте МВА-интерполяцию с данных на свои точки, и вычтите ваши значения в них, получите аналог residuals. А дальше обычный анализ выборки - mean(), rmse, sd(), quantiles(), hist(), curtosis(), ...
как всегда вовремя, как всегда мудро... спасибо!
На самом деле, уважаемый gamm, такая штука:
с послевоеннных лет регулярно наблюдали по времени и пространству довольно обширную акваторию, затем в 93-ем - крайняя съёмка - "кердык" дорогостоящим экспедициям и только эпизодические наблюдения... 2007 г. - МеждГеофизГод - всплеск - что -то по инерции ещё до 2010, сейчас опять всё затухает... каким инструментарием попробовать заполнить временную лакуну с 94 по 2006? - при наличии, повторяю, эпизодических наблюдений за этот период на крайне малых участках исследуемой акватории? Сейчас вот пытаюсь NCEP'овским всякими моделями "жонглировать"? А как грамотно по-"шаманить" не ведаю... Подскажите, что посмотреть, поучить...

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

Re: Сравнение неких модельных полей с фактическими данными

Сообщение gamm » 24 апр 2015, 19:24

nickleb писал(а):каким инструментарием попробовать заполнить временную лакуну с 94 по 2006? - при наличии, повторяю, эпизодических наблюдений за этот период на крайне малых участках исследуемой акватории? Сейчас вот пытаюсь NCEP'овским всякими моделями "жонглировать"? А как грамотно по-"шаманить" не ведаю... Подскажите, что посмотреть, поучить...
единственный подход, который может реально помочь - это использовать Байесовские методы (MCMC - Markov chain Monte-Carlo). Сейчас есть хорошая штука - INLA, которая работает очень быстро и стабильно, мы пробовали (пока прикидку) для восполнения дырявых временных рядов по метеостанциям с построением полного поля на всю территорию. Но это не кнопочный продукт - из нашего опыта, модель нужно строить несколько месяцев, пока заработает + нужен эксперт, который понимает, чего может быть, а чего - нет. Если такой эксперт на вашей стороне есть, то можно пообсуждать задачку - летом, когда посвободнее буду.

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

Re: Сравнение неких модельных полей с фактическими данными

Сообщение nickleb » 25 апр 2015, 09:18

единственный подход, который может реально помочь - это использовать Байесовские методы (MCMC - Markov chain Monte-Carlo). Сейчас есть хорошая штука - INLA, которая работает очень быстро и стабильно, мы пробовали (пока прикидку) для восполнения дырявых временных рядов по метеостанциям с построением полного поля на всю территорию. Но это не кнопочный продукт - из нашего опыта, модель нужно строить несколько месяцев, пока заработает + нужен эксперт, который понимает, чего может быть, а чего - нет. Если такой эксперт на вашей стороне есть, то можно пообсуждать задачку - летом, когда посвободнее буду.
хорошо, gamm

Ответить

Вернуться в «Общий - ПО»

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

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 25 гостей