# #——加载方案,消息= FALSE ---------------------------------------------- 图书馆(LinTInd) # #——消息= FALSE ------------------------------------------------------------ 数据< -paste0(系统。file("extdata"，package = 'LinTInd')，"/CB_UMI") fafile<-paste0(system. txt)file("extdata"，package = 'LinTInd')，"/V3.fasta")file("extdata"，package = 'LinTInd')，"/V3.cutSites")文件(“extdata”,包=“LinTInd”),< -read.table / celltype.tsv)数据(数据,9月=“t \”,头= TRUE) ref < -ReadFasta (fafile) cutsite < -read.table (cutsite col.names = c(“indx”、“开始”、“结束”))celltype < -read.table (celltype头= TRUE, stringsAsFactors = FALSE ) ## ----------------------------------------------------------------------------- 头(数据,3)ref cutsite头(celltype, 3) # #——找到indels并生成数组形式的字符串,message=FALSE--------------- scarinfo<-FindIndel(data=data,scarfull=ref,scar=cutsite,indel.coverage="All"，type="test"，cln=1) scarinfo<-IndelForm(scarinfo,cln=1) ## ----IndelIdents，消息= FALSE ----------------------------------------------- cellsinfo < -IndelIdents (scarinfo method.use =“umi.num cln = 1) # #——IndelPlot ---------------------------------------------------------------- IndelPlot (cellsinfo = cellsinfo) # #——TagProcess --------------------------------------------------------------- 标记< -TagProcess (cellsinfo $ info,细胞= celltype) # #——TagDist ------------------------------------------------------------------ tag_dist = TagDist(标签,方法=“Jaccard”)tag_dist # #——BuildTree ---------------------------------------------------------------- treeinfo < -BuildTree # #——PlotTree(标签) ----------------------------------------------------------------- plotinfo < -PlotTree (treeinfo = treeinfo数据。提取= " TRUE ",注释= " TRUE ") plotinfo $ p # #——sessionInfo回声= TRUE --------------------------------------------------- sessionInfo ()