還是繼續(xù)復(fù)現(xiàn)NC圖表，這幅圖展示的是相關(guān)性分析結(jié)果，可以叫做氣泡熱圖您没，我們之前說過熱圖也可以這么做（熱圖5：ggplot2畫熱圖及個性化修飾）茎刚。相關(guān)性分析這里就不展示了，我們主要是復(fù)現(xiàn)這個圖的做法蒲列，原文作者在supplement里提供了這個圖的作圖數(shù)據(jù)，直接使用搀罢。這個圖的特點是一幅圖里面即顯示了正相關(guān)蝗岖、也顯示了負相關(guān)關(guān)系，以及P值榔至，而且相關(guān)系數(shù)也用氣泡大小表示出來抵赢。還有一個是P值的體現(xiàn)用的是離散型的數(shù)據(jù)。
（Reference：9p21 loss confers a cold tumor immune microenvironment and primary resistance to immune checkpoint therapy）
跟著我們公眾號這么久了，我就不信你不會氣泡圖铅鲤。這個圖的重點在于上下調(diào)的展示划提，這就需要構(gòu)建數(shù)據(jù)，也就是添加新的分組列邢享，方法其實在火山圖系列講過（轉(zhuǎn)錄組不求人系列(十): NCS級別的火山圖鹏往，總有一款適合你！）骇塘，例如添加上下調(diào)基因標(biāo)識伊履，一定要記住，東西是相通的款违，接下來我們從添加組開始唐瀑。
對了，我的復(fù)現(xiàn)結(jié)果如下奠货，和原圖是一樣的介褥，只是基因排列不同：
加載ggplot2，讀入相關(guān)分析結(jié)果數(shù)據(jù)递惋。需要添加兩列柔滔，第一列是分出正調(diào)控和負調(diào)控。第二列是將P值分等級萍虽。最后相關(guān)系數(shù)取絕對值睛廊，因為相關(guān)性強弱和正負沒有關(guān)系，原圖也是用正值表示相關(guān)性的杉编。
更多精彩內(nèi)容請至我的公眾號《KS科研分享與服務(wù)》
