est_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 est_a1b1 estimate shape parameters a1,b1

Usage

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.

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