calcAsymptoticBayesFactorWakefield: Asymptotic Bayes factor

In timflutre/rutilstimflutre: Timothee Flutre's personal R code

Asymptotic Bayes factor

Description

Calculate the asymptotic Bayes factor proposed by Wakefield in Genetic Epidemiology 33:79-86 (2009, http://dx.doi.org/10.1002/gepi.20359). Can also be averaged over a grid of values of W, as done in various papers from Matthew Stephens' lab.

Usage

calcAsymptoticBayesFactorWakefield(
  theta.hat,
  V,
  W = c(0.1, 0.2, 0.4, 0.8, 1.6),
  weights = NULL,
  log10 = TRUE
)

Arguments

theta.hat	vector of MLE(s) of the additive genetic effect(s)
V	vector of the corresponding variance(s) of theta.hat
W	vector of variance(s) of the prior on theta (if several values, the ABF will be averaged over them); the vector of default values comes from the single-SNP R implementation of BLMM by Wen for his 2015 article (see https://github.com/xqwen/blmm)
weights	weights used to average over the grid of W (all equal by default)
log10	return the log10 of the ABF

Value

numeric

Author(s)

Timothee Flutre

timflutre/rutilstimflutre documentation built on Feb. 7, 2024, 8:17 a.m.