With bootstrapping you can obtain the variability of various summary statistics. Recently I came across the Bayesian bootstrap (http://www.sumsar.net/blog/2015/07/easy-bayesian-bootstrap-in-r/). I was curious how large the difference are between the Bayesian and Normal bootstrap. As difference I focus on variability of the mean difference between two hypothetical normal distributions (x = N(0.1, 1) and y = N(0, 1)) with corresponding 95% confidence and credibility intervals (?). I played around with various extremes of the sample sizes (n = 3 or n = 1000). My expectations were that the results from the Bayesian bootstrap would be more conservation (e.g. wider credibility intervals (?)) than that of the normal bootstrap, this seems not to be the case. My question is why?

I am sorry for my ignorance regarding the Baysian bootstrap principle and I am also not certain if the 2.5% and 97.5% can be named credibility intervals. I often like to say that I just obtained the 2.5% and 97.5% of of the bootstrap to avoid using these terms. I think I might do something incredibly wrong here.

Similar questions and discussions