Notebook on nbviewer (original) (raw)

%%R -o counts -o meta

suppressMessages(library(splatter))

params <- newSplatParams()

params <- setParam(params, "nGenes", 5000) params <- setParam(params, "batchCells", c(500,500,500)) params <- setParam(params, "batch.facLoc", 0.5) params <- setParam(params, "batch.facScale", 0.5) params <- setParam(params, "group.prob", c(1/3,1/3,1/3))

sim <- splatSimulate(params, method="groups", verbose=FALSE)

counts = data.frame(counts(sim)) meta = data.frame(colData(sim))