This function samples from the posterior distribution of log(OR) based on the observation of a case genotype <code>post_lor</code> samples from the posterior distribution of log(OR) given genotype.

This function samples from the posterior distribution of log(OR) based on the observation of a case genotype post_lor samples from the posterior distribution of log(OR) given genotype.

post_lor(gt = c(0, 1, 2), a1, b1, p0, nsim)

Arguments

gt	a scalar - takes the values 0 - reference hom, 1 - het, 2 - alternative hom.
a1	a scalar - beta distribution alpha shape parameter for prior on allele frequency in cases
b1	a scalar - beta distribution beta shape parameter for prior on allele frequency in cases
p0	a vector - a set of samples from beta(a0,b0) i.e. prior on control allele frequency
nsim	a scalar - number of log odds ratios to simulate from the posterior distribution

Value

vector of simulated log(OR) from the posterior.

This function samples from the posterior distribution of log(OR) based on the observation of a case genotype post_lor samples from the posterior distribution of log(OR) given genotype.

Arguments

Value

This function samples from the posterior distribution of log(OR) based on the observation of a case genotype `post_lor` samples from the posterior distribution of log(OR) given genotype.