Rescale Phi when sampling in BASiCSSimulate

DESCRIPTION

 Package: splatter
 Type: Package
 Title: Simple Simulation of Single-cell RNA Sequencing Data
-Version: 1.9.6
-Date: 2019-09-17
+Version: 1.9.7
+Date: 2019-09-19
 Author: Luke Zappia
 Authors@R:
     c(person("Luke", "Zappia", role = c("aut", "cre"),

NEWS.md

+## Version 1.9.7 (2019-09-19)
+
+* Rescale when sampling Phi in BASiCSSimulate
+
 ## Version 1.9.6 (2019-09-17)
 
 * Adjust BASiCSSimulate to match development version of BASiCS

R/BASiCS-simulate.R

@@ -72,6 +72,8 @@ BASiCSSimulate <- function(params = newBASiCSParams(), verbose = TRUE, ...) {
                     "spike.means will be sampled.")
             selected <- sample(length(spike.mu), nSpikes, replace = TRUE)
             spike.mu <- spike.mu[selected]
+
+            params <- setParam(params, "spike.mu", spike.mu)
         }
     }
 
@@ -81,6 +83,10 @@ BASiCSSimulate <- function(params = newBASiCSParams(), verbose = TRUE, ...) {
                 "cell.params will be sampled.")
         selected <- sample(nrow(cell.params), nCells, replace = TRUE)
         cell.params <- cell.params[selected, ]
+
+        cell.params$Phi <- (cell.params$Phi / sum(cell.params$Phi)) * nCells
+
+        params <- setParam(params, "cell.params", cell.params)
     }
 
     thetas <- getParam(params, "theta")