Skip to content
Snippets Groups Projects
Commit c88d11c4 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Fix cca hvg intersect

parent f372f6c0
Branches
1 merge request!2Merge develop into master
......@@ -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
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment