This function estimates shape parameters for a prior distribution of control allele frequencies est_a1b1 estimate shape parameters a1,b1

This function estimates shape parameters for a prior distribution of control allele frequencies est_a1b1 estimate shape parameters a1,b1

est_a1b1(a0, b0, target.or, target.prob, nsim)

Arguments

a0	a scalar - beta distribution alpha shape parameter for prior on allele frequency in controls
b0	a scalar - beta distribution beta shape parameter for prior on allele frequency in controls
target.or	a scalar - an odds ratio threshold to compute P(sim.or > target.or)
target.prob	a scalar - a probability that a sampled variant will exceed target.or
nsim	a scalar - number of log odds ratios to simulate from the posterior distribution

Value

list - shape parameters for beta distribution satisfying target.or,target.prob and control prior distribution constraints.