CytoTRACE | 推断拟时序起点

  1. R包下载网站:https://cytotrace.stanford.edu/CytoTRACE_0.3.3.tar.gz
#install
devtools::install_local("yourpath/CytoTRACE_0.3.3.tar.gz")
#load
library(CytoTRACE)

sva/ccaPP/HiClimR/ncdf4这几个包没有的用BiocManager安装一下,安装HiClimR/ncdf4的时候可能会报错(dependency ‘ncdf4’ is not available for package ‘HiClimR’),在linux终端安装netcdf就好了。

sudo apt-get install libnetcdf-dev
  1. 推断拟时序起点,可视化结果储存在本地
####提取表型文件
table(sc$celltype)
phe <- sc$celltype
phe = as.character(phe)
names(phe) <- rownames(sc@meta.data)
 
####提取表达矩阵
mat <- as.matrix(sc@assays$RNA@counts)
mat[1:4,1:4]
results <- CytoTRACE(mat = mat)
 
#可视化
plotCytoGenes(results, numOfGenes = 10)
plotCytoTRACE(results, phenotype = phe)

参考:CytoTRACE (stanford.edu)

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容