Как проверить сходимость в модели ошибок?

Я хочу объяснить сходимость в модели ошибок с помощью команды plot(). Пример вывода показан на следующем рисункевведите здесь описание изображения

Я не уверен, что смогу хорошо прочитать этот вывод, спасибо всем :)


person angeella    schedule 28.01.2016    source источник


Ответы (2)


К сожалению, не похоже, что вы можете подтвердить конвергенцию на рисунке, который вы показываете (EDIT: есть по крайней мере некоторая информация, см. Ниже). Левая часть рисунка представляет собой просто график гусеницы, который фактически просто показывает 95% интервалов распределения для каждого параметра.

Оценка конвергенции — гораздо более тонкий процесс, поскольку существует несколько способов определить, сходится ли ваша модель. Вам нужно будет определить, что ваша модель надлежащим образом исследовала пространство параметров для каждого параметра (через графики трассировки, функция traceplot в библиотеке coda), между и внутри цепочки дисперсии (диагностика Гельмана-Рубина, gelman.diag в библиотеке coda) и автокорреляция в ваших цепочках (autocorr.plot в coda). Существует множество других мер, которые другие предложили для оценки того, сходится ли ваша модель, и просмотр остальной части пакета coda проиллюстрирует это.

Я настоятельно рекомендую вам ознакомиться с учебным пособием WINBUGS в их руководство пользователя (ссылка на pdf), в нем есть раздел, посвященный проверке сходимости модели. Вы хотите убедиться, что трассировочные графики хорошо перемешаны (посмотрите в руководстве, чтобы узнать, что это значит), что ваша диагностика Гельмана Рубина составляет ‹ 1,10 для каждого параметра (общее правило) и что ваши цепочки не слишком коррелированы (это уменьшит ваш эффективный размер выборки в ваших цепочках).

Удачи, и почитайте немного по этому вопросу, это очень поможет вам, если вы интересуетесь байесовским выводом!

Редактировать

Как указал @jacobsocolar, и я полностью пропустил, графики, доступные в этом вопросе, по крайней мере содержат некоторую информацию, указывающую на то, что модель действительно сошлась. Я не видел график R-шляпы в правой части левого графика. Эти значения должны быть меньше 1,1 для каждого параметра, если модель действительно сходилась. Взгляд на приведенный выше график намекает на то, что модель сошлась, но это было бы намного легче увидеть, если бы на отметке 1,1 на графике была вертикальная линия, которой нет.

person mfidino    schedule 05.02.2016
comment
Большое спасибо! теперь я попытаюсь понять код пакета, чтобы провести еще один анализ сходимости. Я понял, что трудно подтвердить конвергенцию только этим графиком, но теперь я могу прочитать этот график. Спасибо :) - person angeella; 05.02.2016
comment
Рядом с графиком гусеницы находится график значений r-hat (psrf). Я согласен с тем, что оценка конвергенции имеет нюансы, но это хорошее начало. В целом, похоже, что модель сошлась, предполагая, что начальные значения были правильно разогнаны и ничего странного не происходит. - person Jacob Socolar; 29.02.2016
comment
Приношу свои извинения, как отметил @JacobSocolar, здесь есть некоторая информация. Я обновил свой ответ, чтобы решить эту проблему. - person mfidino; 06.04.2017

Вашей выходной цифры действительно достаточно, чтобы (начать) оценивать конвергенцию, вопреки ответу M_Fidino. Рядом с графиком гусеницы есть график значений «r-шляпы». Это статистика Гельмана-Рубина — отношение дисперсии между цепями к дисперсии внутри цепочки, и все они ‹ 1,10. рассредоточенный. В остальном я согласен со всем в ответе M_Fidino.

person Jacob Socolar    schedule 29.02.2016