# NOT RUN {
# load data classified previously (see SignacFast)
P <- readRDS("celltypes.rds")
S <- readRDS("pbmcs.rds")
# get edges from default assay from Seurat object
default.assay <- Seurat::DefaultAssay(S)
edges = S@graphs[[which(grepl(paste0(default.assay, "_nn"), names(S@graphs)))]]
# get distance matrix
D = CID.GetDistMat(edges)
# smooth labels
smoothed = CID.smooth(ac = P$CellTypes, dM = D[[1]])
# }
Run the code above in your browser using DataLab