This is an R package for quickly simulating single cell RNAseq read count data. The model is based on data from droplet based methods, for other methods Splatter may be better suited. Sinsynthr shares a similar core model to Splatter but is overall much simpler.
# install.packages("devtools")
devtools::install_github("marchinilab/dropsim")library(dropsim)
# Use default parameters
new_parameters <- dropsim_parameters()
# Simulate counts
simulation <- simulateDGE(new_parameters)
# Normalise Counts matrix
normalised_dge <- normaliseDGE(simulation$counts)
# Do a PCA to check
dge_pca <- prcomp(normalised_dge)
qplot(dge_pca$x[,1], dge_pca$x[,2], colour=rownames(dge_pca$x)) + labs(colour="Group", y="PC2", x="PC1")
For more detailed examples of how to use this package please see the vignette.
The R/ directory contains source code of the functions to simulate the data.
The man/ directory contains the manual pages for the functions, compiled by roxygen.
The vignettes/ directory contains uncompiled Rmarkdown with examples of use.
The tests/ directory contains the unit tests which are carried out by the testthat R package.
Contributions of any size or form are welcome!