### Function: Likelihood Plot Values likelihood <- function(N,Y){ a <- Y + 1 b <- N - Y + 1 dom <- seq(0,1,0.005) val <- dbeta(dom,a,b) return(data.frame('x'=dom, 'y'=val)) }