Can't build metacells

[f6v]

Hi, I'm trying to run the workflow, but it fails at building the metacells with the following error(truncated):

Error in `$<-.data.frame`(`*tmp*`, "sz_cor", value = c(`5S_rRNA` = 0.005, : replacement has 28936 rows, data has 21700
Traceback:

1. build_metacells()
2. mel_build_master_mat(force_new = rebuild_data)
3. mcell_add_gene_stat(all_id, all_id, force = force_new)
4. scm_gene_stat(mat_id, downsample_n = downsample_n)
5. `$<-`(`*tmp*`, "sz_cor", value = c(`5S_rRNA` = 0.005, `7SK` = 0.055, 
 . A1BG = 0.006, `A1BG-AS1` = 0.045, A1CF = 0.001, A2M = 0.255, 
 . `A2M-AS1` = 0.019, A2ML1 = 0.004, A2MP1 = 0.029, A4GALT = 0.061, 
 . AAAS = 0.144, AACS = 0.103, AACSP1 = 0.028, AADAC = -0.002, AADAT = 0.096, 
 . AAED1 = 0.098, AAGAB = 0.098, AAK1 = 0.074, AAMDC = 0.102, AAMP = 0.172, 
 . AAR2 = 0.134, AARD = 0.014, AARS = 0.255, AARS2 = 0.091, AARSD1 = 0.127, 
 . AASDH = 0.084, AASDHPPT = 0.166, AASS = 0.038, AATF = 0.167, 
 . AATK = 0.121, AB019441.29 = 0.42, ABAT = 0.039, `ABC7-42389800N19.1` = 0.047, 
 . `ABC7-43041300I9.1` = 0.03, ABCA1 = 0.11, ABCA10 = 0.022, ABCA11P = 0.011, 

Any idea why that could happen?

[yaniv-el]

Hi,

This is caused by failure of some threads during the calculation of some gene stats, and should be fixed in the metacell package.

I've contacted the current metacell developers so this should be fixed soon. I'll update here when it is fixed.

Best, Yaniv