Add simulateLibSizes function

DESCRIPTION

 Package: splatter
 Type: Package
 Title: Simple Simulation of Single-cell RNA Sequencing Data
-Version: 0.3.4
+Version: 0.3.5
 Date: 2016-10-09
 Author: Luke Zappia
 Authors@R: as.person(c(

R/simulate.R

@@ -71,6 +71,7 @@ splat <- function(params = defaultParams(), method = c("groups", "paths"),
 
     sim <- simulateGeneMeans(sim, params)
     sim <- simulateDE(sim, params)
+    sim <- simulateLibSizes(sim, params)
 
     # Create new SCESet to make sure values are calculated correctly
     sce <- newSCESet(countData = counts(sim),
@@ -126,6 +127,18 @@ simulateDE <- function(sim, params) {
     return(sim)
 }
 
+simulateLibSizes <- function(sim, params) {
+
+    n.cells <- getParams(params, "nCells")
+    lib.loc <- getParams(params, "lib.loc")
+    lib.scale <- getParams(params, "lib.scale")
+
+    exp.lib.size <- rlnorm(n.cells, lib.loc, lib.scale)
+    pData(sim)$ExpLibSize <- exp.lib.size
+
+    return(sim)
+}
+
 #' Get log-normal factors
 #'
 #' Randomly generate multiplication factors from a log-normal distribution.