cell_marker_data <- vroom::vroom('http://bio-bigdata.hrbmu.edu.cn/CellMarker/download/Human_cell_markers.txt') ## instead of `cellName`, users can use other features (e.g. `cancerType`) cells <- cell_marker_data %>% dplyr::select(cellName, geneID) %>% dplyr::mutate(geneID = strsplit(geneID, ', ')) %>% tidyr::unnest()
图片alt
x <- enricher(gene, TERM2GENE = cells) setReadable(x,OrgDb = "org.Hs.eg.db",keyType = "ENTREZID") |> head()
y <- GSEA(geneList, TERM2GENE = cells) head(y)