We need to fix the number of cells that we use as the internal white list presented in
/Solo.out/GeneFull/filtered
the test has shown that the numbers of barcodes in m1 and gene expression are not compatible. The problem might comes from this part:
# Define paths
expression_matrix_dir <- paste0(data_dir, "/matrix.mtx")
expression_barcodes_dir <- paste0(data_dir, "/barcodes.tsv")
expression_features_dir <- paste0(data_dir, "/features.tsv")
filtered_barcodes_dir <- paste0(expression_dir, "../GeneFull/filtered/barcodes.tsv") # For barcode filtration
We need to fix the number of cells that we use as the internal white list presented in
/Solo.out/GeneFull/filteredthe test has shown that the numbers of barcodes in m1 and gene expression are not compatible. The problem might comes from this part: