上一節(jié)介紹了兩組比較中如何展示差異基因數(shù)目队丝，這一節(jié)的內(nèi)容更進(jìn)了一步，我造了一種圖呢堰，可以用來表示具體的差異基因。
想象一下你的單細(xì)胞數(shù)據(jù)同時包含對照組和實驗組凡泣，在做完前期的基礎(chǔ)分析之后枉疼，得到了A B C D等共享的亞群皮假，現(xiàn)在你想知道在A亞群中，實驗組和對照組比較有哪些差異基因骂维。
類似的場景也可以是你的數(shù)據(jù)有多個腫瘤患者惹资，每個病人都有配對的原發(fā)灶、轉(zhuǎn)移灶航闺，現(xiàn)在想每個病人單獨看看原發(fā)褪测、轉(zhuǎn)移的差異基因。
怎么想到這個圖的
說到差異基因潦刃，火山圖大家一定都畫過侮措，但是如果有多個亞群、多個病人乖杠，那么多張火山圖又顯得冗雜了分扎。
再來看看這個圖
每一列是一個亞群，縱軸是log2FC胧洒，有正有負(fù)畏吓。一眼看上去也只能大致看出基因的多少，包含信息并不多卫漫。
如果把這兩種圖結(jié)合起來呢菲饼？就是這樣子的
這張圖怎么看
每個panel表示一個亞群或是一個病人，橫軸表示顯著性列赎，從中間向兩邊都是越來越顯著（整張圖右下角的箭頭就是說的這個意思）宏悦；縱軸表示log2FC，有正有負(fù)粥谬。中間兩條橫向虛線是log2FC的閾值肛根，只有當(dāng)p值和log2FC同時滿足一定閾值才會算作差異基因辫塌。
所以每個panel里面只有右上角漏策、左下角才是差異基因，才會被涂上顏色臼氨。另一方面掺喻，右下角、左上角空出來了储矩，可以把基因label加上去感耙，充分利用空間，這些基因名也是對齊的持隧。
這張圖的難點有以下幾個：
- 圖形的整體構(gòu)思
- 基因名如何對齊
- 分面還是拼圖的選擇（最終選擇了拼圖）
獲取代碼
這個系列都會采取限時公開的方式共享代碼即硼，24小時內(nèi)是免費的。超過這個時間如何獲取屡拨，后臺回復(fù)2022A可知（可能需要你動動小手轉(zhuǎn)發(fā)一下）只酥。
代碼是從count開始的褥实，即使完全不懂單細(xì)胞，也能把圖畫出來裂允，比較友好损离，希望對大家有用。代碼和測試數(shù)據(jù)的網(wǎng)盤鏈接如下：
- 鏈接：https://pan.baidu.com/s/1aWkC9Ddj_RAWp0ecy6juTA
- 提取碼：94r7
