bayesian - 使用 stan_glmer 进行 beta 回归的 rstanarm 错误

Question

stan_glmer()在将 a或stan_glm()函数作为一个家庭进行拟合后mcgv::betar，当我尝试调用posterior_predict它时出现错误。R 说：

exp(eta) 中的错误：数学函数的非数字参数

一个最小的例子：

library(rstanarm)
library(loo)
library(mgcv)

a <- rnorm(100, 0.5, 0.1)
b <- a+rnorm(100, 0.6, 0.01)
d <- data.frame(a=a, b=b)

fit <- stan_glm(a ~ b,
               data = d,
               family=betar,
               chains = 10,
               seed = 1)

posterior_predict(fit)

Answer 1

我在这里找到了答案：https ://discourse.mc-stan.org/t/rstanarm-mgcv-betar-family/2947/4 。这是 rstanarm 中的一个错误。