/BayesianGLM

R scripts for Bayesian generalized linear model in La Manno 2016

Primary LanguageR

BayesianGLM

Bayesian generalized linear regression model for estimation of cell-specific effect for gene expression(La Manno et al.Cell 2016) in R and Rstan.
Easy to understand example of MCMC is here

How to use

Installation

  1. Clone or download repository and create OUTPUT and DATA directory in the repository directory.
  2. Download datafiles(.cef) to DATA directory from the authors repository.
    Note: If you use UNIX, these procedure are automated in init.sh.

Package Requirement

Rstan, ggplot2, pipeR, foreach, reshape2, dplyr, tibble

Option for parallelization doParallel

Execution

Run "runGLM.R"
You can modify DATASET/GENESET/CELLTYPESET in the script to create your own plot.