### Function: Mode of Posterior Beta mode_of_posterior <- function(m,n,N,Y){ a <- Y + (n*m) -1 b <- N - Y + (n*(1-m)) - 1 mode_posterior <- (a-1)/(a+b-2) return(mode_posterior) }