0

目标:在 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
4

0 回答 0