From c88d11c4e6e59094b9e0ed93a93032b9f3a7c5c7 Mon Sep 17 00:00:00 2001 From: "Gervaise H. Henry" <gervaise.henry@utsouthwestern.edu> Date: Mon, 1 Oct 2018 12:02:28 -0500 Subject: [PATCH] Fix cca hvg intersect --- r.scripts/sc-TissueMapper.R | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/r.scripts/sc-TissueMapper.R b/r.scripts/sc-TissueMapper.R index 1f69cc1..ea54e0f 100644 --- a/r.scripts/sc-TissueMapper.R +++ b/r.scripts/sc-TissueMapper.R @@ -1097,9 +1097,9 @@ sc3CCA <- function(sc10x.1,sc10x.2,sc10x.3,nm.1="D17",nm.2="D27",nm.3="D35",cc= genes.hvg.2 <- sc10x.2@var.genes genes.hvg.3 <- sc10x.3@var.genes genes.hvg.Comb <- unique(c(genes.hvg.1,genes.hvg.2,genes.hvg.3)) - genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.1@scale.data)) - genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.2@scale.data)) - genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.3@scale.data)) + genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.1@data)) + genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.2@data)) + genes.hvg.Comb <- intersect(genes.hvg.Comb,rownames(sc10x.3@data)) sc10x.1@meta.data$patient <- nm.1 sc10x.2@meta.data$patient <- nm.2 -- GitLab