ggplot2繪制GO注釋二級(jí)分類統(tǒng)計(jì)圖

>library(ggplot2) #載入ggplot2

>data=read.table("GO.count.txt",,header=T,sep="\t") #讀取數(shù)據(jù)

>dorder = factor(as.integer(rownames(data)),labels=data$Description) #對(duì)數(shù)據(jù)按Description排序

>ggplot(data,aes(x=Description,y=Counts,fill=Class)) #定義X軸妹沙,Y軸的數(shù)據(jù)和顏色填充
+geom_bar(stat="identity",position=position_dodge(0.7),width=0.5,aes(x=dorder)) #定義柱形圖的寬度和間距
+coord_flip() #轉(zhuǎn)換橫縱坐標(biāo)
+scale_y_log10(breaks = c(1,10,100,1000)) #Y軸按log10排序,展示1哗讥,10隅肥,100姆坚,1000刻度
+scale_fill_discrete(name="Ontology") #修改legend的tittle
+theme(panel.background = element_rect(fill = "transparent",colour = NA)) #清楚背景顏色
+xlab("Term") #修改X軸標(biāo)簽

GO_classification_L2.jpg

根據(jù)鏈接的指導(dǎo)返吻,結(jié)合自己的數(shù)據(jù)實(shí)現(xiàn)上圖,代碼略作修改去扣。
用R畫GO注釋二級(jí)分類統(tǒng)計(jì)圖 | 生信菜鳥團(tuán) http://www.bio-info-trainee.com/771.html

這是我自己的數(shù)據(jù)
GOid Description Class Counts TotalGOid TotalGeneAnnotated
GO:0045735 nutrient reservoir activity molecular_function 39 3384 9842
GO:0060089 molecular transducer activity molecular_function 42 3384 9842
GO:0004871 signal transducer activity molecular_function 56 3384 9842
GO:0016209 antioxidant activity molecular_function 64 3384 9842
GO:0098772 molecular function regulator molecular_function 95 3384 9842
GO:0005198 structural molecule activity molecular_function 166 3384 9842
GO:0140110 transcription regulator activity molecular_function 180 3384 9842

GO:0005215 transporter activity molecular_function 478 3384 9842
GO:0003824 catalytic activity molecular_function 3805 3384 9842
GO:0005488 binding molecular_function 5401 3384 9842
GO:0019012 virion cellular_component 1 3384 9842
GO:0099080 supramolecular complex cellular_component 9 3384 9842
GO:0044421 extracellular region part cellular_component 35 3384 9842
GO:0031974 membrane-enclosed lumen cellular_component 54 3384 9842
GO:0005576 extracellular region cellular_component 81 3384 9842
GO:0044422 organelle part cellular_component 284 3384 9842
GO:0032991 macromolecular complex cellular_component 517 3384 9842
GO:0043226 organelle cellular_component 649 3384 9842
GO:0044425 membrane part cellular_component 702 3384 9842
GO:0005623 cell cellular_component 996 3384 9842
GO:0044464 cell part cellular_component 996 3384 9842
GO:0016020 membrane cellular_component 1219 3384 9842
GO:0022610 biological adhesion biological_process 1 3384 9842
GO:0008283 cell proliferation biological_process 1 3384 9842
GO:0098754 detoxification biological_process 1 3384 9842
GO:0019740 nitrogen utilization biological_process 1 3384 9842
GO:0048511 rhythmic process biological_process 1 3384 9842
GO:0002376 immune system process biological_process 5 3384 9842
GO:0015976 carbon utilization biological_process 6 3384 9842
GO:0040007 growth biological_process 13 3384 9842
GO:0051704 multi-organism process biological_process 43 3384 9842
GO:0000003 reproduction biological_process 52 3384 9842
GO:0022414 reproductive process biological_process 52 3384 9842
GO:0032502 developmental process biological_process 59 3384 9842
GO:0032501 multicellular organismal process biological_process 85 3384 9842
GO:0023052 signaling biological_process 183 3384 9842
GO:0071840 cellular component organization or biogenesis biological_process 265 3384 9842
GO:0050896 response to stimulus biological_process 501 3384 9842
GO:0065007 biological regulation biological_process 664 3384 9842
GO:0051179 localization biological_process 756 3384 9842
GO:0009987 cellular process biological_process 3122 3384 9842
GO:0008152 metabolic process biological_process 4225 3384 9842

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末风喇,一起剝皮案震驚了整個(gè)濱河市宁改,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌魂莫,老刑警劉巖还蹲,帶你破解...
    沈念sama閱讀 218,682評(píng)論 6 507
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡谜喊,警方通過查閱死者的電腦和手機(jī)潭兽,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,277評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來锅论,“玉大人讼溺,你說我怎么就攤上這事∽钜祝” “怎么了怒坯?”我有些...
    開封第一講書人閱讀 165,083評(píng)論 0 355
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)藻懒。 經(jīng)常有香客問我剔猿,道長(zhǎng),這世上最難降的妖魔是什么嬉荆? 我笑而不...
    開封第一講書人閱讀 58,763評(píng)論 1 295
  • 正文 為了忘掉前任归敬,我火速辦了婚禮,結(jié)果婚禮上鄙早,老公的妹妹穿的比我還像新娘汪茧。我一直安慰自己,他們只是感情好限番,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,785評(píng)論 6 392
  • 文/花漫 我一把揭開白布舱污。 她就那樣靜靜地躺著,像睡著了一般弥虐。 火紅的嫁衣襯著肌膚如雪扩灯。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,624評(píng)論 1 305
  • 那天霜瘪,我揣著相機(jī)與錄音珠插,去河邊找鬼。 笑死颖对,一個(gè)胖子當(dāng)著我的面吹牛捻撑,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播惜互,決...
    沈念sama閱讀 40,358評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼布讹,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來了训堆?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,261評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤白嘁,失蹤者是張志新(化名)和其女友劉穎坑鱼,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,722評(píng)論 1 315
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡鲁沥,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,900評(píng)論 3 336
  • 正文 我和宋清朗相戀三年呼股,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片画恰。...
    茶點(diǎn)故事閱讀 40,030評(píng)論 1 350
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡彭谁,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出允扇,到底是詐尸還是另有隱情缠局,我是刑警寧澤,帶...
    沈念sama閱讀 35,737評(píng)論 5 346
  • 正文 年R本政府宣布考润,位于F島的核電站狭园,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏糊治。R本人自食惡果不足惜唱矛,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,360評(píng)論 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望井辜。 院中可真熱鬧绎谦,春花似錦、人聲如沸粥脚。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,941評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)阿逃。三九已至铭拧,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間恃锉,已是汗流浹背搀菩。 一陣腳步聲響...
    開封第一講書人閱讀 33,057評(píng)論 1 270
  • 我被黑心中介騙來泰國(guó)打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留破托,地道東北人肪跋。 一個(gè)月前我還...
    沈念sama閱讀 48,237評(píng)論 3 371
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像土砂,于是被迫代替她去往敵國(guó)和親州既。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,976評(píng)論 2 355