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的细胞数和细胞比例
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
