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

Question

在将stan_glmer()作为一个族使用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中的错误。