Как сравнить два гиперпараметра в иерархической модели? - PullRequest
0 голосов
/ 18 декабря 2018

В одной иерархической модели у нас есть два параметра: dnorm (A_mu, 0.25 ^ -2) и dnorm (B_mu, 0.25 ^ -2).В этом случае 0,25 является SD, я использую фиксированное число.A_mu и B_mu представляют среднее значение на уровне группы.После подбора данных по rjags мы получаем распределения для каждого параметра.Так что я просто сравниваю самый высокий задний интервал плотности (ИЧР) A_mu и B_mu?Нужно ли что-то вычислять, используя sd (0.25)?

В другом случае, если sd двух гиперпараметров не фиксирован, например: dnorm (A_mu, A_sd) и dnorm (B_mu, B_sd).Как я могу сравнить два гиперпараметра и принять решение, например, эта группа значительно отличается от другой группы?

1 Ответ

0 голосов
/ 07 января 2019

Помните, что вы получаете последующие распределения для A_mu и B_mu.Это облегчает сравнение, так как вы можете посмотреть 95-процентные доверительные интервалы (CI) для параметров (или выбрать доверительное значение, которое удовлетворяет вашим потребностям).Я полагаю, что JAGS использует выборку Гиббса, поэтому вы сможете получить необработанные сэмплы от авторов для A_mu и B_mu.Затем вы можете спросить «какова вероятность того, что B_mu больше некоторого значения?» , рассчитав процент задних выборок, превышающих это значение.В качестве альтернативы и аналогично частому тестированию гипотез, вы можете спросить, какова вероятность того, что среднее значение B_mu является ничьей с задней стороны A_mu.Таким образом, ключ заключается в том, чтобы просто использовать образцы из вашего апостериорного.Я бы порекомендовал взглянуть на учебник Эндрю Гельмана по BDA3 (глава 4), чтобы получить действительно хороший справочник по этим понятиям.

Несколько вещей, которые следует иметь в виду, прежде чем делать выводы из данных: (1) вы должны всегдапроверьте достоверность своих цепей Маркова, оценивая такие вещи, как автокорреляция (2), попробуйте выполнить последующую прогностическую проверку, чтобы убедиться, что ваша модель хорошо вписывается в данные.Если ваша модель плохо соответствует данным, вы можете получить очень вводящие в заблуждение результаты описанной выше процедуры.

...