search

chuiqin / irGSEA

The integration of single cell rank-based gene set enrichment analysis
Other
106 stars 17 forks source link

Error in ridgeplot #16

Open lkyhfx opened 1 year ago

lkyhfx commented 1 year ago 
Idents(sObj) <- sObj$monaco.main
sObj <- irGSEA.score(object = sObj, assay = "RNA", 
                             slot = "data", seeds = 123, ncores = 1,
                             min.cells = 3, min.feature = 0,
                             custom = F, geneset = NULL, msigdb = T, 
                             species = "Homo sapiens", category = "H",  
                             subcategory = NULL, geneid = "symbol",
                             method = c("UCell"),
                             aucell.MaxRank = NULL, ucell.MaxRank = NULL, 
                             kcdf = 'Gaussian')

ridgeplot <- irGSEA.ridgeplot(object = sObj,
                              method = "UCell",
                              show.geneset = "HALLMARK-INFLAMMATORY-RESPONSE")
ridgeplot

这是我的代码，当使用 seurat 计算的 cluster 作为 idents 时出图是正常的，但是 idents 设定为细胞类型后出图报错

Picking joint bandwidth of 0.00332
Error in `[[<-.data.frame`(`*tmp*`, i, value = c("#5050FFFF", "#CE3D32FF",  : 
  替换数据里有10行，但数据有11