Как извлечь оценки из объекта stanfit

Question

Пусть fit будет stanfit объектом. Затем я могу извлечь оценки параметра AAA следующим образом:

Ожидаемый A Posterior

   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)).Однако конечные точки 95% вероятных интервалов оцениваются не очень точно с настройками по умолчанию для Stan.