This is an R package to compute posteriors from a gene-level hierarchical prior, as described in the manuscript "Bayesian genome- and epigenome-wide association studies with gene level dependence" [1]. Also includes functions to perform two-class testing with shared kernels for methylation array data (see [2]). The file BinaryExample.r gives example code for simulated categorical data. The file GliomaMethylExample.r gives example code for an application to methylation screening data.
The package can be installed in R, directly from GitHub, using the devtools library:
install.packages(devtools)
library(devtools)
install_github("lockEF/BayesianScreening")
[1] Lock, E. F. & Dunson, D. B. (2016). Bayesian genome- and epigenome-wide association studies with gene level dependence. arXiv preprint: http://arxiv.org/abs/1604.08654 .
[2] Lock, E. F. & Dunson, D. B. (2015). Shared kernel Bayesian screening. Biometrika, 102 (4): 829-842.