SCEVAN icon indicating copy to clipboard operation
SCEVAN copied to clipboard

count_mtx choose for integrated Seurat object data

Open Lawwwwo opened this issue 1 year ago • 1 comments

Hi! I am using SCEVAN to analyze an integrated Seurat object. I have no idea which count_mtx data should I choose from my Seurat object. There are three versions: RNA, SCT(transformed version), and integrated as the figure inserted below. Screenshot 2023-04-23 at 2 56 21 PM When I choose RNA's matrix as input, it can run normally. But when I choose integrated it shows errors: The number of elements of the gene set is less than the minimum allowed. Error in hclust(d, method = "ward.D2") : NA/NaN/Inf in foreign function call (arg 10) Calls: <Anonymous> ... classifyTumorCells -> removeSyntheticBaseline -> hclust In addition: Warning messages: 1: In sqrt(count_mtx_proc) : NaNs produced 2: In sqrt(count_mtx_proc + 1) : NaNs produced Execution halted And screen output stops here: [1] " raw data - genes: 5000 cells: 13546" [1] "1) Filter: cells > 200 genes" [1] "low data quality" [1] "2) Filter: genes > 5% of cells" [1] "4935 genes past filtering" [1] "3) Annotations gene coordinates" [1] "found 0 confident non malignant cells" [1] "4797 genes annotated" [1] "4) Filter: genes involved in the cell cycle" [1] "4465 genes past filtering " [1] "5) Filter: cells > 5genes per chromosome " [1] "6) Log Freeman Turkey transformation" [1] "A total of 13546 cells, 4465 genes after preprocessing" [1] "7) Measuring baselines (pure tumor - synthetic normal cells)" And here is the code with using integrated code: Integrationdata_123M_CNV <- SCEVAN::pipelineCNA(Integrationdata_123M@assays[["integrated"]]@data, sample = "BreastCancerIntegration", par_cores = 80, SUBCLONES = FALSE, plotTree = FALSE) Can you help me figure out this problem? And which input should I choose?

Lawwwwo avatar Apr 23 '23 07:04 Lawwwwo

Hi Antonio! I just tried another integration method with my data. The same problem occurs which stops at step 7) Measuring baseline: Screenshot 2023-04-24 at 11 00 09 AM

Lawwwwo avatar Apr 24 '23 03:04 Lawwwwo