bayes.helpers.R

require('dplyr')
require('rethinking')

dmode <- function(x) {
  den <- density(x, kernel=c("gaussian"))
  ( den$x[den$y==max(den$y)] )   
}  

meanhpdi <- function(x, prob=.95) {
  i <- rethinking::HPDI(x, prob=prob)
  data_frame(y=mean(x), ymin=i[1], ymax=i[2] )
}

mmmhpdi <- function(x, prob=.95) {
  i <- rethinking::HPDI(x, prob=prob)
  data_frame(y=mean(x), mode=dmode(x), median=median(x), ymin=i[1], ymax=i[2] )
}