单细胞分析,统计各个cluster的细胞数和细胞比例
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的细胞数和细胞比例
本文标题:单细胞分析,统计各个cluster的细胞数和细胞比例
本文链接:https://www.haomeiwen.com/subject/ctafrktx.html
网友评论