SCP icon indicating copy to clipboard operation
SCP copied to clipboard

Error in RunCellQC(srt = seu_obj)

Open WANGROUOYI opened this issue 1 year ago • 3 comments

May I ask how to solve this problem? Error in RunCellQC(srt = seu_obj) : The db_rate is equal to or greater than 1!

WANGROUOYI avatar Sep 27 '23 14:09 WANGROUOYI

The default doublet rate is determined based on the number of input cells, using the formula db_rate = ncol(srt) / 1000 * 0.01. If the cell count surpasses 100,000, this would typically not be a standard single experiment. It appears that you might be using a merged Seurat object. I recommend installing the latest version of SCP and running RunCellQC with the batch parameter along with db_rate or db_coefficient manually (for instance, db_coefficient is 0.01 for 10X Genomics and 0.005 for Singleron GEXSCOPE).

zhanghao-njmu avatar Oct 18 '23 07:10 zhanghao-njmu

I ran the RunCellQC function for my data, the percent.mito is the same value for all cells. The code: RunCellQC( srt = sce, UMI_threshold = 500, gene_threshold = 400, mito_threshold = 10, return_filtered = F )

chenhy-lab avatar Nov 10 '23 03:11 chenhy-lab

percent.mito <- srt[[paste0(c("percent.mito", sp), collapse = ".")]] <- PercentageFeatureSet(object = srt, assay = assay, pattern = paste0("(", paste0("^", prefix, "-", mito_pattern), ")", collapse = "|"), features = mito_gene)[[1]] percent.ribo <- srt[[paste0(c("percent.ribo", sp), collapse = ".")]] <- PercentageFeatureSet(object = srt, assay = assay, pattern = paste0("(", paste0("^", prefix, "-", ribo_pattern), ")", collapse = "|"), features = ribo_gene)[[1]] percent.genome <- srt[[paste0(c("percent.genome", sp), collapse = ".")]] <- PercentageFeatureSet(object = srt, assay = assay, pattern = paste0("^", prefix))[[1]]

Here [[1]] seems to be the cause of the above error.

chenhy-lab avatar Nov 10 '23 06:11 chenhy-lab