單細胞UMAP圖可視化工具--PieParty

pieparty是一個可以在UMAP圖上展示pie chart的工具,圖如下,

install

直接下載张肾,解壓就能用芭析。

usage:

至少需要三個文件:

1. expression matrix,可以用Seurat生成吞瞪,

2. gene list? 隨意指定

3. tSNE或UMAP坐標馁启,可以用Seurat生成

4. cluster信息(可選)

先走Seurat標準流程:生成以上所需文件

library(Seurat)?

library(dplyr)?

data = Read10X('./filtered_feature_bc_matrix')?

seurat_obj = CreateSeuratObject(data)?

seurat_obj = NormalizeData(seurat_obj) %>% FindVariableFeatures() %>% ScaleData() %>% ? RunPCA() %>% RunUMAP(dims=1:50) %>% FindNeighbors(dims=1:50) %>% ? ? FindClusters(resolution=0.2)

Seurat object 將expression data 導出/寫入csv

write.csv(GetAssayData(seurat_obj, slot="data"), file = "expression_data.csv")

Seurat object 將UMAP/tSNE 信息導出/寫入csv

write.csv(Embeddings(seurat_obj, reduction = "umap"), file = "cell_embeddings.UMAP.csv")

圖1:

python3 PieParty.py -g expression_data.csv -c cell_embeddings.UMAP.csv -l gene_list

直接這樣運行非常耗時,可以選擇下面的方法芍秆,增加cluster信息惯疙,在cluster上畫piechart。

圖2:需要cluster信息

write.csv(seurat_obj@meta.data[,'seurat_clusters',drop=F], file = "cell_clusters.csv")

python3 PieParty.py -g expression_file.csv -c cell_coordinates.csv -l genelist1.csv genelist2.csv -gp True -cf clusters.csv -lc False



暫時用到的地方不多妖啥。

?著作權歸作者所有,轉載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末霉颠,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子荆虱,更是在濱河造成了極大的恐慌蒿偎,老刑警劉巖,帶你破解...
    沈念sama閱讀 217,406評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件怀读,死亡現(xiàn)場離奇詭異诉位,居然都是意外死亡,警方通過查閱死者的電腦和手機菜枷,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,732評論 3 393
  • 文/潘曉璐 我一進店門不从,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人犁跪,你說我怎么就攤上這事椿息。” “怎么了坷衍?”我有些...
    開封第一講書人閱讀 163,711評論 0 353
  • 文/不壞的土叔 我叫張陵寝优,是天一觀的道長。 經(jīng)常有香客問我枫耳,道長乏矾,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,380評論 1 293
  • 正文 為了忘掉前任迁杨,我火速辦了婚禮钻心,結果婚禮上,老公的妹妹穿的比我還像新娘铅协。我一直安慰自己捷沸,他們只是感情好,可當我...
    茶點故事閱讀 67,432評論 6 392
  • 文/花漫 我一把揭開白布狐史。 她就那樣靜靜地躺著痒给,像睡著了一般说墨。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上苍柏,一...
    開封第一講書人閱讀 51,301評論 1 301
  • 那天尼斧,我揣著相機與錄音,去河邊找鬼试吁。 笑死棺棵,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的熄捍。 我是一名探鬼主播烛恤,決...
    沈念sama閱讀 40,145評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼治唤!你這毒婦竟也來了棒动?” 一聲冷哼從身側響起糙申,我...
    開封第一講書人閱讀 39,008評論 0 276
  • 序言:老撾萬榮一對情侶失蹤宾添,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后柜裸,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體缕陕,經(jīng)...
    沈念sama閱讀 45,443評論 1 314
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,649評論 3 334
  • 正文 我和宋清朗相戀三年疙挺,在試婚紗的時候發(fā)現(xiàn)自己被綠了扛邑。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,795評論 1 347
  • 序言:一個原本活蹦亂跳的男人離奇死亡铐然,死狀恐怖蔬崩,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情搀暑,我是刑警寧澤沥阳,帶...
    沈念sama閱讀 35,501評論 5 345
  • 正文 年R本政府宣布,位于F島的核電站自点,受9級特大地震影響桐罕,放射性物質發(fā)生泄漏。R本人自食惡果不足惜桂敛,卻給世界環(huán)境...
    茶點故事閱讀 41,119評論 3 328
  • 文/蒙蒙 一功炮、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧术唬,春花似錦薪伏、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,731評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽博秫。三九已至,卻和暖如春眶掌,著一層夾襖步出監(jiān)牢的瞬間挡育,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,865評論 1 269
  • 我被黑心中介騙來泰國打工朴爬, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留即寒,地道東北人。 一個月前我還...
    沈念sama閱讀 47,899評論 2 370
  • 正文 我出身青樓召噩,卻偏偏與公主長得像母赵,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子具滴,可洞房花燭夜當晚...
    茶點故事閱讀 44,724評論 2 354

推薦閱讀更多精彩內(nèi)容