rstanarm - 如何在 rstanarm 中指定对数正态先验或用户定义的非负先验？

Question

假设我有一个混合效应模型的以下公式： Precipitation ~ s(month,bs="cc")+s(time)+ humidity,random= ~(humidity|year) 并且我知道湿度只会对降水产生积极影响。所以我想为 rstanarm 中的湿度指定一个对数正态先验或其他非负先验。假设湿度的系数是beta。我会指定先验为log(beta) ~ normal(0,10^4)我应该怎么做？我不认为默认的指数分布是适当的先验。

Answer 1

这不在 rstanarm 支持的先验之中。相同的语法适用于 brms，尽管您需要指定先验

my_prior <- prior(lognormal(0, 1e4), coef = "humidity")

但是，以 10^4 作为对数标准差的对数正态先验是荒谬的，因为它对在 64 位计算机上溢出到无穷大的值具有相当大的概率。