opt_a1b1: This function estimates shape parameters for a prior...

In ollyburren/cupcake: PCA for GWAS

Description

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

Usage

opt_a1b1(a0, b0, target.or, target.prob, n.steps)

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
n.steps	an integer - length of search grid to employ - larger gives a more accurate integral estimate

Value

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

ollyburren/cupcake documentation built on July 30, 2020, 9:58 a.m.