SCEVAN icon indicating copy to clipboard operation
SCEVAN copied to clipboard

Error in Step9

Open Drjingluo opened this issue 1 year ago • 0 comments

Hi, I got the list file from my Seurat subject and perfomed the following script. can you help me? bd8e5aa645372348f7b3c915807fc91

results <- SCEVAN::multiSampleComparisonClonalCN(count_mtx_list, analysisName = "all", organism = "human" , par_cores = 48) [1] " raw data - genes: 33691 cells: 1351" [1] "1) Filter: cells > 200 genes" [1] "low data quality" [1] "2) Filter: genes > 5% of cells" [1] "6734 genes past filtering" [1] "3) Annotations gene coordinates" [1] "found 30 confident non malignant cells" [1] "6455 genes annotated" [1] "4) Filter: genes involved in the cell cycle" [1] "6129 genes past filtering " [1] "5) Filter: cells > 5genes per chromosome " [1] "6) Log Freeman Turkey transformation" [1] "A total of 1234 cells, 6129 genes after preprocessing" [1] "7) Measuring baselines (confident normal cells)" [1] "8) Smoothing data" [1] "9) Segmentation (VegaMC)" [1] "10) Adjust baseline" [1] "11) plot heatmap" [1] "found 688 tumor cells" [1] "time classify tumor cells: 1.22660799423854" [1] " raw data - genes: 33691 cells: 3000" [1] "1) Filter: cells > 200 genes" [1] "low data quality" [1] "2) Filter: genes > 5% of cells" [1] "7392 genes past filtering" [1] "3) Annotations gene coordinates" [1] "found 30 confident non malignant cells" [1] "7086 genes annotated" [1] "4) Filter: genes involved in the cell cycle" [1] "6705 genes past filtering " [1] "5) Filter: cells > 5genes per chromosome " [1] "6) Log Freeman Turkey transformation" [1] "A total of 2928 cells, 6705 genes after preprocessing" [1] "7) Measuring baselines (confident normal cells)" [1] "8) Smoothing data" [1] "9) Segmentation (VegaMC)" Error in data.frame(..., check.names = FALSE) : arguments imply differing number of rows: 6705, 5588 In addition: Warning messages: 1: In parallel::mclapply(1:ncol(count_mtx_relat), nonLinSmooth, mc.cores = par_cores) : scheduled cores 1, 5, 15, 18, 22, 25, 26, 32 did not deliver results, all values of the jobs will be affected 2: In matrix(unlist(test.mc), ncol = ncol(count_mtx_relat), byrow = FALSE) : data length [16360200] is not a sub-multiple or multiple of the number of rows [5588]

Drjingluo avatar Nov 13 '22 05:11 Drjingluo