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.