實用干貨 | 輕松繪制Circos圖——環形熱圖-產品資訊-資訊-生物在線

實用干貨 | 輕松繪制Circos圖——環形熱圖

作者:上海阿趣生物科技有限公司 2021-11-25T11:05 (訪問量:9007)

Arial, sans-serif; font-size: 15px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: 1.5px; orphans: 2; text-align: justify; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(241, 244, 245); text-decoration-style: initial; text-decoration-color: initial;">##代碼3

col_fun1 = colorRamp2(c(-2, 0, 2), c("#26B9CB", "#FFFFFF", "#B72865"))##設置熱圖顏色

column_od = hclust(dist(t(mat1)))$order #對列聚類

circos.par(gap.after = c(10))##為添加列名留出空間

circos.heatmap(mat1[, column_od], ##將列聚類后重新排序的矩陣

col = col_fun1, ##設置顏色

dend.side = "inside",##樹狀圖在圈內

rownames.side = "outside",##行名在圈外

dend.track.height = 0.2,

dend.callback = function(dend, m, si) {

# when k = 1, it renders one same color for the whole dendrogram

color_branches(dend, k = 4, col = 2:5)##對樹狀圖進行著色

}

)

circos.track(track.index = 2, ##將列名添加在第二個軌道(就是熱圖所在的環形軌道)

panel.fun = function(x, y) {

if(CELL_META$sector.numeric.index == 1) { # the last sector

cn = colnames(mat1[, column_od])##取得列名

n = length(cn)

circos.text(rep(CELL_META$cell.xlim[2], n) + convert_x(1, "mm"), ##x軸坐標

1:n - convert_y(0.5, "mm"), ##y軸坐標

cn, ##輸入要展示的列名

cex = 0.25, ##列名的大小

adj = c(0, 0.5),

facing = "inside")

}

}, bg.border = NA)


circos.clear()

```


圖2


2> 修改文字的大小(rownames.cex)。字體(rownames.font)和顏色(rownames.col),使用代碼4,結果見圖3。


```

##代碼4

col_fun1 = colorRamp2(c(-2, 0, 2), c("#26B9CB", "#FFFFFF", "#B72865"))##設置熱圖顏色

column_od = hclust(dist(t(mat1)))$order #對列聚類

circos.par(gap.after = c(15))##為添加列名留出空間

circos.heatmap(mat1[, column_od], ##將列聚類后重新排序的矩陣

col = col_fun1, ##設置顏色

dend.side = "inside",##樹狀圖在圈內

rownames.side = "outside",##行名在圈外

rownames.col = 1:nrow(mat1) %% 10 + 1,##行名添加顏色

rownames.cex = runif(nrow(mat1), min = 0.3, max = 2),##行名的大小

rownames.font = 1:nrow(mat1) %% 4 + 1,##行名字體

dend.track.height = 0.2,

dend.callback = function(dend, m, si) {

# when k = 1, it renders one same color for the whole dendrogram

color_branches(dend, k = 4, col = 2:5)##對樹狀圖進行著色

}

)

circos.track(track.index = 2, ##將列名添加在第二個軌道(就是熱圖所在的環形軌道)

panel.fun = function(x, y) {

if(CELL_META$sector.numeric.index == 1) { # the last sector

cn = colnames(mat1[, column_od])##取得列名

n = length(cn)

circos.text(rep(CELL_META$cell.xlim[2], n) + convert_x(1, "mm"), ##x軸坐標

1:n - convert_y(0.5, "mm"), ##y軸坐標

cn, ##輸入要展示的列名

cex = 0.25, ##列名的大小

adj = c(0, 0.5),

facing = "inside")

}

}, bg.border = NA)

circos.clear()

```


圖3


參考文獻

[1] Gu, Z. (2014) circlize implements and enhances circular visualization in R. Bioinformatics. DOI: 10.1093/b

地 址: 嘉定區新培路51號焦點夢想園5層

聯系人: 高小姐

電 話: 400-664-9912

傳 真:

Email:marketing@biotree.cn

相關咨詢

Circulation (IF=38.7)|上海兒童醫學中心張浩/劉一為聯合南京鼓樓醫院王東進團隊破解肥胖心衰患者LVAD療效困局,胰島素增敏劑點亮心肌恢復新希望! (2026-04-25T00:00 瀏覽數:7957)

J. Adv. Res. (IF=13)|血漿蛋白質組學賦能:暴發性心肌炎精準診斷與靶向治療新突破 (暫無發布時間 瀏覽數:7476)

NP(IF=8.1)|西北農林研究團隊解碼黑枸杞“高品質密碼”:光響應轉錄調控網絡如何驅動類黃酮合成與藥用活性? (暫無發布時間 瀏覽數:7501)

LWT(IF=6.6)|四川大學研究團隊通過同位素標記代謝流與脂質組學的“夢幻聯動”,帶你看清酵母的耐鹽策略 (暫無發布時間 瀏覽數:8163)

破微量極限,啟代謝新程|微量NGM 3 Pro,以微量樣本,全景解碼生命代謝 (暫無發布時間 瀏覽數:8179)

JEP (IF=5.4)|福建中醫藥大學研究團隊:多組學揭秘丹瓜方調控GLUD1改善糖尿病代謝機制 (暫無發布時間 瀏覽數:9409)

Food Chem.(IF=9.8)|南農陳暄教授研究團隊多組學揭秘:云南“金花茶”的風味煉金術—匍匐散囊菌如何點“茶”成“香”? (暫無發布時間 瀏覽數:9677)

J.Agric.Food Chem.(IF=6.2)|中國熱科院胡偉/丁澤紅團隊突破性成果!解碼木薯塊根3類蛋白修飾圖譜,解鎖淀粉積累與抗逆調控密碼 (暫無發布時間 瀏覽數:9864)

nat commun(IF=15.7)|北大白玉教授團隊突破靜態局限:動態單細胞代謝組學揭秘腫瘤-巨噬細胞互作新機制 (暫無發布時間 瀏覽數:11428)

nat.commun.(IF=15.7)|復旦大學附屬中山醫院詹成團隊:靶向多胺代謝提升KRAS靶向藥療效 (暫無發布時間 瀏覽數:12678)

ADVERTISEMENT