Hi, I am running the numbat on scRNA-seq data, but I am getting error using run_numbat().
Error
Details of the log.txt
INFO [2025-06-29 15:22:32]
numbat version: 1.5.0
scistreer version: 1.2.0
hahmmr version: 1.0.0
Running under parameters:
t = 1e-05
alpha = 1e-04
gamma = 20
min_cells = 50
init_k = 3
max_cost = 3637.8
n_cut = 0
max_iter = 2
max_nni = 100
min_depth = 0
use_loh = auto
segs_loh = None
call_clonal_loh = FALSE
segs_consensus_fix = None
multi_allelic = TRUE
min_LLR = 5
min_overlap = 0.45
max_entropy = 0.6
skip_nj = FALSE
diploid_chroms = None
ncores = 38
ncores_nni = 38
common_diploid = TRUE
tau = 0.3
check_convergence = FALSE
plot = TRUE
genome = hg38
Input metrics:
12126 cells
INFO [2025-06-29 15:22:34] Mem used: 4.39Gb
INFO [2025-06-29 15:22:44] Approximating initial clusters using smoothed expression ..
INFO [2025-06-29 15:22:46] Mem used: 4.38Gb
INFO [2025-06-29 16:00:36] running hclust...
INFO [2025-06-29 16:19:47] Iteration 1
INFO [2025-06-29 16:19:53] Mem used: 17.2Gb
INFO [2025-06-29 16:20:20] Expression noise level (MSE): low (0.34).
INFO [2025-06-29 16:20:20] Running HMMs on 5 cell groups..
INFO [2025-06-29 16:21:00] diploid regions: 1a,2a,3a,3c,4a,5a,6a,7a,8a,9a,10a,10c,11a,11c,12a,12c,12e,13a,14a,14c,15a,16a,17a,18a,19a,20a,21a,22a
INFO [2025-06-29 16:24:14] Running HMMs on 3 cell groups..
INFO [2025-06-29 16:24:40] diploid regions: 1a,2a,3a,4a,5a,6a,7a,8a,9a,10a,10c,11a,12a,12c,12e,13a,14a,15a,16a,17a,18a,19a,20a,21a,22a
INFO [2025-06-29 16:26:18] Testing for multi-allelic CNVs ..
INFO [2025-06-29 16:26:18] 0 multi-allelic CNVs found:
INFO [2025-06-29 16:26:18] Evaluating CNV per cell ..
INFO [2025-06-29 16:26:21] Mem used: 7.67Gb
INFO [2025-06-29 17:11:51] All cells succeeded
INFO [2025-06-29 17:12:14] Expanding allelic states..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:15] Building phylogeny ..
INFO [2025-06-29 17:12:17] Mem used: 7.71Gb
INFO [2025-06-29 17:12:17] No CNV remains after filtering by entropy in single cells. Consider increasing max_entropy.
Error in h(simpleError(msg, call)) : glue failed in formatter_glue on: 'try-error' chr "Error in optim(fn = function(w) { : non-finite value supplied by optim\n" - attr(, "condition")=List of 2 ..$ message: chr "non-finite value supplied by optim" ..$ call : language optim(fn = function(w) { w = w/sum(w) ... ..- attr(, "class")= chr [1:3] "simpleError" "error" "condition" Raw error message: Expecting '}' Please consider using another log_formatter or skip_formatter on strings with curly braces.
Firstly, I checked my count_mat. And it seems no problem.
Then I checked the lambdas_ref and no NaN is produced.
Details of the lambdas_ref

And I tried to change my ref_internal using ref_internal[ref_internal == 0] <- 1e-8 and ref_internal[ref_internal == 0] <-1e-5, but it doesn't wok.
Could you tell me how to solve this problem? Thank you!
Elbert
Hi, I am running the numbat on scRNA-seq data, but I am getting error using run_numbat().
Error
Details of the log.txt
INFO [2025-06-29 15:22:32]
numbat version: 1.5.0
scistreer version: 1.2.0
hahmmr version: 1.0.0
Running under parameters:
t = 1e-05
alpha = 1e-04
gamma = 20
min_cells = 50
init_k = 3
max_cost = 3637.8
n_cut = 0
max_iter = 2
max_nni = 100
min_depth = 0
use_loh = auto
segs_loh = None
call_clonal_loh = FALSE
segs_consensus_fix = None
multi_allelic = TRUE
min_LLR = 5
min_overlap = 0.45
max_entropy = 0.6
skip_nj = FALSE
diploid_chroms = None
ncores = 38
ncores_nni = 38
common_diploid = TRUE
tau = 0.3
check_convergence = FALSE
plot = TRUE
genome = hg38
Input metrics:
12126 cells
INFO [2025-06-29 15:22:34] Mem used: 4.39Gb
INFO [2025-06-29 15:22:44] Approximating initial clusters using smoothed expression ..
INFO [2025-06-29 15:22:46] Mem used: 4.38Gb
INFO [2025-06-29 16:00:36] running hclust...
INFO [2025-06-29 16:19:47] Iteration 1
INFO [2025-06-29 16:19:53] Mem used: 17.2Gb
INFO [2025-06-29 16:20:20] Expression noise level (MSE): low (0.34).
INFO [2025-06-29 16:20:20] Running HMMs on 5 cell groups..
INFO [2025-06-29 16:21:00] diploid regions: 1a,2a,3a,3c,4a,5a,6a,7a,8a,9a,10a,10c,11a,11c,12a,12c,12e,13a,14a,14c,15a,16a,17a,18a,19a,20a,21a,22a
INFO [2025-06-29 16:24:14] Running HMMs on 3 cell groups..
INFO [2025-06-29 16:24:40] diploid regions: 1a,2a,3a,4a,5a,6a,7a,8a,9a,10a,10c,11a,12a,12c,12e,13a,14a,15a,16a,17a,18a,19a,20a,21a,22a
INFO [2025-06-29 16:26:18] Testing for multi-allelic CNVs ..
INFO [2025-06-29 16:26:18] 0 multi-allelic CNVs found:
INFO [2025-06-29 16:26:18] Evaluating CNV per cell ..
INFO [2025-06-29 16:26:21] Mem used: 7.67Gb
INFO [2025-06-29 17:11:51] All cells succeeded
INFO [2025-06-29 17:12:14] Expanding allelic states..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:14] No multi-allelic CNVs, skipping ..
INFO [2025-06-29 17:12:15] Building phylogeny ..
INFO [2025-06-29 17:12:17] Mem used: 7.71Gb
INFO [2025-06-29 17:12:17] No CNV remains after filtering by entropy in single cells. Consider increasing max_entropy.
Error in h(simpleError(msg, call)) :
gluefailed informatter_glueon: 'try-error' chr "Error in optim(fn = function(w) { : non-finite value supplied by optim\n" - attr(, "condition")=List of 2 ..$ message: chr "non-finite value supplied by optim" ..$ call : language optim(fn = function(w) { w = w/sum(w) ... ..- attr(, "class")= chr [1:3] "simpleError" "error" "condition" Raw error message: Expecting '}' Please consider using anotherlog_formatterorskip_formatteron strings with curly braces.Firstly, I checked my count_mat. And it seems no problem.
Then I checked the lambdas_ref and no NaN is produced.
Details of the lambdas_ref
And I tried to change my ref_internal using ref_internal[ref_internal == 0] <- 1e-8 and ref_internal[ref_internal == 0] <-1e-5, but it doesn't wok.
Could you tell me how to solve this problem? Thank you!
Elbert