单细胞转录组降维聚类分群过滤基因和过滤细胞的区别

文章预览

我们给了大家一个常规的可以批量导入全部的10x单细胞表达量矩阵文件夹的代码： # 参考：https://mp.weixin.qq.com/s/tw7lygmGDAbpzMTx57VvFw dir= 'GSE163558_RAW' samples=list.files( dir  ) #samples=samples[1:2] sceList = lapply(samples, function (pro){     # pro=samples[1]    print(pro)     tmp = Read10X(file.path(dir,pro ) )    print(dim(tmp))    if (length(tmp)== 2 ){     ct = tmp[[ 1 ]]    } else {ct = tmp}   print(dim(ct))   sce =CreateSeuratObject(counts =  ct ,                           project =  gsub( 'GSM[0-9]*_' , '' ,pro)  ,                           min.cells =  5 ,                           min.features =  300  )    return (sce) })  里面有两个很重要的参数，  min.cells = 5,  和  min.features = 300 ，它们都是为了过滤我们的单细胞表达量矩阵，但是前者是过滤基因，后者 ………………………………