knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
The functions in this package focus on analyzing case-control association studies involving a group of genetic variants. In particular, we are intersted in modeling the outcome variable as a function of multivariate genetic profile using Bayesian model uncertainty and variable selection techniques. The package incorporates functions to analyze data sets involving common variants as well as externsions to model rare variants via the Bayesian Risk Index (Quintana and Conti, 2011), as well as haplotypes. Finally, allows the incorporation of biological information to inform the marginal inclusion probabilities via the iBMU (Quintana and Conti (submitted)).
library(devtools) # Master branch install_github("USCbiostats/bvs")
As a first example, below is code to run a rare variant analysis using the "sample" method on a binary outcome.
# load bvs library(bvs) # load the rare variant data set, first column is outcome variable data(RareData) # run 2000 iterations of the MH algorithm (rare = TRUE will construct the BRI) rare_analysis <- bvs(y = RareData$case, x = RareData[, -1], family = "binomial", method = "sample", rare = TRUE, iter = 2000) # summarize results of analysis rare_results <- summary(rare_analysis) # plot of top 25 models and top 10 variants plot(rare_results, num_models = 25, num_snps = 10)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.