接下來兩節(jié)會介紹兩種比較新的圖形枚赡，可能大多數(shù)做單細(xì)胞分析的小伙伴
沒有見過。
今天講的這種圖形是用來展示差異基因數(shù)目的，比如你的單細(xì)胞數(shù)據(jù)同時包含對照組和實驗組买羞，在做完前期的基礎(chǔ)分析之后非春，得到了A B C D等共享的亞群柱徙，現(xiàn)在你想知道在A亞群中缓屠，實驗組和對照組比較有哪些差異基因，差異基因數(shù)目是多少护侮。
下面這種圖就很適合用來展示差異基因數(shù)目：
這張圖是這么看的：左側(cè)圖是Arep組相對于Brep組高表達(dá)的基因敌完，每一列是一個亞群，每一行是一個基因羊初；基因又分為兩個panel滨溉，下層panel是unique基因，即該基因只在一個亞群的Arep VS Brep比較中出現(xiàn)长赞，上層panel是share的晦攒，即該基因至少在兩個亞群的Arep VS Brep比較中出現(xiàn)。圖中數(shù)字表示對應(yīng)的基因數(shù)目得哆。右側(cè)圖是Arep組相對于Brep組低表達(dá)的基因脯颜。
這張圖整體比較簡單，難點主要是橫軸亞群贩据、縱軸基因的排序栋操。我的代碼中提供的排序規(guī)則是我自己的想法，排序方案并不唯一饱亮。
圖形改編自這篇文獻(xiàn)：
- A Single-Cell Transcriptomic Atlas of Human Skin Aging
- https://doi.org/10.1016/j.devcel.2020.11.002
獲取代碼
這個系列都會采取限時公開的方式共享代碼矾芙，24小時內(nèi)是免費的。超過這個時間如何獲取近尚，公粽號后臺回復(fù)2022A可知（可能需要你動動小手轉(zhuǎn)發(fā)一下）蠕啄。
代碼是從count開始的，即使完全不懂單細(xì)胞戈锻，也能把圖畫出來歼跟，比較友好，希望對大家有用格遭。代碼和測試數(shù)據(jù)的網(wǎng)盤鏈接如下：
- 鏈接：https://pan.baidu.com/s/114toYWT2MlpbJilLkgs_ZQ
- 提取碼：hx51
