Nature子刊經典網絡圖在R中的繪製,網絡模塊性可視化

2021-02-13 土壤微生物組

2019年Wagg等發表在Nature Communications上的文章"Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning"是研究微生物群落和生態系統多功能性的典範,通過構建微生物網絡的功能模塊,表明了微生物互作網絡複雜性的提高能夠促進生態系統功能。其中的網絡圖無論是畫法還是配色都值得學習,小編在這裡嘗試用R畫了文章網絡圖,以供大家學習!!

NC原文網絡圖

重現網絡圖

(顏色設置不太好,可自行更改)


數據分析代碼

################ 導入數據

otu<-read.csv("D://PJJ/test_otu.csv",row.names = 1)

design<-read.csv("D://PJJ/test_design.csv")

#### 加載包

library(psych)

library(igraph)

##加載函數 (可後臺獲取)

setwd("D:/R_input")

source("CorrDF.R")

source("cor.mtest.R")

## 進行相關分析

aa <- rcorr(t(otu[1:500,6:14]), type=c("spearman"))

## Create data frame of co-occurring OTUs

net <- CorrDF(aa$r, aa$P)

net$padj <- p.adjust(net$p, method="none")

## 設置閾值

net_padj <- net[which(abs(net$cor) >0.97),]

net_padj <- net_padj[which(abs(net_padj$padj) < 0.01),]

######構建網絡文件

nodeattrib_root_combine <- data.frame(node=union(net_padj$from, net_padj$to))

nodeattrib_root_combine$indicgroup <- 0

all_root_net <- graph_from_data_frame(net_padj, direct=F, vertices=nodeattrib_root_combine)

####可視化

igraph.weight = E(all_root_net)$weight

E(all_root_net)$weight = NA

fc = cluster_fast_greedy(all_root_net,weights =NULL)

##計算模塊性

modularity = modularity(all_root_net,membership(fc))

comps = membership(fc)

##設置模塊顏色

colbar=c("#FCBF26","#B0E2FF","#EE4A25","#BCEE68","#FCBF26","#B0E2FF","#BCEE68","#EE4A25","#FCBF26","#B0E2FF","#EE4A25","#FCBF26","#FCBF26","#B0E2FF","#EE4A25","#BCEE68","#FCBF26","#B0E2FF","#EE4A25","#B0E2FF","#EE4A25","#B0E2FF","#EE4A25","#EE4A25","#B0E2FF","#EE4A25","#B0E2FF","#EE4A25","#EE4A25")

V(all_root_net)$color = colbar[comps] 

#出圖

set.seed(9939)

coords<-layout_(all_root_net,with_fr(niter=99, grid="nogrid"))

plot(all_root_net, vertex.label=NA, edge.width=0.01,vertex.size=6, layout=layout.fruchterman.reingold)

文章概覽

  土壤微生物群高度多樣化,包括多達四分之一的地球生物多樣性。然而,這樣一個多樣化和功能複雜的微生物群如何影響生態系統的功能仍然不清楚。在這裡,我們進行了草地生態系統中的土壤微生物試驗,並觀察到微生物群落的多樣性和微生物網絡的複雜性影響了多個生態系統功能。微生物網絡結構單一、微生物豐富度降低的草地微宇宙,由於支持相同功能的微生物少和較低具有獨特功能的微生物,表現出最低的多功能性。此此外,不同的微生物類群解釋了不同的生態系統功能,指出了微生物群落功能多樣性的意義。這些發現表明真菌和細菌群落內部和之間微生物相互作用對提高生態系統性能的重要性,並表明複雜微生物互作網絡的消失會損害生態系統功能。

Wagg, C. and K. Schlaeppi, et al. (2019). "Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning." Nature Communications 10 (1).

相關焦點

  • 一文學會網絡分析——Co-occurrence網絡圖在R中的實現
    目前生態學領域大家用到的網絡圖多為基於群落數據相關性構建的Co-occurrence網絡圖。此類網絡可以採用R中igraph包構建並實現出圖。當然,除此之外,還有一些非命令行的軟體,例如cytoscape,gephi,pajek等。但我認為,對於R使用者來說,通過R做圖還是最方便的。大致的流程如下圖所示:1)根據觀察,實驗或者相關性推斷來確定物種間的聯繫。
  • 如何繪製漂亮的網絡圖?
    (Soil Biology & Biochemistry, 2018)說到網絡圖的繪製,大家可能首先想到的是Cytoscape,它在生物學相關的領域應用非常廣泛。Gephi最常見是用於社交網絡數據分析,如下,是我手工整理的《天龍八部》中部分人物的關係網絡。我們可以看到段王爺是如何實力搶鏡力壓三大男主的,儘管金庸老爺子心中的男一號可能是段譽。
  • 雙代號網絡圖的繪製
    雙代號網絡圖是目前推廣應用較為廣泛的一種網絡計劃形式,它是以箭線及其兩端節點的編號表示工作的一種網絡圖。(2) 雙代號網絡圖中不允許出現一個以上的起始節點或終點節點。(3)在網絡圖中嚴禁出現循環迴路。右圖中工作C、D、E 形成了閉合迴路,則這個網絡圖肯定是錯誤的。
  • 網絡圖怎麼畫?手把手教你快速繪製高大上的網絡圖
    網絡拓撲圖是由箭頭線條、節點、路線三個因素組成形如網絡結構的圖示類型。網絡工程師繪製網絡圖,通常將網絡節點設備和通信介質進行物理布局,以清晰地展示網絡傳輸方式和途徑。網絡拓撲圖有許多種結構類型,通常有星型、環型、總線型、樹型、網狀、分布式結構、等等。
  • 進度控制知識:網絡圖的繪製規則
    在繪製雙代號網絡圖時,一般應遵循以下基本規則:(1)網絡圖必須按照已定的邏輯關係繪製。由於網絡圖是有向、有序網狀圖形,所以其必須嚴格按照工作之間的邏輯關係繪製,這同時也是為保證工程質量和資源優化配置及合理使用所必需的。(2)網絡圖中嚴禁出現從一個節點出發,順箭頭方向又回到原出發點的循環迴路。
  • Nature又上線2本新子刊!
    你的機會來了——2019年2月,Nature先後推出了兩本新子刊,分別是 Nature Food 和Nature Cancer。Nature雜誌的子刊達到了53本,其中20本為綜述期刊。2019年2月24日,Nature官網正式上線了一個新子刊:Nature Food。這是一本在線期刊,將於2020年1月正式啟動。
  • 網絡分析概述之網絡基礎簡介
    在數學中,「網絡」(networks)通常被稱為「圖」(graphs),關於圖研究的數學領域稱為「圖理論」(graph theory)。網絡圖中的基本元素是節點和邊,構建網絡圖時,將對象稱為「節點」(vertices或nodes),並通常將它們繪製為點;將節點之間的連接稱為「邊」(edges或links),並且通常將它們繪製為點之間的線。如此,得到了網絡圖。
  • 最佳大數據可視化技術
    基本技術如下圖:線圖線圖是最簡單的技術,用於繪製一個變量與另一個變量之間的關係或依存關係。 要繪製兩個變量之間的關係,我們可以簡單地調用plot函數。箱形圖通常用於了解數據中的異常值。可以用於半結構化或非結構化數據的另一種可視化技術是網絡圖。
  • 單代號網絡圖講解
    由一個節點表示一項工作,以箭線表示工作順序的網絡圖稱為單代號網絡圖。單代號網絡圖的邏輯關係容易表達,且不用虛箭線,便於檢查和修改。但不易繪製成時標網絡計劃,使用不直觀。一、單代號網絡圖的基本構成與符號單代號網絡圖中,箭線表示相鄰工作之間的邏輯關係。一個節點表示一項工作,一般用圓圈或矩形表示。
  • Nature子刊:教你零基礎開展微生物組數據分析和可視化
    SDP的獨特功能是其基於模塊,途徑和代謝網絡的功能注釋。MicrobiomeAnalyst使用戶可以輕鬆地可視化這些功能在樣本和研究條件中的分布。它還支持顯式統計檢驗以識別豐富的功能。用戶可以在新陳代謝網絡環境中以交互方式瀏覽結果,以進一步了解功能。
  • Nature子刊:太陽能電池效率或大幅提升!
    論文連結:https://www.nature.com/articles/s41557-019-0297-7現代太陽能電池板採用工作原理基本一樣:一個光子產生一個激子,激子轉化為電流。該研究成果不僅可用於下一代太陽能設備生產,還可推進化學、傳感器和成像中的光催化過程,用以製造藥品、塑料和許多其他類型的消費化學品。
  • 根據下列邏輯關係表繪製的雙代號網絡圖如下圖所示,其存在的錯誤
    根據下列邏輯關係表繪製的雙代號網絡圖如下圖所示,其存在的錯誤   根據下列邏輯關係表繪製的雙代號網絡圖如下圖所示,其存在的錯誤是()
  • 從零開始學Python可視化(二):掌握統計圖形的繪製
    這次我們會嘗試更多常見的統計圖形,比如條形圖、直方圖、餅圖等,我們的目標是掌握工作學習過程中使用最頻繁的圖形技能。1. 柱狀圖(條形圖) —— bar() / barh()柱狀圖又叫條形圖,用於繪製定性(分類)數據的分布特徵,比如不同國家的GDP、不同年齡段的平均體重等。這次,我們以一組學生的身高來演示它的作圖方法。
  • Python數據可視化實例之繪製圖表
    Python數據可視化實例之繪製圖表原創 蟲蟲安全 2018-09-05 17:41:57得利於語言的簡單明了、豐富的數據結構、豐富的類和模塊,Python如今成了數據科學中的香餑餑,成了matlab、R語言之外又一強大的數據分析工具。拋開其他方面的、今天蟲蟲帶大家一起來探索Python在數據可視化方面的應用。
  • 二建考試單代號網絡圖、雙代號網絡圖怎麼學?竅門在這裡!
    答:包餃子的緊前工作有兩個,分別是剁餡和擀餃子皮三、雙代號網絡計劃的繪圖規則1、雙代號網絡圖必須正確表達已定的邏輯關係2、雙代號網絡圖中,嚴禁出現循環迴路>3、雙代號網絡圖中,在節點之間嚴禁出現帶雙向箭頭或無箭頭的連線4、雙代號網絡圖中,嚴禁出現沒有箭頭節點或沒有箭尾節點的箭線5、當雙代號網絡圖的某些節點有多條外向箭線或多條內向箭線時,為使圖形簡潔,可使用母線法繪製
  • Nature子刊封面:哈佛等開發微型手術機器人
    Wood) 基於摺紙結構研發了一款超輕便,高精度的顯微外科手術機器人——mini-RCM,該項研究發表於nature旗下新子刊《nature machine intelligence(自然-機器智能)》,並登上該期刊8月份封面。
  • Nature子刊,今年首個IF就突破12分,明年或可衝18分
    來自WoS截圖今年是Nature Sustainability拿到的第1個IF,作為Nature子刊,首個IF就超12分,確實值得我們關注~期刊基本信息關於此刊的具體情況,大家可以看看:https://www.nature.com/natsustain/
  • R可視化 | 馬賽克圖
    作者:莊閃閃來源:莊閃閃的成長手冊馬賽克圖(mosaic plot),顯示分類數據中一對變量之間的關係,原理類似雙向的100%堆疊式條形圖,但其中所有條形在數值/標尺軸上具有相等長度,並會被劃分成段。可以通過這兩個變量來檢測類別與其子類別之間的關係。主要優點馬賽克圖能按行或按列展示多個類別的比較關係。主要缺點難以閱讀,特別是當含有大量分段的時候。此外,我們也很難準確地對每個分段進行比較,因為它們並非沿著共同基線排列在一起。
  • Nature子刊封面:哈佛大學等研發的微型手術機器人
    Wood) 基於摺紙結構研發了一款超輕便,高精度的顯微外科手術機器人——mini-RCM,該項研究發表於nature旗下新子刊《nature machine intelligence(自然-機器智能)》,並登上該期刊8月份封面。
  • 單代號網絡圖,你還記得多少?
    由一個節點表示一項工作,以箭線表示工作順序的網絡圖稱為單代號網絡圖。單代號網絡圖的邏輯關係容易表達,且不用虛箭線,便於檢查和修改。但不易繪製成時標網絡計劃,使用不直觀。1、單代號網絡圖的基本構成與符號單代號網絡圖中,箭線表示相鄰工作之間的邏輯關係。一個節點表示一項工作,一般用圓圈或矩形表示。