Questions tagged [credible-interval]

Credible intervals are an interval within which an unobserved parameter value falls with a particular probability. They are an important concept in Bayesian statistics.

The purpose of credible intervals is to describe and summarise the uncertainty of statistical parameters.

A common way to sharpen the understanding of credibility is the comparison with the term "confidence".

Credibility vs Confidence

Credible intervals can be confused with confidence intervals. However, while their goal is similar, their statistical definition annd meaning is very different. Confidence is rooted in Bayesian inference and is obtained through a complex algorithm full of rarely-tested assumptions and approximations, credibility is fairly straightforward to compute. Credibility means that there is 95% likelihood that a population parameter lies in the designated interval. Many times this interpretation is falsly attributed to confidence (see e.g. Hoekstra, 2014). Confidence, however, works differently: "If we repeat the experiment infinitely many times, 95% of the experiments will capture the population parameter in their confidence intervals."

Finally, the typical interval level for confidence is 95%. While this was the first choice for credibility for some time several authors started to question this and suggested a default credibility of roughly 90% (Kruschke, 2014; McElreath, 2014, 2018).

References

  • Hoekstra, R., Morey, R. D., Rouder, J. N., & Wagenmakers, E. J. (2014). Robust misinterpretation of confidence intervals. Psychonomic Bulletin & Review, 21(5), 1157-1164.
  • Kruschke, J. (2014). Doing bayesian data analysis: A tutorial with r, jags, and stan. Academic Press.
  • McElreath, R. (2018). Statistical rethinking: A bayesian course with examples in r and stan. Chapman; Hall/CRC.
20 questions
0
votes
1 answer

Extract credible intervals for robust correlations in R

I currently know how to use pbcor from the WRS2 package to extract robust correlations. This function calculates the 95% bootstrap confidence intervals around the estimated robust correlation. For consistency with the rest of my analyses and…
0
votes
2 answers

determine equal-tail credible interval

I have obtained the posterior density for part d: $2 theta^{-1}(1- theta)^{-1}$. How do I plot in R the distribution to find the l and u such that $F_{theta| x} (l) = 0.025$ and $F_{theta| x} (u) = 0.975$? (the equal-tail interval)
siegfried
  • 461
  • 1
  • 5
  • 20
0
votes
0 answers

Coverage probability estimation for Bayesian Credible interval (Normal Distribution)

Bayesian Inference for the Normal Distribution, I use the following r code to obtain the posterior distribution. install.packages(c("mvtnorm","loo","coda"), repos="https://cloud.r-project.org/",dependencies=TRUE) options(repos=c(getOption('repos'),…
score324
  • 687
  • 10
  • 18
0
votes
2 answers

How to calculate a 95 credible region for a 2D joint distribution?

Suppose we have a joint distribution p(x_1,x_2), and we know x_1,x_2,p. Both are discrete, (x_1,x_2) is scatter, its contour could be drawn, marginal as well. I would like to show the area of 95% quantile (a scale of 95% data will be contained) of…
TxWang
  • 57
  • 2
  • 11
0
votes
0 answers

Reduce credible intervals in Causal Impact model

I'm hitting an issue with a causal impact model that I'm building. I'm trying to create a counter factual for daily sales at one store (nseasons = 7). I've included sales for 5 other stores nearby. Eyeballing a lineplot, it appears to me that trends…
1
2