如何从 stanfit 对象中提取估计值答案

【问题标题】：How to extract the estimates from a stanfit object如何从 stanfit 对象中提取估计值
【发布时间】：2019-01-21 01:25:32
【问题描述】：

让fit 成为stanfit 对象。然后我可以通过以下方式提取参数AAA的估计值：

预期的后验

   EAP <- as.data.frame(summary(fit)[[1]])["AAA","mean"]

95% 可信区间

   lower.CI <- as.data.frame(summary(fit)[[1]])["AAA","2.5%"]
   upper.CI <- as.data.frame(summary(fit)[[1]])["AAA","97.5%"]

但我不确定，这是一种规范的方法吗？如果有更简单的代码，请告诉我。

谅解备忘录

EAP

EAP <- get_posterior_mean(fit,par=c("AAA"))

【问题讨论】：

有一个 RStan 小插图详细说明了姿势提取：cran.r-project.org/web/packages/rstan/vignettes/…

标签： stan rstan

【解决方案1】：

平均而言，get_posterior_mean 函数可能更规范一些。对于分位数，我会做类似quantile(extract(fit, pars = "AAA")[[1]], probs = c(0.1, 0.9)) 的事情。但是，使用 Stan 的默认设置无法非常精确地估计 95% 可信区间的端点。

【讨论】：

感谢您的快速和友好的回复，我非常感谢您的每一次回复！
要估计可信区间的精确端点，我必须从默认值更改什么？或者任何网页解释它？
这是probs 的参数quantile。