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

Question

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

Answer 1

这不在rstanarm支持的先例中。同样的语法也适用于brms，尽管您需要指定之前的类似

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

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