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
- Clone or download repository and create OUTPUT and DATA directory in the repository directory.
- Download datafiles(.cef) to DATA directory from the authors repository.
Note: If you use UNIX, these procedure are automated in init.sh.
Rstan, ggplot2, pipeR, foreach, reshape2, dplyr, tibble
Option for parallelization doParallel
Run "runGLM.R"
You can modify DATASET/GENESET/CELLTYPESET in the script to create your own plot.