Added parameter cutoffVersion to function power.sameReadDepth.restric…

…tedDoublets (missed when introducing the parameter)
heiniglab · Jun 25, 2021 · 7deb134 · 7deb134
1 parent f28265a
commit 7deb134
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 1 deletion.
diff --git a/R/power.R b/R/power.R
@@ -658,7 +658,8 @@ power.sameReadDepth.withDoublets<-function(nSamples,nCells,ct.freq,
   } else if (multipletRateGrowth == "constant") {
     multipletFraction<-multipletRate
   } else {
-    stop("No known option for multipletRateGrowth. Use the values 'linear' or 'constant'.")
+    stop(paste("Input", multipletRateGrowth, "No known option for multipletRateGrowth.",
+               "Use the values 'linear' or 'constant'."))
   }
 
   #Check that the number of cells entered does not provide a multiplet rate of >100%
@@ -754,6 +755,7 @@ power.sameReadDepth.restrictedDoublets<-function(nSamples,nCells,ct.freq,
                               mappingEfficiency=0.8,
                               multipletRate=7.67e-06,multipletFactor=1.82,
                               min.UMI.counts=3,perc.indiv.expr=0.5,
+                              cutoffVersion="absolute",
                               nGenes=21000,samplingMethod="quantiles",
                               multipletRateGrowth="linear",
                               sign.threshold=0.05, MTmethod="Bonferroni",
@@ -779,6 +781,7 @@ power.sameReadDepth.restrictedDoublets<-function(nSamples,nCells,ct.freq,
                                           mappingEfficiency,
                                           multipletRate,multipletFactor,
                                           min.UMI.counts,perc.indiv.expr,
+                                          cutoffVersion,
                                           nGenes,samplingMethod,
                                           multipletRateGrowth,
                                           sign.threshold,

diff --git a/man/power.sameReadDepth.restrictedDoublets.Rd b/man/power.sameReadDepth.restrictedDoublets.Rd