# NOT RUN {
   n <- 12         # sample size
   y <- (0:(n-1))  # observations
   a <- b <- 1/3   # the neutral beta prior
   rvbeta(1, shape1=a+y, shape2=b+n-y)
   rvnbeta(1, shape1=y, shape2=n-y)
# }
Run the code above in your browser using DataLab