如何从stanfit对象中提取估计

Question

让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"))

Answer 1

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