火山圖 Volcano plot
在統計學上，火山圖是一種類型的散點圖，被用於在大數據中快速鑒定變化。由於它的形成像火山噴發的樣子，所以被稱為火山圖。和上文講的曼哈頓圖類似。
火山圖基本元素
火山圖也有很多種樣式，在生物學高通量測序結果中，常見有的X和Y軸分別為aboundance 和 fold-change，或p-value和fold-change兩種樣式。如上圖中為p-value versus fold-change的樣式,，先此為例進行圖中基本元素解讀：
- X軸：通常為兩組基因表達或OTU相對豐度比較的變化倍數取2的對數log2(Fold change)。 取log2的原因：由於兩數相比，倍數取值范圍為0 - 正無窮，上調的倍數為1 - 正無窮，比較容易理解；而下調倍數為0-1很難理解。如果取對數可將倍數以1為中心左右對稱，取2的對數是因為經驗上2倍差異比較常用，故取log2變換既可以使上調或下調數據對稱(上、下調的目標同等重要，無偏好)，而且坐標軸1個單位刻度的變化即為2倍差異(方便快速篩選足夠差異的點)；
- Y軸：兩組基因表達或OTU相對豐度相比統計上的顯著性Pvalue值，通常取負對數10變換 -log10(P-value)。取負對數與轉換下調表達數值原理類似，將0-1范圍pvalue數據變為0-正無窮，且將越小越顯著的P-vale轉換為與顯著性正相關，方便觀察顯著差異目標。Pvalue取負對數10轉換后，坐標軸數值2，3對應0.01和0.001的顯著閾值，便於理解。有時也用-ln或-log2對Pvalue值進行轉換。
- 圖中點：代表每個基因或OTU在兩組相比時差異倍數和顯著性的值；通常將顯著變化並且差異倍數大於指定閾值的點規定為顯著差異的結果，常用顏色高亮顯示(圖中綠色點)。有些重點關注的點可以添加箭頭或標簽進一步突出顯示(圖中綠色箭頭)。有時還會有點大小代表相關豐度，點的形狀代表其它一些屬性，如物種分類信息等。
火山圖繪制
常用R語言的ggplot2繪制。接一來的統計和繪圖課會有詳細源代碼和分析實例。
看圖實戰(Result)
示例1. 各比較組間上/下調顯著差異OTU的數量
這篇文章分析了水稻根不同區域的細菌組成，16S分析文章較系統的作品，兩年被引用147次，推薦閱讀
圖2.A 水稻根不同生態位相對於土壤中顯著差異的OTU
1. 圖中元素解釋：不同於上面例子，此圖為Fold change versus Average Abundance樣式，大家注意一下，常見就這兩種。
- X軸標簽為Log10(Average Abundance)，是相對豐度平均值 取10的對數，這里用X軸展示了OTU的豐度，而不是pvalue值。其實作者X標題可能將Log2錯寫為了Log10，因為OTU相對豐度分析常用百分、千分和百萬分數，即使是最大的百萬分數取Log10最大值僅為6，而不會坐標軸刻度至15，這種0-15的刻度是典型的Log2(RPM)百萬分數的數據。
- Y軸標簽為Log10(Fold Change)是指兩組間比較差異倍數經10的對數變換；標簽也可能是錯誤的，因為想讓倍數對稱必須取負對數，而且通常用-Log2，Log10的值是10時就有上百億倍，擴增子不可能有這么大的差異或深度(擴增子差異倍數為了防止分母為零需要加某值，因此也減小了差異倍數，不會有非常大的情況；擴增子測序深度一般只有幾千至幾十萬，不可能有上億的量)。
- 圖中的每個點代表一個OTU，其中圖上部有顏色的點為顯著富集(enriched)的OTU，下部有顏色的點為顯著下降(depleted)的OTU，黑灰色的點為沒有顯著差異的OTU；圖中標記的數字為顯著富集或下降的OTU的數量；
2. 圖表結果：圖中展示了從根際、根表到根內與土壤相比，所有OTU的相對豐度和差異倍數，其中高亮顯示顯著富集或下降差異OTU和數量；
3. 圖表結論和規律：從根際-根表-根內，從外到內，顯著差異的OTU總數量是逐漸上升的，代表變化越來越大；其中根際有較多富集的OTU，而根表和根內則較多為下降的OTU；
4. 經驗和技巧：本圖雖然有X/Y軸標簽描述有待商榷，但許多優點值得學習；
1. 大家經常看到的看火山圖都是豎起的，左、右部分高亮的為差異OTU，而作者將圖像放倒(順時針旋轉90度)，以Y軸零刻度為分界面，上部為富集OTU，下部為下降OTU也很直觀，而且耳目一新；
2. 圖中差異OTU的數量標出，不用讀者自己去估計，直接明了，感覺信息量更豐富；
3. 三組圖並排且變化明顯，不用多説，讀者自然會去比較差異，引導讀者參與分析比較和規律發現。