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

Specify assay and slot for QUSAGE

parent afb7560d
2 merge requests!6Develop,!5Refactor
......@@ -824,7 +824,7 @@ scQuSAGE <- function(sc10x,gs,save=FALSE,type,id,ds=0,nm="pops",print="tsne"){
}
cell.sample <- c(cell.sample,cell)
}
data <- as.data.frame(as.matrix(GetAssayData(sc10x[,colnames(sc10x) %in% cell.sample])))
data <- as.data.frame(as.matrix(GetAssayData(sc10x[,colnames(sc10x) %in% cell.sample],assay="SCT",slot="scale.data")))
labels <- labels[colnames(sc10x) %in% cell.sample]
#groups <- sort(unique(labels))
groups <- paste0("Cluster_",levels(sc10x@active.ident))
......
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