基因组数据展示离不开圈图,科研日常同样离不开圈图,下面就来共同学习经典的圈图吧。
1.数据准备
rm(list = ls())
library(circlize)
set.seed(1234)
data <- data.frame(x=rnorm(432,0,1),
y=runif(432,min = 2,max=6),
fac=sample(letters[1:8],432,replace = T))
head (data)
2.开始画图
attach(data)
circos.initialize(factors = fac,x=x) #图形布局
circos.track(factors=fac, #分类变量
y=y,
panel.fun=function(x,y){
circos.text(CELL_META$xcenter,
CELL_META$cell.ylim[1:5]+uy(25,"mm"),
CELL_META$sector.index,
cex = 1.0,
col = "black")
circos.axis(labels.cex = 0.6)
})
circos.trackPoints(fac,x,y,col=as.numeric(factor(fac)))
circos.trackHist(fac,x,col=as.numeric(factor(fac)))
detach(data)
结果解读:图形被分为8个部分,每部分代表fac中的一种类型,通过这个“甜甜圈”我们可以清晰的看出每个fac中y的频数分布。
备注:本推文学习整理自环状直方图和散点图简书推文,仅供学习参考,如有侵权,请联系删除。
参考链接:
环状直方图和散点图
网友评论