npatient <- 3; target <- 0.3
para.prior = list(alp.prior = target, bet.prior = 1 - target)
result <- gammatable(npatient, target, para.prior)
plot(result)
#This example may cost you a long time to run
npatient <- 30; target <- 0.3
para.prior = list(alp.prior = target, bet.prior = 1 - target)
result <- gammatable(npatient, target, para.prior)
plot(result)
Run the code above in your browser using DataLab