如何从标准拟合对象中提取估计值

Question

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

预期 A 后验

   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 <- 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% 可信区间的端点。