Как указать строгие априоры для фиксированных коэффициентов в rstanarm?

Как получить априорное распределение для коэффициента со средним значением в диапазоне 1e-05 и стандартным отклонением (sd) в диапазоне 1e-06?

Какие распределения выбрать для такой высокой точности (sd ‹= 1e-06) параметров в rstanarm?

1) Студенческая семья

2) Иерархическая семья усадки

3) семья Лапласа

4) Семья Дирихле

5) другие

Я пытался использовать prior = normal(location = 1e-05, scale = 1e-06, autoscale = FALSE), но выборка занимает много времени.


r regression bayesian rstanarm
person Mukund Komati    schedule 28.05.2020    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Я думаю, что prior = normal(location = 1e-05, scale = 1e-06, autoscale = FALSE) - ваш лучший выбор. Тот факт, что для выборки требуется много времени, не является признаком того, что что-то не так с отрисовкой, которую она производит. Но я бы умножил результат на 1 миллион или разделил предиктор на 1 миллион, чтобы prior = normal(location = .1, scale = 1, autoscale = FALSE) было разумным.

person Ben Goodrich    schedule 28.05.2020
comment
Спасибо за предложение, масштабирование значительно улучшило время выборки. Кроме того, что означает Prior_PD = TRUE и QR = TRUE внутри вызова stan_lmer(...)? Как это влияет на расчет? - person Mukund Komati; 29.05.2020
comment
Оба (фактически все) эти аргументы для stan_lmer задокументированы. - person Ben Goodrich; 29.05.2020