美文网首页scrna相关代码
【小脚本】统计各个cluster的细胞数和细胞比例

【小脚本】统计各个cluster的细胞数和细胞比例

作者: Geekero | 来源:发表于2020-06-29 17:23 被阅读0次
Count2Ratio<-function(sam){
  df<-(t(table(sam@meta.data[c('orig.ident', 'seurat_clusters')])))
  df<-as.data.frame(df)
  df<-dcast(df, seurat_clusters~orig.ident)
  write.table(df, file='./count.xls', sep='\t', quote=F, row.names=F)
  
  for(col in colnames(df)){
    if(col=='seurat_clusters'){
      next
    }
    print(col)
    # all_sum = sum(df[col])
    # print(df[col])
    df[col]<-prop.table(df[col])
  }
  write.table(df, file='./ratio.xls', sep='\t', quote=F, row.names=F)
}

相关文章

网友评论

    本文标题:【小脚本】统计各个cluster的细胞数和细胞比例

    本文链接:https://www.haomeiwen.com/subject/iehafktx.html