Rank abundance

擴(kuò)增子分析中邑雅,rank abundance已經(jīng)是必備的一項(xiàng)分析內(nèi)容坑质,它可以從OTU的層面總體反應(yīng)出物種的分布情況(豐度和均勻度)硬梁,如下圖:


image.png

注:橫坐標(biāo):OTU 等級(jí)前硫,“500”代表樣本中按照豐度排列第500 位的OTU;縱坐標(biāo):該等級(jí)OTU 中序列數(shù)的相對(duì)百分含量荧止,即屬于該OTU 的序列數(shù)除以總序列數(shù)屹电,縱坐標(biāo)軸上數(shù)字阶剑,例如“100”代表相對(duì)豐度為100%,“10”代表相對(duì)豐度為10%危号,依次類(lèi)推牧愁。

Rank abundance構(gòu)建方法是統(tǒng)計(jì)單一樣本中,每一個(gè)OTU所含的序列數(shù)外莲,將OTUs按豐度(所含有的序列條數(shù))由大到小等級(jí)排序猪半,再以O(shè)TU等級(jí)為橫坐標(biāo),以每個(gè)OTU中所含的序列數(shù)(也可用OTU中序列數(shù)的相對(duì)百分含量)為縱坐標(biāo)做圖苍狰。

Rank abundance曲線(xiàn)可用來(lái)解釋多樣性的兩個(gè)方面,即物種豐度和物種均勻度烘绽。在水平方向淋昭,物種的豐度由曲線(xiàn)的寬度來(lái)反映,物種的豐度越高安接,曲線(xiàn)在橫軸上的范圍越大翔忽;曲線(xiàn)的形狀(平滑程度)反映了樣本中物種的均度,曲線(xiàn)越平緩盏檐,物種分布越均勻歇式。

Rank abundance計(jì)算過(guò)程:

  1. 獲取每個(gè)樣本中OTU的豐度值,例如從QIIME的otu table中獲取胡野。通俗來(lái)講材失,即搞清楚樣本(如上圖中Group1),中每個(gè)OTU它代表了多少條序列硫豆。
  2. 將每個(gè)樣本中OTU的豐度值按照從大到小順序進(jìn)行排序龙巨。注意:這步最精華的部分,因此整個(gè)分析過(guò)程叫做rank abundacne!
  3. 將每個(gè)樣本歸一化熊响,獲取每個(gè)樣本OTU的相對(duì)豐度旨别。例如樣本1中有5個(gè)OTU,豐度分別為5汗茄,4秸弛,3,3洪碳,5递览,那么排序后為5,5瞳腌,4非迹,3,3纯趋,總豐度即sum(5憎兽,5冷离,4,3纯命,3)=20,歸一化后的相對(duì)豐度為5/20=0.25西剥, 5/20=0.25, ... 。
    注意:樣本的相對(duì)豐度和其它樣本無(wú)關(guān)亿汞,例如樣本1中有5個(gè)OTU瞭空,樣本2中有7個(gè)OTU,那么樣本1的sum值為這5個(gè)OTU豐度的總和疗我,樣本的sum值為7個(gè)OTU豐度的總和咆畏。
  4. 對(duì)于每個(gè)樣本, x軸是OTU的種類(lèi)數(shù)吴裤,y軸是對(duì)應(yīng)的OTU的豐度旧找。例如樣本1中有5個(gè)OTU,豐度分別為0.2麦牺,0.25钮蛛,0.25,0.2剖膳,0.15魏颓,0.15(注意:這里是排好序的相對(duì)豐度)。那么按點(diǎn)(1,0.25),(2,0.25),(3,0.2),(4,0.15),(5,0.15)畫(huà)出梯形圖吱晒, 一般畫(huà)的都是下梯形(和正常的樓梯相似)甸饱。
  5. 一般會(huì)將y軸取對(duì)數(shù),要不然會(huì)相當(dāng)難看仑濒。同時(shí)將所有樣本y軸的刻度調(diào)整到一致柜候。這樣就完成了rank abundance圖的繪制。
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末躏精,一起剝皮案震驚了整個(gè)濱河市渣刷,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌矗烛,老刑警劉巖辅柴,帶你破解...
    沈念sama閱讀 211,194評(píng)論 6 490
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異瞭吃,居然都是意外死亡碌嘀,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,058評(píng)論 2 385
  • 文/潘曉璐 我一進(jìn)店門(mén)歪架,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)股冗,“玉大人,你說(shuō)我怎么就攤上這事和蚪≈棺矗” “怎么了烹棉?”我有些...
    開(kāi)封第一講書(shū)人閱讀 156,780評(píng)論 0 346
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)怯疤。 經(jīng)常有香客問(wèn)我浆洗,道長(zhǎng),這世上最難降的妖魔是什么集峦? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 56,388評(píng)論 1 283
  • 正文 為了忘掉前任伏社,我火速辦了婚禮,結(jié)果婚禮上塔淤,老公的妹妹穿的比我還像新娘摘昌。我一直安慰自己,他們只是感情好高蜂,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,430評(píng)論 5 384
  • 文/花漫 我一把揭開(kāi)白布聪黎。 她就那樣靜靜地躺著,像睡著了一般妨马。 火紅的嫁衣襯著肌膚如雪挺举。 梳的紋絲不亂的頭發(fā)上杀赢,一...
    開(kāi)封第一講書(shū)人閱讀 49,764評(píng)論 1 290
  • 那天烘跺,我揣著相機(jī)與錄音,去河邊找鬼脂崔。 笑死滤淳,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的砌左。 我是一名探鬼主播脖咐,決...
    沈念sama閱讀 38,907評(píng)論 3 406
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼汇歹!你這毒婦竟也來(lái)了屁擅?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 37,679評(píng)論 0 266
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤产弹,失蹤者是張志新(化名)和其女友劉穎派歌,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體痰哨,經(jīng)...
    沈念sama閱讀 44,122評(píng)論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡胶果,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,459評(píng)論 2 325
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了斤斧。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片早抠。...
    茶點(diǎn)故事閱讀 38,605評(píng)論 1 340
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖撬讽,靈堂內(nèi)的尸體忽然破棺而出蕊连,到底是詐尸還是另有隱情悬垃,我是刑警寧澤,帶...
    沈念sama閱讀 34,270評(píng)論 4 329
  • 正文 年R本政府宣布咪奖,位于F島的核電站盗忱,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏羊赵。R本人自食惡果不足惜趟佃,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,867評(píng)論 3 312
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望昧捷。 院中可真熱鬧闲昭,春花似錦、人聲如沸靡挥。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,734評(píng)論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)跋破。三九已至簸淀,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間毒返,已是汗流浹背租幕。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,961評(píng)論 1 265
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留拧簸,地道東北人劲绪。 一個(gè)月前我還...
    沈念sama閱讀 46,297評(píng)論 2 360
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像盆赤,于是被迫代替她去往敵國(guó)和親贾富。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,472評(píng)論 2 348

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