目标:在 stan_glm 模型的一个命令中估计具有置信区间的系数。
尝试使用 fixef() 提取我的 stan_glm 模型的系数 CI 时,出现以下错误。我尝试使用来自bayr 的fixef() 而不是rstanarm,因为它可以方便地提供置信区间。
Error: Can't combine `..1$fixef` <character> and `..3$fixef` <logical>.
这是我尝试过的最简单的模型:
M1 <- stan_glm(variable ~ 1, data = data)
fixef(M1)
我首先认为这可能是库中的冲突,所以我尝试以不同的顺序激活使用 fixef (rstanarm, bayr) 的库,但这没有帮助。目前我使用以下库(按此顺序):
看门人、避风港、data.table、modeest、dplyr、tidyverse、readxl、gridExtra、rstanarm、devtools、mascutils、bayr
关于该错误的来源有什么建议吗?
编辑:
rstanarm::fixef(M1)
有效,但给出没有置信区间的系数。bayr::fixef() 会自动给他们。我试过也试过以下:
coef(M1)
"Error: Can't combine `..1$fixef` <character> and `..3$fixef` <logical>."
posterior_interval(M01) #gives only confidence interval without the coefficient estimate