["ACTIVITY"]]<-CreateAssayObject(counts=activity.matrix)meta<-read.table("data/atac_v1_pbmc_10k_singlecell.csv",sep=",",header=TRUE,row.names=1,stringsAsFactors=FALSE)meta<-meta[colnames(pbmc.atac),]pbmc.atac<-AddMetaData(pbmc.atac,metadata=meta)#质控,过滤掉scATAC-seq数据中总数少于5K的...