Drokhlyansky, Smillie, van Wittenberghe, et al. (Cell, 2020)
This file contains code for detecting heterotypic multiplets in MIRACL-Seq data. Code for other analyses found throughout the paper is available here:
https://github.com/cssmillie/ulcerative_colitis
For other questions, please email Chris Smillie (cssmillie [AT] gmail.com).
Before using this code, you must install the "dnormpar2" R package (included):
R CMD INSTALL dnormpar2.tar.gz
This contains a fast method for calculating the doublet statistics.
source('MIRACL.r')
res = find_doublets(data, ident, sig.use)
doublets = res[[length(res)]]$doublets
Please see code comments for more detailed usage instructions.