聚类教程：cisTopic——来自单细胞表观基因组学数据的顺式调控主题的概率建模

使用cisTopic对scATAC-seq数据聚类：

rm(list = ls())
library(dplyr)
library(patchwork)
library(reticulate)
library(data.table)
library(cisTopic)
mix_atac<-read.table('HT_p0_atac.tsv',sep='\t') #输入文件是cell(列)×peak(行)的count.注意：peak的格式必须是 chr染色体：start-end,不是这种格式要改一下
cisTopicObject <- createcisTopicObject(mix_atac,min.cells = 10