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)
}
【小腳本】統(tǒng)計各個cluster的細胞數(shù)和細胞比例
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
- 文/潘曉璐 我一進店門罚渐，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人驯妄，你說我怎么就攤上這事荷并。” “怎么了青扔？”我有些...
- 文/不壞的土叔 我叫張陵源织，是天一觀的道長翩伪。 經(jīng)常有香客問我，道長谈息，這世上最難降的妖魔是什么缘屹？ 我笑而不...
- 文/花漫 我一把揭開白布盲憎。 她就那樣靜靜地躺著嗅骄，像睡著了一般。 火紅的嫁衣襯著肌膚如雪饼疙。 梳的紋絲不亂的頭發(fā)上溺森，一...
- 文/蒼蘭香墨 我猛地睜開眼卷要，長吁一口氣：“原來是場噩夢啊……” “哼渣聚！你這毒婦竟也來了？” 一聲冷哼從身側響起僧叉，我...
- 正文 年R本政府宣布壹置，位于F島的核電站，受9級特大地震影響表谊，放射性物質發(fā)生泄漏钞护。R本人自食惡果不足惜，卻給世界環(huán)境...
- 文/蒙蒙 一爆办、第九天 我趴在偏房一處隱蔽的房頂上張望难咕。 院中可真熱鬧，春花似錦距辆、人聲如沸余佃。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽爆土。三九已至，卻和暖如春诸蚕，著一層夾襖步出監(jiān)牢的瞬間步势，已是汗流浹背。 一陣腳步聲響...
