vsbuffalo · February 25, 2015 06:45
diff --git a/hpdi_stat.R b/hpdi_stat.R
 hpdi_stat <- function(x, prob=0.95) {
   # stat functions for ggplot that handle HPDI from samples
   # automatically.
   hpdi <- HPDinterval(as.mcmc(x))
   c(y=mean(x), ymin=hpdi[1, 'lower'], ymax=hpdi[1, 'upper'])
 }
	hpdi_stat <- function(x, prob=0.95) {
	# stat functions for ggplot that handle HPDI from samples
	# automatically.
	hpdi <- HPDinterval(as.mcmc(x))
	c(y=mean(x), ymin=hpdi[1, 'lower'], ymax=hpdi[1, 'upper'])
	}