После установки функций stan_glmer()
или stan_glm()
с mcgv::betar
в качестве семейства я получаю сообщение об ошибке при попытке вызвать posterior_predict
на нем. Р говорит:
Ошибка в выражении (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)