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

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

How to extract the estimates from a stanfit object

stan

rstan