Эффективный размер выборки для последующего вывода из выборки MCMC

Вопрос, который вы задаете, отличается от «диагностики сходимости». Допустим, вы выполнили всю диагностику конвергенции (выберите свою любимую (ые)) и теперь готовы начать отбор проб сзади.

Существует два варианта с точки зрения эффективного размера выборки (ESS), вы можете выбрать одномерный ESS или многовариантный ESS. Однофакторный ESS обеспечит эффективный размер выборки для каждого параметра в отдельности, а консервативные методы требуют, чтобы вы выбирали наименьшую оценку. Этот метод игнорирует все взаимные корреляции между компонентами. Это, вероятно, то, что большинство людей использовали некоторое время

Недавно было представлено многомерное определение ESS. Многомерный ESS возвращает одно число для эффективного размера выборки для количеств, которые вы хотите оценить; и делает это путем учета всех взаимных корреляций в процессе. Лично я предпочитаю многомерный ESS. Предположим, вас интересует вектор средних значений апостериорного распределения. MESS определяется следующим образом: $p$ Вот

mESS = n {(\frac{| Λ |}{| Σ |})}^{1 / p} .

$\text{mESS} = n \left(\dfrac{|\Lambda|}{|\Sigma|}\right)^{1/p}.$

- ковариационная структура апостериора (также асимптотическая ковариация в CLT, если у вас были независимые выборки) $\Lambda$
$\Sigma$ $\Lambda$
$p$
$|\cdot|$

$\Lambda$ $\Sigma$ multiESS

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

$\epsilon$ $\epsilon$
$\alpha$
$p$

mESS \geq \frac{2^{2 / p} π}{(p Γ (p / 2))^{2 / p}} \frac{χ_{1 - α, p}^{2}}{ϵ^{2}},

$\text{mESS} \geq \dfrac{2^{2/p} \pi}{(p \Gamma(p/2))^{2/p}} \dfrac{\chi^2_{1-\alpha, p}}{\epsilon^2},$

$\Gamma(\cdot)$ minESS

$p = 20$ $95\%$ $\epsilon = .05$

> minESS(p = 20, alpha = .05, eps = .05)
[1] 8716

Это верно для любой проблемы (в условиях регулярности). Способ, которым этот метод адаптируется от проблемы к проблеме, состоит в том, что медленное перемешивание цепей Маркова занимает больше времени, чтобы достичь этой нижней границы, поскольку mESS будет меньше. Так что теперь вы можете пару раз multiESSпроверить, достигла ли ваша цепь Маркова этой границы; если не пойти и взять больше образцов.

Greenparker
источник

(+1) Отличный ответ. Вы знаете, была ли функция multiESSзакодирована для других языков, таких как MATLAB? (или это будет трудно переопределить?)

lacerbi

Σ

$\Sigma$

Σ

$\Sigma$

@lacerbi Я рад, что вы смогли написать код в Matlab. Если возможно, ответьте на этот комментарий, когда он закончится, чтобы я мог его использовать. Спасибо

Greenparker

Моя реализация multiESS для MATLAB доступна здесь . Это рабочая версия, хотя она потребует дополнительного тестирования (я не знаком с R, иначе я бы сравнил ее с реализацией R).

Lacerbi

Сходимость зависит от нескольких факторов: количества параметров, самой модели, алгоритма выборки, данных ...

Я хотел бы предложить избегать какого-либо общего правила и использовать несколько инструментов диагностики сходимости для определения соответствующего количества выгорания и уменьшения числа итераций в каждом конкретном примере. Смотрите также http://www.johnmyleswhite.com/notebook/2010/08/29/mcmc-diagnostics-in-r-with-the-coda-package/, http://users.stat.umn.edu/~geyer/mcmc/diag.html.

Montecristo
источник

Эффективный размер выборки для последующего вывода из выборки MCMC

Ответы: