Error when using ExampleSimSCData

[micw42]

Hi, I installed the developer version of SCnorm and am trying it on the ExampleSimSCData:

data(ExampleSimSCData)
Conditions = rep(c(1, 2), each= 45)

DataNorm <- SCnorm(Data = ExampleSimSCData, Conditions = Conditions, 
                     PrintProgressPlots = TRUE, FilterCellNum = 10, 
                   reportSF = TRUE)

However, I get this error:

Setting up parallel computation using 95 cores

Gene filter is applied within each condition.

288 genes in condition 1 will not be included in the normalization due to 
             the specified filter criteria.

314 genes in condition 2 will not be included in the normalization due to 
             the specified filter criteria.

A list of these genes can be accessed in output, 
    see vignette for example.

Finding K for Condition 1

Trying K = 1

Warning message in mclapply(X = seq_len(nrow(Grid)), FUN = GetTD, InputData = list(O, :
"scheduled cores 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 18, 19, 40, 41, 45, 51, 52, 53, 57, 58, 59, 64, 65, 69, 73, 74, 75, 79 did not deliver results, all values of the jobs will be affected"
Warning message in mclapply(X = seq_len(nrow(Grid)), FUN = GetTD, InputData = list(O, :
"scheduled cores 78, 77, 67, 66, 62, 55, 54, 46, 44, 39, 37, 32, 31, 27, 23, 21, 17, 16, 12 encountered errors in user code, all values of the jobs will be affected"
Error in PEAK - AllIter: non-numeric argument to binary operator
Traceback:

1. SCnorm(Data = ExampleSimSCData, Conditions = Conditions, PrintProgressPlots = TRUE, 
 .     FilterCellNum = 10, reportSF = TRUE)
2. lapply(seq_along(Levels), function(x) {
 .     normWrapper(Data = DataList[[x]], SeqDepth = SeqDepthList[[x]], 
 .         Slopes = SlopesList[[x]], CondNum = Levels[x], PrintProgressPlots = PrintProgressPlots, 
 .         PropToUse = PropToUse, Tau = Tau, Thresh = Thresh, ditherCounts = ditherCounts)
 . })
3. FUN(X[[i]], ...)
4. normWrapper(Data = DataList[[x]], SeqDepth = SeqDepthList[[x]], 
 .     Slopes = SlopesList[[x]], CondNum = Levels[x], PrintProgressPlots = PrintProgressPlots, 
 .     PropToUse = PropToUse, Tau = Tau, Thresh = Thresh, ditherCounts = ditherCounts)
5. SCnormFit(Data = Data, SeqDepth = SeqDepth, Slopes = Slopes, 
 .     K = i, PropToUse = PropToUse, Tau = Tau, ditherCounts = ditherCounts)
6. Grid[which.min(abs(PEAK - AllIter)), 2]
7. `[.data.frame`(Grid, which.min(abs(PEAK - AllIter)), 2)
8. which.min(abs(PEAK - AllIter))

Can you please help me troubleshoot this? Thank you so much!

[rhondabacher]

Thanks for the message. Are you able to use a more recent version of the package?