实用干货 | 轻松绘制Circos图——环形热图-产品资讯-资讯-生物在线

实用干货 | 轻松绘制Circos图——环形热图

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

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

相关咨询

New Phytol.(IF=8.1)|安农大宋传奎教授团队最新研究刷新认知!甜菜碱竟能“解锁”茶树抗寒密码? (2025-07-30T00:00 浏览数:35)

IF=6.2!中科院茶叶研究所林智团队突破性成果:非靶向特定修饰代谢组学首次系统解析植物糖基化代谢物 (暂无发布时间 浏览数:49)

会议预告 | 百趣生物诚邀您参加2025中国植物生理与植物分子生物学学会 (暂无发布时间 浏览数:2537)

Metabolism (IF=11.9)|中国药科大学团队:线粒体“泄漏”为脑缺血帮凶,揭示神经酰胺介导的mtDNA炎症机制! (暂无发布时间 浏览数:3235)

Food Res Int. (IF=8.0)|合工大徐宝才教授团队破解细胞级抗衰密码:家禽血“变身”超级抗氧化剂 (暂无发布时间 浏览数:2909)

Food Chem. (IF=9.8)重磅发文|上海农科院团队揭示黄桃从“青草味”到“浓郁果香”的风味转化机制 (暂无发布时间 浏览数:4023)

IF=9.2!肠乳杆菌在他汀治疗中的屏障保护与抗纤维化机制 (暂无发布时间 浏览数:5226)

磁珠富集的前分析「陷阱」:血细胞、离心与抗凝剂对血浆蛋白组学的多重干扰 (暂无发布时间 浏览数:5355)

JACS(IF=14.4)|复旦陆豪杰教授团队攻克Tn抗原检测难题:一步酶促策略如何实现糖基化位点的精准Mapping? (暂无发布时间 浏览数:5136)

nature aging(IF=17.0)|浙大黄河教授团队解码造血干细胞 “衰老密码”:代谢图谱锁定尿苷并揭示功能逆转新机制 (暂无发布时间 浏览数:5295)

ADVERTISEMENT